Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Spectrophotometric Study of Charge-Transfer Complexes of Ruxolitinib with Chloranilic Acid and 2,3-Dichloro-5,6-dicyano-1,4-benzoquinone: An Application to the Development of a Green and High-Throughput Microwell Method for Quantification of Ruxolitinib in Its Pharmaceutical Formulations
Next Article in Special Issue
Electronic Properties of Graphene Nano-Parallelograms: A Thermally Assisted Occupation DFT Computational Study
Previous Article in Journal
Effect of Inlet Air Temperature and Quinoa Starch/Gum Arabic Ratio on Nanoencapsulation of Bioactive Compounds from Andean Potato Cultivars by Spray-Drying
Previous Article in Special Issue
Real-Time Extension of TAO-DFT
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 23
10.3390/molecules28237873
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Density Functional Theory Investigation of Temperature-Dependent Properties of Cu-Nitrogen-Doped Graphene as a Cathode Material in Fuel Cell Applications

by
Yashas Balasooriya
1,
Pubudu Samarasekara
2,
Chee Ming Lim
3,
Yuan-Fong Chou Chau
3,
Muhammad Raziq Rahimi Kooh
3,* and
Roshan Thotagamuge
4,*
1
Postgraduate Institute of Science, University of Peradeniya, Peradeniya 20400, Sri Lanka
2
Department of Physics, Faculty of Science, University of Peradeniya, Peradeniya 20400, Sri Lanka
3
Centre for Advanced Material and Energy Sciences, Universiti Brunei Darussalam, Jalan Tungku Link, Gadong BE 1410, Brunei
4
Department of Nano Science Technology, Faculty of Technology, Wayamba University of Sri Lanka, Kuliyapitiya 60200, Sri Lanka
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(23), 7873; https://doi.org/10.3390/molecules28237873
Submission received: 7 November 2023 / Revised: 24 November 2023 / Accepted: 28 November 2023 / Published: 30 November 2023
(This article belongs to the Special Issue Multiconfigurational and DFT Methods Applied to Chemical Systems)
Browse Figures
Graphical abstract
"> Figure 1
<p>Electrostatic potential diagram of the defect structure at all temperatures (red and blue areas represent the negative and positive charges, respectively, and the increasing density of the color represents the increasing charge density).</p> "> Figure 2
<p>Electrostatic potential diagram of the OH + H<sub>2</sub>O intermediate at the temperatures of (<b>a</b>) 298.15 K (<b>b</b>) 353.15 K, and (<b>c</b>) 393.15 K (red and blue areas represent the negative and positive charges, respectively, and increasing density of color represents increasing charge density).</p> "> Figure 3
<p>Top and side displacements of the *OH + H + *OH intermediate around the 289.32 cm<sup>−1</sup> frequency according to vibrational data at 353.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).</p> "> Figure 4
<p>IR spectra of different temperatures for (<b>a</b>) the catalyst (*), (<b>b</b>) *O<sub>2</sub>, (<b>c</b>) OH + H<sub>2</sub>O, and (<b>d</b>) H<sub>2</sub>O + H<sub>2</sub>O.</p> "> Figure 5
<p>Optimized intermediate of OH + H<sub>2</sub>O (*OH + H + *OH) at the temperatures of (<b>a</b>) 298.15 K (<b>b</b>) 353.15 K, and (<b>c</b>) 393.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).</p> "> Figure 6
<p>Optimization process when the H atom is introduced to the first O atom of the OOH intermediate at a temperature of 298.15 K. (<b>a</b>) Initially, the H atom binds to the first O atom of OOH to form H<sub>2</sub>O<sub>2</sub>, (<b>b</b>) in the middle state of optimization, and (<b>c</b>) in the final optimization step where dual OH molecules are formed instead of H<sub>2</sub>O<sub>2</sub> (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).</p> "> Figure 7
<p>Optimization process when a H atom is introduced to the second O atom of the OOH intermediate at the temperature of 298.15 K: (<b>a</b>) initially, the H atom binds with the second O atom of the OOH intermediate to form H<sub>2</sub>O; (<b>b</b>) in the next step, the H atom binds to the second O atom, causing the previously bonded H atom to break its bond and migrate to the first O atom; and (<b>c</b>) in the final optimization step, a dual OH molecule is formed instead of H<sub>2</sub>O (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).</p> "> Figure 8
<p>Top and side view of the optimized ORR steps for the Cu<sub>2</sub>-N<sub>8</sub>/Gr structure at 298.15 K. (<b>a</b>) O<sub>2</sub>, (<b>b</b>) OOH, (<b>c</b>) OH + OH, (<b>d</b>) OH + H + OH, and (<b>e</b>) 2H<sub>2</sub>O molecules binding to Cu<sub>2</sub>-N<sub>8</sub>/Gr (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).</p> "> Figure 9
<p>Free energy diagrams of Cu<sub>2</sub>-N<sub>8</sub>/Gr at different temperatures under various cell potentials in an acidic medium: (<b>a</b>) 298.15 K, (<b>b</b>) 353.15 K, and (<b>c</b>) 393.15 K (note: the OH + H<sub>2</sub>O intermediate is the same as the OH + H +OH intermediate, as explained previously).</p> ">
Versions Notes

Abstract

:
In this study, density functional theory (DFT) was used to investigate the influence of temperature on the performance of a novel Cu-nitrogen-doped graphene Cu2-N8/Gr nanocomposite as a catalyst for the oxygen reduction reaction (ORR) in fuel cell applications. Our DFT calculations, conducted using Gaussian 09w with the 3–21G/B3LYP basis set, focus on the Cu-nitrogen-doped graphene nanocomposite cathode catalyst, exploring its behavior at three distinct temperatures: 298.15 K, 353.15 K, and 393.15 K, under acidic conditions. Our analysis of formation energies indicates that the structural stability of the catalyst remains unaffected as the temperature varies within the potential range of 0–7.21 V. Notably, the stability of the ORR steps experiences a marginal decrease with increasing temperature, with the exception of the intermediate OH + H2O (*OH + H + *OH). Interestingly, the optimization reveals the absence of single OH and H2O intermediates during the reactions. Furthermore, the OH + H2O step is optimized to form the OH + H + OH intermediate, featuring the sharing of a hydrogen atom between dual OH intermediates. Free energy calculations elucidate that the catalyst supports spontaneous ORR at all temperatures. The highest recorded maximum cell potential, 0.69 V, is observed at 393.15 K, while the lowest, 0.61 V, is recorded at 353.15 K. In particular, the Cu2-N8/Gr catalyst structure demonstrates a reduced favorability for the H2O2 generation at all temperatures, resulting in the formation of dual OH intermediates rather than H2O2. In conclusion, at 393.15 K, Cu2-N8/Gr exhibits enhanced catalyst performance compared to 353.15 K and 298.15 K, making it a promising candidate for ORR catalysis in fuel cell applications.

Graphical Abstract">

Graphical Abstract

1. Introduction

