Shun Lu1,2, Matthew Hummel2, Shuai Kang3, Rajesh Pathak4, Wei He5, Xueqiang Qi1,6, Zhengrong Gu2. 1. College of Chemistry and Chemical Engineering, Chongqing University of Technology, Chongqing 400054, China. 2. Department of Agricultural and Biosystems Engineering, South Dakota State University, Brookings, South Dakota 57007, United States. 3. Micro-nano Manufacturing and System Integration Center, Chongqing Institute of Green and Intelligent Technology (CIGIT), Chinese Academy of Sciences, Chongqing 400714, China. 4. Applied Materials Division, Argonne National Laboratory, 9700 Cass Ave, Lemont, Illinois 60439, United States. 5. Department of Electrical Engineering and Computer Science, South Dakota State University, Brookings, South Dakota 57007, United States. 6. School of Chemistry and Chemical Engineering, Chongqing University, Shazhengjie 174, Chongqing 400044, China.
Abstract
Developing efficient and low-cost urea oxidation reaction (UOR) catalysts is a promising but still challenging task for environment and energy conversion technologies such as wastewater remediation and urea electrolysis. In this work, NiO nanoparticles that incorporated graphene as the NiO@Graphene composite were constructed to study the UOR process in terms of density functional theory. The single-atom model, which differed from the previous heterojunction model, was employed for the adsorption/desorption of urea and CO2 in the alkaline media. As demonstrated from the calculated results, NiO@Graphene prefers to adsorb the hydroxyl group than urea in the initial stage due to the stronger adsorption energy of the hydroxyl group. After NiOOH@Graphene was formed in the alkaline electrolyte, it presents excellent desorption energy of CO2 in the rate-determining step. Electronic density difference and the d band center diagram further confirmed that the Ni(III) species is the most favorable site for urea oxidation while facilitating charge transfer between urea and NiO@Graphene. Moreover, graphene provides a large surface for the incorporation of NiO nanoparticles, enhancing the electron transfer between NiOOH and graphene and promoting the mass transport in the alkaline electrolyte. Notably, this work provides theoretical guidance for the electrochemical urea oxidation work.
Developing efficient and low-cost urea oxidation reaction (UOR) catalysts is a promising but still challenging task for environment and energy conversion technologies such as wastewater remediation and urea electrolysis. In this work, NiO nanoparticles that incorporated graphene as the NiO@Graphenecomposite were constructed to study the UOR process in terms of density functional theory. The single-atom model, which differed from the previous heterojunction model, was employed for the adsorption/desorption of urea and CO2 in the alkaline media. As demonstrated from the calculated results, NiO@Graphene prefers to adsorb the hydroxyl group than urea in the initial stage due to the stronger adsorption energy of the hydroxyl group. After NiOOH@Graphene was formed in the alkaline electrolyte, it presents excellent desorption energy of CO2 in the rate-determining step. Electronic density difference and the d band center diagram further confirmed that the Ni(III) species is the most favorable site for urea oxidation while facilitating charge transfer between urea and NiO@Graphene. Moreover, graphene provides a large surface for the incorporation of NiO nanoparticles, enhancing the electron transfer between NiOOH and graphene and promoting the mass transport in the alkaline electrolyte. Notably, this work provides theoretical guidance for the electrochemical urea oxidation work.
The continuous increase in the energy demand needs to pursuit a
clean and renewable energy source because non-renewable energy sources
such as traditional fossil fuels are limited and lead to global warming.[1,2] It is necessary to bridge the gap between academia and the industry
with extensive research and their practical applications. Although
hydrogen and oxygen evolution reactions (HER and OER, respectively)
have been considered as revolutionary fuel cell designs that utilize
water-splitting technology, they suffer from efficiency drawbacks
and can be realized in limited pristine fresh water and with the use
of noble metal catalysts.[3−5] Urea oxidation reaction (UOR)
is a fundamental step in fulfilling the need for practical green energy
because it does not need a high-voltage supply and also does not release
both O2 and H2 gases simultaneously, encountered
during water splitting.[6,7] Furthermore, urea is an abundant
component of human and animal waste, which can result in the production
of problematic ammonia under normal degradation or standard hydrolysis
practices.[8] More importantly, the UOR process
could provide an opportunity for waste disposal and green hydrogen
production.[9]In a typical UOR, urea in an alkali electrolyte is oxidized to
the production of N2 and CO2 at the anode and
H2 on the cathode from water electrolysis.[10] This process is depicted in eqs –3, respectively.The UOR process is slow and inefficient under normal conditions
due to the six-electron transfer process from the anode to the cathode.[11,12] Thus, it is necessary to modify the working electrode using a catalyst.
Nickel-based materials are considered as one of the most promising
groups of materials for catalysis in UOR owing to their low cost,
easy synthesis route, and abundance in nature.[13,14] For example, Tammam and Saleh developed a NiO-modified electrode
for electrocatalytic urea oxidation in the alkaline media and confirmed
that the UOR process is a completely irreversible diffusion-controlled
route.[15] To improve the conductivity of
Ni-based materials while maintaining the compounds’ catalytic
performance, several effective strategies were applied, including
the introduction of conductive support,[16] elemental doping,[17] high valence Ni-based
materials,[18] and defect engineering.[19] Many of these effective strategies were developed
to realize the commercial implementation of Ni-catalyst driven UOR.
Nevertheless, the in-depth theoretical-fundamental understanding on
the UOR was not studied due to its complicated multistep gas adsorption
and desorption.Moreover, one major theoretical drawback found by density functional
theory calculations is the rate-limiting intermediate step of CO2 desorption during the reaction on the anode.[20] In our previous work, we investigated the use of nano-NiO
supported on eggshell membrane-derived carbon for a Ni-catalyzed UOR,
and the periodic heterojunction model was selected to illustrate the
impact between porous carbon and NiO nanoparticles.[16] NiO also possesses many merits, including easy-to-obtain,
low-cost, and exchangeable valence states.[7,16] Meanwhile,
graphene was employed as the alternative porous carbon for simulations.
It is important to add theoretical simulation to understand the UOR
process using different models, especially for the existence of a
single-atom model due to attractive findings.[21,22] The in-depth mechanism of UOR in NiO@Graphene is still not clear.
Consequently, it is useful to employ the single-atom model to investigate
the role of NiO@Graphene in the UOR process.In this work, a single-atom model was built to understand the relationship
of the NiO@Graphenecomposite and its urea oxidation behavior. The
single-atom model, which differs from the previous theoretical model,
was applied to illustrate the influence between graphene and NiO nanoparticles.
Meanwhile, this work also served as an important theoretical supplement
for the previous research. Prior to the investigation of UOR, the
adsorption of the hydroxyl group and urea on NiO@Graphene was compared.
Then, the adsorption of urea and CO2 on NiO and NiOOH with
graphene was calculated and compared. The electron density difference
map was also used to study the electron transfer of the NiO@Graphenecomposite.
Results and Discussion
Competitive Adsorption of the Urea/Hydroxyl Group on NiO@Graphene
In an earlier literature of UOR toward NiO and its derivatives,
most of the theoretical studies have been performed by considering
the adsorption/desorption of CO2 in the gas–liquid
interface, which matched well with the formed *COO species.[16] However, the competitive adsorption of urea
and the hydroxyl group on the surface of the NiO-based electrocatalyst
cannot be ignored before urea oxidation. Thus, to understand the competitive
adsorption, both the adsorption energy of urea and the hydroxyl group
were calculated. As shown in Scheme , the adsorption energy of the hydroxyl group (route
I) on NiO@Graphene is calculated to be −3.49 eV, which is higher
than that of urea (route II). It means that NiOOH@Graphene was formed
in the alkaline media without the disturbance of urea adsorption under
the applied potential.
Scheme 1
Schematic Illustration of the Adsorption Route of Urea/Hydroxyl on
NiO@Graphene
(I) Adsorption of the hydroxyl
group first; (II) adsorption of urea first (gray for C, red for O,
blue for N, white for H, and light blue for Ni atoms).
Schematic Illustration of the Adsorption Route of Urea/Hydroxyl on
NiO@Graphene
(I) Adsorption of the hydroxyl
group first; (II) adsorption of urea first (gray for C, red for O,
blue for N, white for H, and light blue for Ni atoms).
Theoretical Analysis on NiO@Graphene and NiOOH@Graphene
To get an understanding into the electrocatalytic urea oxidation
mechanism of NiO@Graphene, DFT calculations were utilized based on
the single-atom model. Generally, NiO(II) nanoparticles will be oxidized
into NiOOH(III) in the alkaline environment, which is due to the redox
reaction (eq ) occurring
at the NiO(II) nanoparticles, as shown in Scheme .In a typical cyclic
voltammetry diagram of pure NiO in the presence of the alkaline electrolyte,
the oxidation peak around 0.35 V corresponds to the transformation
of Ni(II) to Ni(III), and the reduction peak around 0.15 V corresponds
to the transformation of Ni(III) to Ni(II).[23] The formed Ni(III) species was regarded as the active site for the
UOR process. Before the research on nickel oxide, nickel hydroxide
was first electrooxidized to NiOOH species in alkaline media (eq ), and then urea molecules
adsorbed on the NiOOH species via bridging coordination, whereby a
Ni atom interacts with a C atom (urea). It means that the onset potential
for the UOR has the potential for the formation of NiOOH via Ni electrooxidation.
The next dissociation of urea on NiOOH is multiple processes, producing
a variety of intermediate species, such as the typical reaction pathway
that was proposed as *CO(NH2)2 → *CO(·HNNH2) → *CO(·HNNH) → *CO(·HNN) →
*CO(N·N) → *CO(OH) → *CO(NH2)2 → *CO(OH·OH) → *COO. By analyzing the Gibbs energy
and the resistance of each step (Table S1), they uncovered that the rate-determining step (RDS) is the desorption
of CO2 from NiOOH species. Thus, to simplify the DFT calculation,
the adsorption of CO2 on the given composite was used to
simulate the RDS process of urea oxidation.Graphene was an alternative instead of porous carbon in our previous
work with many advantages of high conductivity and stable support.[16] It is noted that graphene preferentially anchors
NiO nanoparticles via electrostatic attraction due to its electronegativity.[23] Hence, we considered that the connection of
NiO on the given graphene substrate was achieved by the Ni–O–C
bond. The optimized configurations are listed in Figure , and the corresponding adsorption
energies of urea and CO2 molecules are listed in Table and Table S2.
Figure 1
Optimized structure of (a) NiO@Graphene and (b) NiOOH@Graphene.
The adsorption of urea (c) and CO2 (d) on the surface of
NiO@Graphene and the adsorption of urea (e) and CO2 (f)
on the surface of NiOOH@Graphene (red for O, white for H, gray for
C, blue for N, and light blue for Ni atoms).
Table 1
DFT Data for the Adsorption of CO2 and Urea on the Surface of NiO@Graphene and NiOOH@Graphene,
Respectively
species
optimized
energy (eV)
adsorption
energy (eV)
the shortest
distance between CO2 and the catalyst (Å)
the shortest
distance between urea and the catalyst (Å)
NiO@Graphene-Urea
–11619.6959
–1.377778046
Ni–O: 1.920
O–H: 1.744
NiOOH@Graphene-Urea
–12072.8731
–0.77729845
Ni–O: 1.920
O–H: 1.744
NiO@Graphene-CO2
–10569.1883
–0.628792601
Ni–O: 2.040
NiOOH@Graphene-CO2
–11022.4195
–0.082424595
Ni–O: 2.569
Optimized structure of (a) NiO@Graphene and (b) NiOOH@Graphene.
The adsorption of urea (c) and CO2 (d) on the surface of
NiO@Graphene and the adsorption of urea (e) and CO2 (f)
on the surface of NiOOH@Graphene (red for O, white for H, gray for
C, blue for N, and light blue for Ni atoms).A weaker adsorption energy generally corresponds to a more stable
system. NiO nanoparticles reacted with the hydroxyl group in alkaline
to form NiOOH followed by linking the graphene substrate with the
Ni–O–C bond. In the case of the NiO@Graphenecomposite,
Ni(II) species were supposed to be the active sites, and urea molecules
were attracted on the given composite. The most active site of the
NiO@Graphenecomposite adsorption is the Ni(II) species with an adsorption
energy of −1.37 eV (Figure a). This interaction may originate from the interaction
(Ni–O:urea) between Ni 3d and O 2p orbital electrons of the
urea molecule. However, the adsorption energy of NiOOH@Graphene toward
the urea molecule is −0.77 eV. Ni(III) species played as active
sites for urea oxidation, and the lower urea adsorption energy of
NiOOH@Graphene suggests that it is difficult to adsorb ureacompared
to the performance of NiO@Graphene. The possible reason for this phenomenon
is that the Ni–O:urea bond between NiOOH@Graphene and urea
is affected from the around group (i.e., OH–). To
verify this point, the adsorption energy of NiOOH@Graphene over the
hydroxyl group was calculated to be −2.32 eV, which is a strong
interaction. This result further confirms the above explanation and
also illustrates that urea adsorbed on the surface of the catalyst
can accelerate the electrochemical process to some extent.
Figure 2
Adsorption energies of (a) urea and (b) CO2 molecules
on NiO@Graphene and NiOOH@Graphene.
Adsorption energies of (a) urea and (b) CO2 molecules
on NiO@Graphene and NiOOH@Graphene.Moreover, the desorption step of CO2 from the Ni species
is regarded as the rate-determining step for urea oxidation on the
Ni-based electrocatalysts. Therefore, the adsorption energies of CO2 on NiO@Graphene and NiOOH@Graphene were calculated (Figure b). The CO2 adsorption of NiO@Graphene is −0.62 eV, which is higher than
that of NiOOH@Graphene (−0.08 eV). It may be due to the molecular
interaction between CO2 and the catalyst. So, the distance
between CO2 and the catalyst is an important parameter
to explain this phenomenon. Consequently, the shortest distance between
CO2 and NiO@Graphene is 2.040 Å, which is slightly
smaller than that of NiOOH@Graphene (Figures S1 and S2). This result further confirms that NiOOH@Graphene was
easily desorbed on the CO2 molecule during urea oxidation.
Moreover, Ni(III) species play the active sites for efficient urea
oxidation. This result also is further explained by the electron density
difference diagrams.The electron density difference map is calculated to investigate
the electron transfer between the given sample and urea/CO2 molecules intuitively. It can be observed from Figure a,b that electrons transfer
from the center Ni atom and the adjacent O atom of urea to the intermediate
region between NiO@Graphene or NiOOH@Graphene and urea. There is a
significant charge aggregation between Ni and O atoms (Figure a,b), suggesting the potential
formation of Ni–O covalent bonds. This is due to the charge
transfer and redistribution that will lead to the hybridization of
the Ni 3d and O 2p orbitals. Considering the effect of graphene, the
control experiments were also carried out under the same conditions.
Graphene does not have any contributions to urea/CO2 adsorption
(Figure S3). In the viewpoint of electron
configuration, Ni(II) species ([Ar]3d8) with two unsaturated
d orbitals could be filled by O ([He]2s22p4)
well compared to those of Ni(III) species ([Ar]3d7). Generally,
the stronger the hybridization of Ni-O bond, the stronger the adsorption
toward the urea molecule. NiO@Graphene presents the larger overlapping
area of charge density and with the stronger covalent interaction
(Figure S4). As for CO2 adsorption,
the NiO@Graphenecomposite presents a similar adsorption behavior.
In the whole electrochemical urea oxidation process, the adsorption/desorption
of CO2 is regarded as the key descriptor for the UOR in
an alkaline environment. So, the NiOOH@Graphenecomposite performs
a better desorption behavior than NiO@Graphene (Figure S5).
Figure 3
Electron density difference of the urea molecule adsorbed on (a)
NiO@Graphene and (b) NiOOH@Graphene and CO2 adsorbed on
(c) NiO@Graphene and (d) NiOOH@Graphene. The red hooded face means
the enrichment of electrons, while the blue one means the deficiency
of electrons.
Electron density difference of the urea molecule adsorbed on (a)
NiO@Graphene and (b) NiOOH@Graphene and CO2 adsorbed on
(c) NiO@Graphene and (d) NiOOH@Graphene. The red hooded face means
the enrichment of electrons, while the blue one means the deficiency
of electrons.To further confirm this conclusion, the slight electron density
difference in Figure shows little bit more electrons from the Ni atom to the O atom of
urea on the NiOOH@Graphenecomposite than those on the NiO@Graphenecomposite, which indicates that the Ni(III) species presented its
favorable active sites in the key step of UOR. This is kept in line
with the adsorption energy results where the NiO@Graphenecomposite
shows favorable adsorption toward urea, and the NiOOH@Graphenecomposite
shows better CO2 desorption performance in the key reaction.
Figure 4
Slice images of the adsorption of the CO2 molecule on
the surface of (a) NiO@Graphene and (b) NiOOH@Graphene and the corresponding
slice of the electron density difference. The contour around the atoms
represents electron accumulation (red) or electron deletion (blue).
Slice images of the adsorption of the CO2 molecule on
the surface of (a) NiO@Graphene and (b) NiOOH@Graphene and the corresponding
slice of the electron density difference. The contour around the atoms
represents electron accumulation (red) or electron deletion (blue).The above results demonstrate that the NiO@Graphenecomposite can
effectively adsorb urea and then form a Ni–O covalent bond
with urea, which is well related with the adsorption/desorption properties
and surface charge density of molecules. After NiOOH@Graphene is formed
in the alkaline media, Ni 3d orbitals bond to the O 2p orbital of
urea near the Fermi level, suggesting the charge transfer between
Ni and O atoms. Furthermore, based on the DOS diagram in NiO@Graphene
and NiOOH@Graphene, the peaks of NiOOH@Graphene become tightened (Figure ). The peaks of Ni
3d in NiOOH@Graphene shift lower energy near the Fermi level compared
to those of NiO@Graphene. The center of the Ni 3d orbital peak moves
from −1.61 eV (NiO@Graphene) to −2.42 eV (NiOOH@Graphene).
Additionally, more electrons flow from the Ni atom to the O atom to
establish a strong Ni–O bond, resulting in C–Oads with a weakened bond energy, thus facilitating the *COO desorption
from urea. Generally, the Fermi level can expose the ability of electron
transfer on the electrocatalyst surface. The larger the Fermi level,
the higher the electron transfer capability. Compared to the Fermi
level of NiO@Graphene (1.61 eV), the much larger Fermi level of NiOOH@Graphene
(2.42 eV) suggests that the Ni(III) sample can significantly improve
the electron transfer ability of NiOOH@Graphene. In addition, the
lower d band center leads to a weaker adsorption for CO2. Based on the above DFT analysis, we conclude that the active electron
density of NiOOH@Graphene was effectively upshifted in the UOR process.
Figure 5
d density of states of Ni in NiO (dotted red line)
and NiOOH (blue line). The Fermi level is set to zero, and the vertical
lines represent the d band center.
d density of states of Ni in NiO (dotted red line)
and NiOOH (blue line). The Fermi level is set to zero, and the vertical
lines represent the d band center.
Comparison of the Heterojunction Model and Single-Atom Model
In our previous work, the heterojunction model based on NiO nanoparticles
and graphene was constructed to simplify the theoretical work. The
heterojunction structure is a general model to simulate the multicomponent
composite. DFT calculations with the CASTEP package were employed
to reveal the effect of biomass-derived porous carbon on the UOR performance
of the NiO@C nanocomposite and uncover the Ni(III) species working
as active sites for UOR. In this process, the NiO@C nanocomposite
just played the intermediate role. It turned into the NiOOH@C nanocomposite
with the assistance of the alkaline media; then, Ni(III) species in
NiOOH@C acted as active sites for efficient urea oxidation. However,
the above model just qualitatively revealed the influence of the porous
structure on the electronic structure of NiO nanoparticles and the
synergistic effect between NiO nanoparticles and porous carbon. From
the observed SEM and TEM images, it is also applicable for the single-atom
model considering the Ni species playing as active sites and the ideal
carbon substrate in the UOR. Furthermore, DFT calculations with the
Dmol3 package were used in this study to investigate the
role of Ni(III) species in the UOR. Impressively, a different phenomenon
was observed where NiO@Graphene has the favorable adsorption of urea.
It indicates that NiO@Graphene turned into NiOOH@Graphene in the alkaline
electrolyte first, and then urea absorbed on NiO@Graphene.Based
on the above analysis, it can be demonstrated that NiO@Graphene shows
favorable adsorption of the hydroxyl group in the first stage and
then turns into NiOOH@Graphene under the alkalineconditions for efficient
urea oxidation. The presence of Ni(III) species and excellent electrical
conductivity of NiOOH@Graphene show better desorption of CO2. Moreover, benefiting from the excellent conductivity of graphene,
electrons transferred from urea to NiOOH@Graphene through the Ni–O:urea
bond easier. Graphene also provides a function for facilitating alkaline
electrolyte diffusion, ensuring the formation of Ni(III) species and
promoting mass transfer effectively. Such a composite structure has
the above merits to guarantee the stability and efficiency of NiO@Graphene
as an efficient UOR electrocatalyst.
Conclusions
In this work, the single-atom model where NiO nanoparticles were
bonded with graphene as the NiO@Graphenecomposite was constructed
for the electrochemical urea oxidation in terms of theoretical view.
DFT calculations showed that NiO nanoparticles dispersed on graphene
provide strong adsorption with the hydroxyl group; then, NiOOH@Graphene
was formed after NiO@Graphene reacted with the hydroxyl group. Compared
to NiO@Graphene, NiOOH@Graphene presents a higher desorption energy
of CO2 molecules in the key rate-determining step. Notably,
the Ni(III) species in NiOOH@Graphene is the most favorable site for
the urea oxidation reaction. Moreover, NiOOH@Graphene not only guarantees
the stability of NiOOH and graphene but also promotes the electron
transfer between NiOOH and graphene. Benefiting from the coupling
effect between the Ni(III) species and graphene, NiO@Graphene can
reach excellent electrocatalytic urea oxidation theoretically. These
studies provide theoretical guidance that NiO@Graphene played the
intermediate role in the urea oxidation process before Ni(III) species
formed in the alkaline electrolyte. NiOOH@Graphene also facilitates
the desorption of CO2 from the catalyst surface for UOR
catalysis. More experimental investigations based on the NiO nanoparticles
and graphene will be verified in the future.
Computational Methods
To understand the origin of the electrocatalytic urea oxidation
mechanism of the NiO@Graphenecomposite, density functional theory
(DFT) calculations were conducted using the Dmol3 package
with the Perdew–Burke–Ernzerhof (PBE) formulation of
the generalized gradient approximation (GGA) program.[24] The adsorption of urea and CO2 on NiO@Graphene
was investigated compared to that on NiOOH@Graphene. The single-atom
structure was selected in this investigation. The core electrons were
treated by DFT semi-core pseudopotentials.[25] The DNP basis set was chosen as it can provide more precision for
hydrogen-involved calculations. The convergence thresholds for energy
change, maximum force, and maximum displacement are set to be 2 ×
10–5 Hartree, 0.004 Hartree Å–1, and 0.005 Å, respectively. A vacuum layer of 15 Å thick
was employed along the z direction to eliminate the
interactions between different surfaces. In this work, the adsorption
energy of urea or CO2 is an important reference point for
determining the activity and stability of a urea electrocatalyst.
Therefore, the adsorption energies of urea or CO2 over
NiO@Graphene and NiOOH@Graphene are calculated according to eq .where Eads is the adsorption energy, Etotal is the total energy for the adsorption state, Eslab is the energy of the optimized surface of C@NiO or
C@NiOOH, and Eadsorbate is the energy
of a single urea or CO2 molecule. So, a more negative Eads in eq implies that the adsorption is thermodynamically more favorable.