Abstract

The use of propane fuel in high temperature (120°C) polymer electrolyte membrane (PEM) fuel cells that do not require a platinum group metal catalyst is being investigated in our laboratory. Density functional theory (DFT) was used to determine propane adsorption energies, desorption energies, and transition state energies for both dehydrogenation and hydroxylation reactions on a Ni(100) anode catalyst surface. The Boltzmann factor for the hydroxylation of a propyl species to form propanol and its subsequent desorption was compared to that for the dehydrogenation of a propyl species. The large ratio of the respective Boltzmann factors indicated that the formation of a completely reacted product (carbon dioxide) is much more likely than the formation of partially reacted products (alcohols, aldehydes, carboxylic acids, and carbon monoxide). That finding is evidence for the major proportion of the chemical energy of the propane fuel being converted to either electrical or thermal energy in the fuel cell rather than remaining unused when partially reacted species are formed.

1. Introduction

Fuel cells convert the chemical energy of a fuel into electrical energy. Theoretically they can produce more electrical energy from a fuel than either batteries or combustion processes. Propane was the fuel investigated for the fuel cells described in this work. Trucks currently deliver liquefied petroleum gas (LPG) (mostly propane) in rural areas where roads exist. The cost of delivering conventional electrical power in rural areas is approximately an order of magnitude greater than that in urban areas, even though the price charged by utility companies for electrical power is often similar. Therefore, a greater capital cost for fuel cells can be accepted in rural areas than in urban areas. This indicates that a niche market for propane fuel cells in rural areas would be profitable prior to a profitable market in urban areas.

Most fuel cells use either hydrogen or methanol as the fuel. They have several disadvantages. Unfortunately no infrastructure exists for their distribution and storage. Furthermore they are both manufactured from natural gas (primarily methane) using complex reactor systems that have a large capital cost. They also have a large operating cost because 25% of the natural gas is consumed to provide the endothermic heat for the steam reforming reaction. Hydrogen gas requires sophisticated storage systems [1]. None of those disadvantages occur when propane reacts directly at the anode of a fuel cell.

The type of propane fuel cell we are investigating is a high temperature polymer electrolyte membrane (PEM) fuel cell operating near 120°C. It has several advantages over conventional PEM fuel cells that operate near 80°C. The principal one is the elimination of liquid phase water, which eliminates the accumulation of acidic species and the resulting corrosion. With a less corrosive environment, the platinum group metals used as catalysts in conventional PEM fuel cells can be replaced by less expensive metals such as nickel (at the anode) and silver (at the cathode). Unfortunately the proton conductivity of the Nafion electrolyte used in conventional PEM fuel cells diminishes by an order of magnitude when the operating temperature increases from 80 to 120°C. To circumvent this difficulty we have begun an experimental program on the development of a modified zirconium phosphate, ZrP, membrane that appears to be suitable for 120°C operation (Al-Othman et al. [2, 3]).

The primary challenge with fuel cells that use propane and other hydrocarbons as their fuel is that the current densities (reaction rates) are small. During the 1960s an extensive amount of research on direct hydrocarbon fuel cells (DHFCs) including propane fuel cells was performed [46]. Since that time, there have been relatively few experimental investigations. They include hydrocarbons being fed to the anode of solid oxide fuel cells [7] operations at 300°C using 1 mm thick Sn0.9In0.1P2O7 [8], kinetics and mechanisms [911], and the use of polymer electrolytes [1214].

