3.1. Nanocages (C50 and Si50) and their Ni-doped derivatives (Ni-C50 and Ni-Si50)
In this section the properties of nanocages (C50 and Si50) and their Ni-doped derivatives (Ni-C50 and Ni-Si50) are examined. The structures of C50, Si50 and Ni-C50 and Ni-Si50 nanocages are presented in Fig. 1. The energy of adoption (∆Eadoption) of Ni atoms on C50 and Si50 nanocages are calculate as following [41–43] and results are reported in Table 1:
∆Eadoption = ENi−nanocage – Enanocage – ENi (3)
Where the ENi and Enanocage are total energies of C50 and Si50 nanocages and total energy of Ni atom and the ENi−nanocage are the total energies of complexes of Ni atoms with C50 and Si50 nanocages. Results shown that the calculated ∆Eadoption values of Ni-C50 and Ni-Si50 are negative values and so the adoption of Ni atoms on C50 and Si50 nanocages are exothermic processes from thermodynamic viewpoint. Results indicated that the Ni-Si50 nanocage has higher ∆Eadoption values than Ni-C50 nanocage.
The formation energies (∆Eformation) of nanocages (C50 and Si50) and their Ni-doped derivatives (Ni-C50 and Ni-Si50) are calculated as following [44–46] and results are reported in Table 1:
∆Eformation = Enanocage – 50 * EX (4)
∆Eformation = ENi−nanocage – 50 * Enanocage – ENi (5)
Where the ENi−nanocage are the total energies of complexes of Ni atoms with C50 and Si50 nanocages, the ENi and Enanocage are total energies of nanocages and Ni atoms and the EX is the total energies of C and Si atoms. Results shown that the calculated ∆Eformation values of C50, Si50 and Ni-C50 and Ni-Si50 nanocages are negative and valuable amounts.
The ∆Eformation values of Si nanocages are more negative than C nanocages. The derivatives of Ni doped nanocages have higher ∆Eformation values than corresponding nanocages, it can be concluded the Ni atoms can increase the stability of nanocages. The Ni-Si50 has the highest ∆Eformation value than other studied nanoacges. It can be concluded the Ni atoms of metal doped nanocages are catalytic active sites and possible positions to adsorb the CO2 and other possible intermediate species for CO2 reduction reaction.
Table 1. Calculated ∆Eformation and ∆Eadoption values of metal doped nanocages, calculated ∆Eadsorption of possible intermediates of CO2 reduction reactions on metal doped nanocages and the calculated ∆Gactivation, ∆Greaction and overpotential values of reaction steps of possible mechanisms for CO2 reduction reactions on metal doped nanocages.
PW91PW91/6-311+G (2d, 2p)
|
M06-2X/cc-pVQZ
|
Nanocage
|
Eformation
|
Eadoption
|
Eformation
|
Eadoption
|
Ni-C50
|
-4.35
|
-5.65
|
-4.59
|
-5.97
|
Ni-Si50
|
-4.56
|
-5.87
|
-4.82
|
-6.20
|
Model
|
PW91PW91/6-311+G (2d, 2p)
|
M06-2X/cc-pVQZ
|
Nanocage
|
Ni-C50
|
Ni-Si50
|
Ni-C50
|
Ni-Si50
|
Species
|
∆Eadsorption
|
∆Eadsorption
|
∆Eadsorption
|
∆Eadsorption
|
CO2 (a)
|
-0.87
|
-0.94
|
-0.91
|
-0.98
|
CO2 (b)
|
-0.79
|
-0.85
|
-0.82
|
-0.89
|
CO
|
-2.47
|
-2.66
|
-2.58
|
-2.78
|
HCOOH
|
-1.41
|
-1.52
|
-1.47
|
-1.59
|
HCHO
|
-1.76
|
-1.90
|
-1.84
|
-1.98
|
CH3OH
|
-1.52
|
-1.64
|
-1.59
|
-1.71
|
CH4
|
-0.75
|
-0.81
|
-0.78
|
-0.84
|
PW91PW91/6-311+G (2d, 2p)
|
Nanocages
|
Ni-C50
|
Ni-Si50
|
Reaction steps for CO2 reduction
|
∆Gbarrier
|
∆Greaction
|
∆Gbarrier
|
∆Greaction
|
nanocage-*CO2 → nanocage-*COOH
|
0.15
|
-0.25
|
0.14
|
-0.27
|
nanocage-*COOH → nanocage-*CO
|
0.23
|
-0.26
|
0.21
|
-0.28
|
nanocage-*CO → nanocage-*CHO
|
0.12
|
0.24
|
0.11
|
0.26
|
nanocage-*CHO → nanocage-*CH2O
|
0.24
|
0.21
|
0.22
|
0.23
|
nanocage-*CH2O → nanocage-*CH3O
|
0.11
|
-0.38
|
0.10
|
-0.41
|
nanocage-*CH3O → nanocage* + CH3OH
|
0.35
|
-0.58
|
0.32
|
-0.62
|
nanocage-*CH3O → nanocage-*O + CH4
|
0.25
|
-1.21
|
0.23
|
-1.30
|
nanocage-*COOH → nanocage*-HCOOH
|
0.41
|
0.17
|
0.38
|
0.18
|
nanocage-*CO2 → nanocage-*OCHO
|
0.23
|
-0.32
|
0.21
|
-0.34
|
nanocage-*OCHO → nanocage*-HCOOH
|
0.33
|
-0.07
|
0.31
|
-0.08
|
nanocage-*HCOOH → nanocage-*CHO
|
0.45
|
-0.11
|
0.42
|
-0.12
|
nanocage-*CHO → nanocage-*HCHO
|
0.32
|
-0.56
|
0.30
|
-0.60
|
M06-2X/cc-pVQZ
|
Nanocages
|
Ni-C50
|
Ni-Si50
|
Reaction steps for CO2 reduction
|
∆Gbarrier
|
∆Greaction
|
∆Gbarrier
|
∆Greaction
|
nanocage-*CO2 → nanocage-*COOH
|
0.14
|
-0.27
|
0.13
|
-0.29
|
nanocage-*COOH → nanocage-*CO
|
0.21
|
-0.28
|
0.20
|
-0.30
|
nanocage-*CO → nanocage-*CHO
|
0.11
|
0.26
|
0.10
|
0.28
|
nanocage-*CHO → nanocage-*CH2O
|
0.22
|
0.23
|
0.20
|
0.25
|
nanocage-*CH2O → nanocage-*CH3O
|
0.10
|
-0.41
|
0.09
|
-0.44
|
nanocage-*CH3O → nanocage* + CH3OH
|
0.32
|
-0.63
|
0.30
|
-0.68
|
nanocage-*CH3O → nanocage-*O + CH4
|
0.23
|
-1.32
|
0.21
|
-1.42
|
nanocage-*COOH → nanocage*-HCOOH
|
0.38
|
0.18
|
0.35
|
0.20
|
nanocage-*CO2 → nanocage-*OCHO
|
0.21
|
-0.35
|
0.20
|
-0.37
|
nanocage-*OCHO → nanocage*-HCOOH
|
0.30
|
-0.08
|
0.28
|
-0.08
|
nanocage-*HCOOH → nanocage-*CHO
|
0.41
|
-0.12
|
0.38
|
-0.13
|
nanocage-*CHO → nanocage-*HCHO
|
0.29
|
-0.61
|
0.27
|
-0.66
|
Overpotential
|
PW91PW91/6-311+G (2d, 2p)
|
M06-2X/cc-pVQZ
|
nanocages
|
Ni-C50
|
Ni-Si50
|
Ni-C50
|
Ni-Si50
|
reactions to produce the CO
|
0.32
|
0.30
|
0.30
|
0.28
|
reactions to produce the HCOOH
|
0.26
|
0.24
|
0.25
|
0.23
|
reactions to produce the HCHO
|
0.29
|
0.27
|
0.28
|
0.26
|
reactions to produce the CH3OH
|
0.23
|
0.21
|
0.22
|
0.20
|
reactions to produce the CH4
|
0.21
|
0.19
|
0.20
|
0.18
|
|
|
|
|
|
|
|
|
3.2. Reduction reaction of CO2 to CO, CH4, HCOOH, HCHO and CH3OH species
The adsorption of CO2 molecule on surfaces of Ni-C50 and Ni-Si50 nanocages as first step of CO2 reduction reaction is examined. The structures of complexes of possible configuration of CO2 on Ni-C50 and Ni-Si50 nanocages are presented in Fig. 1. The calculated ∆Eadsorption values of CO2 molecule on Ni-C50 and Ni-Si50 nanocages are reported in Table 1.
Results shown that the Ni atoms of Ni-C50 and Ni-Si50 nanocages is acceptable sites to adsorb the CO2 molecules. Results indicated that the CO2 molecules can adsorb to Ni atoms through O and C atoms as reported in Fig. 1. The complexes of CO2 with Ni doped nanocages via C atoms have higher ∆Eadsorption values than corresponding complexes via O atoms, significantly. The Ni-Si50 nanocage has higher ∆Eadsorption values than Ni-C50 nanocage [47–50].
In this section the possible mechanisms and pathways for CO2 reduction reaction to produce the CO, CH4, HCOOH, HCHO and CH3OH species on surfaces of Ni-C50 and Ni-Si50 nanocages are examined. The structures of important intermediates of reaction steps for CO2 reduction reaction on surfaces of Ni-C50 and Ni-Si50 nanocages are presented in Fig. 1. The calculated ΔGreaction values and free barrier energy values of CO2 reduction reaction steps by possible mechanisms to produce the CO, CH4, HCOOH, HCHO and CH3OH species are reported in Table 1.
After adsorption of CO2 molecule on Ni-C50 and Ni-Si50 nanocages the next step of CO2 reduction reaction is the protonation of surfaces to produce the nanocage-*COOH (nanocage-*CO2 + H+ + e– → nanocage-*COOH) and nanocage-*OCHO (nanocage-*CO2 + H+ + e– → nanocage-*OCHO). The free energy barrier of nanocage-*OCHO in Table 1 are higher than corresponding values of nanocage-*COOH on surfaces of Ni-C50 and Ni-Si50 nanocages. The formation of nanocage-*OCHO has more negative ΔGreaction values than nanocage-*COOH creation on surfaces of Ni-C50 and Ni-Si50 nanocages in Table 1. The next steps of CO2 reduction reaction can be processed by reduction of nanocage-*COOH and nanocage-*OCHO to produce the CO, HCOOH, HCHO, CH3OH and CH4 species.
The nanocage-*CO creation is done from hydrogenation of nanocage-*COOH (nanocage-*COOH → nanocage-*CO + H2O) on surfaces of Ni-C50 and Ni-Si50 nanocages. The desorption of CO from nanocage-*CO is need the high free energy and so the efficiency of CO production on surfaces of Ni-C50 and Ni-Si50 nanocages is low as reported in Table 1. Results shown that the nanocage-*COOH can hydrogenated on surfaces of Ni-C50 and Ni-Si50 nanocages and the HCOOH is produced (nanocage-*COOH → nanocage + HCOOH). The HCOOH is also can be created from the hydrogenation of nanocage-*OCHO on surfaces of Ni-C50 and Ni-Si50 nanocages (nanocage-*OCHO → nanocage + HCOOH).
Results demonstrated that for production the HCOOH specie the pathway of nanocage-*CO2 → nanocage-*OCHO → nanocage + HCOOH is more favored than other studied pathway from thermodynamic view point. Results indicate that production of CO and HCOOH on Ni-Si50 nanocage have more negative ΔGreaction values and lower free barrier energy values than Ni-C50 nanocage, significantly as reported in Table 1. Results shown that the nanocage-*COOH → nanocage-*CO + H2O reaction is the rate-limiting step on surfaces of Ni-C50 and Ni-Si50 nanocages.
The CH3OH creation from CO2 reduction reaction on surfaces of Ni-C50 and Ni-Si50 nanocages can be processed by two difference pathways as follow: nanocage-*CO2 → nanocage-*COOH → nanocage-*CO → nanocage-*CHO → nanocage-*CH2O → nanocage-*CH3O → nanocage + CH3OH and the nanocage-*CO2 → nanocage-*OCHO → nanocage + HCOOH → nanocage-*CHO → nanocage-*CH2O → nanocage-*CH3O → nanocage + CH3OH. Results indicated that the nanocage-*CO → nanocage-*CHO is the rate limiting step for CH3OH creation from CO2 reduction reaction on Ni-C50 and Ni-Si50 nanocages. Results shown that the ΔGreaction values of these six reaction steps for CH3OH creation on Ni-Si50 nanocage are more negative than Ni-C50 nanocage as reported in Table 1. The reaction steps of CH3OH creation on Ni-Si50 nanocage have lower free barrier energy values than Ni-C50 nanocage, significantly.
The production of HCHO from CO2 reduction reaction on surfaces of Ni-C50 and Ni-Si50 nanocages can be processed as following by two different pathways: nanocage-*CO2 → nanocage-*COOH → nanocage-*CO → nanocage-*CHO → nanocage + HCHO and other mechanism is nanocage-*CO2 → nanocage-*OCHO → nanocage-*HCOOH → nanocage-*CHO → nanocage + HCHO mechanisms. Results shown that the nanocage-*CO → nanocage-*CHO is the rate-limiting step of production of HCHO on surfaces of Ni-C50 and Ni-Si50 nanocages. Results demonstrated that for production the HCHO specie the pathway of nanocage-*CO2 → nanocage-*COOH → nanocage-*CO → nanocage-*CHO → nanocage + HCHO is more favored than other studied pathway from thermodynamic view point.
The CO2 molecule on surfaces of Ni-C50 and Ni-Si50 nanocages can process the eight reduction reaction steps and finally the CH4 is created as following: nanocage-*CO2 → nanocage-*COOH → nanocage-*CO → nanocage-*CHO → nanocage-*CH2O → nanocage-*CH3O → nanocage-*O + CH4. Results shown that rate limiting step for CH4 production on Ni-C50 and Ni-Si50 nanocages is the nanocage-*CO → nanocage-*CHO. The Ni-Si50 nanocage has more negative ΔGreaction values than Ni-C50 nanocage to process these eight reduction reaction and to create the CH4, significantly.
Results shown that the overpotential for CH4 and CH3OH production are lower than HCOOH and HCHO creation on Ni-C50 and Ni-Si50 nanocages. Results shown that overpotential values to produce the CO, CH4, HCOOH, HCHO and CH3OH species on Ni-Si50 nanocage are lower than corresponding overpotential values on Ni-C50 nanocage. Results indicated that the overpotential of CO2 reduction reaction on Ni-C50 and Ni-Si50 nanocages are lower than corresponding values on various metal catalysts, significantly. It can be concluded that the Ni-C50 and Ni-Si50 nanocages can catalyze the CO2 reduction reaction to produce the CO, CH4, HCOOH, HCHO and CH3OH species with high efficiency.