Experimental and Process Modelling Investigation of the Hydrogen Generation from Formic Acid Decomposition Using a Pd/Zn Catalyst

: The use of hydrogen as a renewable fuel has attracted great attention in recent years. The decomposition of formic acid under mild conditions was investigated using a 2%Pd 6 Zn 4 catalyst in a batch reactor. The results showed that the conversion of formic acid increases with reaction temperature and with the formic acid concentration. A process-simulation model was developed to predict the decomposition of formic acid using 2%Pd 6 Zn 4 in a batch reactor. The model demonstrated very good validation with the experimental work. Further comparisons between the 2%Pd 6 Zn 4 catalyst and a commercial Pd/C catalyst were carried out. It was found that the 2%Pd 6 Zn 4 demonstrated signiﬁcantly higher conversions when compared with the commercial catalyst. catalyst. The results show that, during continuous-ﬂow operation, pore contamination and poisoning caused by formate ions lead to catalyst deactivation. Although these factors cause extensive deactivation in the plug-ﬂow mode, by operating in a continuous stirred tank reactor to reduce the steady-state formic acid concentration, encouraging results can be obtained. Therefore, under mild conditions and without stoichiometric additives, the system can operate continuously for more than 2500 cycles without loss of activity.


Introduction
The use of conventional fossil fuels for energy production has led to serious global consequences for the environment and climate. Renewable fuels have gained increasing attention in recent years to combat the current climate issues. Formic acid (HCOOH) has proven immensely popular as a favourable hydrogen storage/generation material due to its high gravimetric and volumetric hydrogen capacity, ease of handling, non-toxicity, stability at room temperature, and abundant supply from the conversion of biomass and carbon dioxide (CO 2 ) [1]. The decomposition of formic acid reaction can be conducted in a liquid phase (at temperatures lower than 373 K) as well as in a gas phase (at excess temperatures), resulting in the dehydrogenation or dehydration pathways [2]. The former produces CO 2 and hydrogen (H 2 ) as products, and the latter generates carbon monoxide (CO) and water (H 2 O) products. Performing the reaction at high temperatures (greater than 423 K) may lead to lower selectivity in the hydrogen production, because of the reverse water-gas shift reaction. Therefore, selective catalysts for formic acid dehydrogenation under these conditions are in high demand [3].
Formic acid decomposition for H 2 production has recently been investigated using Pd/ZnO catalysts, which are typically used for methanol steam reforming. Bulushev et al. [3] investigated the catalytic properties of PdZn/ZnO in formic acid decomposition. In addition, these catalysts were compared with Pt/ZnO and Pd/Al 2 O 3 catalysts as well as ZnO support. The results show that the measured catalyst activity corresponds to the following sequence: Pd/Al 2 O 3 ≥ Pd/ZnO > Pt/ZnO > ZnO. The Pd/ZnO catalyst shows the highest selectivity to hydrogen (up to 99.3%). This is attributed to the PdZn alloy formation during the catalyst reduction pretreatment. Increasing the pretreatment temperature from 573 K to 773 K resulted in a significant increase in the mean size of the PdZn (PtZn) nanoparticles. However, the activity of the catalyst remains unchanged, though the hydrogen selectivity increases.
Yao et al. [35] investigated Pd-supported catalysts for enhancing the decomposition of formic acid. The porous carbon support was synthesised using zeolitic imidazole frameworks (ZIFs) as a precursor. Furthermore, the porous carbon was produced by the one-step method. By doping the Zn into the predetermined bimetallic ZnCo-ZIFs, the operation of the catalyst was enhanced. The greatest catalytic activity of the prepared catalyst (Pd/Co@CN-2), with a 2.6 nm average diameter of PdCo, was observed when the Zn/Co molar ratio was 2 and demonstrated an initial turnover frequency value (TOF) as high as 2302 h −1 at 30 • C.
Ding et al. [36] investigated the synthesis of a hybrid material of Pd nanoparticles encapsulated within carbon nanotubes (CNTs) (Pd-CNTs-in). The results show that the prepared Pd-CNTs-in catalyst exhibits particularly high formic acid decomposition activity and durability at room temperature. The TOF value was found to be 1135 h −1 in the first 10 min, and there was no significant decrease in the continuous three-fold recovery study. The conclusion reached is that this work can provide a promising strategy for manufacturing cost-effective and highly active Pd-based catalysts for the dehydrogenation of formic acid.
The investigation of formic acid decomposition in batch reactors is often limited by their ability to generate a continuous stream of H 2 for fuel cell applications. Moreover, these reactors do not permit a detailed study of the catalyst lifetime to give an indication of its long-term performance. Therefore, the study of formic acid decomposition in continuous-flow reactors is essential to the comprehensive understanding of this reaction. The decomposition of formic acid has been explored using homogeneous and heterogeneous catalysts in continuous stirred tank reactors (CSTR) [37][38][39].
Caiti et al. [40] investigated the production of hydrogen via the low-temperature (less than 110 • C), additive-free decomposition of formic acid over a heterogeneous Pd/C catalyst. The results show that, during continuous-flow operation, pore contamination and poisoning caused by formate ions lead to catalyst deactivation. Although these factors cause extensive deactivation in the plug-flow mode, by operating in a continuous stirred tank reactor to reduce the steady-state formic acid concentration, encouraging results can be obtained. Therefore, under mild conditions and without stoichiometric additives, the system can operate continuously for more than 2500 cycles without loss of activity.
The decomposition of formic acid has often been investigated on an experimental basis. The utilisation of software for conducting detailed numerical studies can facilitate an understanding of parameter optimisation for the decomposition of formic acid reaction. Furthermore, theoretical studies are a valuable tool to assess the validity of the experimental findings [41,42]. Hafeez and Sanchez et al. [43] investigated the liquid-phase decomposition of formic acid using a Pd/C catalyst in a packed bed microreactor. A comprehensive heterogeneous computational fluid dynamic (CFD) model was developed to validate the experimental data. The results showed that the CFD model demonstrated very good validation with the experimental work. Furthermore, the novel model could successfully depict the deactivation of the Pd/C catalyst based on the poisoning species CO. The model developed can be used to effectively predict the decomposition of formic acid in microreactors for potential fuel cell applications.
The process-simulation modelling of formic acid decomposition has not been implemented for the current reaction system, which enhances the uniqueness of the current work. Process-simulation or flowsheet modelling can offer a solid foundation based on mass and energy balances, mass and heat transfer, reaction-system modelling, and phase equilibria to simulate the decomposition of formic acid. We utilised the capabilities of Aspen Plus to develop a simple innovative reaction system by integrating the phase-equilibrium exercises. It has been observed that both CFD and process-simulation modelling have the ability to predict different reactions with similar results [44]. The production of H 2 from the decomposition of formic acid has been investigated using a 2%Pd 6 Zn 4 catalyst in a batch reactor. A comprehensive process-simulation model has been developed based on the batch stirred reactor to validate the experimental findings. Further studies are performed to assess the characterisation and activity of the catalyst.