Fuel cell reactions can be examined at different levels of detail. The overall propane reaction in a propane high temperature PEM fuel cell is shown as The overall half-cell reaction for a propane fuel cell anode is shown as where (g) is gas phase, (ZrP) is Zr(HPO4)2H2O–glycerol electrolyte, and (Ni) is nickel metal anode electrocatalyst. Nickel was chosen for the anode electrocatalyst, because it was known to be an effective nonplatinum catalyst in alkaline fuel cells [15]. The mechanism by which the propane reaction in (2) is expected to occur is shown in Figure 1. The reaction mechanism in Figure 1 is analogous to the reaction mechanism for methane that was reported by Bagotzky et al. [16]. Propane in the gas phase is first adsorbed on the nickel metal electrocatalyst. Then it is dehydrogenated to form an adsorbed propyl species. The investigation reported here primarily concerns the subsequent reactions of the adsorbed propyl species. It could react with a hydroxyl group to form propanol (an alcohol), C3H7OH(Ni), followed by desorption to the gas phase. Alternatively a hydrogen atom could be removed to form an adsorbed allyl species, C3H6(Ni). An appropriate combination of subsequent dehydrogenation and hydroxylation reactions could lead to complete reaction and the formation of CO2. At the anode there is only one electrochemical reaction, H−−> H+ + e-. Figure 1 shows that the majority of the elementary fuel cell reactions on a propane fuel cell nickel anode electrocatalyst do not involve charge transfer.

The general purpose of this work was to investigate the reactions in the network shown in Figure 1. One specific objective was to determine whether or not the hydrocarbon reaction intermediates remain on the catalyst surface where a series of dehydrogenation reactions can combine with hydroxylation reactions to ultimately form carbon dioxide, the completely reacted product. The undesirable alternative would be for the hydrocarbon surface species to retain their hydrogen atoms and to form partially reacted products such as alcohols, aldehydes, and carboxylic acids. For fuel cells complete reaction is desired since some of the fuel’s chemical energy remains unused when partially reacted products are formed.

Catalyst composition appears to be the factor that determines if the products are partially reacted or completely reacted. Several experiments with a platinum anode have shown that CO2 was the predominant product formed in specific situations. The only product reported was CO2 when the hydrocarbon fuel was (a) ethylene being fed to a sulphuric acid fuel cell [17], (b) propane being fed to a phosphoric acid fuel cell [18], and (c) methane being fed to a PEMFC [19]. When a PEMFC with a platinum anode was being fed propane, the major product was CO2 at cell potentials less than 0.1 V and CO when the cell potential was greater than 0.4 V [20]. In contrast with a platinum-ruthenium anode, when methane was being fed to a PEMFC, all of the possible partial reaction products, CH3OH, CH2O, COOH, CO, and CO2, were formed [19]. This work was performed to obtain an indication of whether complete or partial reaction products would be obtained when using a nickel catalyst.

2. Methodology and Computational Details

SIESTA software based on Kohn-Sham density functional theory was used to perform quantum chemical computations (Soler et al. [21]). The generalized gradient approximation (GGA) method was used as the exchange correlation functional type, with Becke, Lee, Yang, and Parr (BLYP) parameterization. The default basis set, a double-ζ polarization set composed of a compact and a diffused orbital basis, was used in these calculations. The tolerance of the density matrix was set to 10−3. This value sets the maximum allowable difference between the output and the input on each element of the density matrix in a self-consistent field (SCF) cycle. A (/  -points) Monkhorst-Pack -point mesh was used. The convergence as a function of the number of -points was carefully monitored. Increasing the -point mesh from 32 -points to 48 -points changed the adsorption energies by an insignificant amount (~4 × 10−5 eV or ~8 × 10−3 kJ/mol). In some calculations, the coordinates of atoms were allowed to relax (change position) to determine the geometry having the minimum energy. When the difference in energy between successive calculations was less than 10−3 eV and when the maximum atomic force was less than 0.01 eV/Å, the convergence criteria were attained and the atomic geometry was optimized.

