For the metal carbonyl species, only the most stable ones are listed, other optimized geometries and feasible spin multiplicities of the reaction can be obtained from the Supporting Information (Table 1 and Fig S1-S4). Three typical adsorption sites are considered, including top sites, bridge sites and a facet site. The calculated adsorption energies are in agreement with previous experiment and calculations [33, 34].
Table 1
Vibrational frequency, total energy ET and stepwise CO adsorption energy Ead (kJ/mol) for the Fe3(CO)n (n = 0–12).
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
11Fe3
|
73.94
|
-3791.41
|
-
|
5Fe3(CO)7
|
12.79
|
-4585.20
|
-132.61
|
9Fe3(CO)1
|
46.02
|
-3904.81
|
-168.79
|
3Fe3(CO)8
|
30.41
|
-4698.60
|
-144.87
|
9Fe3(CO)2
|
18.97
|
-4018.22
|
-187.03
|
1Fe3(CO)9
|
22.92
|
-4811.99
|
-136.06
|
9Fe3(CO)3
|
20.79
|
-4131.62
|
-166.55
|
1Fe3(CO)10
|
20.29
|
-4925.38
|
-152.60
|
7Fe3 (CO)4
|
29.24
|
-4245.02
|
-167.80
|
1Fe3(CO)11
|
7.79
|
-5038.75
|
-66.52
|
7Fe3 (CO)5
|
16.39
|
-4358.42
|
-152.62
|
1Fe3(CO)12
|
32.44
|
-5152.14
|
-139.11
|
5Fe3(CO)6
|
31.15
|
-4471.81
|
-140.83
|
CO
|
2126.25
|
-113.34
|
-
|
3.1 Fe3(CO)n (n = 1–12)
The geometrical parameters of the Fe3(CO)n (n = 1–12) cluster systems are shown in Fig. 1. Vibrational frequency, total energies ET, and stepwise CO adsorption energy for the Fe3 systems are collected in Table 1. According to our study and previous experimental study, the Fe3 electronic ground state is found to be 11A'. Table 1 indicates the following, all CO adsorption steps are exothermic, from n = 1–12 for Fe13, the most CO adsorption energy (-187.63 kJ/mol) is found for n = 2 in Fe3 system. According to research, the most stable geometries tend to occur on one Fe atom to form Fe2-Fe(CO)4 structure, and then on the second Fe atom until it is saturated, last, up to the last structure Fe3(CO)12. According to the calculated results, the bond length of CO molecule is 1.139 Å. It is clear from Fig. 1, after CO molecules adsorbed on these Fe3 clusters, the elongation of C-O are extended with the range from 1.3–4.5%, it can think that CO has been activated. The most adsorption energy structure of Fe3 carbonyl species is 9Fe3(CO)2 with Cs symmetry, with the most stepwise CO adsorption energy by 187.03 kJ/mol.
3.2 Fe2Cu (CO)n (n = 1–11)
According to the theoretical calculation, the geometries and parameters of Fe2Cu (CO)n (n = 1–11) are listed in the Fig. 2 and Table 2. On the basis of the theoretical studies, the ground state of Fe2Cu is 8A', which is slower than 6A' and 10A'' by 30.6 and 46.5 kJ/mol, respectively. We first study CO adsorption behavior on both Fe and Cu tops. In contrast, the corresponding adsorption energy value on one Fe side of Fe2Cu is 188.02 kJ/mol, and the corresponding adsorption energy value on Cu is 121.99 kJ/mol, which is slower than that on Fe atom by 66.03 kJ/mol. The most stable structure of Fe2Cu(CO)n (n = 1–11) is found that the CO coordination tend to occur on one Fe atom, later, on the next Fe atom to constitute Fe(CO)1−3-Cu-Fe2(CO)4, last, up to the saturated CO coordination structure. In the carbonyl species, the 6Fe2Cu(CO)1 with the most Eads by 188.02 kJ/mol, in which, Fe-CO and C-O distances are 1.768 and 1.178 Å, respectively. In 2Fe2Cu(CO)11 complex, the energy of a structure with top site adsorption of CO (doublet state with C2v symmetry) is only 6.71 kJ/mol.
Table 2
Vibrational frequency, total energy ET and stepwise CO adsorption energy Ead (kJ/mol) for the Fe2Cu (CO)n (n = 0–11).
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
8Fe2Cu
|
73.25
|
-4168.25
|
-
|
2Fe2Cu(CO)6
|
27.45
|
-4848.66
|
-146.82
|
6Fe2Cu(CO)1
|
46.49
|
-4281.66
|
-188.02
|
2Fe2Cu(CO)7
|
7.77
|
-4962.05
|
-120.35
|
6Fe2Cu(CO)2
|
37.55
|
-4395.06
|
-163.76
|
2Fe2Cu(CO)8
|
19.98
|
-5075.43
|
-128.27
|
6Fe2Cu(CO)3
|
31.10
|
-4508.46
|
-161.69
|
2Fe2Cu(CO)9
|
20.30
|
-5188.81
|
-114.97
|
4Fe2Cu(CO)4
|
33.98
|
-4621.87
|
-185.43
|
2Fe2Cu(CO)10
|
17.99
|
-5302.17
|
-41.33
|
4Fe2Cu(CO)5
|
22.35
|
-4735.27
|
-155.02
|
2Fe2Cu(CO)11
|
15.38
|
-5415.51
|
-6.71
|
3.3 FeCu2(CO)n (n = 1–10)
Figure 3 shows the most stable geometries of FeCu2(CO)n (n = 1–10). There are several possible ways to evaluate the feature of the FeCu2(CO)n (n = 1–10) structures predicted in Table 3. The BPW91 calculations showed that the FeCu2 ground state is the quintet state, which is slower than the triplet and septet states are 20.74 and 61.99 kJ/mol respectively. Analysis of the entire geometries of FeCu2(CO)n (n = 1–4) shows that the isomer with CO adsorbed on Fe atom have more stable compared to those with CO combine with Cu atoms. For the FeCu2(CO)6, it is found that the CO coordination prefer to take place on different Cu atoms, which has lower energy than that of CO coordination on only one Cu atom. In the carbonyl species, the 1FeCu2(CO)4 has the most Eads by 190.58 kJ/mol.
Table 3
Vibrational frequency, total energy ET and stepwise CO adsorption energy Ead (kJ/mol) for the FeCu2(CO)n (n = 1–10).
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
5FeCu2
|
110.47
|
-4545.09
|
-
|
1FeCu2(CO)6
|
18.47
|
-5225.50
|
-118.18
|
3FeCu2(CO)1
|
49.61
|
-4658.51
|
-189.67
|
1FeCu2(CO)7
|
15.42
|
-5338.85
|
-34.52
|
3FeCu2(CO)2
|
31.99
|
-4771.91
|
-170.23
|
1FeCu2(CO)8
|
14.07
|
-5452.20
|
-44.43
|
1FeCu2(CO)3
|
33.93
|
-4885.32
|
-181.31
|
1FeCu2(CO)9
|
3.01
|
-5565.55
|
-30.41
|
1FeCu2(CO)4
|
14.87
|
-4998.73
|
-190.58
|
1FeCu2(CO)10
|
13.92
|
-5678.90
|
-34.58
|
1FeCu2(CO)5
|
18.27
|
-5112.11
|
-122.57
|
|
|
|
|
3.4 Cu3(CO)n (n = 1–8)
On the basis of calculation, the Cu3 ground state is the doublet state, which is lower than the quartet state 166.96 kJ/mol in energy. For Cu3(CO)n (n = 1–8), most stable structures and the corresponding geometrical parameters are shown in Fig. 4 and Table 4. In the beginning of the chemisorption process, we study CO adsorption behavior on difference sites of Cu3 cluster. For the Cu3CO, the lowest energy state of the complex is a doublet with C2v symmetry, in which CO is adsorbed on the top adsorption site. In the Cu3(CO)2 complex, a CO molecule is adsorbed on the one of the other two Cu atoms of Cu2-CuCO, the lowest energy state of the complex is a doublet with C2v symmetry. For the Cu3(CO)5, the stepwise CO adsorption tend to adsorbed on the top site of last Cu atom of Cu3(CO)4. In the Cu3 carbonyl species, the 2Cu3(CO)1 has the most Eads by 138.93 kJ/mol.
Table 4
Vibrational frequency, total energy ET and stepwise CO adsorption energy Ead (kJ/mol) for the Cu3(CO)n (n = 1–8).
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
Species
|
NIMAG
(cm− 1)
|
ET (hartree)
|
Ead
(kJ/mol)
|
2Cu3
|
105.60
|
-4921.96
|
-
|
2Cu3 (CO)5
|
18.21
|
-5488.83
|
-68.03
|
2Cu3(CO)1
|
38.18
|
-5035.35
|
-138.93
|
2Cu3(CO)6
|
13.00
|
-5602.19
|
-56.90
|
2Cu3(CO)2
|
30.30
|
-5148.72
|
-100.32
|
2Cu3(CO)7
|
6.22
|
-5715.55
|
-64.98
|
2Cu3(CO)3
|
15.38
|
-5262.09
|
-78.72
|
2Cu3(CO)8
|
18.19
|
-5828.90
|
-33.72
|
2Cu3 (CO)4
|
19.54
|
-5375.46
|
-91.61
|
|
|
|
|
3.4 Adsorption energy and molecular orbital properties
The change curve of CO adsorption energy on FexCuy (x + y = 3) are collected in Fig. 5. Overall, the Cu3 system has lower CO adsorption capacity on ground state, while the Fe3 system has a relatively high value of stepwise CO adsorption energy. For the bimetallic clusters, it has been found that the stepwise CO adsorption energy of Fe2Cu is higher than that on Fe3 and Cu3 with n = 1, 4, 5 and 6, that the stepwise CO adsorption energy of FeCu2 is higher than that on Fe3 and Cu3 with n = 1, 3 and 4, respectively. For bimetallic systems with n ≥ 7, the stepwise CO adsorption energy are between Fe3 and Cu3. The largest stepwise CO coordination energy (-190.58 kJ· mol− 1) is found for FeCu2(CO)4.
Frontier molecular orbitals of FexCuy (x + y = 3) cluster systems are collected in Fig. 6 and Supporting Information (Tables S2), where HOMO is the highest occupied molecular orbital of FexCuy (x + y = 3) carbonyl complexes, ΔE is the energy level difference between the HOMO of carbonyl complexes and the LUMO of CO. In the case of the CO adsorbs on the FexCuy (x + y = 3) carbonyl complexes, CO molecule obtains electrons from the HOMO of clusters. Based on the study, the molecular reaction can happen only if the ΔE is less than 0.2205 a.u. (579 kJ/mol) [35]. The LUMO energy of free CO is -0.3855 a.u., while the HOMO energies of FexCuy (x + y = 3) carbonyl complexes are all less than − 0.1650 a.u., it means the transfer of electrons from the HOMO of FexCuy (x + y = 3) carbonyl complexes to the LUMO of CO is possible. In brief, the smaller of the ΔE is, the easier of the electron transfer to CO would be, also, the reactions of CO with bimetallic carbonyl complexes are easier than pure clusters system, except Fe3 cluster with first CO molecular.
Frontier molecular orbital theory predicts that the smaller the energy gap (Eg) between the HOMO and LUMO of reactant involved is, the more reactive activity of the molecular would be [36]. The Eg of FexCuy (x + y = 3) carbonyl complexes are showed in Fig. 7 and Supporting Information (Tables S2). Collectively, we can obtain that bimetallic carbonyl complexes have lower value of Eg than that of corresponding monometallic carbonyl complexes, except Fe2Cu and FeCu2(CO)10. This study also gets a good relation of contacting the values of the adsorption energy with Eg of carbonyl complexes, observing in this case that the smaller Eg is, on the contrary, the greater adsorption energy is. The lowest Eg of all carbonyl species is 2Fe2Cu(CO)8 with the value of 2.63 kJ/mol.