Catalyst Preparation
The 2%Pd 6 Zn 4 @HHT catalyst was prepared using the galvanic replacement procedure employed by Liao et al. [43], modified for the Pd-Zn system. This typical method consists of a two-step synthesis. In the first part, Zn(0) was obtained using a reducing agent (NaBH 4 ) from ZnCl 2 salt. The second step involves the spontaneous displacement reduction between (PdCl 4 ) 2and Zn. First, the support (HHT-CNFs) was added in milli-Q water under sonication conditions for 40 min. Subsequently, the mixture was degassed in a vacuum and purged in N 2 for 1 h. After that, 20 mg of NaBH 4 and 3.8 mg of Zn as ZnCl 2 were added. The mixture was then shaken by sonication for 5 min and placed overnight in a vacuum under vigorous stirring. In the second part, 6.2 mg of Pd as K 2 PdCl 4 was added. The mixture was then vigorously stirred in vacuum for 20 min. After that, the catalyst was filtered, washed with deionised water, and dried in the oven for 48 h.

Catalytic Tests
The liquid-phase formic acid dehydrogenation was performed using a two-necked round-bottom flask furnished with a magnetic stirrer and reflux condenser and placed in a water bath. Usually, 10 mL of formic acid (FA) (aqueous) was loaded into the reactor. Different concentrations of the substrate solution were utilised: 0. to establish the stirring rate were previously shown by our group [4]. The catalyst was also tested at various temperatures in order to study the kinetic behaviour of the catalyst, 30 • C, 40 • C, 50 • C, and 60 • C. Stability tests were performed using the procedure previously discussed, reusing the catalyst without any further purification procedure.
High-performance liquid chromatography (HPLC) was implemented to assess the conversion of formic acid. A H+ chromatographic column (Alltech OA-10,308, 300 mm 7.8 mm) was utilised with a UV detector established at 210 nm. The liquid samples of 200 µL were obtained after certain time intervals and subsequently diluted to 5 mL with a solvent of H 2 SO 4 solution (0.005 N). This isocratic effluent flow was maintained at a flow rate of 0.5 mL/min. The volume of gases, CO and CO 2 , was monitored and recorded using an online micro-gas chromatograph (Agilent 3000A) every 2 h. This equipment consisted of a molecular sieve module, an OV-1 module (stationary phase of polydimethylsiloxilane), and a TCD detector. The measurement of CO and CO 2 was achieved from calibration curves taken from commercial standards.