Fuel cells operating at low temperatures have attracted considerable attention due to their enhanced durability, reliability, and efficiency with low weight, making them more portable compared to their counterparts. These cells typically operate at temperatures below 473.15 K, often around 373.15 K, and commonly rely on carbon-supported catalysts to enhance the oxygen reduction reaction (ORR). However, these carbon supports are susceptible to rapid corrosion at higher temperatures, resulting in catalyst agglomeration, material durability, and subsequent reductions in efficiency [1].
Walch et al. revealed the effectiveness of the DFT method in studying ORR performance and catalyst properties [2]. Commercially available low-temperature fuel cells predominantly employ platinum group metals (PGMs), such as platinum (Pt), palladium (Pd), and iridium (Ir), as catalysts. These costly materials effectively facilitate hydrogen ionization at the anode [3,4] and provide an efficient ORR surface at the cathode [5]. However, as the carbon support undergoes oxidation, the catalyst layered upon it tends to agglomerate, resulting in reduced fuel cell efficiency over short operational cycles. Catalyst materials made of transition-metal-based perovskite oxides have also been used as effective catalysts compared to PGM. However, the process of creating those materials is relatively costly due to the need for La, Sr-like alkaline-earth/rare-earth materials, and transition metals [6]. To address this issue and lower the costs, non-platinum group metals (non-PGM) have been explored as potential catalysts. These alternatives, typically composed of nitrogen atoms doped onto graphite layers [7,8,9] or transition metal layers [10], offer advantages such as ease of fabrication, high electron conductivity, and cost-effectiveness [11,12,13].
Low-temperature fuel cells, such as proton exchange membrane (PEM) fuel cells, have shown improved efficiency with increasing temperature [14]. A temperature range of 353.15–403.15 K has been identified as particularly promising for low-temperature fuel cells [15]. Some studies have reported significant performance enhancements in PEM fuel cells at a temperature of 393.15 K [15,16]. Raising the temperature has the benefits of increased voltage and efficiency, but it also introduces the drawback of reduced durability [14,16,17].
The output voltage, as per the Nernst equation, is directly proportional to the temperature, with a higher temperature resulting in enhanced charge particle kinetics and an increased fuel cell voltage [18]. However, the lowest performance of the fuel cell is observed at around 343.15 K [19]. As Guan et al. explained, the reason for the lowest performance at 343.15 K and the highest performance at 393.15 K may be a phase change that occurs with increasing temperature. Furthermore, as temperature increases, the oxygen evolution reaction (OER) takes place, leading to the reconstruction of the catalyst surface [20]. To achieve a voltage increase, the enhancement in voltage must exceed the voltage loss due to the negative thermodynamic correlation between the open-circuit voltage and temperature [18,21].
The rate of catalyst decay is influenced by several factors, including temperature, ORR, hydrogen oxidation reaction (HOR), high potential, and the pH of the medium [14]. The formation of carbon monoxide (CO) at the cathode is due to the reaction between water and the carbon support of the catalyst that can cover the catalyst, and this significantly reduces the electrochemical reaction [22]. This process takes place mainly at temperatures below 373.15 K [23]. Above 373.15 K, the water generated from the reaction turns into vapor and is rapidly removed. This dehydration process causes a reduction in CO generation from the catalyst [24].
Numerous non-PGM catalysts, such as Fe/N/Gr, Co/N/Gr, and Mn/N/Gr, have been extensively studied, and some have even reached commercialization [25,26]. For instance, Mn/N/Gr catalysts have been demonstrated to be unstable for potential above 0.53 V [27], while Co/N/Gr-based catalysts tend to produce hydrogen peroxide (H2O2) during ORR, which reduces the current density and can damage the fuel cell membrane [28]. In contrast, Fe-based non-PGM catalysts have been shown to be promising, especially with dual metal structures such as Fe2-N6/Gr and Fe2-N8/Gr [29,30]. Increasing the Fe content to 1.0 wt % can significantly reduce H2O2 generation [31]. However, a drawback of Fe-based catalysts is the formation of Fe(OH)2 during ORR, which can lead to the removal of Fe atoms from the catalyst layer [32]. Furthermore, Fe oxidation can contribute to the catalyst layer degradation [33].
Previous research on copper-based non-PGM single metal catalysts has shown them to have good spontaneous ORR [34,35]. In our previous work, we performed a DFT investigation on Cu-based novel duel metal non-PGM catalyst structures (Cu2-N6/Gr and Cu2-N8/Gr) and a single metal structure (Cu-N4/Gr) [36], and our study showed that all three structures exhibit higher stability than the previously studied Cu-N2/Gr structure. The study also showed that the possibility of H2O2 formation is much less with dual metal structures compared to single metal structures. It is noted that the study was carried out at an initial temperature of 298.15 K and pressure of 1 atm without applying van der Waals force correction. In conclusion, the Cu2-N8/Gr structure showed promising results, confirming its potential for high cell potential with spontaneous ORR and the absence of H2O2 formation. Additionally, Cu-N4/Gr was also shown to have a high cell potential with spontaneous ORR but has a higher probability of forming H2O2.
In this study, DFT calculations were performed to determine the effect of temperature on the Cu2-N8/Gr catalyst in an acidic medium. To avoid the limitations of previous studies, a dispersion correction was applied to the DFT calculation. Salam et al., Amrit et al., and Ferraris et al. have reported that the performance of low-temperature fuel cells changes dramatically at temperatures around 298.15, 353.15, and 393.15 K due to changes in the particle kinetics with temperature using experimental methods, as previously stated [14,15,19]. Therefore, for this study, 298.15, 353.15, and 393.15 K temperature values were selected to investigate the impact of temperature on the fuel cell performance of Cu2-N8/Gr as a cathode catalyst in an acidic medium. At these three temperature values, the structural stability, binding ability of ORR steps, the break-free ability of H2O, evidence of H2O2 generation, and maximum cell potential were predicted as assessments of objectives.

2. Results and Discussion

2.1. Stability and Reactivity

Stability and reactivity were predicted using formation energy data [28,35], and the HOMO–LUMO energy calculation [37,38,39,40] is elaborated in Equation (6) and Equations (15)–(18), respectively. Under open circuit conditions, the Cu2-N8/Gr structure provided formation energy values of −28.856, −28.865, and −28.864 eV for the three different temperatures (298.15, 353.15, and 393.15 K) at zero potential (U = 0), as presented in Table 1. According to Equation (6), an increase in the potential value leads to an enhancement of the formation energy, resulting in a decrease in the structure’s stability. Minor changes in formation energy, at the same potential, generated due to thermal correction were small in comparison to the optimized energy. Negative values in the formation energies indicate that the energy gained favors structural stability. Notably, the formation energy of the structure at all the selected temperatures remained negative for external potentials ranging between 0 and 7.21 V. According to free energy data, the highest maximum cell potential of 0.69 V is observed under an external potential of 7.21 V, which favors this situation. Table 1 illustrates the results showing that the catalyst structure remains stable at all three temperatures as long as the external potential values do not exceed 7.21 V.
In the case of the dual metal structure, Cu2-N8/Gr involves multiple copper atoms that coordinate with eight nitrogen atoms. This configuration results in a significantly more extensive bonding network. It enhances coordination stability when compared to single metal structures, such as Cu-N4/Gr and Cu-N2/Gr, in which only one metal atom coordinates with a smaller number of nitrogen atoms. This increased stability is compatible with findings reported by Yang et al. [41]. Additionally, the dual metal structure is likely to facilitate a more efficient distribution and delocalization of electrons due to the presence of multiple Cu and N atoms. The electrostatic potential (ESP) diagram (Figure 1), which supports the electron distribution theory, for all three temperatures, shows no significant differences in electronic structure. The red area represents the distribution of negative charge (electrons), and the blue and white areas represent positive charge and depleted charge density, respectively.
HOMO–LUMO energy calculations, as explained in Equations (15)–(18), were carried out to investigate the stability and reactivity of the Cu2-N8/Gr and *O2, *H2O, and OH + H2O intermediates [37,38,39,40]. Specifically, the energy gap (Eg), chemical hardness (η), chemical potential (μ), and electrophilicity index (ω) were calculated for these structures. Here, the *O2 and 2*H2O intermediates were considered for the initial and final ORR steps, respectively. These intermediates are important when compared to other ORR steps due to the breaking of the O2 bond with the catalyst that initiates the ORR, and 2H2O is released from the catalyst to conclude the ORR. The OH + H2O intermediate was also considered due to its abnormality. Within all the temperature ranges considered, the HOMO–LUMO calculations remained constant, except for the OH + H2O intermediate, indicating that the catalyst structure did not change for the given temperatures.
Chemical hardness is a measure that reflects the stability and reactivity of a structure, where stable structures tend to be less reactive, and highly reactive structures are typically less stable [9,38,42]. Table 2 shows that the OH + H2O intermediate’s stability increases in the following temperature order: 353.15 K < 298.15 K < 393.15 K. This indicates that at 393.15 K, the intermediate becomes more stable, while at 353.15 K it becomes highly reactive compared to the other temperatures. The ESP diagram in Figure 2 reveals that as the temperature increases from 298.15 to 353.15 K, the electron charge density decreases around the Cu-N structure. The pointed green area in Figure 2b indicates the potential depletion of the electron density. This causes the lower stability of the OH + H2O intermediate at 353.15 K. However, when the temperature reaches 393.15 K, this depletion area is significantly reduced. Furthermore, at 393.15 K, the electron distribution becomes more symmetric (Figure 2c) around the Cu-N structure, as indicated by the similarly sized green area on the opposite side of the blue area. The high stability of the OH + H2O intermediate at 393.15 K is attributed to this symmetric electron charge distribution. According to Watanabe et al., vibrational motion can significantly change the electron density [43], and as the temperature increases, the vibration frequency and direction of the atoms change rapidly. This is attributed to variations in the ESP at different temperatures. To confirm the vibration effects, the vibration modes at different frequencies were examined at a temperature of 353.15 K. A frequency around 289.32 cm−1 indicates the displacement of H atoms, as shown in Figure 3. This mode also has a significantly high infrared value of 25.30 compared to other frequencies. Other vibrational modes were considered but then neglected because their vibration direction did not align with the charge distribution in the ESP diagram. Here, the H atoms represent the most positive charge in the ESP diagram. Therefore, the long displacement vector of the H atom (left H atom) compared to the abnormal electron depletion area is highlighted in Figure 2b. The displacement vector of the hydrogen atom located on the right-hand side in Figure 3 exhibits a smaller magnitude in comparison to the hydrogen atom situated on the left-hand side, an observation consistent with the electrostatic potential (ESP) depicted in Figure 2b. The middle H atom, with the smallest displacement vector, creates a bottleneck in the dumbbell-shaped electron-depletion area in the ESP diagram.
The chemical reactivity of the structure increases as the chemical potential decreases. Chattaraj et al. [38] established a direct connection between Mulliken charge negativity and chemical potential. For all temperatures, the chemical potentials of the O2 and 2H2O intermediates, and the Cu2-N8/Gr structure remain unchanged. The chemical potentials of O2 and 2H2O, which are −3.615 and −3.322 eV, respectively, indicate that O2 has higher reactivity than H2O. Therefore, the optimized 2H2O intermediate is more stable, suggesting that it can easily detach from the catalyst structure due to its weak bond strength and long bond length. The chemical potential of the OH + H2O intermediate decreases in the following order for different temperatures: 298.15 K > 393.15 K > 353.15 K. This order aligns with the high reactivity of the OH + H2O intermediate observed in the chemical hardness data.
The electrophilicity indices of O2 and 2H2O at all temperatures are 8.854 and 8.048 eV, respectively. Since the electrophilicity index represents the ability to attract electrons from the surroundings to form a stable structure, the O2 intermediate is more reactive than 2H2O, as explained in the chemical potential data. For the OH + H2O intermediate, the lowest and highest electrophilicity indices are observed at temperatures of 393.15 and 353.15 K, respectively, which is in agreement with the chemical hardness data. The order of reactivity for the OH + H2O intermediate is as follows for different temperatures: 353.15 K > 298.15 K > 393.15 K, which confirms the chemical potential data, indicating its highest reactivity at 353.15 K.
The infrared (IR) spectra, illustrated in Figure 4a–d, were examined for the catalyst (*), *O2, *OH + H2O, and *H2O + H2O. Notably, the peaks remained unchanged with variations in temperature, except for *OH + H2O, as depicted in Figure 4c. This consistent behavior in the IR spectrum substantiates stability across temperature variations and aligns with the results obtained from other calculations, such as HOMO–LUMO energy calculations assessing stability and reactivity.
Moreover, Figure 4c highlights that at temperatures of 298.15 K and 393.15 K, similar peaks are observed. However, at 353.15 K, shifts or variations in peaks are evident, indicating distinct behaviors of the *OH + H2O intermediate at this temperature. Simulation results suggest vibrational activity around the frequency of 289.32 cm−1 for the OH + H2O intermediate, while higher peaks signify vibrations of C and H atoms at the catalyst layer’s edge.
Furthermore, at 353.15 K, certain peaks of the OH + H2O intermediate exhibit a shift towards lower wavenumbers within the range of 1200–1600 cm−1. This observation implies thermal denaturation, as explained by Panick et al. [44], or a phase change associated with temperature, as elucidated by Guan et al. [20]. These findings contribute valuable insights into the temperature-dependent dynamics of the studied OH + H2O intermediate.

