Low Nickel, Ceria Zirconia-Based Micro-Tubular Solid Oxide Fuel Cell: A Study of Composition and Oxidation Using Hydrogen and Methane Fuel

: The study examines the effect of using low nickel (Ni) with high ceria (CeO 2 ) anode content towards the oxidation of H 2 and CH 4 fuel by evaluating the activation energy of the ohmic process and charge transfer process. Using a micro-tubular solid oxide fuel cell (MT-SOFC), the anodes are made up of 50% YSZ with varying NiO:CeO 2 percentages from 0% NiO, 50% CeO 2 to 50% NiO, 0% CeO 2 . The performance is measured based on maximum power density (MPD), electrochemical impedance spectroscopy (EIS) and activation energy, E a of the ohmic (Rohm) and charge transfer (R ct ) processes. We found that by lowering the Ni content to lower than 50% NiO, anode conductivity will drop by 7-fold. An anode containing 37.5% NiO, 12.5% CeO 2 yield MPD of 41.1 and 2.9 mW cm − 2 when tested on H 2 and CH 4 fuels thus have the lowest Ni content without an abrupt negative effect on the MPD and EIS. The signiﬁcant effect of conductivity drops on MPD and EIS are observed to occur at 25% NiO, 25% CeO 2 and lower NiO content. However, anode content of 25% NiO, 25% CeO 2 has the lowest E a for R ct (29.74 kJ mol − 1 ) for operation in CH 4 , making it the best anode composition to oxidize CH 4 . As a conclusion, an anode containing 25% NiO:25% CeO 2 :50% YSZ and 37.5% NiO:12.5% CeO 2 :50% YSZ shows promising results in becoming the low Ni anode for coking-tolerant SOFC.


Introduction
The emission of greenhouse gases, mainly carbon dioxide, from fossil fuel burning has led to climate changes where the global average temperature has been rising by 0.08 • C per decade since 1880 [1]. This emission is mainly attributed to electricity generation by coal burning, involving multiple steps of energy conversion. The stored chemical energy in and complete cells of anode, electrolyte and cathode can be fabricated as a triple-layer hollow fiber (TL-HF).
In anodes, the charge transfer process involves a reaction between oxide ion and fuel that leads to oxidation of the fuel and release of electrons, which generates electricity. The study of the charge transfer process is important since the process depends greatly on the catalytic properties of the anode material [10]. Similar to every other chemical process, reactions between fuel and oxide ions on the anode surface require the reactant to overcome an energy barrier, termed activation energy or E a . Lower E a will result in more reactants having minimum energy requirements and hence more reactions occurring, thus increasing the current density of the fuel cell. Investigating the E a value of the charge transfer process can be achieved via an Arrhenius plot of the R ct value at different operating temperatures.
The objective of this study is to investigate the effect of CeO 2 and Ni content in the anode of MT-SOFC on that of the oxidation of H 2 and CH 4 fuel by using evaluation of the activation energy of the ohmic process and charge transfer process. Unlike other studies, CeO 2 was fabricated together with the NiO-YSZ anode to produce the Ni-CeO 2 -YSZ anode with a YSZ electrolyte. The two layers were fabricated together in the form of dual-layer hollow fiber using a co-extrusion-based phase inversion process. The MT-SOFC uses pure YSZ electrolytes, and the anode is composed of 50% YSZ, with the remaining 50% made up of mixtures of Ni and CeO 2 . The Ni and CeO 2 content were varied to investigate the effect of the Ni and CeO 2 content on the physical structure such as morphology, gas permeability and mechanical strength. The effect was also investigated in terms of anode conductivity, power density and impedance spectroscopy in H 2 and CH 4 fuel at varying temperatures. Finally, analysis of activation energy, E a , by means of an Arrhenius plot was performed on the ohmic resistance, R ohm , and charge transfer resistance, R ct , from the EIS of the SOFC obtained at temperatures from 700 to 900 • C on H 2 and CH 4 fuel. This is to examine the effect of temperature on the ohmic and charge transfer process of the varying Ni and CeO 2 content. The produced MT-SOFC was expected to increase the performance of MT-SOFC on operating with CH 4 fuel while using a low content of Ni. This work will enable MT-SOFC to operate stably using methane fuel, which can later be extended to hydrocarbon-based renewable fuel such as biogas as well.

Fabrication of DL-HF
The fabrication of anode/electrolyte DL-HF was performed as mentioned in [13,14]. Briefly, the suspension of the inner and outer layer of DL-HF was prepared separately. The outer layer made up the electrolyte layer entirely and consists of YSZ, while the inner layer formed the anode layer and was made up of mixtures of NiO, CeO 2 and YSZ as the ceramic component, as shown in Table 1. Five compositions of anodes were prepared with varying NiO and CeO 2 content, named Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. The electrolyte layer was made of 100YSZ as the ceramic component. The suspension was made up of 65% ceramic, 6.5% polyethylethersulfone (PESf) as polymer binder, 0.12% Arlacel P135 as surfactant and 28.38% N-Methyl-2-pyrrolidone (NMP) as solvent. The ceramic components were mixed with the solvent and dispersant in a planetary ball milling machine (Magna Value), then ball milled for 24 h, followed by the addition of polymer binder and further ball milled for 24 h. The suspensions were pumped through an orifice opening called spinneret, as shown in Figure A1 in Appendix A.
The flow rates used were 3, 10 and 10 mL min −1 each for the outer layer, inner layer and water acting as an internal coagulant into the water bath, which acts as external coagulant at an air gap of 10 cm. Phase inversion follows in which solvent/non-solvent exchange occurred, leaving a solidified polymer phase with ceramic powders. The produced material, called DL-HFs, was soaked in the water and air-dried for 24 h. Then, the produced DL-HF was calcined at 600 • C for 1 h (ramp up rate = 5 • C min −1 ) followed by sintering at 1450 • C for 8 h (ramp up rate 10 • C min −1 , ramp down rate 5 • C min −1 ) using a tube furnace (Magna Value). Both calcination and sintering were performed using air as the atmosphere. These ceramic DL-HFs produced were subjected to characterization analysis.

