Material properties and process modeling with density functional theory (DFT) is an accurate method to facilitate the study and the design of materials computationally for the development of different electrochemical technologies such as fuel cells, solar cells, and batteries, among others, mainly to achieve alternative ways for energy conversion and storage. Considering the relevance of DFT in the development of these alternative technologies for energy generation and storage, in this chapter, the application of DFT to study catalytic alloys and their reactivity processes to develop polymer membrane fuel cells (PEMFCs) is presented. In this sense, firstly, a brief review of the application of DFT to develop catalysts for PEMFCs and the relation with the concept of chemical descriptors is presented. Secondly, the main chemical descriptors for this task are presented and discussed. Finally, a summary of the main findings of the modeling with DFT and chemical descriptors approach of catalytic alloys for PEMFCs is presented and analyzed.
- chemical descriptors
- catalytic alloys
The rational design of materials is an important and challenging task, which implies the in-depth understanding of determining factors for the desired properties of the materials. These properties emerge from the material chemistry (atoms involved and its electronic structure) and material structure (organization at nano and microscale) [1, 2], for which the understanding at atomistic level is key and could enable the tailoring of materials that improve the efficiency of almost any technology. Particularly, catalyst development is crucial for efficient chemical production and clean energy generation and storage technologies such as photoelectrochemical cells, fuel cells, and batteries, which are required for an environmentally responsible economic growth .
In the last few decades, density functional theory (DFT) has been successfully used to study heterogeneous catalysis process and the properties of the catalysts, achieving an understanding of reaction mechanisms and the determining factors in the catalytic activity of the materials. For example, it has been found that the electronic structure of the catalytic surface determines its properties [4, 5, 6]. Also, some general relations between these factors and the catalytic activity have been developed and used for the search and design of catalytic materials for different applications [7, 8, 9, 10, 11, 12, 13, 14].
In the next sections, the application of DFT to study catalytic materials and their reactivity to develop polymer membrane fuel cells (PEMFCs) is presented. The first section is devoted to review briefly and in a general way, the different works, where DFT is applied for studying catalysts for the electrodes of PEMFCs and some determining factors on catalysis found in some of these works are drawn. Secondly, the main factors used for the development of catalysts through the concept of chemical descriptors are presented and discussed. Finally, a summary of the main findings of the modeling with DFT and chemical descriptors approach of catalytic alloys for PEMFCs is presented and analyzed.
2. Using DFT to study heterogeneous catalysis in PEMFCs
Fuel cells are devices that convert chemical energy of supplied fuel and oxidant into electrical energy. Their operation involves two electrodes, an anode and a cathode, separated by an electrolyte. The anode fuel, such as hydrogen, methanol, and ethanol, which have been proposed for PEMFCs, is oxidized and donates electrons. Electrons travel through the extern circuit, where they perform work and reach the cathode, where oxidant, generally oxygen, is reduced combining with electrons. The electrolyte, which in the case of PEMFCs is a conducting polymer, conducts electricity through the movement of ions, completing the circuit. A schematic representation of a unit cell, illustrating the operation principle explained above is shown in Figure 1.
The PEMFCs require the improvement of their performance and the decreasing cost. A principal factor for their performance and its cost are the electrodes, where Pt catalysts are commonly used. So, it required the development of low-cost catalysts, mainly for the oxygen reduction reaction (ORR) at the cathode and hydrogen, methanol or ethanol oxidation reaction (EOR) at anode. The catalysts should reduce the Pt content increasing the rates of the corresponding reactions. Looking for better catalysts is an important challenge, which should be addressed with a comprehensive approach that involves the understanding of the reasons for which a material is an efficient catalyst for a given reaction; in this way, materials can be designed with a “scientific” base.
In this sense, the use of simulations to study catalysis is a great tool, particularly, with the use of atomistic models based on DFT and molecular dynamics. DFT is suited to study adsorbate-catalyst interactions, which allow determining the potential energy surface and thus the energy barrier of the chemical transformations, also allows calculating the electronic structure of materials, which can help to understand bond-breaking and bond-making process. Also DFT, through thermodynamics approach, facilitates to describe reaction process and energetics of intermediates, reactants, and products under realistic conditions. Additionally, the information from atomistic scale can be used to compute reaction rates. In the next paragraphs, some works about catalysts for PEMFCs are presented.
Heterogeneous catalytic reactions in PEMFCs involve adsorption of reactants, different bond breaking and/or bond making processes and some electron transfer processes, and desorption of products. The reaction follows the path with minimal energy to go from reactants to product configuration, so this path depends, essentially, on the energetics of the possible elementary steps. So, many studies of catalysts applying DFT is focused in considering each possible reaction step to calculate energetics of reactants and products of each step from which energy barrier for each step can be calculated with different methods such as nudged elastic band [15, 16, 17], synchronous transit-guided quasi Newton (STQN) methods , among others. With this method, a free-energy diagram for different process can be drawn, an example for dissociative mechanism of oxygen on Pt is shown in Figure 2.
The above described approach was used to study EOR or the ethanol decomposition in other processes on Pt [19, 20], Rh , Pd , Au , and some Pt alloys as PtRh, PtRe, PtRu, PtSn, and PtPd . In some of these works, only some steps or the complete reaction network was considered. For example, it was calculated by Alcalá et al., the transition states for C-C and C-O bond cleavage on Pt using constrained optimization, where the length of the bond representing the reaction coordinate is constrained; they identified the steps with lowest activation energy for these bond cleavage reactions. The evaluation of complete reaction network is a formidable task due to the higher quantity of possible intermediate derivatives, but in the ethanol oxidation, reaction network is explored thoroughly identifying key steps for the selectivity on ethanol oxidation and the intermediates, which can lead to C-C bond cleavage in some surfaces.
Also, the mechanistic details of ORR on Pt and Pt-based alloys were investigated using DFT to calculate the intermediate energies and to estimate reaction and barrier energy of each step considered in the reaction pathway [25, 26, 27, 28, 29, 30, 31, 32, 33]. Different pathways are considered such as the O2 dissociation mechanism, and the associative mechanism implying the addition of H to the O2-containing species. Also, it considered two different ways for hydrogen transfer reaction within each of these pathways: the direct addition of H from the electrolyte (Eley-Rideal type) or the H adsorption and subsequent transfer to O2-containing species (Langmuir-Hinshelwood type) . For example, in , the dissociative mechanism finding that the high stability of oxygen and hydroxyl can imply kinetic limitations in proton/electron transfer to adsorbed oxygen or hydroxyl was studied. In other study, different mechanisms were considered and it is found that under determined conditions, the determining step is the O2 dissociation, so the associative mechanism is preferred .
One of the challenges in the study of electrochemical reactions is considering the effects of the solvation of molecules, the electric field in the interface electrolyte-electrode, and the other adsorbed molecules, nonetheless, some developments have successfully addressed those issues. It is widely accepted that the water layers between the polymer and the catalysts are very important in ORR. It was intended to study the effects of solvation using nine water molecules and one hydronium ion in an ab initio molecular dynamic study . Also, in other studies, adsorption energies of intermediates were calculated adding water adsorbed to the surface or water on top of the adsorbed intermediates [30, 35]. The solvent can be also described in averaged manner as dielectric continuum characterized by a dielectric constant, which is used to estimate free energy of solute-solvent interactions [36, 37, 38]; this approach is named implicit solvation method.
As it is envisioned above, significant progress for the accurate description of electrochemical reactions have been achieved. However, it also can be noted that the complete kinetic description of a given system is quite a demanding task. But, different works show that enough knowledge about catalytic activity could be gained studying some properties of the system at the atomic scale, which can be calculated in a simpler way. These properties are named descriptors, that is, an observable that can be directly linked to a fundamental property [10, 14]; in this case, the catalytic activity. For example, Nørskov et al. show that the oxygen binding energy on metal electrodes is related to the turnover frequency of the oxygen reduction reaction .
ORR involves two key steps, the H transfer to adsorbed O2 to form adsorbed OOH before the O-O bond breaking and the H transfer to adsorbed OH to form H2O. The activation energy of each step is related with the energy change in each process, which in turn is related with stability of adsorbed OOH and OH, respectively. These adsorbed species will be more stable if O binding to the surface is stronger; therefore, the O binding energy can be related with the reaction rates of these two key steps. The reasoning for this relation is as follows: If O binding is stronger, it is difficult to destabilize the Pt-OH bond to form H2O; if it is weak, it is difficult to stabilize the OOH to allow the O-O bond to break, so the O binding energy should be moderated. Precisely, in different works, bimetallic Pt and Pd alloys promising as catalysts for ORR looking for materials with moderated oxygen binding energy were identified [39, 40].
The descriptor approach also was used to search materials for hydrogen oxidation reaction on anode of PEMFC. For this reaction, the binding energy of H atom is a key, and according to Sabatier principle, optimal catalytic activity can be achieved on surfaces with moderated binding energies [9, 41]. Sabatier principle states that the bond between the catalyst and the key adsorbate should not be too weak, such that the reactants do not bind, and are not too strong, since this leads to catalyst poisoning . Similarly, in other work, the CO binding energy is considered along the H binding energy to find active catalysts tolerant to CO poisoning .
Other property that can be used as a descriptor of catalytic activity is the d-band center of the transition metal surfaces, that is, the average energy of the d states. As is proposed by the d-band model [6, 43], the d-band center is correlated with the binding energy of adsorbates to the surface. The d-band model is fulfilled for bonding of several atoms and molecules, then, as is shown in several works, tuning the electronic structure of a metal by the addition of alloying atom can lead to a better catalyst for ORR [44, 45, 46] or catalysts more tolerant to CO . Also, it identified other electronic structure factors, such as transition metal, e.g., band filling and covalency, that can be used as descriptors for ORR on metal oxide catalysts .
It is worth to mention that the review presented is very brief and there are other extensive reviews that should be consulted [48, 49]. Besides, from the works reviewed above, it should be remarked that it seems that one key property for catalytic activity is the adsorption energy of key intermediate molecules, which can be used to search better catalysts for a determined reaction. Next section reviews the principal concepts in catalysis used for the identification of appropriate descriptor, in this sense scaling relations between activation energies and binding energies, Sabatier principle and d-band model are reviewed.
3. Development of catalysts using chemical descriptors
In the last section, different properties of a material that can be related with its catalytic activity were observed. For instance, binding energies of atoms and molecules can be related with activation energies, so with reaction rates. Sabatier principle shows that a determined value of the adsorption energy of intermediates will lead to higher catalytic activity. Finally, it is mentioned that modification of the electronic structure changes the catalytic activity.
3.1. Correlations between adsorption and activation energies and Sabatier principle
The central issue in catalysis is the effect of the catalyst on the rate of a chemical reaction, the rate constant of an elementary reaction can be calculated using an Arrhenius expression in terms of a pre-factor () and an activation energy ()as follows:
The activation energy is an important quantity, which varies with the change of catalyst, so reviewing the origin of activation energy is useful to understand the effect of catalyst on the reaction rate.
As it is illustrated in Figure 3, any chemical reaction can be described as a transition between two local minima on the potential energy surface of the system as a function of the spatial coordinates of all the involved atoms. The reaction path is defined as the minimum energy path, which connects the configuration of the reactants with minimum energy (R) to the configuration of the products with minimum energy (P) along the reaction path; there is a configuration with the highest energy, which is called transition state (TS). The difference between the TS and R energies is the activation energy () [50, 51, 52]. The configurations of R, P, and TS in a catalytic process involve adsorbed molecules, then the energy of these configurations, and hence activation energy, depends on each surface.
Also, it is worth to remark that as transition states have the same basic physics as adsorbed species, then the transition state energy () is correlated with adsorption energies (). In fact, linear correlations between transition state and adsorption energies have been found [53, 54, 55, 56]. This type of relationship is a well-established approach in the understanding of trends in chemical reactions that dates back to Brønsted in 1928, also Evans and Polanyi in 1938 , who stated a linear relationship between activation energies and reaction energies, named Brønsted-Evans-Polanyi (BEP) relationship.
The linear relationship between of a chemical reaction and a particular relevant for that reaction can be understood thinking that TS is an absorbed specie with a configuration relatively close to the adsorption configuration of the molecule X with which is correlated. Then, it is reasonable to express the as the sum of two contributions, the first is part of adsorption energy of the molecule X (), which accounts, mainly, for the interaction energy between the surface and the molecule. The second is the quantity of energy (), which accounts for the reorganization of the molecule X to transform into the transition state molecule. So, the general form of linear relationship is:
It is common to use as the energy of products configuration. In this case, parameter is a value between 0 and 1, which indicates the position of the transition state along the reaction coordinate. If alpha is close to 1, the transition state has a configuration similar to the products (late transition state) , whereas is a part of the energy involved in the reverse reaction (energy to transform products into transition state molecule) .
Using DFT calculations to calculate transition state energies and product energies, it has been possible to establish transition state scaling relations with sufficient statistical significance. This linear relationship provides a straightforward way to study the effect of the surface on the activity for a given reaction, as it enables the estimation of activation energies from adsorption energies, which could be calculated in a simpler way using computational calculations.
Transition state scaling relationship implies that the activation energies are lower in surfaces, where the products are more stable (more negative adsorption energies of products) and hence the reaction rate is higher. However, the very strong adsorption of products leads to poison the surface, since the desorption of products is more difficult, hence the rate of desorption is lower. On the other hand, if adsorption of products is very weak, the reaction rate would be lower and the desorption rate higher. Therefore, the adsorption energy of products should not be too weak, such that the reaction does not happens, and not too strong since this leads to catalyst poisoning . The above mentioned is precisely the Sabatier principle, a valuable statement for catalyst design.
It allows an explanation of volcano curves, which show the catalytic activity for a given elementary reaction as function of the adsorption energy. As shown in Figure 4, at low values of adsorption energy (more positive values), the reaction is slow because the rate of the reaction is rate-limiting, whereas at high values of adsorption energy, desorption becomes the rate-limiting step. Then, the intermediate values of adsorption energies are required in order to obtain the highest activity (top of the volcano) [42, 58, 59, 60]. The volcano plot in the figure was obtained from a theoretical calculation, where the rate constants for temperature value of 298 K are obtained from Eq. (1) (), the activation energy for the reaction is obtained from Eqs. (2) and (3) () with and equal to 1, and equal to −1.2. The activation energy for desorption process is the negative of .
As mentioned above, a key issue in catalysis is chemisorption, the bond formed between surface and adsorbate can be described in terms of electronic structure of the surface. For transition metals, it is described in terms of d-bands of metals, as is proposed in the “d-band model,” which is presented in the next subsection.
3.2. d-band model
The electronic structure of transition metals involves a broad s-band and narrow d-band. The interaction of the adsorbate valence electrons with s-band electrons can be assumed to be similar for all the different transition metals; hence, the principal differences in the interaction emerge from d-band electrons [4, 5, 43]. The interaction between the adsorbate and surface d-electrons can be represented with molecular orbitals theory, as illustrated in Figure 5, the interaction causes the formation of new shared states, which can be classified into bonding and antibonding states; the antibonding states are higher in energy and its filling leads to the weakening of the bond between adsorbate and surface [4, 5, 43, 61].
If more of antibonding states are empty, then the surface bonds strongly the molecule. Since electrons fill up all states located below the Fermi energy, in order to estimate the strength of the bond is enough to know how many antibonding states are higher in energy than Fermi energy (
The average is a representative value for a group the values or a statistical distribution is the average. In this case, energy average of d-states can be used to indicate the energy of all metal d-states, and it is calculated using Eq. (4) . The difference between this value and Fermi energy can be used to compare in an effective way the d-states energy of different metals. This value is named d-band center (), if it is higher (less negative), then more antibonding states would be above Fermi energy and the bond would be stronger. This is the called d-band model proposed by Hammer and Nørskov [6, 43].
The concepts above presented represent the basic theoretical instruments used to investigate catalytic chemical reactions and to provide satisfactory descriptions with a reduced set of properties compared to complete kinetic description.
3.3. Summarizing the descriptors approach
As remarked in the above sections, the descriptor approach is based on the linking of catalytic activity with a limited set of properties, named descriptors. The main step of this approach is to find a suitable reactivity descriptor and the other is to calculate it. As it was remarked, the adsorption energies of key intermediates can be used as descriptors, the identification of particular adsorption energies that are descriptors can be achieved considering the underlying reaction mechanism in order to identify the rate limiting step, also a Sabatier analysis and the known BEP relations can be used to identify the key adsorption energy. Also, an adsorption energy descriptor can be found through an educated guess, resulting in the knowledge on similar reactions [4, 5, 14].
Since adsorbate-surface interaction is a phenomenon originated by electronic interaction, the electronic structure of the systems can be related with adsorption energies, and some characteristics of this electronic structure can be used as descriptor, one of them is the d-band center. But additional research is required to identify more reliable properties, which can be used in other systems different to transition metals. The next section shows the main results of the modeling with chemical descriptors approach of catalytic carbon-carbon bond cleavage in the ethanol oxidation, reaction important for PEMFCs.
4. Analysis of C-C bond cleavage using chemical descriptors
The catalytic cleavage of the C-C bond in ethanol is fundamental for improving the proton exchange membrane fuel cells fueled with ethanol, whose efficiency is still far from obtained with hydrogen fuel cells, which can be partly attributed to the incomplete oxidation of ethanol to , so energy from ethanol is not totally exploited [62, 63, 64, 65, 66, 67]. Catalysts actually used only promotes partial ethanol oxidation to acetic acid, the limiting step for total oxidation to is, precisely, the C-C bond cleavage. Some works suggest that addition of metals as Ni , Rh , Ir , and Re  to Pt helps in C-C breakage at temperatures below 100°C. So, it is of interest to model this reaction step on Pt-based catalysts in order to clarify the effect of different co-catalysts on the C-C bond cleavage and to understand the conditions under which the C-C bond can be broken. It can be achieved using the descriptor approach above presented.
4.1. Reactivity descriptor for C-C bond cleavage
The C-C cleavage in molecules with more hydrogen atoms attached to the carbon atoms, as ethanol molecule (CH3CH2OH) or acetaldehyde (CH3CHO), has higher activation energies [19, 21, 22]. It is attributed to the high directionality of C-C bond σ-orbital in these molecules, which is constrained along the bond axis, and the several substituents (hydrogen and oxygen atoms) on both carbon atoms constraining sterically this bond, which makes the interaction of the metal with this bond more difficult than with C-H single bond . Probably, breaking of C-H bonds prior to C-C bond cleavage is required, thus the metal could interact with π-orbitals of double C-C bond, which are oriented sideway and have less steric constrains . In fact, it is found that in hydrocarbon reaction on metal surfaces, molecules dehydrogenate before the cleaving of C-C bond .
Experimental works on ethanol oxidation reaction [62, 74, 75, 76] shows that ethanol oxidation involves dehydrogenation reactions, forming acetyl CH3CO; this molecule could be further dehydrogenated forming CH2CO and CHCO. In all the molecules, the C-C bond cleavage is more feasible, so it is worth to focus on them. The reaction for C-C bond cleavage in these CHXCO molecules is written as follows
Linear relationship between products energy and transition state energy has been previously developed for this reaction [19, 22, 53, 77]. So, the activation energy would be related with the interaction of the surface with CHX and CO molecules. It is expected that the trends in the reactivity do not changes with the molecule CHXCO chosen, since the adsorption energy of the CHX molecules is expected to be similar , that is, if rhodium has higher adsorption of CH than Ru, it would have higher adsorption energy of CH2 and CH3 than Ru. So, in this work, to study the C-C bond cleavage in the simplest molecule CHCO, also this step, according to previous works [19, 22], has the lowest activation energy for the C-C bond cleavage on different transition metals. The linear relationship , which will be used, is shown in Eq. (6).
and are, respectively, the transition state and product energies respect to the energy of reactants in gas phase and the clean slab. So, is the energy involved to go from CHCO and the surface far apart to the CH and CO adsorbed on the surface that is expressed in Eq. (7)
Activation energy is expressed as in Eq. (2), in this case, is the energy of adsorbed respect to the energy of reactants (CHCO) in gas phase and the clean slab, that is, the adsorption energy of CHCO. So, the final expression for activation energy, which states the relationship between catalytic activity and adsorption energy, is:
Considering that CO2, which desorbs from surface easily, can be readily obtained from CO oxidation by OH, which is available in alkaline media and in acidic media can be obtained from water activation facilitated by oxyphilic atoms as Sn, Ru, and Mo [20, 75, 78, 79]. So, it seems reasonable to expect that in the appropriated conditions (alkaline media or oxyphilic metal as co-catalyst), the surfaces with higher CO formation also have higher CO2 formation. Therefore, the study of the C-C bond cleavage step can give enough information about which catalysts can helps to obtain higher CO2. The reactivity of the surfaces for C-C bond cleavage is studied through DFT modeling of adsorption of CH, CO, and CHCO.
The adsorption energies of these molecules were calculated for different Pt3M alloys (M = Sn, Re, Ru, Rh, and Ni). Since the alloys with this composition or lower content of M are commonly proposed as a convenient composition to improve the catalytic activity to ethanol oxidation in experimental studies [68, 69, 71, 80, 81, 82].
4.2. Calculation of reactivity descriptors
The adsorption energies of CH, CO, and CHCO are calculated using the slab model to represent the (111) facet of the face-centered cubic metals (Pt3M) and DFT as is implemented in Quantum Espresso . It used plane waves, exchange correlation energy described by Perdew and Wang (PW91) , and ultra-soft pseudopotentials . Convergence was tested with respect to total energy for the Brillouin sampling and for the kinetic energy cut-off in the plane wave expansion, so 3 × 3 × 1 Monkhorst-Pack k-mesh  and the kinetic energy cut-off for the plane waves of 544 eV was used.
For slab construction, the lattice parameters obtained from geometric optimization with DFT was used. The effect of the size of the supercell has also been carefully tested for the adsorption energy. Increasing the supercell size from 2 × 2, which we have used, to 4 × 4 changes the adsorption energy by only 0.008 eV, so the description with 2 × 2 is considered good enough considering the computational cost of the simulation of large systems. Five layers of atoms were employed for representing the portion of solid and empty space of 9.00 Å was left above of the atomic surface to avoid self-interaction of the slab. With this empty space, the adsorption energies changed just 0.02 eV compared to the adsorption energies obtained using an empty space of 15.00 Å, indicating that the top of one slab has essentially no effect on the bottom of the next.
The adsorption energy was calculated as total energy of the slab with the i-species adsorbed minus the total energy of clean slab and the i-species in gas phase .
Adsorption of , , and were tested on several sites, atop sites (A-Pt and A-M), bridge sites (B-Pt and B-M), and hollow sites with and without an atom directly below in the subsequent layer, which are called HCP and FCC hollow sites (HCP-Pt, HCP-M, FCC-Pt, and FCC-M).
The most stable adsorption configuration for CO adsorption on Pt3M is at atop site with molecular axis of CO normal to the surface and C atom bonded to the surface in agreement with the experimental findings for Pt-based catalysts . The CO adsorption energies on A-Pt and A-M for the different surfaces are shown in Figure 6. It is found that CO adsorption energy on A-Pt sites in the alloys is reduced respect to pure Pt (dotted line) and it is more reduced when the metal atom adsorbs more strongly the CO molecule. Also it is remarkable that adsorption energy on tin atom is only −0.08 eV. It is in agreement with other works, which points out that Pt-CO binding is weakened in the PtM alloy [88, 89, 90], whereas the M-CO bond is strengthened with the presence of Pt and CO adsorption on tin is not stable [91, 92, 93].
The most stable adsorption configuration for CH adsorption on Pt3M is at hollow site (both hollow sites have similar adsorption energy) with molecular axis of CH normal to the surface and C atom bonded to an atom on the surface in agreement with other theoretical studies [94, 95, 96]. It is noted in Figure 6, a trend similar to CO adsorption, adsorption on sites with only Pt in the alloy is weaker than pure Pt, whereas the sites with metal atoms have higher or similar adsorption energy than pure Pt with the exception of Pt3Sn, which is significantly lower.
The most stable adsorption configuration for CHCO on the alloys (except for Pt3Sn) and pure Pt is shown in
Adsorption energies of former configuration with CH and CO parts bonded to the surface ranges from −3.66 to −3.98 and C-C bond length is stretched by 0.20 ± 0.01 Ǻ, whereas adsorption energies of last configuration ranges from −2.68 to −3.24 eV and C-C bond length is only stretched slightly by 0.08 ± 0.01 Ǻ. So, the ensembles with only Pt atoms in Pt3M (with exception of Pt3Rh) have adsorption energies significantly lower than pure Pt and do not stretch significantly the C-C bond, whereas the ensembles with M atom have adsorption energies similar to Pt and stretch significantly the C-C bond. Indicating that the addition of M to Pt decreases the ability of Pt for interacting with CHCO and hence part of the alloy is inactive, except for Pt3Rh.
Summarizing, the adsorption of CH, CO, and CHCO shows that the addition of M (M = Re, Rh, Ni, Ru, or Sn) to Pt leads to obtain two active sites with different reactivities. The Pt atoms in all these alloys form bonds with the described molecules weaker than Pt atoms in pure Pt, whereas the Re, Rh, and Ru atoms in the alloys have a stronger interaction with CH, CO, and CHCO than pure Pt. This difference explains that adsorption energies of CO on metal atoms in Pt3M (M = Rh, Re, and Ru) are higher than Pt (111), whereas the adsorption energies of CH on sites with M atoms in Pt3M (M = Rh, Re, and Ru) are similar to Pt (111), since CH interacts with two Pt and one M atom at hollow site. Then, the stronger interaction of mentioned M atoms should compensate the weaker interaction of Pt in these alloys.
Besides, Ni and Pt atoms in Pt3Ni have slightly weaker interaction than pure Pt (111), so the adsorption energies at different sites in this alloy are slightly lower than pure Pt. On the other hand, Sn and Pt atoms in Pt3Sn have a weaker interaction with the molecules than pure Pt (111), suggesting that Pt3Sn has low catalytic activity toward C-C bond cleavage, confirming previous experimental and theoretical works that affirms that PtSn alloys have lower catalytic activity for C-C bond cleavage, but they increase the catalytic activity for EOR by promoting bifunctional mechanism and the pathway toward acetic acid [70, 98]. Other remarkable finding is the two configurations for CHCO adsorption, one of them destabilizes the molecule binding it strongly and stretching significantly the C-C bond (around to 0.20 Ǻ) and the other is found in metals inactive, it suggests that part of alloy is inactive, except for Pt3Rh.
The adsorption energies obtained can be related with activation energy for C-C bond cleavage through Eq. (8), only the active sites are considered since the linear relationship between ETS and Eads only have sense for them. As is can be expected from Eq. (8), the surfaces with higher adsorption energies of CH and CO have lower activation energy, that is the sites with Re, Rh, and Ru atoms, it can be observed in Figure 8, which shows the product energies Ep (see Eq. (7)) and activation energies for the different active sites.
These results indicate that Pt3Sn is the catalyst with the lowest activity for C-C bond cleavage, since in this alloy Sn has a weak interaction with all molecules studied and Pt atoms have a weaker interaction than Pt atoms in pure Pt (111). Also, it is expected that Pt3Rh alloy has higher reactivity for C-C bond cleavage since all sites, sites with only Pt, and sites with Rh atoms have similar or higher reactivity than pure Pt. Whereas Pt3Re and Pt3Ru have a nonsignificant improvement of the reactivity with respect to pure Pt, since they have one site with higher reactivity (site with M atom) and other with lesser reactivity (site with only Pt atoms); between these ones, Pt3Re has the most reactive site. Pt3Ni has slightly weaker interaction with molecules than Pt. Consequently, the ranking from high to low reactivity for C-C bond cleavage is Pt3Rh, Pt3Re, Pt3Ru, Pt, Pt3Ni, and Pt3Sn.
This trend agrees with previous experimental results, where it is reported that Pt-based alloys with rhodium (PtRh and PtSnRh) show a production of CO2 higher than Pt and Pt alloys [69, 99]. In contrast, addition of Rhenium to Pt improves performance for C-C cleavage, but this improvement is not much higher, it is reported that Re addition to Pt improves the acetaldehyde decomposition, increasing the produced CO (product of C-C bond cleavage) from 10 to 18% . On the contrary, it is reported that Sn addition to Pt decreases the CO2 product respect to Pt [70, 98]. But the PtSn alloys are recognized as the better catalyst for EOR, where the enhancement of the activity with this alloy is attributed to the promotion of bifunctional mechanism promoting the acetic acid production.
Due to mentioned issues, different authors propose materials with Pt, Sn, and Rh as effective catalysts in splitting the C-C bond and oxidizing intermediates molecules of EOR via bifunctional mechanism. This mix in the appropriate composition will provide high overall conversion of ethanol and higher selectivity to CO2 than Pt. However, the selectivity can be considered low and the total conversion to CO2 remains elusive. It can be attributed to the presence of oxidants as OH, which should promote the oxidation of intermediates to acetic acid and decreases the C-C bond cleavage rate . But the unavailability of oxidants limits the conversion of CO to CO2, poisoning the surface. Finding a better catalyst for EOR requires considering, in addition of C-C cleavage, the step of oxidation to produce acetic acid and the dehydrogenation step.
Kaxiras E. Atomic and Electronic Structure of Solids. New York, United States of America, Cambridge: Cambridge University Press; 2003
Callister WD, Rethwisch DG. Materials Science and Engineering: An Introduction. 7th ed. John Wiley & Sons, Inc; 2007
Seh ZW, Kibsgaard J, Dickens CF, Chorkendorff I, Nørskov JK, Jaramillo TF. Combining theory and experiment in electrocatalysis: Insights into materials design. Science. 2017; 355(6321):eaad4998
Nørskov JK, Bligaard T, Rossmeisl J, Christensen CH. Towards the computational design of solid catalysts. Nature Chemistry. 2009; 1(1):37-46
Nørskov JK, Abild-Pedersen F, Studt F, Bligaard T. Density functional theory in surface chemistry and catalysis. Proceedings of the National Academy of Sciences of the United States of America. 2011; 108(3):937-943
Hammer B, Nørskov JK. Electronic factors determining the reactivity of metal surfaces. Surface Science. 1995; 343(3):211-220
Linic S, Jankowiak J, Barteau MA. Selectivity driven design of bimetallic ethylene epoxidation catalysts from first principles. Journal of Catalysis. 2004; 224(2):489-493
Greeley J, Mavrikakis M. Alloy catalysts designed from first principles. Nature Materials. 2004; 3(11):810-815
Greeley J, Jaramillo TF, Bonde J, Chorkendorff I, Nørskov JK. Computational high-throughput screening of electrocatalytic materials for hydrogen evolution. Nature Materials. 2006; 5(11):909-913
Andersson M, Bligaard T, Kustov A, Larsen K, Greeley J, Johannessen T, et al. Toward computational screening in heterogeneous catalysis: Pareto-optimal methanation catalysts. Journal of Catalysis. 2006; 239(2):501-506
Strasser P, Fan Q, Devenney M, Weinberg WH, Liu P, Nørskov JK. High throughput experimental and theoretical predictive screening of materials–A comparative study of search strategies for new fuel cell anode catalysts. The Journal of Physical Chemistry. B. 2003; 107(40):11013-11021
Grabow LC, Hvolbæk B, Nørskov JK. Understanding trends in catalytic activity: The effect of adsorbate-adsorbate interactions for CO oxidation over transition metals. Topics in Catalysis. 2010; 53(5-6):298-310
Suntivich J, Gasteiger HA, Yabuuchi N, Nakanishi H, Goodenough JB, Shao-Horn Y. Design principles for oxygen-reduction activity on perovskite oxide catalysts for fuel cells and metal–air batteries. Nature Chemistry. 2011; 3(8):647-647
Grabow LC. Computational catalyst screening. In: Asthagiri A, Janik MJ, editors. Computational Catalysis. Cambridge, UK: Royal Society of Chemistry; 2013. pp. 1-58
Jónsson H, Mills G, Jacobsen KW. Nudged elastic band method for finding minimum energy paths of transitions. In: Classical and Quantum Dynamics in Condensed Phase Simulations. Singapore: World Scientific; 1998. pp. 385-404
Henkelman G, Uberuaga BP, Jónsson H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. The Journal of Chemical Physics. 2000; 113(22):9901
Henkelman G, Jónsson H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. The Journal of Chemical Physics. 2000; 113(22):9978-9985
Peng C, Bernhard Schlegel H. Combining synchronous transit and quasi-Newton methods to find transition states. Israel Journal of Chemistry. 1993; 33(4):449-454
Alcalá R, Mavrikakis M, Dumesic JA. DFT studies for cleavage of C-C and C-O bonds in surface species derived from ethanol on Pt(111). Journal of Catalysis. 2003; 218:178-190
Wang H, Liu Z. Comprehensive mechanism and structure-sensitivity of ethanol oxidation on platinum: New transition-state searching method for resolving the complex reaction. Network. 2008; 10:10996-11004
Choi Y, Liu P. Understanding of ethanol decomposition on Rh(111) from density functional theory and kinetic Monte Carlo simulations. Catalysis Today. 2011; 165(1):64-70
Li M, Guo W, Jiang R, Zhao L, Shan H. Decomposition of ethanol on Pd(111): A density functional theory study. Langmuir. 2010; 26(3):1879-1888
Ishimoto T, Kazuno H, Kishida T, Koyama M. Theoretical study on oxidation reaction mechanism on Au catalyst in direct alkaline fuel cell. Solid State Ionics. 2014; 262:328-331
Xu Z-F, Wang Y. Effects of alloyed metal on the catalysis activity of Pt for ethanol partial oxidation: Adsorption and dehydrogenation on Pt(3)M (M=Pt, Ru, Sn, Re, Rh, and Pd). The Journal of Physical Chemistry. C, Nanomaterials and Interfaces. 2011; 115(42):20565-20571
Pérez-Mendoza AE, Castañeda-Ramírez S, Ribadeneira RE, Perez Mendoza AE, Castaneda Ramirez S, Ribadeneira Paz RE. Multiscale modeling of a proton exchange membrane fuel cell: Atomistic oxygen reduction reaction model. ECS Transactions. 2015; 66(28):7-18
Sidik RA, Anderson AB. Density functional theory study of O2 electroreduction when bonded to a Pt dual site. Journal of Electroanalytical Chemistry. 2002; 528(1-2):69-76
JA K, Jerkiewicz G, Jacob T. Theoretical investigations of the oxygen reduction reaction on Pt(111). Chemphyschem. 2010; 11(13):2779-2794
Eichler A, Mittendorfer F, Hafner J. Precursor-mediated adsorption of oxygen on the (111) surfaces of platinum-group metals. Physical Review B. 2000; 62(7):4744-4755
Eichler A, Hafner J. Molecular precursors in the dissociative adsorption of O2 on Pt(111). Physical Review Letters. 1997; 79(22):4481-4484
Nørskov JK, Rossmeisl J, Logadottir A, Lindqvist L, Lyngby D, Jo H. Origin of the overpotential for oxygen reduction at a fuel-cell cathode. The Journal of Physical Chemistry. B. 2004: 108;17886-17892
Walch S, Dhanda A, Aryanpour M, Pitsch H. Mechanism of molecular oxygen reduction at the cathode of a PEM fuel cell: Non-electrochemical reactions on catalytic Pt particles. Journal of Physical Chemistry C. 2008; 112(22):8464-8475
Balbuena PB, Altomare D, Agapito L, Seminario JM. Theoretical analysis of oxygen adsorption on Pt-based clusters alloyed with Co, Ni, or Cr embedded in a Pt matrix. The Journal of Physical Chemistry. B. 2003; 107(49):13671-13680
Panchenko A, Koper MTM, Shubina TE, Mitchell SJ, Roduner E. Ab initio calculations of intermediates of oxygen reduction on low-index platinum surfaces. Journal of the Electrochemical Society. 2004; 151(12):A2016
Jinnouchi R, Okazaki K. New insight into microscale transport phenomena in pefc by quantum md. Microscale Thermophysical Engineering. 2003; 7(1):15-31
Liu S, White MG, Liu P. Mechanism of oxygen reduction reaction on Pt(111) in alkaline solution: Importance of chemisorbed water on surface. Journal of Physical Chemistry C. 2016; 120(28):15288-15298
Sha Y, Yu TH, Liu Y, Merinov BV, Goddard WA. Theoretical study of solvent effects on the platinum-catalyzed oxygen reduction reaction. Journal of Physical Chemistry Letters. 2010; 1(5):856-861
Jinnouchi R, Kodama K, Hatanaka T, Morimoto Y. First principles based mean field model for oxygen reduction reaction. Physical Chemistry Chemical Physics. 2011; 13(47):21070
Mathew K, Sundararaman R, Letchworth-Weaver K, Arias TA, Hennig RG. Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways. The Journal of Chemical Physics. 2014; 140(8):084106-1-084106-8
Greeley J, Stephens IEL, Bondarenko AS, Johansson TP, Hansen HA, Jaramillo TF, et al. Alloys of platinum and early transition metals as oxygen reduction electrocatalysts. Nature Chemistry. 2009; 1(7):552-556
Ou L. Design of Pd-based bimetallic catalysts for ORR: A DFT calculation study. Journal of Chemistry. 2015; 2015:1-11
Greeley J, Mavrikakis M. Near-surface alloys for hydrogen fuel cell applications. Catalysis Today. 2006; 111(1-2):52-58
Knozinger H, Deutschmann O. Heterogeneous catalysis and solid catalysts. In: Ullmann’s Encyclopedia of Industrial Chemistry. Weinheim, Germany: Wiley-VCH Verlag GmbH & Co. KGaA; 2009
Hammer B, Nørskov JK. Theoretical surface science and catalysis: Calculations and concepts. Advances in Catalysis. 2000; 45
Stamenkovic V, Mun BS, Mayrhofer KJJ, Ross PN, Markovic NM, Rossmeisl J, et al. Changing the activity of electrocatalysts for oxygen reduction by tuning the surface electronic structure. Angewandte Chemie, International Edition. 2006; 45(18):2897-2901
Stamenkovic VR, Mun BS, Arenz M, Mayrhofer KJJ, Lucas CA, Wang G, et al. Trends in electrocatalysis on extended and nanoscale Pt-bimetallic alloy surfaces. Nature Materials. 2007; 6(3):241-247
Bandarenka AS, Koper MTM. Structural and electronic effects in heterogeneous electrocatalysis: Toward a rational design of electrocatalysts. Journal of Catalysis. 2013; 308:11-24
Bligaard T, Nørskov JK. Ligand effects in heterogeneous catalysis and electrochemistry. Electrochimica Acta. 2007; 52(18):5512-5516
Lile JRD, Zhou S. Theoretical modeling of the PEMFC catalyst layer: A review of atomistic methods. Electrochimica Acta. 2015; 177:4-20
Shi Z. Application of first principles methods in the study of fuel cell air-cathode electrocatalysis. In: Zhang J, editor. PEM Fuel Cell Electrocatalysts and Catalyst Layers. London: Springer; 2008. p. 289-329
Groß A, Gross A. Theoretical Surface Science. A Microscopic Perspective. Berlin, Heidelberg: Springer Berlin Heidelberg; 2009
Norskov JK, Studt F, Abild-Pedersen F, Bligaard T. Fundamental Concepts in Heterogeneous Catalysis. Hoboken, United States of America: John Wiley & Sons, Inc; 2014. 208 p
Eyring H. The theory of absolute reaction rates. Transactions of the Faraday Society. 1938; 34:41-48
Wang S, Temel B, Shen J, Jones G, Grabow LC, Studt F, et al. Universal Brønsted-Evans-Polanyi relations for C–C, C–O, C–N, N–O, N–N, and O–O dissociation reactions. Catalysis Letters. 2011; 141(3):370-373
Wang S, Vorotnikov V, Sutton JE, Vlachos DG. Brønsted–Evans–Polanyi and transition state scaling relations of furan derivatives on Pd(111) and their relation to those of small molecules. ACS Catalysis. 2014; 4(2):604-612
Nørskov JK. Universality in heterogeneous catalysis. Journal of Catalysis. 2002; 209(2):275-278
Michaelides A, Liu ZP, Zhang CJ, Alavi A, King DA, Hu P. Identification of general linear relationships between activation energies and enthalpy changes for dissociation reactions at surfaces. Journal of the American Chemical Society. 2003; 125(13):3704-3705
Evans MG, Polanyi M. Inertia and driving force of chemical reactions. Transactions of the Faraday Society. 1938; 34:11
Bligaard T, Nørskov JKK, Dahl S, Matthiesen J, Christensen CHH, Sehested J. The Brønsted-Evans-Polanyi relation and the volcano curve in heterogeneous catalysis. Journal of Catalysis. 2004; 224(1):206-217
Nørskov JK, Bligaard T, Hvolbaek B, Abild-Pedersen F, Chorkendorff I, Christensen CH. The nature of the active site in heterogeneous metal catalysis. Chemical Society Reviews. 2008; 37(10):2163-2171
Logadottir A, Rod TH, Nørskov JK, Hammer B, Dahl S, Jacobsen CJH. The Brønsted–Evans–Polanyi relation and the volcano plot for ammonia synthesis over transition metal catalysts. Journal of Catalysis. 2001; 197(2):229-231
Nilsson A, Pettersson LGM, Hammer B, Bligaard T, Christensen CH, Norskov JK. The electronic structure effect in heterogeneous catalysis. Catalysis Letters. 2005; 100(3-4):111-114
García-rodríguez S, Herranz T, Rojas S. New and future developments in catalysis. In: Electroctalysts for the Electrooxidation of Ethanol. Amsterdam, The Netherlands: Elsevier; 2013. pp. 33-67
Lamy C. The direct ethanol fuel cell: A challenge to convert bioethanol cleanly into electric energy. In: Barbaro P, Bianchini C, editors. Catalysis for Sustainable Energy Production. Weinheim, Germany: WILEY-VCH Verlag GmbH & Co. KGaA; 2009. pp. 1-46
Kutz RB, Braunschweig B, Mukherjee P, Behrens RL, Dlott DD, Wieckowski A. Reaction pathways of ethanol electrooxidation on polycrystalline platinum catalysts in acidic electrolytes. Journal of Catalysis. 2011; 278(2):181-188
Badwal SPS, Giddey S, Kulkarni A, Goel J, Basu S. Direct ethanol fuel cells for transport and stationary applications – A comprehensive review. Applied EnergyAppl Energy. 2015; 145:80-103
Zhao TS, Li YS, Shen SY. Anion-exchange membrane direct ethanol fuel cells: Status and perspective. Frontiers of Energy and Power Engineering in China. 2010; 4(4):443-458
Kamarudin MZF, Kamarudin SK, Masdar MS, Daud WRW. Review: Direct ethanol fuel cells. International Journal of Hydrogen Energy. 2012; 38:1-16
Beyhan S, Coutanceau C, Léger JM, Napporn TW, Kadirgan F. Promising anode candidates for direct ethanol fuel cell: Carbon supported PtSn-based trimetallic catalysts prepared by Bönnemann method. International Journal of Hydrogen Energy. 2013; 38:6830-6841
Kowal A, Li M, Shao M, Sasaki K, Vukmirovic MB, Zhang J, et al. Ternary Pt/Rh/SnO2 electrocatalysts for oxidizing ethanol to CO2. Nature Materials. 2009; 8(4):325-330
Li M, Cullen DA, Sasaki K, Marinkovic NS, More K, Adzic RR. Ternary electrocatalysts for oxidizing ethanol to carbon dioxide: Making ir capable of splitting C-C bond. Journal of the American Chemical Society. 2013; 135(1):132-141
Tayal J, Rawat B, Basu S. Effect of addition of rhenium to Pt-based anode catalysts in electro-oxidation of ethanol in direct ethanol PEM fuel cell. International Journal of Hydrogen Energy. 2012; 37(5):4597-4605
Murakami M, Ito Y. Cleavage of carbon—carbon single bonds by transition metals. In: Activation of Unreactive Bonds and Organic Synthesis. Berlin Heidelberg: Springer; 1999
Mueller JE, Van Duin ACT, Goddard WA. Application of the ReaxFF reactive force field to reactive dynamics of hydrocarbon chemisorption and decomposition. Journal of Physical Chemistry C. 2010; 114(12):5675-5685
Shao MH, Adzic RR. Electrooxidation of ethanol on a Pt electrode in acid solutions: In situ ATR-SEIRAS study. Electrochimica Acta. 2005; 50(12):2415-2422
Vigier F, Coutanceau C, Hahn F, Belgsir EMM, Lamy C. On the mechanism of ethanol electro-oxidation on Pt and PtSn catalysts: Electrochemical and in situ IR reflectance spectroscopy studies. Journal of Electroanalytical Chemistry. 2004; 563(1):81-89
Shen SY, Zhao TS, Wu QX. Product analysis of the ethanol oxidation reaction on palladium-based catalysts in an anion-exchange membrane fuel cell environment. International Journal of Hydrogen Energy. 2012; 37(1):575-582
Ferrin P, Simonetti D, Kandoi S, Kunkes E, Dumesic JA, Nørskov JK, et al. Modeling ethanol decomposition on transition metals: A combined application of scaling and Brønsted-Evans-Polanyi relations. Journal of the American Chemical Society. 2009; 131(16):5809-5815
Antolini E, Gonzalez ERR. Effect of synthesis method and structural characteristics of Pt-Sn fuel cell catalysts on the electro-oxidation of CH3OH and CH3CH2OH in acid medium. Catalysis Today. 2011; 160(1):28-38
Gong X-Q, Hu P, Raval R. The catalytic role of water in CO oxidation. The Journal of Chemical Physics. 2003; 119(12):6324-6334
Guo Y, Zheng Y, Huang M. Enhanced activity of PtSn/C anodic electrocatalyst prepared by formic acid reduction for direct ethanol fuel cells. Electrochimica Acta. 2008; 53:3102-3108
Li H, Sun G, Cao L, Jiang L, Xin Q. Comparison of different promotion effect of PtRu/C and PtSn/C electrocatalysts for ethanol electro-oxidation. Electrochimica Acta. 2007; 52(24):6622-6629
Vigier F, Coutanceau C, Perrard A, Belgsir EMM, Lamy C. Development of anode catalysts for a direct ethanol fuel cell. Journal of Applied Electrochemistry. 2004; 34(4):439-446
Giannozzi P, Baroni S. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. Journal of Physics: Condensed Matter. 2009; 36:395502-1-395502-36
Perdew JP. Unified theory of exchange and correlation beyond the local density approximation. In: Ziesche P, Eschrig H, editors. Electronic Structure of Solids ‘91. Berlin: Akademie Verlag; 1991. pp. 11-20
Vanderbilt D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Physical Review B. 1990; 41(11):7892-7895
Pack JD, Monkhorst HJ. “Special points for Brillouin-zone integrations”-a reply. Physical Review B. 1977; 16(4):1748-1749
Van Hove MA. The adsorption sites of CO and NO molecules on metal surfaces. Israel Journal of Chemistry. 1998; 38(4):349-352
Liao MS, Cabrera CR, Ishikawa Y. Theoretical study of CO adsorption on Pt, Ru and Pt-M (M = Ru, Sn, Ge) clusters. Surface Science. 2000; 445(2-3):267-282
Su HY, Bao XH, Li WX. Modulating the reactivity of Ni-containing Pt(111)-skin catalysts by density functional theory calculations. The Journal of Chemical Physics. 2008; 128(19):1-7
Ramstad A, Strisland F, Raaen S, Borg A, Berg C. CO and O2 adsorption on the Re/Pt(111) surface studied by photoemission and thermal desorption. Surface Science. 1999; 440(1-2):290-300
Shubina TE, Koper MTM. Quantum-chemical calculations of CO and OH interacting with bimetallic surfaces. Electrochimica Acta. 2002; 47(22-23):3621-3628
Dupont C, Jugnet Y, Loffreda D. Theoretical evidence of PtSn alloy efficiency for CO oxidation. Journal of the American Chemical Society. 2006; 128(28):9129-9136
Liu P, Logadottir A, Nørskov JK. Modeling the electro-oxidation of CO and H2/CO on Pt, Ru, PtRu and Pt3Sn. Electrochimica Acta. 2003; 48(25-26):3731-3742
Kua J, Faglioni F, W a G. Thermochemistry for hydrocarbon intermediates chemisorbed on metal surfaces: CH n- m(CH3) mwith n= 1, 2, 3 on Pt, Ir, Os, Pd, Rh, and Ru. Journal of the American Chemical Society. 2000; 122(5):2309-2321
Hahn K, Mavrikakis M. Atomic and molecular adsorption on Re(0001). Topics in Catalysis. 2014; 57(1-4):54-68
Zhang M, Yang K, Zhang X, Yu Y. Effect of Ni(111) surface alloying by Pt on partial oxidation of methane to syngas: A DFT study. Surface Science. 2014; 630:236-243
Courtois J, Du W, Wong E, Teng X, Deskins NA. Screening iridium-based bimetallic alloys as catalysts for direct ethanol fuel cells. Applied Catalysis A: General. 2014; 483:85-96
Alcala R, Shabaker JW, Huber GW, MA S-C, JA D. Experimental and DFT studies of the conversion of ethanol and acetic acid on PtSn-based catalysts. The Journal of Physical Chemistry. B. 2005; 109(6):2074-2085
Bach Delpeuch A, Maillard F, Chatenet M, Soudant P, Cremers C. Ethanol oxidation reaction (EOR) investigation on Pt/C, Rh/C, and Pt-based bi- and tri-metallic electrocatalysts: A DEMS and in situ FTIR study. Applied Catalysis B: Environmental. 2016; 181:672-680
Ciftci A, Ligthart DAJM, Sen AO, van Hoof AJF, Friedrich H, Hensen EJM. Pt-Re synergy in aqueous-phase reforming of glycerol and the water-gas shift reaction. Journal of Catalysis. 2014; 311:88-101
Kavanagh R, Cao XM, Lin WF, Hardacre C, Hu P. Origin of low CO2 selectivity on platinum in the direct ethanol fuel cell. Angewandte Chemie, International Edition. 2012; 51(7):1572-1575