The pseudopotential for nickel was generated using the ATOM program of SIESTA. The Perdew Burke Ernzerhof (PBE) [22] exchange correlation was used for the generation of the electronic configuration. The improved Troullier-Martins (tm2) method was used to generate the pseudopotential file for nickel in its nonpolarized ground state electron configuration. The pseudopotential input file for a metal required that a particular core radius, , should be specified. It was obtained by trial and error. If the initial value was too small (1.5 for metals used in this study), the pseudopotential could not be generated. After a value of that was slightly greater than the minimum acceptable value had been obtained, the software would generate an exact value for . That value (3.34 bohr for nickel) was used to generate the pseudopotential. The pseudopotential output file was used for all subsequent computations that included nickel atoms. To confirm the validity of our pseudopotential generation method and the resulting pseudopotentials, energy calculations were performed for a Ni unit cell using our custom generated Ni pseudopotential. The lattice constant (LC) that had the minimum total energy for the Ni unit cell was 0.370 nm. There is a difference of 4.9% when the 3.70 nm value is compared to the experimental LC for Ni (0.352 nm) that was measured by X-ray diffraction (XRD; Haynes [23]). This error is smaller than the 6% error that was reported when using the PBE functional (Perdew et al. [24]).

DFT calculations were performed on a system of periodically repeated entities, defined in relation to nickel crystal unit cells. Each entity consisted of one propane molecule, a slab of nickel atoms, and a vacuum space. The slab consisted of one layer of crystal unit cells (two layers of Ni atoms) located on an - plane. The cells were arranged 3 cells 3 cells with 4 atoms per cell ( = 36 atoms per slab). The (100) surface was used for all nickel slabs. The vacuum space had a thickness of nine empty unit cells and was located above the nickel slab in the direction. The vacuum space occupied the distance from the upper surface of a nickel slab within one periodic entity to the bottom surface of the next nickel slab above it in the next periodic entity. Ideally the vacuum layer can prevent or minimize interactions between the periodic entities in the direction perpendicular to the surface. The choice of slabs with 36 nickel atoms per slab is typical of the number of metal atoms in many DFT slab calculations [2528].

The slab surface should be large enough, in relation to the surface species, to avoid both (a) interactions between the surface species on one slab and the nickel atoms in neighbouring slabs and (b) adsorbate-adsorbate interactions. Because there were 18 atoms in the surface of our slabs, the fractional surface coverage of propane was less than 0.2. It has been shown by Grabow et al. [29] that adsorbate-adsorbate interactions can influence the binding energy by as much as 1 eV when the fractional surface coverage increases from 0.2 to 1. Our low surface coverage minimized absorbate-adsorbate contributions to the adsorption energy.

To confirm that the slabs were thick enough to predict propane adsorption on Ni(100) the thickness of the slabs was increased from 2 atoms to 4 atoms (increasing the number of atoms per slab from 36 to 72). The increased slab thickness only changed the value of the propane adsorption energy by 7.3% (~0.04 eV or ~4 kJ/mol). This relatively small variation indicates that calculations using slabs having two layers of atoms are sufficiently accurate to qualitatively predict trends in propane adsorption energy.

Another single adsorption energy calculation for propane was performed on a 2-layer Ni(100) slab in which the bottom layer was fixed at the positions of bulk nickel and the top layer was allowed to relax. When compared to the calculation in which the atoms in both layers were allowed to relax, the change in energy was also very small (4 × 10−4 eV ~ 0.04 kJ/mol).

The adsorption energies () of the adsorbates (initially adsorbed reactants and subsequent adsorbed intermediate species) on the metal slabs were calculated according to the following equation: where is the total energy for the absorbate plus metal slab, is the total energy for the configuration of the adsorbate plus metal slab that has the minimum energy, is the energy of an isolated reactant molecule in the gas phase, is the energy for a metal slab, and is the energy for adsorption of the species on the surface of the metal slab (or heat of adsorption).

The transition state energy, , for each reaction was obtained by calculating as a function of either the C–H bond length or the O–H bond length. For each calculation either the C–H bond length or the O–H bond length was maintained constant, while all other bond lengths and all bond angles in the adsorbed species were relaxed to obtain the adsorbate configuration having the minimum energy for that particular C–H or O–H bond length. Subsequently the resulting energies that had been calculated at each C–H or O–H bond length were compared as a function of their bond lengths. The configuration of the adsorbed species having the maximum energy, , was the transition state (TS).