XPS
The XPS measurements were recorded using a thermo Scientific K-alpha+ spectrometer. The sample analysis was achieved using a monochromatic Al X-ray source running at 72 W, with the signal recorded over an oval-shape area of 600 × 400 µ. The data were recorded at 150 eV for survey scans and 40 eV for high resolution (HR) scans with 1 eV and 0.1 eV step sizes, respectively. The analysis of the data was performed using CASAXPS (v2.3.17 PR1.1), using Scofield sensitivity factors and an energy exponent of −0.6.

TEM
The transmission electron microscopy (TEM) experiments were completed using a double Cs aberration-corrected FEI Titan3 Themis 60-300 microscope. The microscope was equipped with a monochromator, an X-FEG gun, and a high-efficiency XEDS ChemiSTEM. This comprised four windowless SDD detectors. HR-STEM imaging was completed at 200 kV, applying a high-angle annular dark-field (HAADF) detector with a camera length of 11.5 cm. The HAADF-STEM procedure is responsive to the atomic number of the elements, in which the intensity is roughly proportional to the square of the atomic number (Z2), and it provides the opportunity to differentiate the small nanoparticles carried on light supports. The XEDS mappings were achieved using a beam current of 200 pA and a dwell time per pixel of 128 µs. The elemental maps were filtered using a Gaussian blur of 0.8 (Velox software) to enhance the quality of the images. Due to the STEM-HAADF images of the catalysts, the diameters of more than 200 randomly chosen metal particles were measured and the resultant metal particle size distributions (PSD) were defined. As a result, the average particle diameter (d) was computed corresponding to the following expression: d = ∑n i d i /∑n i , where n i ≥ 200. Equally, the total metal dispersion was evaluated according to D = Ns/Nt, where Ns is the total number of surface-metal atoms and Nt is total number of atoms in the metal particle. For the particle size calculations, the ImageJ software was used.

Modelling Methodology
The decomposition of formic acid reaction was modelled using the software Aspen Plus V11. The model selected to accurately predict the behaviour of the vapour-liquid equilibrium of the reaction system was the non-random two-liquid (NRTL) model [44]. The process-modelling studies were performed in a batch reactor with a 100 mL volume and a temperature range of 30-60 • C. Characteristically, 10 mL of the required concentration of formic acid was deposited into the batch reactor. The catalyst used was 2%Pd 6 Zn 4 , the formic acid/metal molar ratio was 2000:1, and the stirring rate was 1400 rpm.
The most significant feature of the NRTL method [45] is its applicability to multicomponent systems, presuming that the local compositions can be depicted by a relationship similar to that acquired for binary systems [46]: where x ji and x ki are the mole fractions of components j and k, respectively, in the neighbourhood of molecule i; g ji and g ki are the Gibbs energies of interaction for the j + i and k + i pairs; and α ij and α ik are characteristic constants of the corresponding binary mixtures. The activity coefficients can be expressed as follows: where and The reactor energy balance, coupled with the material balance, can determine the heating or cooling requirements. The energy balance given for the batch stirred reactor can be expressed as follows: where T is the reaction temperature, r A is the rate of reaction, V is the volume, ∆H RX is the constant heat of reaction, U is the overall heat transfer coefficient, A is the batch reactor heat exchange area, T a is the ambient temperature, N i is the number of moles, and C pi is the average heat capacity of species i. The decomposition of formic acid can take place via two possible pathways dependent upon the catalyst, reaction temperature, and reactant concentrations. The dehydrogenation of formic acid is shown in Equation (1), and the dehydration reaction is shown by Equation (2). The former reaction is slightly exothermic, while the latter reaction is slightly endothermic. The type of catalyst used can determine which pathway is favoured for the decomposition of formic acid [47,48].
The rate of reaction based on the dehydrogenation of formic acid can be expressed using the Power-Law method: where r is the rate of reaction, k is the kinetic rate constant, C is the concentration of formic acid, and n is the order of reaction [5]. The estimated activation value of 22 kJ/mol was obtained for the catalyst used in this study (Pd/Zn catalyst), and the estimation was based on the experimental data of this study. The order of the reaction was accepted to be 0.2, the same as for Pd/C [5], which reproduced the experimental results satisfactorily.