2.2. ORR Steps and Binding Energy

In this study, the ORR steps did not indicate deviations at different temperatures. However, when the intermediate 2OH is bonded with an H atom to proceed to the next step, some abnormalities are observed. The expected outcome after *2OH− + H+ was *OH with a separated H2O molecule. However, after optimization, the *OH + H + *OH intermediate was formed, where both OH molecules shared a single H atom (Figure 5). With the exception of *OH + H + *OH (intermediates after *2OH− + H+), the other intermediates did not indicate any deviations in terms of bond lengths, angles, and binding strength with different temperatures (Table 3).
In the first step of ORR, the oxygen molecule exhibits a significant binding strength value compared with that reported by Kattel et al. and Bhatt et al. [28,35], and it is slightly mitigated with increasing temperatures. This trend has been observed for all other intermediates, as recorded in Table 3, except for the OH + H2O intermediate (*OH + H + *OH). Initially, the *OH + H + *OH intermediate shows a reduction in its binding strength that mirrors the behavior of the other intermediates. However, at 393.15 K, its binding strength increases compared to the lower temperature of 353.15 K.
Binding energy calculations reveal that the highest binding energy for the OH intermediate was at −6.45 eV, which is observed only at an initial temperature of 298.15 K. At other temperatures, the binding energy significantly decreases to −0.25 eV at 353.15 K and −0.24 eV at 393.15 K. These binding energy values at higher temperatures are not favorable for the ORR, as they represent weaker binding strength for all intermediates. In a typical case, the OH intermediate is generated after the water molecules are rendered and moved away, and subsequently, an H atom attaches itself to the O intermediate, as shown in the ORR reaction Equations (11)–(13). The optimized reactions show that the H atom is attached to the 2OH intermediate, instead of forming a water molecule and OH intermediate. This results in the formation of the OH + H + OH intermediate, as previously explained. Therefore, the weakest bond strength of the OH intermediate is not considered a step in the ORR.
In the final step of the ORR, it is expected that the water molecule will form and easily break away from the catalyst. The calculations confirmed that H2O has the lowest binding strength, which favors the described situation. Furthermore, the longest bond length of 2.144 Å between the Cu and O atoms of the catalyst intermediate confirms the weak bond strength. However, according to the optimized ORR steps, the possibility of a single H2O intermediate forming at the end of the reaction is not significant. Due to the presence of the OH + H + OH intermediate, the possibility of the formation of a 2H2O intermediate is substantial at the end of ORR. The binding energies of the 2H2O intermediates are −2.12, −2.11, and −2.10 eV at temperatures of 298.15, 353.15, and 393.15 K, respectively. When compared to the other ORR steps, the 2H2O intermediate possesses the weakest bond strength and the longest bond length at all temperatures, favoring the situation.
As the OOH intermediate emerges during the ORR (reaction Equation (2)), there is a significant possibility of H2O2 generation in the middle of the reaction. To investigate this, the H atom is introduced to the first oxygen atom to determine whether H2O2 is formed. Figure 6 illustrates this reaction at an initial temperature of 298.15 K. As reported in our previous study [36], under the initial temperature and the other two temperatures, the formation of the 2OH intermediate occurs instead of H2O2 molecules. This phenomenon can be attributed to the electronegativity of nitrogen (N) being higher than that of copper (Cu). In the Cu2-N8 structure, where one Cu atom is attached to four N atoms, this arrangement is responsible for the electron cloud depletion around the Cu atoms. Consequently, this can lead to the electrons around the oxygen of the OOH intermediate being highly attracted to Cu atoms. The high binding strengths of the 2OH intermediate, which are −10.97, −8.51, and −8.50 eV at the temperatures of 298.15, 353.15, and 393.15 K, respectively, support this explanation. When a H atom is attached to the first oxygen atom, it can draw electrons away from the first oxygen atom. This effect weakens the O-O bond of the OOH intermediate and leads to the generation of the 2OH intermediate.
According to the studies by Xiao et al. [34] and Noh et al. [45], the H atom could bind to the terminal O atom (Figure 7a) of the OOH intermediate to form an H2O molecule. In this study, the possibility of this reaction was determined via the optimization process. For all three temperatures, the final optimization resulted in the formation of a dual OH intermediate (Figure 7c and Figure 8c) instead of the formation of the water molecule. Figure 7 depicts the reaction at a temperature of 298.15 K. When the introduced H atom attempts to bond with the second O atom of OOH to form H2O, the previously bonded H atom breaks its bond with the second O atom and moves toward the first O atom (Figure 7b). During this process, a water molecule forms in the initial optimization steps, and then the Cu atom near the second O atom attracts it to establish a strong bond. This phenomenon is the same as that explained in previous H2O2 generation by the OOH intermediate. Due to electron depletion around the terminal O of *OOH, the O atom could not bear two H atoms as H2O. This could cause the initially bonded H atom to break its bond with the second O atom.
Figure 8 illustrates the optimized ORR steps, which differ from previous studies [28]. The changes observed in the 2OH and OH + H + HO intermediates, as opposed to O + H2O and OH +H2O in the ORR, may be attributed to the dual metal catalyst structure, as explained before. The following reaction Equations (1)–(5) depict the optimized ORR steps from the current study in an acidic medium. The asterisk (*) indicates the defect’s configuration.
* + O2 + 4H+ + 4e−→ *O2 + 4H+ + 4e−
*O2 + 4H+ + 4e− → *OOH + 3H+ + 3e−
*OOH + 3H+ + 3e−→ *2OH + 2H+ + 2e−
*2OH + 2H+ + 2e− → *(OH + H + HO) + H+ + e−
*(OH + H + HO) + H+ + e− → * + 2H2O