The transition state energies were obtained using the following equation: Variations in the electrical potential of the slabs were not considered (zero electric field).

3. Results and Discussion

A dehydrogenation mechanism for the transformation of propane, in the gas phase, to an adsorbed propyl species is shown in Figure 2. The formula –(Ni–Ni)–SLAB represents a 36-atom slab of Ni atoms in two layers in which the Ni atoms in each layer have fourfold symmetry. The adsorbed propane (Figure 3(a)) is oriented parallel to the Ni slab surface. That is, a plane drawn though the three carbon atoms in the adsorbed propane (Figure 3(a)) is perpendicular to the Ni(100) surface. The central carbon atom of the propane has on-top coordination with the nickel atom below it. Relative to the central carbon atom, the two propane end carbon atoms are symmetrical and are elevated from the slab surface. The hydrogen atoms that are bonded to the propane central carbon atom are also coordinated to the same Ni atom as the propane central carbon atom. In Figure 3 the C–H bond distance between the central carbon atom and its hydrogen atoms is 0.110 nm compared to 0.112 nm when propane is in the gas phase.

The adsorption energy of propane on Ni, 0.534 eV (51.2 kJ/mole), is also shown in Figure 3. Other computed values for propane adsorption on metals include values of 33–38 kJ/mol on Pt, shown by Nykänen and Honkala [30], a value of 55 kJ/mole on Pd, shown by Antony et al. [31] and 39.8–42.1 on Pd, shown by Kao and Madix [32]. Experimental values for propane adsorption on metals include a calorimetric measurement of 32 kJ/mol on a Mo film by Smutek and Černý [33], and temperature programmed desorption measurements on Pd, shown of 41.5 kJ/mol, shown by Kao and Madix [32] and 45 kJ/mol, shown by Antony et al. [31]. Nykänen and Honkala [30] noted that there is considerable uncertainty in heat measurements made by temperature programmed desorption. Furthermore the adsorption energy on Ni would be expected to be different than those on Mo or Pd. Nevertheless, the propane adsorption energy on Ni reported here is of the same magnitude as the values obtained for other metals either by computation or by measurement.

The transition state (TS) between an adsorbed propane molecule and an adsorbed propyl species is shown in Figure 3(b). An important feature of the TS is that both the propyl species and its separated hydrogen atom remain configured to the same Ni atom in the slab (structure A3 in Figure 2). The TS structure is neither symmetrical nor parallel to the Ni slab surface. A plane drawn through the three carbon atoms of the TS structure (Figure 3(b)) intersects with the surface of the Ni(100) slab at an angle that is different from 90°. The carbon-hydrogen distance for the hydrogen being removed from the central propane carbon atom is shown at a distance of 0.170 nm in the TS. That is substantially larger than the 0.110 nm C–H distance in the adsorbed propane.

Figure 4 provides the energies of the four species shown in Figure 2. The extent of the dehydrogenation reaction is indicted by the C–H bond distance between the propane central carbon atom and the hydrogen atom being removed. The energy change upon propane adsorption on the Ni(100) slab is 0.531 ev (51.2 kJ/mole). The barrier height between an adsorbed propane molecule and its TS to an adsorbed propyl species is 0.734 eV (70.8 kJ/mole). The energy of reaction from an adsorbed propane molecule to an adsorbed propyl species and an adsorbed hydrogen atom is 0.163 eV (15.7 kJ/mole). The transition state energy = 0.734 eV we obtained by varying the C–H bond distance is a close approximation to the true transition state energy that would have been obtained from a saddle point calculation.

A mechanism for the dissociation of water is shown in Figure 5. The structures of the water derived adsorbed species are shown in Figure 6. The top view of an adsorbed water molecule suggests that the threefold nature of the three atoms in the water molecule (one oxygen atom and two hydrogen atoms) appears to be coordinated on top of a single Ni atom in the Ni(100) surface. The O–H bond distance in adsorbed water, 0.100 nm, is different than that in gas phase water, 0.145 nm.