Results and Discussion
The 2%Pd 6 Zn 4 @HHT catalyst was prepared using the galvanic replacement procedure employed by Liao et al. [49], modified for the Pd-Zn system. High-heat-treated carbon nanofibers (HHT-CNFs) were utilised as support for the NPs.

TEM Results
The PdZn catalyst was characterised by TEM ( Figure 1) to explore the bulk properties of the system (bulk zinc palladium atomic ratio, particle size distribution, and average particle size). The catalysts displayed a mean particle size of 2.1 nm, with a narrow particle size distribution ( Figure 1b) and with particles well dispersed on the surface of the carbon support. On the other hand, the Pd-monometallic counterpart (1% Pd/HHT) showed a higher average particle size of 3.9 nm [4]. STEM-XEDS analysis confirmed the presence of Pd-Zn bimetallic particles with an average molar ratio of 86:14 lower than the nominal one (60:40) (Figure 1a). Figure 2 displays STEM images representing the d spacing.
obtained for the catalyst used in this study (Pd/Zn catalyst), and the estimation was based on the experimental data of this study. The order of the reaction was accepted to be 0.2, the same as for Pd/C [5], which reproduced the experimental results satisfactorily.

Results and Discussion
The 2%Pd6Zn4@HHT catalyst was prepared using the galvanic replacement procedure employed by Liao et al. [49], modified for the Pd-Zn system. High-heat-treated carbon nanofibers (HHT-CNFs) were utilised as support for the NPs.

TEM Results
The PdZn catalyst was characterised by TEM ( Figure 1) to explore the bulk properties of the system (bulk zinc palladium atomic ratio, particle size distribution, and average particle size). The catalysts displayed a mean particle size of 2.1 nm, with a narrow particle size distribution ( Figure 1b) and with particles well dispersed on the surface of the carbon support. On the other hand, the Pd-monometallic counterpart (1% Pd/HHT) showed a higher average particle size of 3.9 nm [4]. STEM-XEDS analysis confirmed the presence of Pd-Zn bimetallic particles with an average molar ratio of 86:14 lower than the nominal one (60:40) (Figure 1a). Figure 2 displays STEM images representing the d spacing.

XPS Results
X-ray photoelectron spectroscopy (XPS) was performed on the material, obtaining information on the surface properties of the PdZn catalyst. From the XPS survey spectra, both metals were detected. The Pd/Zn atomic ratio detected by XPS is higher than the nominal one but similar to the EDX atomic ratio (Table 1), which indicates an enrichment of Pd on both surface and bulk. In the high-resolution spectra of the Pd 3d region (Figure 3), both metallic Pd and Pd oxide were detected with relative amounts of 9% and 91%, respectively. The novel catalyst shown in this study possesses a high Pd(II) content with respect to the monometallic 1% Pd/HHT, which has a Pd(0) content of 52% [4], suggesting an oxidation effect due to the presence of Zn. By analysing the high-resolution spectra of the Zn 2p region, only ZnO was detected.

XPS Results
X-ray photoelectron spectroscopy (XPS) was performed on the material, obtaining information on the surface properties of the PdZn catalyst. From the XPS survey spectra, both metals were detected. The Pd/Zn atomic ratio detected by XPS is higher than the nominal one but similar to the EDX atomic ratio (Table 1), which indicates an enrichment of Pd on both surface and bulk. In the high-resolution spectra of the Pd 3d region (Figure 3), both metallic Pd and Pd oxide were detected with relative amounts of 9% and 91%, respectively. The novel catalyst shown in this study possesses a high Pd(II) content with respect to the monometallic 1% Pd/HHT, which has a Pd(0) content of 52% [4], suggesting an oxidation effect due to the presence of Zn. By analysing the high-resolution spectra of the Zn 2p region, only ZnO was detected.