2.3. Reaction Spontaneity

All free energy calculations were performed for the ground state of the ORR steps, as shown in Figure 9. After conducting free energy calculations for all three temperatures, unlike in our previous studies [36], it was found that the maximum cell potential was greater than 0.60 V. At the initial temperature of 298.15 K, all reaction coordinates exhibited a downhill process (spontaneous process) until the cell potential reached 0.66 V (maximum cell potential). Beyond the maximum cell potential, the last step of the ORR became an uphill process (Figure 9a), indicating an increase in the energy barrier between the reaction coordinates. Further enhancement of cell potential is attributed to the initiation of the uphill process from the 2OH intermediate step. The difference in maximum cell potentials between the previous and current findings may be attributed to the inclusion of a dispersion correction for van der Waals forces in the Gaussian calculations of the current study.
By increasing the temperature to 353.15 K, the maximum cell potential decreased to 0.61 V (Figure 9b). Beyond 0.61 V, the chemical reactions become uphill processes from the 2OH intermediate step to the OH + H2O intermediate step. This adverse trend intensifies with higher temperatures. When the temperature was raised to 393.15 K, the maximum cell potential reached the highest value of 0.69 V (Figure 9c). For cell potentials higher than 0.69 V, the OH + H2O intermediate step acted as a barrier to the downhill process of the ORR. According to the free energy data, further increases in the cell potential augmented the uphill process between the 2OH and OH + H2O intermediates. The order of maximum cell potential at different temperatures was as follows: 353.15 K < 298.15 K < 393.15 K, indicating that temperature and maximum cell potential were not directly proportional. The decrease in maximum cell voltage at around 353.15 K suggests a decline in cell performance, which aligns with the findings of Ferraris et al. [19], who reported the lowest cell potential around the operating temperature of 343.15 K, supporting our observation. Furthermore, other investigations [15,16] suggest that at 393.15 K PEM fuel cell performance improved dramatically in congruence with cell voltage escalation due to high temperature.
In the free energy calculations, the observed up–down behavior of the maximum cell potential with temperature is mainly caused by OH + H2O intermediate step. The ESP diagram of the OH + H2O intermediate (Figure 2c) delineates that at 393.15 K, this intermediate gains high stability due to a dense and symmetric electron density distribution. This stability, along with higher electron mobility at higher temperatures, could contribute to a higher maximum cell potential observed at 395.15K compared to 353.15 K, a conclusion supported by HOMO–LUMO energy calculations [23,24].
The free energy analysis of the first and last steps in the ORR revealed a 15 eV gap, showing a departure from the typical values observed in conventional ORR catalysts. In a work by Kattel et al. [28], the free energy of O2 was determined through the reaction O2 + 2H2 → 2H2O, resulting in a free energy change of 4.92 eV. This value contrasts with most other density functional theory (DFT) calculations, where such values typically fall below 8 eV [27,46]. While a high free energy gap might suggest reduced effectiveness of the catalyst, insights from Santisouk et al. [47] challenge this notion. Binding energy, bond length, and HOMO–LUMO energy calculations indicate the effectiveness of the catalyst with intermediates. Notably, the introduction of dispersion correction (GD3BJ) in the current study contributes significantly to this observed change. In contrast, a previous study [36] of ours in the same research field did not exhibit such a substantial free energy gap, with the O2 to 2H2O free energy change hovering around 6 eV. This prompts consideration of alternative dispersion correction methods that could potentially enhance the accuracy of the free energy results.
Furthermore, given the novelty of the catalyst’s structure, additional studies are warranted to validate its effectiveness. It is worth noting that some ORR free energy calculations in the literature indicate differences between the first and last ORR steps ranging from 8 to 14 eV [34,47,48] emphasizing the variability in observed outcomes across different studies.
Table 4 presents the binding energy, formation energy, and bond length of molecules on the catalyst surface at the initial temperature for various catalyst types. All the catalysts under investigation are non-platinum group metal (non-PGM) types, featuring transition metals and nitrogen doping. In terms of formation energy, our current study yields the most promising results, with the cobalt (Co)-based catalyst displaying the least promising outcome. Specifically, the titanium (Ti)-based catalyst exhibits the strongest binding capability for O2 and H2O intermediates, suggesting its limited effectiveness as a catalyst material. Conversely, the Co-based catalyst demonstrates the weakest binding capability for the O2 intermediate, indicating its diminished efficacy as a catalyst. In the context of H2O2 formation at the conclusion of the ORR, our current study reveals no H2O2 formation, whereas all other catalysts exhibit such potential. Notably, the Co-based catalyst, characterized by its weakest binding strength and longest bond length, implies the likelihood of H2O2 formation and easy detachment from the catalyst layer. Furthermore, all dual molecules in our current study exhibit dual bonding, contributing to their strong attachment to the catalyst. These findings provide valuable insights into the diverse binding characteristics of the catalysts, shedding light on their potential efficacy in the context of the ORR.
Nonetheless, it is important to acknowledge several limitations of this study. Firstly, the use of a small 3–21G basis set may be considered a limitation. However, this choice was made to compare the results with previous investigations by Bhatt et al. [35], as well as to maintain consistency with our earlier study [36] and to compare with the Cu-N2/Gr results. Additionally, insufficient computational power meant that the study could not use a higher basis set. Furthermore, this study did not incorporate periodic boundary conditions (PBC), which distinguishes it from the approach employed by Xiao et al. [34]. The utilization of different DFT simulation software in these referenced studies precludes direct comparison with our findings. It should also be noted that the preference for Gaussian software over VASP or ab initio molecular dynamics simulation software was influenced by financial constraints and a lack of specialists. In this study, a single layer of graphene was employed to examine the impact of the defect structures on the surface of the catalyst in the context of ORR, a methodology similar to those of Bhatt et al. [35] and Kattel et al. [28]. In future studies related to this research, the use of the PBC approach will be explored to facilitate more direct comparisons with the current study and potentially expand the scope of investigation.

3. Methodology