Characterization
Unreduced ceramic DL-HFs were scanned using SEM-EDX (Hitachi TM3030) to examine the morphology and the distribution of elements in the sample. The samples were carefully snapped in the middle and mounted onto a sample holder to examine the cross section view of the ceramic DL-HFs. Images and EDX elemental mapping were taken at a magnification of 60× for the overall cross section shape, 500× to examine the layering of the electrolyte and anode and 1500× at the anode to examine the microstructure. The samples were then crushed and tested with X-ray Diffraction (XRD) using Rigaku Smart Lab to examine the crystallinity structure of the Ni-CeO 2 -YSZ material. Similarly, the structure of the reduced Mat_C DL-HF was scanned using XRD to examine the reducibility of Ni and the effect of reduction on the reducibility of CeO 2 .
Gas permeation tests were performed on the unreduced and reduced anode layer as performed by [13,14]. To perform this, the sample was mounted on 3/8" tube-thread male connectors, as shown in Figure A2 in Appendix A using epoxy resin. The sample was fitted to a permeation cell where the flow rate of N 2 gas flowing through the anode was measured using a bubble flow meter at varying inlet pressures.
Gas permeability was determined from Equation (1), where P is the value of gas permeability (mol m −2 s −1 Pa −1 ). Q, D i , D o , L and P are each gas permeation rate (mol s −1 ), inner diameter (m), outer diameter (m), length of hollow fiber (m) and pressure difference (P).
Mechanical strength of the ceramic DL-HF was measured via the Modulus of Rupture (MOR) method using mechanical testing equipment INSTRON 8032, as performed by [13,14]. This is to examine the effect of the sintering temperature and anode composition on the strength of the fiber. The sample was placed between two supports at a distance of Sustainability 2021, 13, 13789 5 of 26 5 cm each, as shown in Figure A3 in Appendix A. The load was applied by the machine, and the bending strength was calculated using Equation (2) below.
where N is the load value when the break occurs.

Conductivity Test
The electrical conductivity of the single-layer anode was measured using a 2-point direct-current (DC) technique using a potentiostat/galvanostat (Autolab PGSTAT 30, Eco-Chemie, The Netherlands). Potential between 0 and 1 V was applied across WE/S and RE/CE, and the resultant current was measured. The voltage-current plot was made, where the slope was taken as resistance, ρ a . The conductivities were calculated according to Equation (3).
where σ = conductivities, S cm −1 , L = length of fiber, cm, ρ a = resistance, A = area. A = π (D o 2 − D i 2 )/4. The plot of conductivities against percentage by volume, % volume Ni were plotted. The value of percentage by volume were obtained by dividing the mass with density for each component.

Performance Test
Finally, for the performance test, the ceramic DL-HF was coated with a cathode via a paint brush technique consisting of 3 layers where the first and second layers were 50% LSM and 50% YSZ, and the third layer was 100% LSM. The cathode material was prepared in suspension made of equal weight propylene glycol. The cell is sintered similarly as before but at 1250 • C sintering temperature.
To test for electrochemical performance, completed cells consisting of an anode, electrolyte and cathode were mounted onto the setup, as shown in Figure A4 in Appendix A. To perform this, silver wire was passed through the lumen or the anode side of the completed cells while another silver wire was wrapped around the cathode, as shown in Figure A4. The completed cells were mounted onto two alumina tubes at each end and sealed with a ceramic sealant (Aremco Ceramabond 552).
The anode collecting wire was connected to the WE/S electrode of PGSTAT as in the conductivity test, while the cathode collecting wire was connected to the RE/CE electrode. Fuel gas was passed through the lumen of the completed cells, while air was used as the oxidant on the cathode side. The cells were tested using H 2 and CH 4 gas at 10 mL min −1 at 700, 750, 800, 850 and 900 • C. At each fuel and temperature, maximum power density (MPD) was obtained by testing the cell for the current-voltage (I-V) curve. The completed cells were also tested for electrochemical impedance spectroscopy (EIS). From the EIS result, the ohmic resistance (R ohm ) was obtained by measuring the distance on the real axis from 0 to the left side of the first semicircle on the Nyquist plot, and the charge transfer resistance (R ct ) was measured as the diameter of the semicircle on the real axis. These values were subsequently plotted in the Arrhenius plot to investigate the activation energy (R a ) of the ohmic and charge transfer process.