Batch Reactor Model Validation
A formic acid decomposition reaction was used as a model reaction for hydrogen production to evaluate the catalytic performance of the material. This reaction occurs in two different pathways: dehydrogenation and dehydration. The stirring rate (1400 rpm) was optimised in previous studies in order to ensure the chemical kinetic regime [4]. The experimental results for the decomposition of formic acid were compared to the processmodelling results to assess their validity. Figure 4 depicts the conversion of formic acid at varying reaction temperatures of 30-60 °C. It can be observed that the conversion increases with increasing temperatures, as expected. Furthermore, there is good validation between the experimental and theoretical results. Usually, temperatures greater than 60 °C are not investigated due to the requirement of mild operating conditions in portable devices, which utilise formic acid fuel cells. According to the Arrhenius expression, = , an apparent activation energy of approximately 22 kJ/mol was obtained for this reaction with a 2%Pd6Zn4 catalyst for a batch reactor configuration. The catalytic deactivation is stronger at higher temperatures; therefore, higher reaction temperatures promote the dehydration pathway and CO formation, causing the catalyst to deactivate faster.

Batch Reactor Model Validation
A formic acid decomposition reaction was used as a model reaction for hydrogen production to evaluate the catalytic performance of the material. This reaction occurs in two different pathways: dehydrogenation and dehydration. The stirring rate (1400 rpm) was optimised in previous studies in order to ensure the chemical kinetic regime [4]. The experimental results for the decomposition of formic acid were compared to the processmodelling results to assess their validity. Figure 4 depicts the conversion of formic acid at varying reaction temperatures of 30-60 • C. It can be observed that the conversion increases with increasing temperatures, as expected. Furthermore, there is good validation between the experimental and theoretical results. Usually, temperatures greater than 60 • C are not investigated due to the requirement of mild operating conditions in portable devices, which utilise formic acid fuel cells. According to the Arrhenius expression, k = Ae − Ea RT , an apparent activation energy of approximately 22 kJ/mol was obtained for this reaction with a 2%Pd 6 Zn 4 catalyst for a batch reactor configuration. The catalytic deactivation is stronger at higher temperatures; therefore, higher reaction temperatures promote the dehydration pathway and CO formation, causing the catalyst to deactivate faster.
Kinetic profiles of the catalyst were obtained at different concentrations (0.2-0.5 M). Figure 5 demonstrates the effect of formic acid concentration on the conversion. The varying concentrations of formic acid used were 0.2-0.5 M at a reaction temperature of 30 • C and a substrate/metal molar ratio of 2000. The kinetic regime and these optimal reaction conditions were previously obtained by Sanchez et al. [5]. The results show that higher concentrations lead to higher conversions. In addition, the process-modelling studies were performed at constant reaction conditions with the experiment, and it can be concluded that there is good agreement between the experimental and theoretical data for the batch-reactor system.  Figure 5 demonstrates the effect of formic acid concentration on the conversion. The varying concentrations of formic acid used were 0.2-0.5 M at a reaction temperature of 30 °C and a substrate/metal molar ratio of 2000. The kinetic regime and these optimal reaction conditions were previously obtained by Sanchez et al. [5]. The results show that higher concentrations lead to higher conversions. In addition, the process-modelling studies were performed at constant reaction conditions with the experiment, and it can be concluded that there is good agreement between the experimental and theoretical data for the batch-reactor system.

Pd/Zn vs. Pd/C Catalyst
The experimental and process-modelling results obtained using the novel 2%Pd6Zn4 catalyst were compared with a commercial 5 wt% Pd/C catalyst using a batch reactor [5]. In addition, the catalytic performance of the material presented in this work was compared with published data on 1% Pd/HHT [4,17]. The decomposition of formic acid under mild conditions using a Pd/C and 1% Pd/HHT catalyst occurred under similar conditions to those of the current work. The study was conducted using a two-necked 100 mL roundbottom flask, which was immersed in an oil bath alongside a reflux condenser and a magnetic stirrer. The reaction took place at 30 °C, and 0.5 M of aqueous formic acid solution was placed into the reactor. The substrate/metal molar ratio was 2000:1, and the stirring rate was 1400 rpm [4,5,17].   Figure 5 demonstrates the effect of formic acid concentration on the conversion. The varying concentrations of formic acid used were 0.2-0.5 M at a reaction temperature of 30 °C and a substrate/metal molar ratio of 2000. The kinetic regime and these optimal reaction conditions were previously obtained by Sanchez et al. [5]. The results show that higher concentrations lead to higher conversions. In addition, the process-modelling studies were performed at constant reaction conditions with the experiment, and it can be concluded that there is good agreement between the experimental and theoretical data for the batch-reactor system.