The DFT calculations were performed using the Gaussian 09w software, applying the B3LYP method with a 3–21G basis set. The GD3BJ van der Waals correction (dispersion correction) was applied for every calculation. The form of the defect structures was visualized with GaussView 6.0. All atoms in each structure were relaxed using an optimization function in Gaussview 6 with the previously stated basis set and method. A single pristine graphene layer was created as the base, and from it, the Cu2-N8/Gr structure was developed to investigate the ORR under various temperatures (298.15, 353.15, and 393.15 K). The pristine graphene layer design at the initial temperature has C-C and C-H bond lengths of 1.43 and 1.08 Å, respectively, and a C-C-C angle of 119.92˚ which is compatible with the studies by Bhatt et al. and Hernandez et al. [35,49]. The structure that features the pyridinic nitrogen was constructed. Following the initial optimization, structures with a 1-unit imaginary frequency (=1) were reoptimized by using a stability function in GaussView 6. Individual optimization of the molecules involved in the ORR steps was carried out, with subsequent measurement and comparison of bond lengths and bond angles with previous studies by Xiao et al., Bhatt et al., and Balasooriya et al. [34,35,36] on this subject.
The stability of the catalyst structure with the cathode potential at different temperatures was confirmed. Stability was inspected using the formation energy (ΔE) defined as [28,35],
Δ E = E g r a p h e n e + ( C u 2 − N 8 ) + y μ C − E g r a p h e n e + x μ N + E M
Here, E g r a p h e n e + ( C u 2 − N 8 ) is the energy of the optimized Cu2-N8 graphene layer with the thermal correction to energy. μ C and μ N are the chemical potential of carbon, which is defined as the total energy per carbon atom for defect-free graphene, and the chemical potential of nitrogen, defined as half of the total energy of the N2 molecule with the thermal correction to energy, respectively. x represents the number of nitrogen atoms added, and y represents the carbon atoms removed during the defect formation. M is the metal atom (Cu) doped in the system. E g r a p h e n e is the energy of the optimized pristine graphene layer with the thermal correction to energy. E M   is the total energy of Mn+ defined as
E M n + = E M − n e U
where E M is the total energy with the thermal correction to energy of isolated M (M = Cu) in the gas phase and n, e, and U are the number of electron transfer (+2), electron charge, and the external potential, respectively.
Binding energies (BEs) were calculated for the ORR steps in an acidic medium, as defined by [28,34,35],
B E = E d e f e c t + m o l e c u l e − ( E d e f e c t + E m o l e c u l e )
Here, E d e f e c t + m o l e c u l e is the total energy with the thermal correction to the energy of molecules adsorbed by the defective graphene. E d e f e c t is the total energy with the thermal correction to energy of the defective graphene configuration, and E m o l e c u l e is the energy of isolated molecule species (O2, H2O, 2H2O, OOH, OH, 2OH, H2O2) with the thermal correction to energy. Negative signed binding energies indicate energy released during the bonding process which is more favorable to the process. The formation of H2O2 during ORR is also considered.
The ORR steps of an ordinary catalyst in an H+ medium are defined as follows [28], where * indicates the defect’s configuration.
* + O2 + 4H+ + 4e−→ *O2 + 4H+ + 4e−
*O2 + 4H+ + 4e− → *OOH + 3H+ + 3e−
*OOH + 3H+ + 3e−→ *O + H2O + 2H+ + 2e−
*O + H2O + 2H+ + 2e− → *OH + H2O + H+ + e−
*OH + H2O + H+ + e−→ * + 2H2O
The free energies were calculated for each ORR step of all defect configurations as defined by [28],
Δ G = Δ E + Δ Z P E − T Δ S + Δ G U + Δ G p H + Δ G f i e l d
Δ E represents the energy calculated from DFT related to the relevant reaction step with the thermal correction to energies, while Δ Z P E is the correction obtained from DFT calculations, accounting for the zero-point energy. The value of entropy (S) was obtained from the DFT calculations for the three temperatures (T = 298.15, 353.15, and 393.15 K) as a variable in this study, depending on temperature (T). In the context of this research, ΔGU = −eU where U and e are the electrode potential and the charge transferred, respectively. Additionally, Δ G p H = k B T × l n 10 × p H , where kB is the Boltzmann’s constant, and T is the temperature. Since this study considered acidic media, pH = 0. Therefore, normally the Δ G p H value is considered to be zero. Δ G f i e l d is the contribution of the interaction of adsorbate with the local electric field in the electric double layer formed in the vicinity of the cathode, which is negligible according to Nørskov et al. [50]. Corrections of zero-point energy and entropy values were obtained after frequency calculations. Free energy vs. reaction coordinate graphs were plotted with different potentials (U) for all temperature values.
The energy gap (Eg), chemical hardness (η), chemical potential (μ), and electrophilicity index (ω) were calculated using Koopman’s principle for optimized structures as defined by [37,38,39,40,51,52],
E g = E L U M O − E H O M O
η = I − A 2
μ = − ( I + A ) 2
ω = μ 2 2 η
Here, I and A are the ionization potential (≅−EHOMO) and the electron affinity (≅−ELUMO), respectively.

4. Conclusions

Leveraging the computational prowess of density functional theory, insights were gleaned into the structural dynamics of the Cu2-N8/Gr motif across three pivotal temperature gradients: 298.15, 353.15, and 393.15 K. Within these defined thermal windows, the motif exhibited an inherent predisposition toward spontaneous oxygen reduction reactions (ORRs), with a peak cell potential of 0.69 V recorded at an elevated threshold of 393.15 K. A discernible shift in the maximum cell potential was observed as temperatures progressed from 298.15 to 353.15 K, marking a transition from 0.66 to 0.61 V. Formation energy evaluations underscored the motif’s unwavering stability across a potential (U) spectrum ranging from 0 to 7.21 V, irrespective of the temperature in question. The maximum cell potential of the motif within the potential range of stability delineates better structural stability. Delving deeper into chemical metrics—spanning chemical hardness, potential, and electrophilicity—the motif showcased consistent ORR trajectories across the temperature continuum, with the singular exception of the OH + H2O (manifesting as OH + H +OH) intermediate. Intricate vibrational and electrostatic potential assessments pinpointed this intermediate as the pivotal agent that mediates cell voltage variations across disparate temperatures. Despite a subtle attenuation in ORR stepwise stability with surging temperatures, the overarching cell potential intriguingly reached its apogee at 393.15 K, indicating that catalyst performance had not been reduced. In fact, the HOMO–LUMO evaluations and IR spectrum reiterate this assertion, highlighting the motif’s unwavering stability across escalating temperatures. At the zenith of 393.15 K, the motif not only boasts unparalleled stability but also outperforms others in terms of cell potential, solidifying its credentials as a prime catalyst contender. The binding energy values slightly escalate with increasing temperature, except for the OH + H2O intermediate for the ORR steps. The binding dynamics of this intermediate witnessed a decline as temperatures transitioned from 298.15 K to 353.15 K, followed by an upswing at 393.15 K. Singular OH and H2O intermediates, due to the optimization dynamics, were found to be incongruent with conventional ORR pathways. Free energy assessments across the temperature spectrum showcased the motif’s spontaneous reactivity, consistently exceeding the 0.6 V cell potential threshold. Significantly, at 393.15 K, this potential reached an unparalleled 0.69 V, the pinnacle at all examined temperatures. To resolve the higher gap between the first and last ORR steps of the free energy data, further investigations should be carried out with different dispersion corrections and basis sets. Therefore, considering stability, binding strength, maximum cell potential, and spontaneity, the catalyst shows the best performance at 393.15 K, while its performance at 353.15 K appears sub-optimal, a sentiment echoed in the existing literature [15,16,19]. Similar to our previous work, this study also delineates a low probability of forming H2O2 at all temperatures in DFT optimization. Consequently, adherence to the quartet-electron pathways across all temperature spectra potentially preserves current density uniformity. In summary, the Cu2-N8/Gr catalyst, when assayed at 393.15 K, emerges as the gold standard, melding superior cell potential with minimized H2O2 formation.

Author Contributions

Y.B.: conceptualization, methodology, software, simulation, writing—original draft preparation, investigation, and visualization; P.S.: conceptualization, supervision, methodology, and writing—review and editing; C.M.L.: writing—review and editing; Y.-F.C.C.: writing—review and editing; M.R.R.K.: supervision, software, writing—review and editing, and funding acquisition; R.T.: conceptualization, supervision, methodology, writing—review and editing, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universiti Brunei Darussalam (Grant No.: UBD/RSCH/1.9/FICBF(b)2022/017), and the APC is funded by the corresponding author.

Data Availability Statement

All data used to support the findings of this study are included in the article.

Acknowledgments