Morphology and Properties
SEM imaging and EDX mapping of the cross section of the ceramic DL-HFs are shown in Figure 1. Based on Figure 1i,ii, the finger-like voids were present in all of the DL-HFs. The finger-like voids have been observed and reported in all of the studies, producing hollow fibers using a phase inversion/co-extrusion method and Ni-GDC and Ni-YSZ [8,15]. Finger-like formation occurs due to the phase inversion process between the water-soluble phase in the internal and external coagulant phase and the solvent-polymer phase in the suspension. The movement of the solvent NMP from suspension into the water caused the remaining ceramic polymer to solidify. The affinity of the NMP towards water and ceramicpolymer determines the structure of the solid ceramic-polymer formed. Since different ceramic components are known to have different affinities towards NMP, different ceramic materials will cause different formations, sometimes with undesirable effects. However, in this study, it was found that a finger-like formation with a sponge-like structure is wellformed in the NiO-CeO 2 -YSZ anode and in the electrolyte. This geometry was formed since both the anode and electrolyte layers were in contact with water during the phase inversion process.
Sustainability 2021, 13, x FOR PEER REVIEW finger-like voids have been observed and reported in all of the studies, producing holl bers using a phase inversion/co-extrusion method and Ni-GDC and Ni-YSZ [8,15]. Fing formation occurs due to the phase inversion process between the water-soluble phase internal and external coagulant phase and the solvent-polymer phase in the suspensio movement of the solvent NMP from suspension into the water caused the remaining ce polymer to solidify. The affinity of the NMP towards water and ceramic-polymer deter the structure of the solid ceramic-polymer formed. Since different ceramic componen known to have different affinities towards NMP, different ceramic materials will cause ent formations, sometimes with undesirable effects. However, in this study, it was foun a finger-like formation with a sponge-like structure is well-formed in the NiO-CeO2-Y ode and in the electrolyte. This geometry was formed since both the anode and elect layers were in contact with water during the phase inversion process. The finger-like formation created a large porose structure with a diameter of a 1 micrometer, while the sponge-like structure gave submicron porosity to the stru Based on Figure 1i, the average DL-HF obtained have diameters ranging from 1.23 t The finger-like formation created a large porose structure with a diameter of around 1 micrometer, while the sponge-like structure gave submicron porosity to the structure. Based on Figure 1i, the average DL-HF obtained have diameters ranging from 1.23 to 1.36 mm with lumen from 0.868 to 1.01 mm. The electrolyte thicknesses ranged from 23.4 to 27 µm, which was preferable since the reduced thickness ensures lower resistance towards the flow of oxygen ions [16].
High magnification images of 1500× were taken for all of the fibers, as shown in Figure 1iii, to examine the morphologies at the grain level. Large chunks of matter were found to present increasingly in size from Mat_A to Mat_E in the anode layer. EDX mapping, as shown in Figure 1iv, shows that the large chunks consist mainly of Y and Zr elements. Thus, it was found that the YSZ had sintered and agglomerated into larger chunks as the content of NiO increased (or CeO 2 decreased). Based on XRD results in the next section, YSZ was found to form a new phase with CeO 2 forming multiple ceriazirconia phases. In low NiO content or high CeO 2 content, the YSZ formed a uniform ceria-zirconia structure without large chunks of YSZ. However, in high NiO content or low CeO 2 content, the YSZ agglomerated together to form a larger YSZ structure.
An X-ray diffraction study on unreduced Mat_A, Mat_B, Mat_C, Mat_D and Mat_E is shown in Figure 2a. Based on Figure 2a, in Mat_A (50% YSZ, 50% CeO 2 anode), the XRD result shows a major peak of YSZ and CeO 2 , whilst the minor peaks of the phases of ceria-zirconia (CeZrO 4 ) are present as found using the database search. The major sharp peaks of YSZ were expected to come from the electrolyte and anode layers. Little of the CeO 2 content, however, was found to enter the lattice of the ZrO 2 in YSZ to create the ceria-zirconia compound. This event was observed when zirconia and ceria were present and sintered at high temperatures in another study [17]. In Mat_E, however, the NiO and YSZ mixture was not found to produce another phase, as shown in Figure 2a (50% YSZ, 50% NiO anode). It is interesting to note that two types of ceria-zirconia phases were detected in each Mat_A, Mat_B, Mat_C and Mat_D, which are CeZrO4 and Zr0.82Ce0.12O2. The CeZrO4 found in Mat_A and Mat_B contain higher molar ratios of Ce compared to Zr0.82Ce0.12O2 found in Mat_C and Mat_D. Mat_A to Mat_C were also found to contain minority peaks of Ce(III) in the form of Ce2O3, as found by the database search. Such event was observed in a few It is interesting to note that two types of ceria-zirconia phases were detected in each Mat_A, Mat_B, Mat_C and Mat_D, which are CeZrO 4 and Zr 0.82 Ce 0.12 O 2 . The CeZrO 4 found in Mat_A and Mat_B contain higher molar ratios of Ce compared to Zr 0.82 Ce 0.12 O 2 found in Mat_C and Mat_D. Mat_A to Mat_C were also found to contain minority peaks of Ce(III) in the form of Ce 2 O 3 , as found by the database search. Such event was observed in a few studies, where sintering of CeO 2 at short lengths of time tends to cause low oxygen content in the ceria compound, where the oxide from the CeO 2 escaped, leaving more oxygen-deficient Ce 2 O 3 behind [18]. The presence of Zr is also known to create the Ce(IV)/Ce(III) pair in the CeO 2 compound [19]. The database search also indicates that the Y 2 O 3 phase was present, showing the leaching of Y 2 O 3 out of the YSZ phase. Since the ceria-zirconia phase was also detected, this result might suggest that the CeO 2 has replaced the yttria in the YSZ material, which is also observed in other studies [20].
Two regions of interest are shown in Table 2. The region at 2θ of 47.5 • corresponds to phase (2 0 2) of CeO 2 , while 50.2 • corresponds to the (2 0 2) of ZrO 2 in the YSZ used. Intermediate peaks found between these two peaks are 48.2 • and 48.5 • , corresponding to the phase (4 4 0) and (1 1 2) of CeZrO 4 and Zr 0.82 Ce 0.12 O 2 each. Thus, when the amount of CeO 2 increases in a mixture of CeO 2 and YSZ, sintering will produce a Ce x Zr 1−x O 2 phase with a large x value, increasing the d-spacing in the process. A similar trend can be seen in 2θ of 56.3 • and 59.6 • . No insertion of NiO into the ZrO 2 or CeO 2 was found in any sample. This shows that, in the presence of the NiO-CeO 2 -ZrO 2 system, the interacting materials are only between CeO 2 and ZrO 2 , which made the system behave equally to CeO 2 -ZrO 2 .  Figure 2b. Only results of Mat_C are shown here as it contains similar amounts of NiO and CeO 2 , and it can be used to represent the result of other SOFC since it contains all the ceria-zirconia species present in Mat_A and Mat_E. Upon reduction with H 2 gas, Mat_C produced Ni metal together with decreasing peaks of CeO 2 and increased peaks of Ce 2 O 3 . The increased amount shows that not only the NiO reduction occurred but the reduction of Ce(IV) to Ce(III) also took place in the sample. This phenomenon is well-known to occur in the operation of CeO 2 in a highly reducing environment, where the usage of electrolytes is unfavorable to the n-type conductivity gained in the Ce(IV)-Ce(III) mixture [21]. However, the existence of the Ce(IV)-Ce(III) mixture was suggested by few studies to actually have higher catalytic activity towards oxidizing hydrocarbon compared to Ce(IV) alone [22,23]. The increased electronic conduction and different oxidation state might give more electron and oxide ion mobility to the mixture, causing it to perform a better oxidation reaction. Figure 3 shows the results for the mechanical testing. The bending strength is shown in Figure 3a for Mat_A to Mat_E from sintering temperatures of 1350 to 1550 • C. This temperature range was chosen based on the optimal sintering temperature of CeO 2 and YSZ [24,25]. The bending strength of Mat_A was shown to be the highest in all temperature ranges compared to the other DL-HF.  The bending strength of Mat_A reached the optimal point at sintering temperatures from 1450 to 1500 °C before declining, while the bending strength of Mat_B, Mat_C and Mat_D achieved optimal strength at a temperature of 1450 °C before significantly dropping at 1500 °C. Mat_E achieved optimal bending strength at 1400 °C, before declining slowly. Mat_E with 50% NiO-50% YSZ, nonetheless, has relatively similar bending strength from 1350 to 1550 °C, which is expected since YSZ is known to exhibit almost similar density over 1350 °C sintering temperature [25]. The rate of shrinkage of YSZ also increases as the sintering temperature increases, which might further explain this observation [25].
Mat_A (50% CeO2, 50% YSZ) was found to achieve optimal strength at a slightly higher temperature compared to Mat_E (50% NiO, 50% YSZ). From an XRD perspective, the diffusion of CeO2 into the ZrO2 lattice producing the ceria-zirconia compound was found to produce a ceramic with a higher strength compared to YSZ alone in Mat_A (no diffusion of NiO was observed).
Low mechanical strength in the mixture of the CeO2-Ni anode (Mat_B, Mat_C, Mat_D) was observed at all temperatures. At the low temperature of 1350 °C, mechanical strength decreases from Mat_A to Mat_C and increases back to Mat_E. For sintering temperatures from 1350 to 1550 °C, the lowest mechanical strengths are all in Mat_C. The low mechanical strength may possibly be attributed to the presence of different types of ceriazirconia compounds found in this membrane (and also in Mat_B and Mat_D) compared The bending strength of Mat_A reached the optimal point at sintering temperatures from 1450 to 1500 • C before declining, while the bending strength of Mat_B, Mat_C and Mat_D achieved optimal strength at a temperature of 1450 • C before significantly dropping at 1500 • C. Mat_E achieved optimal bending strength at 1400 • C, before declining slowly. Mat_E with 50% NiO-50% YSZ, nonetheless, has relatively similar bending strength from 1350 to 1550 • C, which is expected since YSZ is known to exhibit almost similar density over 1350 • C sintering temperature [25]. The rate of shrinkage of YSZ also increases as the sintering temperature increases, which might further explain this observation [25].
Mat_A (50% CeO 2 , 50% YSZ) was found to achieve optimal strength at a slightly higher temperature compared to Mat_E (50% NiO, 50% YSZ). From an XRD perspective, the diffusion of CeO 2 into the ZrO 2 lattice producing the ceria-zirconia compound was found to produce a ceramic with a higher strength compared to YSZ alone in Mat_A (no diffusion of NiO was observed).
Low mechanical strength in the mixture of the CeO 2 -Ni anode (Mat_B, Mat_C, Mat_D) was observed at all temperatures. At the low temperature of 1350 • C, mechanical strength decreases from Mat_A to Mat_C and increases back to Mat_E. For sintering temperatures from 1350 to 1550 • C, the lowest mechanical strengths are all in Mat_C. The low mechanical strength may possibly be attributed to the presence of different types of ceria-zirconia compounds found in this membrane (and also in Mat_B and Mat_D) compared to Mat_A. Since Mat_E did not have these Zr 0.82 Ce 0.12 O 2 compounds, it is possible that the compound actually reduces the mechanical strength of the ceramic DL-HF. The same trend was also found in Mat_D, as this sample contains these Zr 0.82 Ce 0.12 O 2 compounds as well. Higher mechanical strength in Mat_A compared to Mat_E, however, indicates that CeZrO 4 actually gives larger mechanical strength. The authors of [24] noted that the increase in gadoliniumdoped ceria (GDC) content in YSZ induces the shift in oxide peak plane in the XRD content, signifying possible dissolution of CeO 2 into the ZrO 2 [26]. The dissolution may create a denser CeO 2 -ZrO 2 structure as compared to YSZ [20]. A study on the GDC-YSZ system shows that as the YSZ component increases, even though the coefficient of thermal expansion (CTE) increases, the chemical expansion is actually higher in the GDC-YSZ mixture compared to GDC or YSZ alone [27]. NiO-YSZ was shown to have a CTE from 10.3 × 10 −6 (pure YSZ) to less than 14 × 10 −6 K −1 (60% NiO:40% YSZ) [28], which is different from the CTE from 13.2 × 10 −6 to 14.7 × 10 −6 K −1 in the ceria-zirconia system [29]. Low mechanical strength experienced in the mixture might be the result of the CTE mismatch between ceria-zirconia species and NiO-YSZ material.
Gas permeability across the anode is given in Figure 3b. Gas permeability is the lowest in Mat_A and increased from Mat_A to Mat_E. This property suggests that the higher the CeO 2 content, the lesser the gas permeation. As observed in Figure 1ii, the low permeability is also caused by the lesser presence of finger-like formations in the anode of ceramic DL-HF [30]. The lower gas permeability in higher CeO 2 content combined with the higher mechanical strength observed from Mat_C to Mat_A signifies that the presence of CeO 2 or ceria-zirconia does not only yield a harder anode but also a less porous one. This might be explained by the smaller finger-like projection and more sponge-like porous structure obtained from Mat_E to Mat_A, as observed in Figure 1ii. The lack of a large porous structure in Mat_A resulted in lower gas permeability. Furthermore, at sintering temperatures above 1350 • C, the rate of shrinkage in YSZ is larger than the rate of shrinkage in the ceria-zirconia structure [31].
Gas permeability was found to increase after reduction by H 2 gas. This is expected as when the NiO was reduced to Ni, volume shrinkage occurred, causing enlargement of pores and hence higher gas permeation. However, it is also to be noted that even though Mat_A contains no NiO, reduction elevates the gas permeability as well. This is caused by a reduction of Ce(IV) to Ce(III) or from CeO 2 to Ce 2 O 3 , as observed in XRD, which similarly leads to volume shrinkage and higher gas permeation. Similar trends of gas permeability are also observed in the gas tightness test, as shown in Figure 3c. We expect that the explanation of the behavior of the gas tightness check of the electrolyte is similar to the gas permeation in the anode. Since the thickness of the electrolyte was the same in all samples, the gas tightness of the electrolyte was found to be affected by the gas permeation rate through the anode in this study.
The anode conductivity obtained is shown in Figure 4. The conductivity of the Nicontaining anode shows percolation was achieved at above 30% Ni by volume, which is Mat_E (50% NiO by weight). Lowering the content of NiO to below 50% NiO will lead the conductivity to drop by a magnitude of six, as observed in Mat_D. This is similarly reported in early studies [32]. It is important to note that the conductivity obtained was lower than that reported in the literature, which is due to different methodologies. The conductivity, as reported in other studies, was measured by pressuring the sharp probe onto the anode surface and is only applicable since the geometry of the SOFC is planar. Here, the same methodology cannot be used due to the tubular geometry, and less pressure is able to be applied to the lumen of the tubular SOFC. Nonetheless, a large jump in conductivities is observed as in other studies, showing the good reliability of our method. Anode conductivity is crucial in the performance of SOFC, and low anode conductivity has been shown to correlate with lower performance. Low conductivity is expected in this study due to the low content of Ni. However, since Ni is also known to cause coking, a negative effect of lowering the Ni content was offset with the addition of CeO2, which was known to provide conductivity upon reduction. It is interesting to note that the conductivity of Mat_A is slightly higher than the conductivity of Mat_B. A similar trend was also observed in the GDC-YSZ system, where it was found that the conductivity of pure GDC and YSZ is higher than the mixture of GDC-YSZ [20]. Another explanation for the higher conductivity is possibly due to the denser packing structure in Mat_A compared to Mat_B, as shown in the SEM imaging in Figure 1. A denser structure allows for a better conductance path of electrons and hence higher conductivity [33].