Pd/Zn vs. Pd/C Catalyst
The experimental and process-modelling results obtained using the novel 2%Pd6Zn4 catalyst were compared with a commercial 5 wt% Pd/C catalyst using a batch reactor [5]. In addition, the catalytic performance of the material presented in this work was compared with published data on 1% Pd/HHT [4,17]. The decomposition of formic acid under mild conditions using a Pd/C and 1% Pd/HHT catalyst occurred under similar conditions to those of the current work. The study was conducted using a two-necked 100 mL roundbottom flask, which was immersed in an oil bath alongside a reflux condenser and a magnetic stirrer. The reaction took place at 30 °C, and 0.5 M of aqueous formic acid solution was placed into the reactor. The substrate/metal molar ratio was 2000:1, and the stirring rate was 1400 rpm [4,5,17].

Pd/Zn vs. Pd/C Catalyst
The experimental and process-modelling results obtained using the novel 2%Pd 6 Zn 4 catalyst were compared with a commercial 5 wt.% Pd/C catalyst using a batch reactor [5]. In addition, the catalytic performance of the material presented in this work was compared with published data on 1% Pd/HHT [4,17]. The decomposition of formic acid under mild conditions using a Pd/C and 1% Pd/HHT catalyst occurred under similar conditions to those of the current work. The study was conducted using a two-necked 100 mL round-bottom flask, which was immersed in an oil bath alongside a reflux condenser and a magnetic stirrer. The reaction took place at 30 • C, and 0.5 M of aqueous formic acid solution was placed into the reactor. The substrate/metal molar ratio was 2000:1, and the stirring rate was 1400 rpm [4,5,17].
The process-simulation studies took place under constant reaction conditions for the comparison between the two catalysts to assess their performance for the reaction. From Table 2, it can be observed that, under constant reaction conditions, the performance of the novel Pd/Zn catalyst is significantly higher than that of the commercial Pd/C catalyst and of the monometallic counterpart, underlining the role of Zn in the catalytic process. As a result, the novel Pd/Zn catalyst could be utilised to enhance the decomposition of formic acid reaction. Furthermore, there is good agreement (less than 5% difference) between the process-simulation modelling and experimental data. Therefore, the theoretical models created in this study have the ability to successfully predict the decomposition of formic acid under mild conditions. A gas analysis was performed on the system using an online micro-GC in order to calculate CO/CO 2 ratio, and a selectivity of 76% for the dehydrogenation reaction was obtained.

Conclusions
The decomposition of formic acid has been investigated using a novel 2%Pd 6 Zn 4 catalyst in a batch reactor. The studies were conducted experimentally, and a processsimulation model was developed to validate the results. It was found that the conversion of formic acid increases with temperatures in the range of 30-60 • C. Temperatures greater than 60 • C are not typically investigated due to the requirement of mild operating conditions for fuel cell applications. Furthermore, the formic acid conversion can be enhanced by increasing the concentration of formic acid in the batch system. For all of the parameter studies performed, the catalyst displayed no loss of activity.
The experimental and modelling data obtained for the 2%Pd 6 Zn 4 catalyst were compared with a commercial Pd/C catalyst. It was found that the novel catalyst demonstrated significantly higher conversions than the commercial catalyst under constant reaction conditions. Furthermore, the theoretical models developed to predict the decomposition of formic acid for both catalysts depicted sound validation with the experimental data, highlighting the importance of utilising computational software to further facilitate the design of experimental studies. The process-simulation-modelling methodology offered in this study can potentially offer shorter computational times and requires less computing power when compared with CFD modelling. Future research could be directed towards implementing the catalyst in a continuous flow system to allow for constant generation of hydrogen for fuel cell applications.