The authors would like to thank the Government of Brunei Darussalam, the Universiti Brunei Darussalam, and the Wayamba University of Sri Lanka for their continuous support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antolini, E. Carbon supports for low-temperature fuel cell catalysts. Appl. Catal. B Environ. 2009, 88, 1–24. [Google Scholar] [CrossRef]
  2. Stephen Walch, A.D. Masoud Aryanpour, and Heinz Pitsch. Mechanism of Molecular Oxygen Reduction at the Cathode of a PEM Fuel Cell Non-Electrochemical Reactions on Catalytic Pt Particles. J. Phys. Chem. 2008, 112, 8464–8475. [Google Scholar]
  3. Roland, U.; Braunschweig, T.; Roessner, F. On the nature of spilt-over hydrogen. J. Mol. Catal. A Chem. 1997, 127, 61–84. [Google Scholar] [CrossRef]
  4. Curtis Conner, W.; Pajonk, G.M.; Teichner, S.J. Spillover of Sorbed Species. In Advances in Catalysis; Academic Press: Cambridge, MA, USA, 1986; pp. 1–79. [Google Scholar]
  5. Lin, R.; Cai, X.; Hao, Z.; Pu, H.; Yan, H. Rapid microwave-assisted solvothermal synthesis of shape-controlled Pt-Ni alloy nanoparticles for PEMFC. Electrochim. Acta 2018, 283, 764–771. [Google Scholar] [CrossRef]
  6. Guan, D.; Xu, H.; Zhang, Q.; Huang, Y.C.; Shi, C.; Chang, Y.C.; Xu, X.; Tang, J.; Gu, Y.; Pao, C.W.; et al. Identifying a Universal Activity Descriptor and a Unifying Mechanism Concept on Perovskite Oxides for Green Hydrogen Production. Adv. Mater. 2023, 35, 2305074. [Google Scholar] [CrossRef] [PubMed]
  7. Sidik, R.A.; Anderson, A.B.; Subramanian, N.P.; Kumaraguru, S.P.; Popov, B.N. O2 Reduction on Graphite and Nitrogen-Doped Graphite Experiment and Theory. J. Phys. Chem 2006, 110, 1787–1793. [Google Scholar] [CrossRef] [PubMed]
  8. Lee, K.R.; Lee, K.U.; Lee, J.W.; Ahn, B.T.; Woo, S.I. Electrochemical oxygen reduction on nitrogen doped graphene sheets in acid media. Electrochem. Commun. 2010, 12, 1052–1055. [Google Scholar] [CrossRef]
  9. Zhang, L.; Xia, Z. Mechanisms of Oxygen Reduction Reaction on Nitrogen-Doped Graphene for Fuel Cells. J. Phys. Chem. C 2011, 115, 11170–11176. [Google Scholar] [CrossRef]
  10. Rao, C.V.; Cabrera, C.R.; Ishikawa, Y. In Search of the Active Site in Nitrogen-Doped Carbon Nanotube Electrodes for the Oxygen Reduction Reaction. J. Phys. Chem. Lett. 2010, 1, 2622–2627. [Google Scholar] [CrossRef]
  11. Lin, R.; Che, L.; Shen, D.; Cai, X. High durability of Pt-Ni-Ir/C ternary catalyst of PEMFC by stepwise reduction synthesis. Electrochim. Acta 2020, 330, 135251. [Google Scholar] [CrossRef]
  12. Beermann, V.; Holtz, M.E.; Padgett, E.; de Araujo, J.F.; Muller, D.A.; Strasser, P. Real-time imaging of activation and degradation of carbon supported octahedral Pt–Ni alloy fuel cell catalysts at the nanoscale using in situ electrochemical liquid cell STEM. Energy Environ. Sci. 2019, 12, 2476–2485. [Google Scholar] [CrossRef]
  13. Pinheiro, A.L.N.; Oliveira-Neto, A.; de Souza, E.C.; Perez, J.; Paganin, V.A.; Ticianelli, E.A.; Gonzalez, E.R. Electrocatalysis-on-Noble-Metal-and-Noble-Metal-Alloys-Dispersed-on-High-Surface-Area-Carbon. J. New Mater. Electrochem. Syst. 2003, 6, 1–8. [Google Scholar]
  14. Salam, M.A.; Habib, M.S.; Ahmed, K.; Uddin, M.S.; Hossain, T.; Papri, N.; Arefin, P. Effect of temperature on the performance factors and durability of proton exchange membrane of hydrogen fuel cell: A narrative review. Mater. Sci. Res. India (Online) 2020, 17, 179–191. [Google Scholar] [CrossRef]
  15. Chandan, A.; Hattenberger, M.; El-Kharouf, A.; Du, S.; Dhir, A.; Self, V.; Pollet, B.G.; Ingram, A.; Bujalski, W. High temperature (HT) polymer electrolyte membrane fuel cells (PEMFC)–A review. J. Power Sources 2013, 231, 264–278. [Google Scholar] [CrossRef]
  16. Esfeh, H.K.; Hamid, M.K.A. Temperature effect on proton exchange membrane fuel cell performance Part II: Parametric study. Energy Procedia 2014, 61, 2617–2620. [Google Scholar] [CrossRef]
  17. Jourdani, M.; Mounir, H.; Marjani, A.E. Latest trends and challenges in proton exchange membrane fuel cell (PEMFC). Open Fuels Energy Sci. J. 2017, 10, 96–105. [Google Scholar] [CrossRef]
  18. Li, C.; Liu, Y.; Xu, B.; Ma, Z. Finite time thermodynamic optimization of an irreversible proton exchange membrane fuel cell for vehicle use. Processes 2019, 7, 419. [Google Scholar] [CrossRef]
  19. Ferraris, A.; Messana, A.; Airale, A.G.; Sisca, L.; de Carvalho Pinheiro, H.; Zevola, F.; Carello, M. Nafion® tubing humidification system for polymer electrolyte membrane fuel cells. Energies 2019, 12, 1773. [Google Scholar] [CrossRef]
  20. Guan, D.; Shi, C.; Xu, H.; Gu, Y.; Zhong, J.; Sha, Y.; Hu, Z.; Ni, M.; Shao, Z. Simultaneously mastering operando strain and reconstruction effects via phase-segregation strategy for enhanced oxygen-evolving electrocatalysis. J. Energy Chem. 2023, 82, 572–580. [Google Scholar] [CrossRef]
  21. Khotseng, L. Fuel cell thermodynamics. In Pem Fuel Cell Testing and Diagnosis; Elsevier: Amsterdam, The Netherlands, 2019; pp. 25–75. [Google Scholar]
  22. Maass, S.; Finsterwalder, F.; Frank, G.; Hartmann, R.; Merten, C. Carbon support oxidation in PEM fuel cell cathodes. J. Power Sources 2008, 176, 444–451. [Google Scholar] [CrossRef]
  23. Zhang, J.; Zhang, H.; Wu, J.; Zhang, J. The effects of temperature on PEM fuel cell kinetics and performance. Pem. Fuel Cell Test. Diagn. 2013, 121–141. [Google Scholar]
  24. Zeis, R. Materials and characterization techniques for high-temperature polymer electrolyte membrane fuel cells. Beilstein J. Nanotechnol. 2015, 6, 68–83. [Google Scholar] [CrossRef]
  25. Orellana, W. Catalytic Properties of Transition Metal–N4 Moieties in Graphene for the Oxygen Reduction Reaction: Evidence of Spin-Dependent Mechanisms. J. Phys. Chem. C 2013, 117, 9812–9818. [Google Scholar] [CrossRef]
  26. Khotseng, L. Oxygen Reduction Reaction. In Electrocatalysts for Fuel Cells and Hydrogen Evolution—Theory to Design; BoD–Books on Demand: Norderstedt, Germany, 2018. [Google Scholar]
  27. Liu, K.; Qiao, Z.; Hwang, S.; Liu, Z.; Zhang, H.; Su, D.; Xu, H.; Wu, G.; Wang, G. Mn- and N- doped carbon as promising catalysts for oxygen reduction reaction: Theoretical prediction and experimental validation. Appl. Catal. B Environ. 2019, 243, 195–203. [Google Scholar] [CrossRef]
  28. Kattel, S.; Atanassov, P.; Kiefer, B. Catalytic activity of Co-N(x)/C electrocatalysts for oxygen reduction reaction: A density functional theory study. Phys. Chem. Chem. Phys. 2013, 15, 148–153. [Google Scholar] [CrossRef]
  29. Chen, C.; Sun, M.; Wang, K.; Li, Y. Dual-metal single-atomic catalyst: The challenge in synthesis, characterization, and mechanistic investigation for electrocatalysis. SmartMat 2022, 3, 533–564. [Google Scholar] [CrossRef]
  30. Zhang, N.; Zhou, T.; Ge, J.; Lin, Y.; Du, Z.; Zhong, C.A.; Wang, W.; Jiao, Q.; Yuan, R.; Tian, Y.; et al. High-Density Planar-like Fe2N6 Structure Catalyzes Efficient Oxygen Reduction. Matter 2020, 3, 509–521. [Google Scholar] [CrossRef]
  31. Masa, J.; Zhao, A.; Xia, W.; Muhler, M.; Schuhmann, W. Metal-free catalysts for oxygen reduction in alkaline electrolytes: Influence of the presence of Co, Fe, Mn and Ni inclusions. Electrochim. Acta 2014, 128, 271–278. [Google Scholar] [CrossRef]
  32. Na Yang, L.P.; Li, L.; Jing, L.; Qiang, L.; Minhua, S.; Zidong, W. Theoretically probing the possible degradation mechanisms of an FeNC catalyst during the oxygen reduction reaction. Chem. Sci. 2021, 12, 12476–12484. [Google Scholar] [CrossRef]
  33. Xie, X.; He, C.; Li, B.; He, Y.; Cullen, D.A.; Wegener, E.C.; Kropf, A.J.; Martinez, U.; Cheng, Y.; Engelhard, M.H.; et al. Performance enhancement and degradation mechanism identification of a single-atom Co–N–C catalyst for proton exchange membrane fuel cells. Nat. Catal. 2020, 3, 1044–1054. [Google Scholar] [CrossRef]
  34. Xiao, Y.; Zhang, W. DFT analysis elementary reaction steps of catalytic activity for ORR on metal-, nitrogen- co-doped graphite embedded structure. SN Appl. Sci. 2020, 2, 194. [Google Scholar] [CrossRef]
  35. Bhatt, M.D.; Lee, G.; Lee, J.S. Density Functional Theory (DFT) Calculations for Oxygen Reduction Reaction Mechanisms on Metal-, Nitrogen- co-doped Graphene (M-N2-G (M = Ti, Cu, Mo, Nb and Ru)) Electrocatalysts. Electrochim. Acta 2017, 228, 619–627. [Google Scholar] [CrossRef]
  36. Balasooriya, Y.; Samarasekara, P.; Lim, C.M.; Chau, Y.-F.C.; Kooh, M.R.R.; Thotagamuge, R. Cu-Nitrogen doped graphene (Cu–N/Gr) nanocomposite as cathode catalyst in fuel cells–DFT study. Heliyon 2023, 9, e15989. [Google Scholar] [CrossRef] [PubMed]
  37. Glossman-Mitnik, D. Computational Study of the Chemical Reactivity Properties of the Rhodamine B Molecule. Procedia Comput. Sci. 2013, 18, 816–825. [Google Scholar] [CrossRef]
  38. Chattaraj, P.K.; Giri, S. Electrophilicity index within a conceptual DFT framework. Annu. Rep. Sect. "C" (Phys. Chem.) 2009, 105, 13–39. [Google Scholar] [CrossRef]
  39. Parr, R.G.; Szentpaly, L.v.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  40. Syaahiran, A.; Lim, C.M.; Kooh, M.R.R.; Mahadi, A.H.; Chau, Y.-F.C.; Thotagamuge, R. A theoretical insight of Cr dopant in tungsten oxide for gas sensor application. Mater. Today Commun. 2021, 28, 102508. [Google Scholar] [CrossRef]
  41. Yang, X.; Priest, C.; Hou, Y.; Wu, G. Atomically dispersed dual-metal-site PGM-free electrocatalysts for oxygen reduction reaction: Opportunities and challenges. SusMat 2022, 2, 569–590. [Google Scholar] [CrossRef]
  42. Syaahiran, M.A.; Mahadi, A.H.; Lim, C.M.; Kooh, M.R.R.; Chau, Y.C.; Chiang, H.P.; Thotagamuge, R. Theoretical Study of CO Adsorption Interactions with Cr-Doped Tungsten Oxide/Graphene Composites for Gas Sensor Application. ACS Omega 2022, 7, 528–539. [Google Scholar] [CrossRef]
  43. Watanabe, N.; Yamazaki, M.; Takahashi, M. Vibrational effects on valence electron momentum distributions of ethylene. J. Chem. Phys. 2012, 137, 114301. [Google Scholar] [CrossRef] [PubMed]
  44. Panick, G.; Malessa, R.; Winter, R. Differences between the pressure-and temperature-induced denaturation and aggregation of β-lactoglobulin A, B, and AB monitored by FT-IR spectroscopy and small-angle X-ray scattering. Biochemistry 1999, 38, 6512–6519. [Google Scholar] [CrossRef] [PubMed]
  45. Noh, S.H.; Kwak, D.H.; Seo, M.H.; Ohsaka, T.; Han, B. First principles study of oxygen reduction reaction mechanisms on N-doped graphene with a transition metal support. Electrochim. Acta 2014, 140, 225–231. [Google Scholar] [CrossRef]
  46. Szakacs, C.E.; Lefevre, M.; Kramm, U.I.; Dodelet, J.P.; Vidal, F. A density functional theory study of catalytic sites for oxygen reduction in Fe/N/C catalysts used in H(2)/O(2) fuel cells. Phys. Chem. Chem. Phys. 2014, 16, 13654–13661. [Google Scholar] [CrossRef] [PubMed]
  47. Santisouk, S.; Sengdala, P.; Jiang, X.; Xue, X.X.; Chen, K.Q.; Feng, Y. Tuning the Electrocatalytic Properties of Black and Gray Arsenene by Introducing Heteroatoms. ACS Omega 2021, 6, 13124–13133. [Google Scholar] [CrossRef] [PubMed]
  48. Hwang, Y.; Yun, K.-H.; Chung, Y.-C. Carbon-free and two-dimensional cathode structure based on silicene for lithium–oxygen batteries: A first-principles calculation. J. Power Sources 2015, 275, 32–37. [Google Scholar] [CrossRef]
  49. Hernandez Rosas, J.J.; Ramirez Gutierrez, R.E.; Escobedo-Morales, A.; Chigo Anota, E. First principles calculations of the electronic and chemical properties of graphene, graphane, and graphene oxide. J Mol Model 2011, 17, 1133–1139. [Google Scholar] [CrossRef] [PubMed]
  50. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist LR, K.J.; Kitchin, J.R.; Bligaard, T.; Jonsson, H. Logadottir, and L. Lindqvist. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  51. Brahim, N.S.; Thotagamuge, R.; Kooh, M.R.R.; Lim, C.M.; Syaahiran, M.A.; Usman, A.; Shahri, N.N.M.; Chou Chau, Y.-F.; Chou Chao, C.-T.; Chiang, H.-P. Enhanced CO Gas Sensing with DFT Optimized PbS Loading on ZnO and CrZnO Nanocomposites. Sustainability 2022, 14, 13978. [Google Scholar] [CrossRef]
  52. Kumara, N.; Kooh, M.R.R.; Lim, A.; Petra, M.I.; Voo, N.Y.; Lim, C.M.; Ekanayake, P. DFT/TDDFT and experimental studies of natural pigments extracted from black tea waste for DSSC application. Int. J. Photoenergy 2013, 2013. [Google Scholar] [CrossRef]