Fuel Cell Performance
Fuel cell performance testing was performed at 700, 750, 800, 850 and 900 °C. The temperature of 700 °C was selected as the lowest temperature for testing since the YSZ electrolyte was found to have lower performance below this point [2]. High operating temperatures are expected to yield better performance for SOFC operating in hydrocarbon fuel since higher temperature ensures operation in the hydrocarbon reforming area [34]. Nonetheless, in this setup, we would like to test the effect of CeO2 as an oxidizing catalyst in SOFC without depending much on the reforming process; thus, 900 °C was selected as the highest operating temperature to be tested.

Power Density Analysis (I-V Diagram)
The current-voltage (I-V) diagram and power measurements are given in Figures A5  and A6, and Electrochemical Impedance Spectroscopy (EIS) is given in Figures A7 and A8 in Appendix A for each H2 and CH4 fuel.
For I-V measurements, generally, when run in H2 fuel, the MPD increased as the operating temperature increased, which is typical in YSZ-type cells for this operating temperature range. The MPD for SOFC running on H2 and CH4 is summarized in Figure 5. The MPD was found to unexpectedly drop from Mat_A to Mat_B and rise back again from Mat_B to Mat_E. The increment observed from Mat_B to Mat_E signifies that Ni presence  Anode conductivity is crucial in the performance of SOFC, and low anode conductivity has been shown to correlate with lower performance. Low conductivity is expected in this study due to the low content of Ni. However, since Ni is also known to cause coking, a negative effect of lowering the Ni content was offset with the addition of CeO 2 , which was known to provide conductivity upon reduction. It is interesting to note that the conductivity of Mat_A is slightly higher than the conductivity of Mat_B. A similar trend was also observed in the GDC-YSZ system, where it was found that the conductivity of pure GDC and YSZ is higher than the mixture of GDC-YSZ [20]. Another explanation for the higher conductivity is possibly due to the denser packing structure in Mat_A compared to Mat_B, as shown in the SEM imaging in Figure 1. A denser structure allows for a better conductance path of electrons and hence higher conductivity [33].