The transition state (TS) from adsorbed water to the combination of an adsorbed hydroxyl species and an adsorbed hydrogen atom is shown in Figure 6(b). The O–H bond distance in the TS has increased to 0.160 nm which is a substantial change when compared to the 0.100 nm O–H bond distance in adsorbed water.

The energies of the four species in Figure 5 are shown in Figure 7, as a function of the dissociation reaction, as measured by the O–H bond distance. The energy change upon water adsorption on the Ni(100) surface is 0.798 ev (77.0 kJ/mole). The barrier height between an adsorbed water molecule and its TS to an adsorbed hydroxyl radical is 0.655 eV (63.2 kJ/mole). There are numerous reports in the literature on the dissociation of water on metals. For example, a DFT study [34] of 6 different metals, Os, Ru, P, Ir, Rh, and Pt, reported barrier heights in the range of 0.343 to 1.25 eV. Another example is a study by Pozzo et al. [35] comparing the dissociation of water on Rh(111) and Ni(111) surfaces. They reported barrier heights of 0.92 and 0.89 for the water dissociation reaction to OH and H on Rh(111) and Ni(111), respectively. The energy change during reaction of an adsorbed water molecule to form an adsorbed hydroxyl species and an adsorbed hydrogen atom is 0.386 eV (37.2 kJ/mole). The barrier height for the combination of formation of hydroxyl species and hydrogen atoms (Figure 7) is approximately 10% less than that for the combination of the formation of propyl species and hydrogen atoms (Figure 4). If a comparison of two different single barrier heights was the only consideration, it would suggest that it might be easier to form hydroxyl species than propyl species and that hydroxyl species might be more abundant on the Ni(100) surface than propyl species.

A mechanism for the second dehydrogenation reaction, conversion of adsorbed propyl species to adsorbed allyl species, is shown in Figure 8. The structures of the adsorbed propyl species and the TS from an adsorbed propyl species to an adsorbed allyl species are shown in Figure 9. The central carbon atom of the adsorbed propyl radical in Figure 9(a) is on top of one of the Ni atoms in the Ni(100) slab. A plane drawn through the three carbon atoms of the propyl species intersects the surface of the intersects the Ni(100) surface at an angle that is much different than 90°, indicating that the propyl species is not situated symmetrically on the Ni(100) surface.

The transition state (TS) from an adsorbed propyl species to an adsorbed allyl species is shown in Figure 9(b). The central carbon atom of the TS is on top of a Ni atom in the surface. The hydrogen atom that is being removed from the central carbon atom is bridge bonded with two Ni atoms that are adjacent to the Ni atom on which the allyl central carbon atom is on top. The distance between the central carbon atom in the TS and the separated hydrogen in the TS has increased to approximately 0.180 nm which is a substantial change when compared to the C–H bond distance in an adsorbed propyl species.

The energies of the three species in the transition from propyl to allyl, in Figure 8, are shown in Figure 10, as a function of the dehydrogenation reaction, as measured by the C–H bond distance. There is a difference in energy of approximately 16 eV between the combination of the propyl species and slab that is represented in Figures 8 and 10 and the combination of propyl radical, hydrogen radical, and slab represented in Figures 2 and 4. The barrier height between an adsorbed propyl species and the transition state (TS) to an adsorbed allyl species and a hydrogen atom, is 1.22 eV (~118 kJ/mole).

A mechanism for the synthesis reaction, between an adsorbed propyl species and an adsorbed hydroxyl species to form propanol, is shown in Figure 11. The structures of the adsorbed propyl and hydroxyl species are shown in Figure 12(a). The central carbon atom of the adsorbed propyl species in Figure 9(a) is on top of one of the Ni atoms in the Ni(100) surface. A plane drawn through the three carbon atoms of the propyl radical intersects the surface of the Ni(100) slab at an angle that is much different than 90°, indicating that the propyl radical is not situated symmetrically on the Ni(100) surface. The C–O bond distance between the central carbon atom in the propyl species and the oxygen atom in the hydroxyl species is 0.4 nm.