Molecules 28 07873 g001
Figure 1. Electrostatic potential diagram of the defect structure at all temperatures (red and blue areas represent the negative and positive charges, respectively, and the increasing density of the color represents the increasing charge density).
Figure 1. Electrostatic potential diagram of the defect structure at all temperatures (red and blue areas represent the negative and positive charges, respectively, and the increasing density of the color represents the increasing charge density).
Molecules 28 07873 g001
Molecules 28 07873 g002
Figure 2. Electrostatic potential diagram of the OH + H2O intermediate at the temperatures of (a) 298.15 K (b) 353.15 K, and (c) 393.15 K (red and blue areas represent the negative and positive charges, respectively, and increasing density of color represents increasing charge density).
Figure 2. Electrostatic potential diagram of the OH + H2O intermediate at the temperatures of (a) 298.15 K (b) 353.15 K, and (c) 393.15 K (red and blue areas represent the negative and positive charges, respectively, and increasing density of color represents increasing charge density).
Molecules 28 07873 g002
Molecules 28 07873 g003
Figure 3. Top and side displacements of the *OH + H + *OH intermediate around the 289.32 cm−1 frequency according to vibrational data at 353.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Figure 3. Top and side displacements of the *OH + H + *OH intermediate around the 289.32 cm−1 frequency according to vibrational data at 353.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Molecules 28 07873 g003
Molecules 28 07873 g004
Figure 4. IR spectra of different temperatures for (a) the catalyst (*), (b) *O2, (c) OH + H2O, and (d) H2O + H2O.
Figure 4. IR spectra of different temperatures for (a) the catalyst (*), (b) *O2, (c) OH + H2O, and (d) H2O + H2O.
Molecules 28 07873 g004
Molecules 28 07873 g005
Figure 5. Optimized intermediate of OH + H2O (*OH + H + *OH) at the temperatures of (a) 298.15 K (b) 353.15 K, and (c) 393.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Figure 5. Optimized intermediate of OH + H2O (*OH + H + *OH) at the temperatures of (a) 298.15 K (b) 353.15 K, and (c) 393.15 K (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Molecules 28 07873 g005
Molecules 28 07873 g006
Figure 6. Optimization process when the H atom is introduced to the first O atom of the OOH intermediate at a temperature of 298.15 K. (a) Initially, the H atom binds to the first O atom of OOH to form H2O2, (b) in the middle state of optimization, and (c) in the final optimization step where dual OH molecules are formed instead of H2O2 (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Figure 6. Optimization process when the H atom is introduced to the first O atom of the OOH intermediate at a temperature of 298.15 K. (a) Initially, the H atom binds to the first O atom of OOH to form H2O2, (b) in the middle state of optimization, and (c) in the final optimization step where dual OH molecules are formed instead of H2O2 (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Molecules 28 07873 g006
Molecules 28 07873 g007
Figure 7. Optimization process when a H atom is introduced to the second O atom of the OOH intermediate at the temperature of 298.15 K: (a) initially, the H atom binds with the second O atom of the OOH intermediate to form H2O; (b) in the next step, the H atom binds to the second O atom, causing the previously bonded H atom to break its bond and migrate to the first O atom; and (c) in the final optimization step, a dual OH molecule is formed instead of H2O (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Figure 7. Optimization process when a H atom is introduced to the second O atom of the OOH intermediate at the temperature of 298.15 K: (a) initially, the H atom binds with the second O atom of the OOH intermediate to form H2O; (b) in the next step, the H atom binds to the second O atom, causing the previously bonded H atom to break its bond and migrate to the first O atom; and (c) in the final optimization step, a dual OH molecule is formed instead of H2O (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Molecules 28 07873 g007
Molecules 28 07873 g008
Figure 8. Top and side view of the optimized ORR steps for the Cu2-N8/Gr structure at 298.15 K. (a) O2, (b) OOH, (c) OH + OH, (d) OH + H + OH, and (e) 2H2O molecules binding to Cu2-N8/Gr (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Figure 8. Top and side view of the optimized ORR steps for the Cu2-N8/Gr structure at 298.15 K. (a) O2, (b) OOH, (c) OH + OH, (d) OH + H + OH, and (e) 2H2O molecules binding to Cu2-N8/Gr (white, grey, blue, orange, and red spheres represent hydrogen, carbon, nitrogen, copper, and oxygen atoms, respectively).
Molecules 28 07873 g008
Molecules 28 07873 g009
Figure 9. Free energy diagrams of Cu2-N8/Gr at different temperatures under various cell potentials in an acidic medium: (a) 298.15 K, (b) 353.15 K, and (c) 393.15 K (note: the OH + H2O intermediate is the same as the OH + H +OH intermediate, as explained previously).
Figure 9. Free energy diagrams of Cu2-N8/Gr at different temperatures under various cell potentials in an acidic medium: (a) 298.15 K, (b) 353.15 K, and (c) 393.15 K (note: the OH + H2O intermediate is the same as the OH + H +OH intermediate, as explained previously).
Molecules 28 07873 g009
Table 1. Formation energies of the Cu2-N8/Gr structure determined at the temperatures (T) of 298.15, 353.15, and 393.15 K under various external potentials (U).
Table 1. Formation energies of the Cu2-N8/Gr structure determined at the temperatures (T) of 298.15, 353.15, and 393.15 K under various external potentials (U).
T/KU/VFormation Energy/eV
298.150−28.856
1−24.856
4−12.856
6−4.856
7.21−0.016
7.220.024
353.150−28.865
1−24.865
4−12.865
6−4.865
7.21−0.025
7.220.015
393.150−28.864
1−24.864
4−12.864
6−4.864
7.21−0.024
7.220.016
Table 2. Energy gap (Eg), chemical hardness (η), chemical potential (μ), and electrophilicity index (ω) values of the O2, OH + H2O, 2H2O intermediates, and defect structures for different temperatures (T).
Table 2. Energy gap (Eg), chemical hardness (η), chemical potential (μ), and electrophilicity index (ω) values of the O2, OH + H2O, 2H2O intermediates, and defect structures for different temperatures (T).
T/KIntermediate/StructureEg/eVη/eVμ/eVω/eV
298.15 Cu2-N8/Gr0.9920.496−3.4696.063
O20.7380.369−3.6158.854
OH + H2O0.9470.473−3.3585.953
2H2O0.6860.343−3.3228.048
353.15 Cu2-N8/Gr0.9920.496−3.4696.063
O20.7380.369−3.6158.854
OH + H2O0.7770.388−3.4187.520
2H2O0.6860.343−3.3228.048
393.15 Cu2-N8/Gr0.9920.496−3.4696.063
O20.7380.369−3.6158.854
OH + H2O1.0390.519−3.4125.605
2H2O0.6860.343−3.3228.048
Table 3. The binding energies (BEs) of the Cu2-N8/Gr structure at different temperatures and the shortest distance (d) between Cu-O and O-O atoms in angstroms (Å).
Table 3. The binding energies (BEs) of the Cu2-N8/Gr structure at different temperatures and the shortest distance (d) between Cu-O and O-O atoms in angstroms (Å).
Temperature/KMoleculeBE/eVdO-O/ÅdCu-O/Å
298.15O2−3.001.4711.893
OOH−9.601.5651.887
OH + OH−10.972.6801.864, 1.809
OH + H2O−5.352.4411.960, 1.939
OH−6.45-1.832
H2O−1.11-2.144
2H2O−2.122.5191.997, 2.909
353.15O2−2.991.4711.893
OOH−9.591.5651.887
OH + OH−8.512.6801.864, 1.809
OH + H2O−5.25-1.943, 1.946
OH−0.25-1.837
H2O−1.09-2.144
2H2O−2.112.5191.997, 2.909
393.15O2−2.981.4711.893
OOH−9.581.5651.887
OH + OH−8.502.6801.864, 1.809
OH + H2O−5.332.4371.946
OH−0.24-1.837
H2O−1.08-2.144
2H2O−2.102.5191.996, 2.908
Table 4. Comparison of binding energy (BE), bond length of molecule to catalyst (dM-O), and formation energy at zero potential (EF(U = 0)) of different catalysts in this study with the published studies, at initial temperature. M-O is defined as M = metal atom in catalyst surface and O = oxygen atom of molecule.
Table 4. Comparison of binding energy (BE), bond length of molecule to catalyst (dM-O), and formation energy at zero potential (EF(U = 0)) of different catalysts in this study with the published studies, at initial temperature. M-O is defined as M = metal atom in catalyst surface and O = oxygen atom of molecule.
CatalystMoleculeBE/eVdM-O/ÅEF(U = 0)/eVReferences
Cu-N2/GrO2−2.902.047−5.68
OOH−1.811.957
O−5.141.890Bhatt et al. [35]
OH−2.681.908
H2O−1.082.203
H2O2−1.492.236
Cu-N4/GrO2−2.761.98−13.22
OOH−1.701.87
O−5.151.88Balasooriya et al. [36]
OH−1.991.84
H2O−1.232.08
H2O2−0.922.18
Co-N4/GrO2−0.671.93−3.56
OOH−1.021.87
O−3.181.71Kattel et al. [28]
OH−2.441.88
H2O2−0.063.36
Cu2-N8/GrO2−3.001.893−28.856Current study
OOH−9.601.887
OH + OH−10.971.864, 1.809
OH + H2O−5.351.960, 1.939
OH−6.451.832
H2O−1.112.144
2H2O−2.121.997, 2.909
H2O2None formingNone forming
Ti-N4/GrO2−8.951.838−3.72
OOH−7.991.829
O−12.061.644Bhatt et al. [35]
OH−7.991.771
H2O−4.142.067
H2O2−9.131.829
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Balasooriya, Y.; Samarasekara, P.; Lim, C.M.; Chau, Y.-F.C.; Kooh, M.R.R.; Thotagamuge, R. Density Functional Theory Investigation of Temperature-Dependent Properties of Cu-Nitrogen-Doped Graphene as a Cathode Material in Fuel Cell Applications. Molecules 2023, 28, 7873. https://doi.org/10.3390/molecules28237873

AMA Style

Balasooriya Y, Samarasekara P, Lim CM, Chau Y-FC, Kooh MRR, Thotagamuge R. Density Functional Theory Investigation of Temperature-Dependent Properties of Cu-Nitrogen-Doped Graphene as a Cathode Material in Fuel Cell Applications. Molecules. 2023; 28(23):7873. https://doi.org/10.3390/molecules28237873

Chicago/Turabian Style

Balasooriya, Yashas, Pubudu Samarasekara, Chee Ming Lim, Yuan-Fong Chou Chau, Muhammad Raziq Rahimi Kooh, and Roshan Thotagamuge. 2023. "Density Functional Theory Investigation of Temperature-Dependent Properties of Cu-Nitrogen-Doped Graphene as a Cathode Material in Fuel Cell Applications" Molecules 28, no. 23: 7873. https://doi.org/10.3390/molecules28237873

Article Metrics

Back to TopTop