Fuel Cell Performance
Fuel cell performance testing was performed at 700, 750, 800, 850 and 900 • C. The temperature of 700 • C was selected as the lowest temperature for testing since the YSZ electrolyte was found to have lower performance below this point [2]. High operating temperatures are expected to yield better performance for SOFC operating in hydrocarbon fuel since higher temperature ensures operation in the hydrocarbon reforming area [34]. Nonetheless, in this setup, we would like to test the effect of CeO 2 as an oxidizing catalyst in SOFC without depending much on the reforming process; thus, 900 • C was selected as the highest operating temperature to be tested.

Power Density Analysis (I-V Diagram)
The current-voltage (I-V) diagram and power measurements are given in Figures A5 and A6, and Electrochemical Impedance Spectroscopy (EIS) is given in Figures A7 and A8 in Appendix A for each H 2 and CH 4 fuel.
For I-V measurements, generally, when run in H 2 fuel, the MPD increased as the operating temperature increased, which is typical in YSZ-type cells for this operating temperature range. The MPD for SOFC running on H 2 and CH 4 is summarized in Figure 5. The MPD was found to unexpectedly drop from Mat_A to Mat_B and rise back again from Mat_B to Mat_E. The increment observed from Mat_B to Mat_E signifies that Ni presence is crucial in SOFC operating in H 2 fuel, in which it acts as a catalyst and conductor. Nonetheless, Mat_A was found to have a higher MPD than Mat_B despite having no Ni content at all. The low catalytic activity of CeO 2 in oxidizing H 2 was well understood and explained by Density Functional Theory (DFT) modeling, which was shown to be inferior compared to Ni [35]. Therefore, this high MPD might be possibly explained by the presence of Ce(III) and possibly the ceria-zirconia compound as well, where both were shown to be able to function as an oxidation catalyst [36,37]. Apart from that, the presence of a large amount of CeO 2 in Mat_A compared to Mat_B might yield more Ce(IV)-Ce(III), increasing the anode conductivity and giving a slightly better performance [21]. This slightly higher conductivity can be observed in Figure 4.
is crucial in SOFC operating in H2 fuel, in which it acts as a catalyst and conductor. Nonetheless, Mat_A was found to have a higher MPD than Mat_B despite having no Ni content at all. The low catalytic activity of CeO2 in oxidizing H2 was well understood and explained by Density Functional Theory (DFT) modeling, which was shown to be inferior compared to Ni [35]. Therefore, this high MPD might be possibly explained by the presence of Ce(III) and possibly the ceria-zirconia compound as well, where both were shown to be able to function as an oxidation catalyst [36,37]. Apart from that, the presence of a large amount of CeO2 in Mat_A compared to Mat_B might yield more Ce(IV)-Ce(III), increasing the anode conductivity and giving a slightly better performance [21]. This slightly higher conductivity can be observed in Figure 4. It is interesting to note that even though Mat_D has conductivity below the threshold value, the MPD values of Mat_D (41 and 2.9 mW cm −2 at 900 °C in H2 and CH4 each) are relatively near to that of Mat_E (185 and 95 mW cm −2 at 900 °C in H2 and CH4 each) compared to that of Mat_C (6.7 and 2.0 mW cm −2 at 900 °C in H2 and CH4 each). Since Mat_D contains 12.5% less NiO content than Mat_E, Mat_D can be considered as a promising low-nickel-content anode to be used.

Electrochemical Impedance Spectroscopy (EIS)
When run on hydrogen fuel, generally the impedance for all of the samples decreases as the temperature increase. This is a typical characteristic of ohmic resistance due to resistance towards oxide ion movement through the electrolyte [20]. As ionic conductivity of YSZ increases as the temperature increases, the ohmic resistance drops as the operating temperature increases, as ionic conductivity is a temperature-activated process [38].
The impedance spectroscopy also shows characteristics of gas transfer limitations in SOFC with higher CeO2 content (Mat_A to Mat_C). This is observed with an increasing value of impedance at low frequency with the characteristic of 45 o inclination, described as the Warburg process [39]. However, in Mat_A to Mat_C, the peak at low frequency started to form a loop at higher temperatures, suggesting the reduction in gas transfer resistance. Another characteristic is that the impedance seems to consist of two loops that converge to form one loop at higher temperatures. As shown in Figure A9 in Appendix A, the first loop is often described as the polarization resistance, (Rct), associated with the accumulation of charge at the electrode due to the difficulty of charge transfer to reactants. It is also associated with the resistance towards the transfer of the electron from the cathode surface to the oxygen, or transfer of electron from the oxide ion to the anode, which It is interesting to note that even though Mat_D has conductivity below the threshold value, the MPD values of Mat_D (41 and 2.9 mW cm −2 at 900 • C in H 2 and CH 4 each) are relatively near to that of Mat_E (185 and 95 mW cm −2 at 900 • C in H 2 and CH 4 each) compared to that of Mat_C (6.7 and 2.0 mW cm −2 at 900 • C in H 2 and CH 4 each). Since Mat_D contains 12.5% less NiO content than Mat_E, Mat_D can be considered as a promising low-nickel-content anode to be used.