The transition state (TS) during the reaction of an adsorbed propyl species and an adsorbed hydroxyl species to form an adsorbed propanol molecule is shown in Figure 12(b). The C–O bond distance in the TS has decreased to 0.215 nm from 0.4 nm for the C–O distance when the two species were in their most stable independent states. It should be noted that both the carbon atom in the propyl species and the oxygen atom in the hydroxyl species are coordinated to the same Ni atom in the TS, prior to the synthesis reaction and the formation of a C–O bond. It appears that they need to be adsorbed on the same atom in order to be geometrically close enough for the bond formation reaction to occur.

The energies of the four species (propyl, hydroxyl, adsorbed propanol, and propanol in the gas phase) are shown in Figure 13, as a function of the C–O distance. The energy of reaction from the combination of an adsorbed propyl species plus an adsorbed hydroxyl species to form an adsorbed propanol molecule is 0.507 eV (48.9 kJ/mole). The barrier height for formation of propanol from combination of an adsorbed propyl species and an adsorbed hydroxyl species is 1.59 eV (~153 kJ/mole). In contrast the barrier height for the formation of allyl species (1.22 eV), discussed above, is approximately 25% less than that of formation of propanol (1.59 eV). Furthermore, to desorb propanol from the Ni surface into the gas phase, an additional amount of energy 0.819 eV (~79.0 J/mole) must be provided.

The Boltzmann factor, exp(/RT), represents the probability that molecules have sufficient energy to attain the transition state, where is the barrier height for the reaction. It can be used to estimate the ratio of two reaction rates. The transformation of a propyl species to an allyl species [with barrier height ] is the desired reaction because the species stays on the surface and therefore retains the possibility of eventually being converted to CO2, the complete reaction product. The reaction of a propyl radical with a hydroxyl species to form propanol in the gas phase [with barrier height ] is the undesired reaction, because propanol is a partially reacted product that retains chemical energy that could have been converted to electrical energy by further reaction to carbon dioxide.

The barrier height for the transformation of adsorbed propyl species to adsorbed allyl species () = 1.22 eV can be compared with the barrier height for the transformation of the adsorbed propyl species to gas phase propanol. Two reaction steps are used to transform propyl species to gas phase propanol molecules. First, there is the reaction of an adsorbed propyl species with a hydroxyl radical to form adsorbed propanol, = 1.59 eV. Then there is the desorption of propanol into the gas phase, = 0.82 eV. Because two steps are required for one transformation and one step for the other, a comparison of two different single barrier heights will not describe the phenomenon. The combination exp[/RT] = exp(−1.59/RT) * exp (−0.82/RT) giving = 2.41 eV represents the product of two phenomena.

The ratio of the two reaction rates is The Ratio calculated at 150°C (423 K) was 1.4 * 1014. That result indicates that propanol is an unlikely product. We are not aware of any experimental fuel cell results on the selectivity for CO2 on a nickel anode. However, platinum is in the same column of the periodic table as nickel, and platinum is frequently used as a fuel cell electrode. These computed results for an Ni(100) surface are consistent with experimental results on platinum catalysts that found carbon dioxide as the predominant product. This study indicates that carbon dioxide may also be the predominant product on nickel and that the formation of substantial amounts of partially oxygenated products is unlikely.

4. Conclusion

The DFT results reported here indicate that the reaction of an adsorbed propyl species to form an adsorbed allyl species on a nickel (100) anode catalyst surface is much more likely than its reaction with a surface hydroxyl species to form gas phase propanol. That suggests that the reaction products from a propane fuel cell having an anode with a nickel electrocatalyst will likely be completely reacted (CO2) rather than partially reacted (alcohols, aldehydes, carboxylic acids, and CO), thereby converting essentially all the chemical energy in the propane to electrical and thermal energy.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.