Electrochemical Impedance Spectroscopy (EIS)
When run on hydrogen fuel, generally the impedance for all of the samples decreases as the temperature increase. This is a typical characteristic of ohmic resistance due to resistance towards oxide ion movement through the electrolyte [20]. As ionic conductivity of YSZ increases as the temperature increases, the ohmic resistance drops as the operating temperature increases, as ionic conductivity is a temperature-activated process [38].
The impedance spectroscopy also shows characteristics of gas transfer limitations in SOFC with higher CeO 2 content (Mat_A to Mat_C). This is observed with an increasing value of impedance at low frequency with the characteristic of 45 o inclination, described as the Warburg process [39]. However, in Mat_A to Mat_C, the peak at low frequency started to form a loop at higher temperatures, suggesting the reduction in gas transfer resistance. Another characteristic is that the impedance seems to consist of two loops that converge to form one loop at higher temperatures. As shown in Figure A9 in Appendix A, the first loop is often described as the polarization resistance, (R ct ), associated with the accumulation of charge at the electrode due to the difficulty of charge transfer to reactants. It is also associated with the resistance towards the transfer of the electron from the cathode surface to the oxygen, or transfer of electron from the oxide ion to the anode, which occurs due to the catalytic limitation of the electrode. The second loop is often called concentration polarization, R conc resistance, which occurs due to mass transport limitation in the electrode associated with morphology of the electrode [39].
The value of R ohm for each EIS was measured and plotted according to temperatures, and subsequently, an Arrhenius plot was prepared from the data to obtain the activation energy, E a , for the ohmic process, as shown in Figure 6a,b,e,f. Similar data processing was applied to R ct and is shown in Figure 6c,d,g,h. The activation energy, E a , obtained for the ohmic process and charge transfer process is shown in Figure 7.
occurs due to the catalytic limitation of the electrode. The second loop is often called concentration polarization, Rconc resistance, which occurs due to mass transport limitation in the electrode associated with morphology of the electrode [39].
The value of Rohm for each EIS was measured and plotted according to temperatures, and subsequently, an Arrhenius plot was prepared from the data to obtain the activation energy, Ea, for the ohmic process, as shown in Figure 6a,b,e,f. Similar data processing was applied to Rct and is shown in Figure 6c,d,g,h. The activation energy, Ea, obtained for the ohmic process and charge transfer process is shown in Figure 7.  Ohmic resistance for H 2 (a) and CH 4 (c) and charge transfer resistance for H 2 (e) and CH 4 (g). Arrhenius plot for ohmic resistance for H 2 (b) and CH 4 (d), and Arrhenius plot for charge transfer resistance for H 2 (f) and CH 4 (h).

R ohm for Operation in H 2 and CH 4
The R ohm value was found to decrease as the Ni content increased from Mat_A to Mat_E in both H 2 and CH 4 fuel. The ohmic resistance is understood to be the process involving the movement of O 2− ions across the cathode, electrolyte and anode, where a lower R ohm is favored for a higher MPD of the cell [20]. This is reflected in the reading of MPD shown in Figure 5. For all SOFC, R ohm decreases as temperature increases, showing typical temperature-activated processes of the O 2− movement. The trend in R ohm follows exactly the trend seen in conductivity, where higher conductivity was found to yield lower R ohm , despite having similar content of YSZ in the cathode, electrolyte and anode. This shows that anode conductivity plays a crucial role in R ohm . Anode conductivity may have made it easier for the electron to move through the anode. The oxidation of O 2at the anode might have been impeded if the anode conductivity was low simply because the electron cannot leave O 2− or the oxidation reaction because there is no conduction path, leaving the O 2transport impeded as well [40].
It is interesting to note that even though Mat_D has conductivity below a threshold value, the R ohm values of Mat_D (6.5 and 10.0 Ω cm 2 at 900 • C in H 2 and CH 4 each) are relatively near to that of Mat_E (4.0 and 2.5 Ω cm 2 at 900 • C in H 2 and CH 4 each) compared to that of Mat_C (60.1 and 70.2 Ω cm 2 at 900 • C in H 2 and CH 4 each). Since Mat_D contains 12.5% less NiO content than Mat_E, Mat_D can be considered as a promising low-nickel-content anode to be used.

R ct for Operation in H 2 and CH 4
R ct involves the transfer of charge where oxygen is reduced at the cathode or the fuel is oxidized at the anode [31]. However, since the cathode thickness is the same for all of the SOFCs, the value of R ct here reflects the oxidation of fuel at the anode.
R ct was found to be lower as Ni content increased in H 2 . This is expected since the oxidation of H 2 is catalyzed by the Ni metal, hence increasing the Ni content resulted in lower R ct values. The trend is actually similar to the trend of anode conductivity as well. As explained in the R ohm section, there might be another factor affecting the R ct when operating in H 2 , which is anode conductivity. Here, the high conductivity and high catalytic rate both led to the reduction of R ct as the Ni content increased. R ct in CH 4 , on the other hand, increased from Mat_A to Mat_C, then decreased rapidly to Mat_D and Mat_E. The charge transfer process in CH 4 oxidation is catalyzed both by Ni and at lower rates by CeO 2 [35]. However, as explained in the R ohm section, R ct may actually depend on the value of R ohm and hence anode conductivity as well, since lower anode conductivity may make oxidation harder to proceed since the electron cannot leave the oxidation reaction [41]. With the lack of Ni (Mat_A and Mat_B), oxidation is catalyzed by CeO 2 , but since the anode conductivity is low, the rate is lower, resulting in higher R ct . In the presence of Ni (Mat_D and Mat_E), it can be achieved with Ni at a higher rate and, coupled with high conductivity, push the R ct value even lower. We believe that the Mat_C mixture provided higher R ct , possibly due to lesser catalytic activity due to less CeO 2 and Ni, coupled with a lower-than-threshold value of anode conductivity.
It needs to be noted that the anode composition of Mat_D has almost similar Rct (4.5 and 1.9 Ω cm 2 at 900 • C in H 2 and CH 4 each) to Mat_E (1.5 and 0.2 Ω cm 2 at 900 • C in H 2 and CH 4 each) compared to Mat_C (7.2 and 34.9 Ω cm 2 at 900 • C in H 2 and CH 4 each), showing 37.5% Ni, 12.5% CeO 2 and 50% YSZ to be the suitable low-Ni anode for use.

E a for R ohm for Operation in H 2 and CH 4
The value of E a for R ohm shows the energy barrier for the ohmic process or the resistance towards the flow of oxide ions in the electrolyte [14,20]. Greater changes in R ohm that result from changes in temperature yield higher E a , signifying that higher temperature will supply more energy for the reactant to overcome the energy barrier. In the case of R ohm , higher E a indicates a higher energy barrier for ionic transportation of mainly the O 2− ion to move through the electrolyte. Higher E a is undesirable since the process usually requires a high operating temperature. O 2− ion to move through the electrolyte. Higher Ea is undesirable since the process u requires a high operating temperature. The value of Ea for Rohm and Rct for each H2 and CH4 is shown in Figure 7. Brief Ea for Rohm for H2 was found to be 16.5, 23.3, 9.7, 10.3 and 8.1 kJ mol −1 for Mat_A, M Mat_C, Mat_D and Mat_E. Ea for Rohm for operation in H2 generally follows the trendline as in conductivity, where Mat_A is lower than Mat_B but becomes lower composition approaches Mat_E. A possible explanation for this observation is that ment of the oxide ion through the electrolyte is a multistep process, which partly co of the process of electron transport movement, as illustrated in Figure 8a,b. For exa for the O 2− ion to move, O 2− production via reduction at the cathode and O 2− oxida the anode involves the transfer of an electron from and to the oxygen, and thus, s the flow of electron be resisted at the anode (or cathode), the flow of oxide ions w be affected across the electrolyte [20].
The Ea for H2 Rohm is also generally lower than for CH4 (except for Mat_B), sh the effect of fuel on the resistance to O 2− movement. This is known since the prese H2 in a CeO2-containing anode is known to reduce the CeO2 to Ce2O3, which off alternate O 2-movement pathway other than through the YSZ, which has lower E increases the O 2− mobility and hence the flow of the O 2− ion [40]. However, in the ope in H2, the anode conductivity seems to be the limiting factor in the Rohm, while in CH anode conductivity is not the limiting factor [31].
The Ea for Rohm for CH4 was found to be 24.4, 20.6, 14.5, 25.1 and 21.9 kJ mo Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. Compared to H2, Ea for Rohm for opera CH4, on the other hand, drops from Mat_A to Mat_B, increased to Mat_D and dec slightly to Mat_E. The possible explanation for this phenomenon is similar to tha for H2, where the movement of oxide ion is a multistep process consisting of an el conduction process and O 2− ion generation process, etc., as shown in Figure 8c.
In a sample with high Ni content (Mat_D to Mat_E), the drop in Ea is actua result of the increase in the anode conductivity, giving lower Ea for electron condu and hence Rohm. The Ea for Mat_A to Mat_C also has a downward trend due to a s effect. However, as the CeO2 content decreases from Mat_A to Mat_E, the conducta O 2− also decreases, giving a higher Ea. This can be observed in higher Ea in Mat_ Mat_E compared to Mat_B and Mat_C. From Mat_A to Mat_C with low Ni conte with high CeO2 content, the low Ea for Rohm is actually caused by the presence o amounts of CeO2, Ce(IV)/Ce(III) pairs and ceria-zirconia compounds, which are k to be able to conduct O 2− [40]. The value of E a for R ohm and R ct for each H 2 and CH 4 is shown in Figure 7. Briefly, the E a for R ohm for H 2 was found to be 16.5, 23.3, 9.7, 10.3 and 8.1 kJ mol −1 for Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. E a for R ohm for operation in H 2 generally follows the same trendline as in conductivity, where Mat_A is lower than Mat_B but becomes lower as the composition approaches Mat_E. A possible explanation for this observation is that movement of the oxide ion through the electrolyte is a multistep process, which partly consists of the process of electron transport movement, as illustrated in Figure 8a,b. For example, for the O 2− ion to move, O 2− production via reduction at the cathode and O 2− oxidation at the anode involves the transfer of an electron from and to the oxygen, and thus, should the flow of electron be resisted at the anode (or cathode), the flow of oxide ions will also be affected across the electrolyte [20].
The E a for H 2 R ohm is also generally lower than for CH 4 (except for Mat_B), showing the effect of fuel on the resistance to O 2− movement. This is known since the presence of H 2 in a CeO 2 -containing anode is known to reduce the CeO 2 to Ce 2 O 3 , which offers an alternate O 2movement pathway other than through the YSZ, which has lower E a . This increases the O 2− mobility and hence the flow of the O 2− ion [40]. However, in the operation in H 2 , the anode conductivity seems to be the limiting factor in the R ohm , while in CH 4 , the anode conductivity is not the limiting factor [31].
The E a for R ohm for CH 4 was found to be 24.4, 20.6, 14.5, 25.1 and 21.9 kJ mol −1 for Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. Compared to H 2 , E a for R ohm for operation in CH 4 , on the other hand, drops from Mat_A to Mat_B, increased to Mat_D and decreased slightly to Mat_E. The possible explanation for this phenomenon is similar to that of R ct for H 2 , where the movement of oxide ion is a multistep process consisting of an electron conduction process and O 2− ion generation process, etc., as shown in Figure 8c.
In a sample with high Ni content (Mat_D to Mat_E), the drop in E a is actually the result of the increase in the anode conductivity, giving lower E a for electron conductance and hence Rohm. The Ea for Mat_A to Mat_C also has a downward trend due to a similar effect. However, as the CeO 2 content decreases from Mat_A to Mat_E, the conductance of O 2− also decreases, giving a higher E a . This can be observed in higher Ea in Mat_D and Mat_E compared to Mat_B and Mat_C. From Mat_A to Mat_C with low Ni content but with high CeO 2 content, the low E a for R ohm is actually caused by the presence of high amounts of CeO 2 , Ce(IV)/Ce(III) pairs and ceria-zirconia compounds, which are known to be able to conduct O 2− [40]. Sustainability 2021, 13, x FOR PEER REVIEW 17 of 29 Figure 8. Illustration of processes involved in SOFC and related Ea.

Ea for Rct for Operation in H2
The Ea for Rct for H2 was found to be 29.2, 36.6, 27.9, 27.1 and 22.7 kJ mol −1 for each Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. The trend of Ea for the charge transfer process in H2 follows a similar pattern as Rohm and anode conductivity. For H2, Ea increases from Mat_A to Mat_B and decreases back to Mat_E, achieving a maximum at Mat_B. While the pattern suggests anode conductivity as the rate-limiting step, we believe that the charge transfer process of H2 is also a multiple-step process, as shown in Figure 8b, but the other process has lower Ea than the electron conduction process. This is because the oxidation of H2 may follow a different route depending on the content of ceria-zirconia, the Ce(IV)/Ce(III) pair and Ni [19]. Lowering the concentration of both CeO2 or NiO (as in Mat_B to Mat_D) results in higher Ea. Both Ni and CeO2 are known to be able to adsorb

E a for R ct for Operation in H 2
The E a for R ct for H 2 was found to be 29.2, 36.6, 27.9, 27.1 and 22.7 kJ mol −1 for each Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. The trend of E a for the charge transfer process in H 2 follows a similar pattern as R ohm and anode conductivity. For H 2 , E a increases from Mat_A to Mat_B and decreases back to Mat_E, achieving a maximum at Mat_B. While the pattern suggests anode conductivity as the rate-limiting step, we believe that the charge transfer process of H 2 is also a multiple-step process, as shown in Figure 8b, but the other process has lower E a than the electron conduction process. This is because the oxidation of H 2 may follow a different route depending on the content of ceria-zirconia, the Ce(IV)/Ce(III) pair and Ni [19]. Lowering the concentration of both CeO 2 or NiO (as in Mat_B to Mat_D) results in higher Ea. Both Ni and CeO 2 are known to be able to adsorb H 2 gas onto their surface, where the rate of adsorption of H 2 onto Ni is greater than the rate of adsorption on CeO 2 [42]. However, unlike Ni that requires ZrO 2 (in YSZ) to oxidize the adsorbed hydrogen, CeO 2 can oxidize the H 2 by directly turning it into Ce 2 O 3 . Thus, from Mat_A to Mat_B where the concentration of Ni is still too low to catalyze H 2 oxidation, the oxidation is performed by CeO 2 , and the increment of E a observed from Mat_A to Mat_B is the result of the lowering of CeO 2 content. However, in Mat_C to Mat_E, the presence of higher content of Ni led to a rapid reduction of E a for R ct of H 2 . This effect is obvious in the charge transfer process of CH 4 .

E a for R ct for Operation in CH 4
The E a for R ct for CH 4 was found to be 40.8, 35.9, 29.7, 39.0 and 64.0 kJ mol −1 for Mat_A, Mat_B, Mat_C, Mat_D and Mat_E. For oxidation of CH 4 , the trend is the reverse of the trend seen in H 2, where the E a decreased from Mat_A to Mat_C and increased back to Mat_E. This trend is also almost similarly observed by [17], which found the presence of Ni on ceria-zirconia increased oxidation of methane [19]. The performance is high when small amounts of Zr are present in the Ni-ceria-zirconia compound, where the rate of methane oxidation was found to be Ni-ZrO 2 < Ni-Ce 0. 17  A summary of the effect of Ni and CeO 2 content on electronic conductivity, O 2− conductivity, H 2 oxidation and CH 4 oxidation was made based on observations from the I-V curve, MPD, EIS and E a and is shown in Table 3. Also based on these observations, a summary of E a for R ohm and R ct when operating in H 2 and CH 4 fuel is shown in Table 4.  Based on the anode conductivity reading, lowering the Ni content to lower than 50% NiO (Mat_E) in the anode was shown to produce a detrimental effect on the conductivity, where the anode conductivity drastically falls more than 10 7 times. However, based on the MPD value, large losses in the value were actually observed when the Ni content lowered to below 37.5% NiO (Mat_D), which sees the abrupt losses in MPD. Similarly, based on EIS data in H 2 and CH 4 fuel, a large value of impedance is observed in lower than 37.5% NiO. However, lowering of Ni coupled with an increasing amount of CeO 2 yield benefits in terms of reduction of E a for R ct for operation in CH 4 , especially when the content of NiO-CeO 2 -YSZ is balanced at 25%:25%:50% (Mat_C). Having relatively similar E a for R ohm for operation under H 2 and CH 4 with the lowest E a for R ct for CH 4 , this anode content may be concluded to be the best to function as an anode for operation in CH 4 . The property exhibited by Mat_C and Mat_D shows great potential in becoming a coking-tolerant anode for SOFC. Further study will be conducted to investigate the long-term performance of this anode in CH 4 fuel.
To date, there is no known study that fabricates CeO 2 directly into the Ni-YSZ anode of MT-SOFC via phase-inversion techniques as exhibited in this study, and hence no direct comparison of performance can be made. However, this study shows that it is possible to use low Ni content in an SOFC anode to limit the coking of an anode when used with hydrocarbon fuel in MT-SOFC. Low content Ni of 37.5% NiO with 12.5% CeO 2 (Mat_D) is shown here to obtain comparable performance with an anode containing 50% NiO (Mat_E) when compared in terms of MPD and impedance. MPD obtained by Mat_D is relatively low compared to traditional SOFC with 60% NiO and 40% YSZ but almost similar to another study that aims to reduce carbon deposition by using CeO 2 or ceria-zirconia anodes, as shown in Table 5. It is important to note that the value of MPD may be low due to different fuel mixes used since this study uses pure CH 4 while other studies used mixtures of CH 4 , H 2 and CO 2 . The later fuel mixture indicates the reforming process in other studies. Furthermore, other studies use mainly planar geometry, which involves more preparation steps and more costs compared to the micro-tubular geometry in this study. The low MPD value is also caused by the thin electrolyte (~20 µm) used in this study that led to low gas tightness [14]. Low gas tightness may lead to leakage of oxygen from the cathode into the anode and fuel into the cathode, leading to the oxidation of fuel at the cathode and the reduction in oxygen in the anode, reducing the power output.
The low MPD problem in this SOFC can be mitigated by fabricating a thicker electrolyte to enhance the gas tightness. It can also be mitigated with a parallel stacking of MT-SOFC to increase current production or stacking in series to increase the voltage of the power sources if the MT-SOFC is to be commercialized [13]. Given the ease of fabrication and the smaller size of MT-SOFC compared to planar SOFC, the study provides great insight into fabricating low-nickel coking-tolerant SOFC. Furthermore, this type of SOFC can be used with renewable energy sources that contain hydrocarbon, such as biogas, or to extract energy from methane-based energy storage technology to ensure the most environmentally friendly application possible [46,47].

Conclusions
From this study, we found that lowering the Ni content to lower than 50% NiO (Mat_E) in the anode caused conductivity to drop by a magnitude of six, below the conductivity percolation threshold. However, the effect of anode conductivity on the MPD and EIS reading is only prominent when it reaches below 37.5% NiO. MPD for low Ni content, which are Mat_A, Mat_B and Mat_C (0% to 25% NiO), is generally low, while MPD for Mat_D and Mat_E (37.5% to 50% NiO) is high. In terms of E a for R ct (which shows the catalytic properties of anodes), Ea are low for Mat_A to Mat_D compared to Mat_E. This resulted in Mat_D (37.5% NiO, 12.5% CeO 2 ) showing the lowest amount of Ni without abrupt negative effect on the MPD and impedance. Nonetheless, Mat_C was found to have the lowest E a for R ct for operation in CH 4 , making it the best content to oxidize CH 4, but at a low MPD value.
By using Ea for Rohm and Rct as a method for investigation, we found a few differences between the oxidation of H 2 and CH 4 in the Ni-CeO 2 -YSZ anode. They are:

1.
In H 2 operation, E a for R ohm and R ct is limited by anode conductivity, where higher anode conductivity will provide lower E a . This implies that anode conductivity is a rate-limiting step in H 2 operation. 2.
In CH 4 operation, E a for R ohm is limited by anode conductivity, oxide ion conductivity and fuel oxidation rate. 3.
In CH 4 operation, lower E a for R ct or CH 4 oxidation is observed in the Ni-CeO 2 -YSZ anode. The presence of Ni in ceria-zirconia increases the oxidation of methane.
Thus, based on this study, anode compositions containing 25% NiO:25% CeO 2 :50% YSZ (Mat_C) and 37.5% NiO:12.5% CeO 2 :50% YSZ (Mat_D) show promising results in becoming the low-Ni anode for coking-tolerant SOFC operation under CH 4 fuel.    . Figure A3. Schematic of the triple point bending test.   . Figure A3. Schematic of the triple point bending test. Figure A3. Schematic of the triple point bending test.    Figure A6. I-V and power measurement for Mat_A to Mat_E from temperatures 700 to 900 °C running in CH4 at 10 mL min −1 . Figure A6. I-V and power measurement for Mat_A to Mat_E from temperatures 700 to 900 • C running in CH 4 at 10 mL min −1 . Sustainability 2021, 13, x FOR PEER REVIEW 25 of 29 Figure A7. EIS for Mat_A to Mat_E from temperatures 700 to 900 °C running in H2 gas at 10 mL min −1 . Figure A7. EIS for Mat_A to Mat_E from temperatures 700 to 900 • C running in H 2 gas at 10 mL min −1 . Sustainability 2021, 13, x FOR PEER REVIEW 26 of 29 Figure A8. EIS for Mat_A to Mat_E from temperatures 700 to 900 °C running in CH4 gas at 10 mL min −1 . Figure A8. EIS for Mat_A to Mat_E from temperatures 700 to 900 • C running in CH 4 gas at 10 mL min −1 . Figure A9. Equivalent circuit and measurement of R ohm and R ct .