Bao-Chen Qiana,
Xiao Wanga,
Qi Wanga,
Xiao-Qing Zhu*b and
Guang-Bin Shen*a
aCollege of Medical Engineering, Jining Medical University, Jining, Shandong 272000, P. R. China. E-mail: gbshen@mail.jnmc.edu.cn
bThe State Key Laboratory of Elemento-Organic Chemistry, Department of Chemistry, Nankai University, Tianjin 300071, China. E-mail: xqzhu@nankai.edu.cn
First published on 2nd January 2024
N-heterocycles are important chemical hydrogen-storage materials, and the acceptorless dehydrogenation and hydrogenation of N-heterocycles as organic hydrogen carriers have been widely studied, with the main focus on the catalyst synthesis and design, investigation of the redox mechanisms, and extension of substrate scope. In this work, the Gibbs free energies of the dehydrogenation of pre-aromatic N-heterocycles (YH2) and the hydrogenation of aromatic N-heterocycles (Y), i.e., ΔGH2R(YH2) and ΔGH2A(Y), were derived by constructing thermodynamic cycles using Hess' law. The thermodynamic abilities for the acceptorless dehydrogenation and hydrogenation of 78 pre-aromatic N-heterocycles (YH2) and related 78 aromatic N-heterocycles (Y) were well evaluated and discussed in acetonitrile. Moreover, the applications of the two thermodynamic parameters in identifying pre-aromatic N-heterocycles possessing reversible dehydrogenation and hydrogenation properties and the selection of the pre-aromatic N-heterocyclic hydrogen reductants in catalytic hydrogenation were considered and are discussed in detail. Undoubtedly, this work focuses on two new thermodynamic parameters of pre-aromatic and aromatic N-heterocycles, namely ΔGH2R(YH2) and ΔGH2A(Y), which are important supplements to our previous work to offer precise insights into the chemical hydrogen storage and hydrogenation reactions of pre-aromatic and aromatic N-heterocycles.
Since the chemical processes for the acceptorless dehydrogenation and hydrogenation of N-heterocycles involve H2 release and acceptance, therefore, the thermodynamics of the dehydrogenation and hydrogenation of N-heterocycles are important thermodynamic parameters to evaluate the hydrogen-storage abilities and hydrogenation abilities of N-heterocycles, especially for pre-aromatic N-heterocycles possessing reversible dehydrogenation and hydrogenation abilities.3–8 This work follows on from our previous study into the thermodynamics of hydrogen transfer for amines.20 This present work focuses on thermodynamic evaluations of the acceptorless dehydrogenation and hydrogenation of a special category of N-heterocycles, namely, pre-aromatic and aromatic N-heterocycles, in acetonitrile.
In our previous research work, we computed the Gibbs free energies of 84 amines (YH2), including 78 pre-aromatic N-heterocycles and 6 general amines, releasing hydrides and their activation free energies for hydride self-exchange reactions, i.e., ΔGH−R(YH2) and ΔG‡(YH2/YH+), using density functional theory (DFT) calculations in acetonitrile.20a Earlier in 2023, thermodynamic evaluations were conducted of 84 amines releasing two hydrogen ions (H− + H+) and the corresponding 84 imines accepting two hydrogen ions (H− + H+), including ΔGH−PR(YH2) and ΔGH−PA(Y).20b Combined with the above thermodynamic data, in this work, the thermodynamic values upon H2 release and acceptance for the dehydrogenation and hydrogenation of pre-aromatic and aromatic N-heterocycles (Scheme 1) were naturally derived via constructing thermodynamic cycles using Hess' law,21 and the thermodynamic abilities of the pre-aromatic N-heterocycles as chemical hydrogen-storage materials were well evaluated and compared. Moreover, application of the thermodynamic data for identifying pre-aromatic N-heterocycles possessing reversible dehydrogenation and hydrogenation abilities3–8 were investigated and discussed in detail. Undoubtedly, this work focuses on new thermodynamic parameters, which is an important supplement to our previous work to offer precise insights into the chemical hydrogen storage and hydrogenation reactions of pre-aromatic N-heterocycles.1–19
Scheme 1 Chemical equations for the acceptorless dehydrogenation and hydrogenation of pre-aromatic and aromatic N-heterocycles. |
Scheme 2 Chemical structures of 84 amines (YH2), including the 78 pre-aromatic N-heterocycles (1H2–78H2) and 6 general amines (79H2–84H2) investigated in this work. |
To obtain the Gibbs free energy of YH2 releasing H2, the constructed thermodynamic cycle based on the processes of YH2 releasing two hydrogen ions or H2 in acetonitrile is shown in Scheme 3. As can be seen from Scheme 3, Step 1 is the chemical process of YH2 releasing H2, YH2 → Y + H2, and the thermodynamic driving force of YH2 dehydrogenation could be described by the Gibbs free energy of YH2 releasing H2, ΔGH2R(YH2). Step 2 is the chemical process of YH2 releasing two hydrogen ions, YH2 → Y + H− + H+, and the corresponding thermodynamic driving force was defined as the Gibbs free energy of YH2 releasing two hydrogen ions, ΔGH−PR(YH2).20 Step 3 is the chemical process of H− reacting with H+ to form H2, H− + H+ → H2, and the related thermodynamic driving force was defined as the Gibbs free energy of one molar H− reacting with one molar H+ to generate one molar H2, ΔGPA(H−). The ΔGPA(H−) value was reported as −76.0 kcal mol−1 in acetonitrile.22–24 Therefore, the value of ΔGH2R(YH2) could be calculated using eqn (1) in Table 1 by establishing a thermodynamic cycle according to Hess' law (Scheme 3),21 where ΔGH2R(YH2) = ΔGH−PR(YH2) + ΔGPA(H−) (eqn (1)). For eqn (1), ΔGH−PR(YH2) values are available from our previous work and are displayed in the fourth column of Table 2. Since the ΔGH−R(YH2) values were computed using the DFT method with a precision of 1.1 kcal mol−1 in our previous work,20a and the pKa values of YH+ were predicted using XGBoost with a 0.87 pKa error,20a the ΔGH2R(YH2) values were derived from ΔGH−R(YH2) and pKa(YH+), whereby ΔGH2R(YH2) = ΔGH−R(YH2) + 1.37pKa(YH+), and therefore, the precision of the ΔGH2R(YH2) values could be estimated within 2.3 kcal mol−1, which was a suitable precision to give practical guidance on the chemical hydrogen storage and hydrogenation reactions, etc.
Scheme 3 Constructed thermodynamic cycle based on the processes of YH2 releasing two hydrogen ions or H2 in acetonitrile. |
Chemical processes | Thermodynamic parameters | Sources or computed equations | Equation X | |
---|---|---|---|---|
Step 1 | YH2 → Y + H2 | ΔGH2R(YH2) | ΔGH2R(YH2) = ΔGH−PR(YH2) + ΔGPA(H−) | 1 |
Step 2 | YH2 → Y + H− + H+ | ΔGH−PR(YH2) | Ref. 20 | — |
Step 3 | H− + H+ → H2 | ΔGPA(H−) | −76.0 kcal mol−1 (ref. 22–24) | — |
Step 4 | Y + H2 → YH2 | ΔGH2A(Y) | ΔGH2A(Y) = −ΔGH2R(YH2) | 2 |
Compounds | Structures | G | ΔGH−PR(YH2)a | ΔGH2R(YH2)b |
---|---|---|---|---|
−ΔGH−PA(Y)a | −ΔGH2A(Y)c | |||
a ΔGH−PR(YH2) and ΔGH−PA(Y) values are from ref. 20.b ΔGH2R(YH2) = ΔGH−PR(YH2) + ΔGPA(H−).c ΔGH2A(Y) = −ΔGH2R(YH2). | ||||
1H2 | — | 68.0 | −8.0 | |
2H2 | CH3 | 68.8 | −7.2 | |
3H2 | tBu | 70.6 | −5.4 | |
4H2 | Ph | 66.5 | −9.5 | |
5H2 | NH2 | 63.5 | −12.5 | |
6H2 | CHO | 75.6 | −0.4 | |
7H2 | CN | 75.2 | −0.8 | |
8H2 | NO | 78.9 | 2.9 | |
9H2 | NO2 | 71.7 | −4.3 | |
10H2 | CH3 | 69.7 | −6.3 | |
11H2 | tBu | 71.8 | −4.2 | |
12H2 | Ph | 70.2 | −5.8 | |
13H2 | NH2 | 67.8 | −8.2 | |
14H2 | CHO | 80.6 | 4.6 | |
15H2 | CN | 78.7 | 2.7 | |
16H2 | NO | 87.3 | 11.3 | |
17H2 | NO2 | 85.7 | 9.7 | |
18H2 | CH3 | 67.4 | −8.6 | |
19H2 | tBu | 72.0 | −4.0 | |
20H2 | Ph | 64.8 | −11.2 | |
21H2 | NH2 | 60.2 | −15.8 | |
22H2 | CHO | 72.3 | −3.7 | |
23H2 | CN | 74.1 | −1.9 | |
24H2 | NO | 76.6 | 0.6 | |
25H2 | NO2 | 84.6 | 8.6 | |
26H2 | — | 83.1 | 7.1 | |
27H2 | — | 83.3 | 7.3 | |
28H2 | — | 83.4 | 7.4 | |
29H2 | — | 77.9 | 1.9 | |
30H2 | CH3 | 75.8 | −0.2 | |
31H2 | tBu | 79.5 | 3.5 | |
32H2 | Ph | 75.9 | −0.1 | |
33H2 | NH2 | 69.6 | −6.4 | |
34H2 | CHO | 86.0 | 10.0 | |
35H2 | CN | 86.4 | 10.4 | |
36H2 | NO | 91.0 | 15.0 | |
37H2 | NO2 | 81.8 | 5.8 | |
38H2 | CH3 | 76.7 | 0.7 | |
39H2 | tBu | 81.1 | 5.1 | |
40H2 | Ph | 78.8 | 2.8 | |
41H2 | NH2 | 76.8 | 0.8 | |
42H2 | CHO | 89.4 | 13.4 | |
43H2 | CN | 87.3 | 11.3 | |
44H2 | NO | 98.0 | 22.0 | |
45H2 | NO2 | 93.2 | 17.2 | |
46H2 | CH3 | 73.5 | −2.5 | |
47H2 | tBu | 84.9 | 8.9 | |
48H2 | Ph | 74.9 | −1.1 | |
49H2 | CH3 | 76.2 | 0.2 | |
50H2 | tBu | 80.1 | 4.1 | |
51H2 | Ph | 75.9 | −0.1 | |
52H2 | NH2 | 74.7 | −1.3 | |
53H2 | CHO | 79.0 | 3.0 | |
54H2 | CN | 80.6 | 4.6 | |
55H2 | NO | 80.4 | 4.4 | |
56H2 | NO2 | 81.2 | 5.2 | |
57H2 | — | 67.6 | −8.4 | |
58H2 | H | 87.5 | 11.5 | |
59H2 | CH3 | 87.5 | 11.5 | |
60H2 | Ph | 87.3 | 11.3 | |
61H2 | — | 76.3 | 0.3 | |
62H2 | — | 87.6 | 11.6 | |
63H2 | — | 64.5 | −11.5 | |
64H2 | — | 73.7 | −2.3 | |
65H2 | — | 74.6 | −1.4 | |
66H2 | — | 77.6 | 1.6 | |
67H2 | — | 79.6 | 3.6 | |
68H2 | — | 78.4 | 2.4 | |
69H2 | — | 84.1 | 8.1 | |
70H2 | — | 76.6 | 0.6 | |
71H2 | — | 80.8 | 4.8 | |
72H2 | 77.0 | 1.0 | ||
73H2 | — | 84.5 | 8.5 | |
74H2 | — | 68.6 | −7.4 | |
75H2 | — | 85.5 | 9.5 | |
76H2 | — | 78.2 | 2.2 | |
77H2 | — | 82.7 | 6.7 | |
78H2 | — | 76.9 | 0.9 | |
79H2 | — | 86.3 | 10.3 | |
80H2 | — | 87.6 | 11.6 | |
81H2 | — | 91.2 | 15.2 | |
82H2 | — | 93.8 | 17.8 | |
83H2 | — | 93.7 | 17.7 | |
84H2 | — | 94.3 | 18.3 |
Step 4 in Scheme 2 is the chemical process of Y accepting H2 to offer YH2, Y + H2 → YH2, and the corresponding thermodynamic driving force could be described as the Gibbs free energy of Y accepting H2 to generate YH2, ΔGH2A(Y). Since the Y hydrogenation and YH2 dehydrogenation are reverse chemical reactions, the ΔGH2A(Y) value was the opposite of the ΔGH2R(YH2) value,21 i.e., ΔGH2A(Y) = −ΔGH2R(YH2) (eqn (2) in Table 1). Herein, the chemical processes, thermodynamic parameters, and data sources or computed equations of step 1–step 4 for YH2 dehydrogenation and Y hydrogenation in acetonitrile are presented in Table 1, meanwhile, the ΔGH−PR(YH2), ΔGH−PA(Y), ΔGH2R(YH2), and ΔGH2A(Y) values of the 84 amines (YH2) dehydrogenation and their relevant imines (Y) hydrogenation in acetonitrile are listed in Table 2.
Scheme 4 Classifications of the 84 amines (YH2) and imines (Y) and the corresponding 4 groups of YH2 and Y, respectively, considered in this work. |
According to the definition of ΔGH2R(YH2), if the ΔGH2R(YH2) value is more negative than 0, ΔGH2R(YH2) < 0, the chemical process of YH2 releasing H2 is thermodynamically favorable, and YH2 is recognized as a thermodynamically excellent H2 donor.21 In contrast, if the ΔGH2R(YH2) value is greater than 0, ΔGH2R(YH2) > 0, the chemical process of YH2 releasing H2 is thermodynamically unfavorable, and YH2 is not a thermodynamically feasible H2 donor.
As can be seen from Table 2, the ΔGH2R(YH2) scale of the 78 pre-aromatic N-heterocycles (1H2–78H2) investigated in this work ranged from −15.8 kcal mol−1 to 22.0 kcal mol−1, which spanned a very wide thermodynamic range of 37.8 kcal mol−1. Moreover, among the 78 pre-aromatic N-heterocycles (1H2–78H2), 21H2 from YIH2 was thermodynamically the best hydrogen donor or carrier (−15.8 kcal mol−1), even better than HCO2H (−5.9 kcal mol−1),25 while 44H2 from YIH2 was thermodynamically the worst hydrogen donor or carrier (22.0 kcal mol−1). For more refined thermodynamic analysis, the distribution of ΔGH2R(YH2) values for the 78 pre-aromatic N-heterocycles (1H2–78H2) in every 5 kcal mol−1 is clearly shown in Fig. 1 with YH2 amounts as the ordinate and ΔGH2R(YH2) ranges as the abscissa.
Fig. 1 Distributions of ΔGH2R(YH2) ranges for the 78 pre-aromatic N-heterocycles in every 5 kcal mol−1 from −20 kcal mol−1 to 25 kcal mol−1. |
From Fig. 1, several interesting conclusions could be drawn as follows. (1) The distribution of YH2 amounts exhibited an excellent normal distribution, and ΔGH2R(YH2) values of pre-aromatic N-heterocycles ranging from 0 to 5 kcal mol−1 were the most common (29.5%). (2) It was found that the ΔGH2R(YH2) values of 30 YH2 were more negative than 0, meaning that 30 YH2 were thermodynamically feasible H2 donors, and belonged to potential chemical hydrogen-storage materials. While the ΔGH2R(YH2) values of 48 YH2 were greater than 0, indicating that the 48 YH2 were thermodynamically unfeasible H2 donors. (3) It was also discovered that 75 ΔGH2R(YH2) values ranged from −15 kcal mol−1 to 15 kcal mol−1, indicating that ∼96% of the thermodynamic driving forces for the pre-aromatic N-heterocycles dehydrogenations were between −15 kcal mol−1 to 15 kcal mol−1, while less than 4% of the thermodynamic driving forces for the pre-aromatic N-heterocycles dehydrogenation were more negative than −15 kcal mol−1 or greater than 15 kcal mol−1. Accordingly, for an unknown pre-aromatic N-heterocycle, the thermodynamic driving force for YH2 dehydrogenation is generally between −15 kcal mol−1 to 15 kcal mol−1 (∼96% possibilities).
To clearly reveal the relationship between the structural features and thermodynamic abilities, the ΔGH2R(YH2) scales of 4 groups of YH2, along with the ΔGH2R(XH2) values of HCO2H, H2, 67H2, 78H2, and HEH2 (Hantzsch ester, 27H2) releasing H2 in acetonitrile are shown in Scheme 5,25–27 because HCO2H,28 H2,29–33 67H2,34,35 78H2,36,37 and HEH2 (27H2)38–40 are the most frequently-used hydrogen reductants (denoted as XH2) used in the hydrogenation reactions.
Scheme 5 ΔGH2R(YH2) scales of 4 groups of YH2 (YIH2–YIVH2), along with the ΔGH2R(XH2) values of HCO2H, H2, 67H2, 78H2, and HEH2 releasing H2 in acetonitrile (kcal mol−1). |
From Scheme 5, the following valuable conclusions could be made: (1) the ΔGH2R(YIH2) scale covered from −15.8 kcal mol−1 to 22.0 kcal mol−1, which spanned the widest thermodynamic range by 37.8 kcal mol−1 among the 4 groups of amines YIH2–YIVH2; (2) according to the ΔGH2R(YH2) scales of YIIH2 (−11.5 to 8.5 kcal mol−1), YIIIH2 (−7.4 to 9.5 kcal mol−1), and YIVH2 (10.3–18.3 kcal mol−1), the dehydrogenation abilities decreased in the order of YIIH2 ≈ YIIIH2 > YIVH2; (3) in view of the relations between sets, it was found that the ΔGH2R(YH2) values displayed the following regular pattern of {ΔGH2R(YIIH2) ∪ ΔGH2R(YIIIH2) ∪ ΔGH2R(YIVH2)} ⊆ ΔGH2R(YIH2); (4) for the pre-aromatic N-heterocycles, the ΔGH2R(YH2) scales ranged from −15.8 to 22.0 kcal mol−1 for YIH2, from −11.5 to 8.5 kcal mol−1 for YIIH2, and from −7.4 to 9.5 kcal mol−1 for YIIIH2, respectively. Since the ΔGH2R(YH2) scales of YIH2, YIIH2, and YIIIH2 crossed negative and positive values (−15.8 to 22.0 kcal mol−1), it was indicated that not all the acceptorless dehydrogenation of pre-aromatic N-heterocycles is a thermodynamically uphill or downhill process under ambient conditions,3–19 and not all pre-aromatic N-heterocycles are thermodynamically feasible to serve as chemical hydrogen-storage materials. ΔGH2R(YH2) is absolutely an important thermodynamic parameter to guide chemists to discover more potentially excellent chemical hydrogen-storage materials.
Initially, we tried to explain the effects of substituents on ΔGH2R(YH2), but we failed to draw a meaningful conclusion. For example, an electron-withdrawing group could decrease the hydride-donating ability of YH2, ΔGH−R(YH2), while it increased the proton-donating ability of YH+, ΔGPR(YH+). Since ΔGH2R(YH2) is derived from ΔGH−R(YH2) and ΔGPR(YH+) essentially, i.e., ΔGH2R(YH2) = ΔGH−R(YH2) + ΔGPR(YH+) + ΔGPA(H−), for YH2 possessing an electron-withdrawing group, whether the ΔGH2R(YH2) value increases or decreases depends on the D-values of ΔGH−R(YH2) decreasing and ΔGPR(YH+) increasing, i.e., ΔΔGH−R(YH2) and ΔΔGPR(YH+). In addition, the steric factors from the substituents will have a significant effect on reaction kinetics, and a slight effect on the reaction thermodynamics.
As for the ΔGH2A(Y) values, we can safely draw some interesting conclusions to inform the hydrogenation reactions of imines, including aromatic N-heterocycles (1–78) and general imines (79–84), from Scheme 6. Specifically: (1) for aromatic N-heterocycles, the ΔGH2A(Y) scales ranged from −22.0 kcal mol−1 to 15.8 kcal mol−1 for YI, from −8.5 to 11.5 kcal mol−1 for YII, and from −9.5 kcal mol−1 to 7.4 kcal mol−1 for YIII, respectively. The ΔGH2A(Y) scale of aromatic N-heterocycles (YI–YIII) ranged from −22.0 kcal mol−1 to 15.8 kcal mol−1, which covered a very large range of 37.8 kcal mol−1. The ΔGH2A(Y) scales indicated that the H2 (0 kcal mol−1),22–24 even the great hydrogen-reductant HCO2H (−5.9 kcal mol−1),25 could not hydrogenate all the aromatic imines (YI, YII, and YIII) to offer pre-aromatic N-heterocycles (YIH2, YIIH2, and YIIIH2). It also could be deduced that some pre-aromatic N-heterocycles from YIH2, YIIH2, and YIIIH2 were thermodynamically better hydrogen reductants than H2 in hydrogenation reactions. In practice, 27H2 (HEH2),38–40 67H2,34,35 73H2,41 74H2,42 and 78H2 (ref. 36 and 37) have already been extensively researched in hydrogenating various unsaturated compounds, including aldehydes ketones, alkenes, imines, and heterocycles;36–42 (2) in contrast, since YIV were general imines without aromatic structures, all the ΔGH2A(YIV) values (−18.3 to −10.3 kcal mol−1) were all greater than the ΔGH2R(H2) values (0 kcal mol−1), meaning that H2 could hydrogenate the related general imines (YIV) in organic synthesis under suitable catalytic conditions, which has been proved by many published studies in the literature.43–47 Similarly, HCO2H (ΔGH2R(HCO2H) = −5.9 kcal mol−1),28 78H2 (ΔGH2R(78H2) = 0.9 kcal mol−1),36,37 67H2 (ΔGH2R(67H2) = 3.6 kcal mol−1),34,35 and HEH2 (ΔGH2R(HEH2) = 7.3 kcal mol−1)38–40 could be applied to hydrogenate YIV to prepare general amines YIVH2 in organic synthesis from a thermodynamics viewpoint.
Scheme 6 ΔGH2A(Y) scales of 4 groups of Y (YI–YIV), along with the ΔGH2R(XH2) values of HCO2H, H2, 67H2, 78H2, and HEH2 in acetonitrile (kcal mol−1). |
Theoretically, according to the reaction rule, for the reversible dehydrogenation and hydrogenation of pre-aromatic and aromatic N-heterocycles, if the Gibbs free energy of a YH2 dehydrogenation (YH2 → Y + H2) is more negative than 0, ΔGH2R(YH2) < 0, then the YH2 dehydrogenation reaction is judged as thermodynamically favorable; while if the Gibbs free energy of the corresponding Y hydrogenation reaction (Y + H2 → YH2) is greater than 0, ΔGH2A(Y) > 0, then the hydrogenation reaction of Y is thermodynamically unfavorable. In contrast, if the Gibbs free energy of a YH2 dehydrogenation (YH2 → Y + H2) is greater than 0, ΔGH2R(YH2) > 0, and the dehydrogenation reaction is considered as thermodynamically unfavorable; while if the Gibbs free energy of the corresponding Y hydrogenation (Y + H2 → YH2) is more negative than 0, ΔGH2A(Y) < 0, and the related Y hydrogenation reaction is regarded as thermodynamically favorable. Therefore, it is very curious why the reversible dehydrogenation and hydrogenation of pre-aromatic and aromatic N-heterocycles could happen, even using the same metal–organic catalyst.3–9
Examining the previous literature, it is believed that the H2 plays a very important role in regulating the reversible dehydrogenation and hydrogenation thermodynamics.3–9 In YH2 dehydrogenation, the H2 is released from the reaction system, which greatly affects the dehydrogenation equilibrium.3–9 Therefore, the Gibbs free energy of YH2 dehydrogenation could decrease to make the reaction happen if the ΔGH2R(YH2) value is not too much greater than 0. Furthermore, in Y hydrogenation, the high H2 pressure would prompt more H2 to dissolve in the reaction solvent, which would have a strong influence on the hydrogenation equilibrium.3–9 Consequently, the Gibbs free energy of Y hydrogenation could decrease to make the reaction occur if the ΔGH2A(Y) is not too much greater than 0. Herein, this leads to the other key question: what are the restrictions on the ΔGH2R(YH2) value for the reversible dehydrogenation and hydrogenation of a pre-aromatic and aromatic N-heterocycle?3
Investigating previous work, the Gibbs solvation energy of H2 was estimated as 3.4 kcal mol−1 in acetonitrile.48 That is, the Gibbs free energy of H2 release from acetonitrile solution was −3.4 kcal mol−1, and the H2 release from the dehydrogenation system could decrease the ΔGH2R(YH2) by 3.4 kcal mol−1. Therefore, it is reasonable to deduce that if the ΔGH2R(YH2) value of a pre-aromatic N-heterocycle is greater than −3.5 kcal mol−1 and more negative than 3.5 kcal mol−1, −3.5 kcal mol−1 < ΔGH2R(YH2) < 3.5 kcal mol−1, the pre-aromatic N-heterocycle could be considered as a potentially reversible chemical organic hydrogen material, which could release H2 and be regenerated by H2 too. According to the above judgment criterion, the ΔGH2R(YH2) values of general amines (YIVH2) (10.3–18.3 kcal mol−1) are much greater than 3.5 kcal mol−1, and so YIVH2 could not be designed as reversible dehydrogenation and hydrogenation materials, which is proved by published works.3–9 Moreover, the ΔGH2R(YH2) scale of 28 pre-aromatic N-heterocycles from YIH2, YIIH2, and YIIIH2 ranged from −3.5 kcal mol−1 to 3.5 kcal mol−1, while the corresponding 28 pre-aromatic N-heterocycles were identified as reversible dehydrogenation and hydrogenation materials. The 28 pre-aromatic N-heterocycles included 6H2 (−0.4 kcal mol−1), 7H2 (−0.8 kcal mol−1), 8H2 (2.9 kcal mol−1), 15H2 (2.7 kcal mol−1), 23H2 (−1.9 kcal mol−1), 24H2 (0.6 kcal mol−1), 29H2 (1.9 kcal mol−1), 30H2 (−0.2 kcal mol−1), 31H2 (3.5 kcal mol−1), 32H2 (−0.1 kcal mol−1), 38H2 (0.7 kcal mol−1), 40H2 (2.8 kcal mol−1), 41H2 (0.8 kcal mol−1), 46H2 (−2.5 kcal mol−1), 48H2 (−1.1 kcal mol−1), 49H2 (0.2 kcal mol−1), 51H2 (−0.1 kcal mol−1), 52H2 (−1.3 kcal mol−1), 53H2 (3.0 kcal mol−1), 61H2 (0.3 kcal mol−1), 64H2 (−2.3 kcal mol−1), 65H2 (−1.4 kcal mol−1), 66H2 (1.6 kcal mol−1), 68H2 (2.4 kcal mol−1), 70H2 (0.6 kcal mol−1), 72H2 (1.0 kcal mol−1), 76H2 (2.2 kcal mol−1), and 78H2 (0.9 kcal mol−1), which need further validation and support in experimental work. Examining the chemical structures of the investigated N-heterocycles in previous literature,3–19 they were not exactly the same pre-aromatic N-heterocycle structure, and most cases involved two H2 molecules release and acceptance from N-heterocycles, which could not provide direct experimental data to validate the thermodynamic model.
If an unsaturated substrate (Sub) is hydrogenated to a reduced unsaturated substrate (SubH2), Sub + H2 → SubH2, whose thermodynamic driving force is defined as the Gibbs free energy of an unsaturated substrate (Sub) accepting H2 to afford SubH2 (equal to the opposite of Gibbs free energy of SubH2 releasing H2 to give Sub, SubH2 → Sub + H2), ΔGH2A(Sub) = −ΔGH2R(SubH2), then the larger the thermodynamic driving force for the pre-aromatic N-heterocycle reductant dehydrogenation, the more favorable the Sub hydrogenation by YH2, i.e., Sub + YH2 → SubH2 + Y.21
Herein, the applications of thermodynamic parameters on choosing suitable catalytic amounts of pre-aromatic N-heterocycle reductants in hydrogenation reactions are displayed in Scheme 8 to aid a clear discussion.
Scheme 8 Applications of the thermodynamic parameters on choosing suitable catalytic amounts of N-heterocyclic hydrogen reductants in hydrogenation reactions. |
It should be noted that not all unsaturated substrates hydrogenation reactions could be designed as organic hydrogen-reductant catalytic reactions. For unsaturated substrates (Sub), if the ΔGH2R(SubH2) values are more negative than −3.5 kcal mol−1, due to the regeneration of aromatic N-heterocycles being impossible (ΔGH2R(YH2) < −3.5 kcal mol−1), the hydrogenation reactions could not be designed as pre-aromatic N-heterocyclic reductant catalyzed reactions. In contrast, if the ΔGH2R(SubH2) values of unsaturated substrates are greater than 3.5 kcal mol−1, whose thermodynamic driving force scale of SubH2 releasing H2 is denoted as ΔGH2R(SubH2)min ∼ ΔGH2R(SubH2)max, then pre-aromatic N-heterocycles with ΔGH2R(YH2) values more negative than ΔGH2R(SubH2)min and greater than −3.5 kcal mol−1 can be recognized as potential catalytic pre-aromatic N-heterocycle reductants, that is, −3.5 kcal mol−1 < ΔGH2R(YH2) < ΔGH2R(SubH2)min. Specially, if 3.5 kcal mol−1 < ΔGH2R(YH2) < ΔGH2R(SubH2)min, then the corresponding pre-aromatic N-heterocycles could be designed as thermodynamically good catalytic hydrogen reductants. If 0 < ΔGH2R(YH2) < 3.5 kcal mol−1, the corresponding pre-aromatic N-heterocycles could be designed as thermodynamically better catalytic hydrogen reductants. If −3.5 kcal mol−1 < ΔGH2R(YH2) < 0 kcal mol−1, the corresponding pre-aromatic N-heterocycles could be designed as the thermodynamically best catalytic hydrogen reductants. If ΔGH2R(YH2) < −3.5 kcal mol−1, the corresponding pre-aromatic N-heterocycles could be considered thermodynamically very excellent reductants, but they could not be regenerated by H2 and designed as catalytic hydrogen reductants.
Thus, it can be seen that the two new thermodynamic parameters ΔGH2R(YH2) and ΔGH2A(Y) can not only help us evaluate and identify good chemical hydrogen-storage materials, but also choose appropriate hydrogen reductants, especially for catalytic amounts of N-heterocyclic hydrogen reductants during the hydrogenation of unsaturated substrates.
(1) ΔGH2R(YH2) is an important thermodynamic parameter to guide chemists to discover more potentially excellent hydrogen carriers. Not all pre-aromatic N-heterocycles are thermodynamically feasible chemical hydrogen-storage materials. The ΔGH2R(YH2) scale of the considered 78 pre-aromatic N-heterocycles ranged from −15.8 kcal mol−1 to 22.0 kcal mol−1, with 30 ΔGH2R(YH2) values being more negative than 0, meaning that the related 30 YH2 were identified as thermodynamically feasible H2 donors, and belong to potential chemical hydrogen-storage materials. The ΔGH2R(YIH2) scale (−15.8 to 22.0 kcal mol−1) spanned the widest thermodynamic range by 37.8 kcal mol−1 among YIH2–YIVH2, and the dehydrogenation abilities of YIIH2–YIVH2 decreased in the order of YIIH2 ≈ YIIIH2 > YIVH2.
(2) As for the thermodynamic abilities of Y hydrogenation, the ΔGH2A(Y) data indicated that for H2, even the great hydrogen-reductant HCO2H (−5.9 kcal mol−1), could not hydrogenate all the aromatic imines (YI, YII, and YIII) to form pre-aromatic N-heterocycles. YIV are general imines without aromatic structures, and all their ΔGH2A(YIV) values (−18.3 to −10.3 kcal mol−1) were greater than ΔGH2R(H2) (0 kcal mol−1), meaning that H2 could hydrogenate the related general imines (YIV) in organic synthesis under suitable catalytic conditions.
(3) H2 plays an important role in regulating reversible dehydrogenation and hydrogenation thermodynamics. The thermodynamic features in the reversible dehydrogenation and hydrogenation of pre-aromatic and aromatic N-heterocycles were clarified such that if −3.5 kcal mol−1 < ΔGH2R(YH2) < 3.5 kcal mol−1, the pre-aromatic N-heterocycle is considered as a thermodynamically potential reversible chemical organic hydrogen material, which could release H2 and be regenerated by H2 too.
(4) The application of thermodynamic parameters to the selection of pre-aromatic N-heterocyclic hydrogen reductants in catalytic hydrogenation was exhibited in this work. If the ΔGH2R(SubH2) values of unsaturated substrates (SubH2) are greater than 3.5 kcal mol−1, pre-aromatic N-heterocycles with −3.5 kcal mol−1 < ΔGH2R(YH2) < ΔGH2R(SubH2)min could be recognized as thermodynamically potential catalytic reductants. Specially, if 3.5 kcal mol−1 < ΔGH2R(YH2) < ΔGH2R(SubH2)min, the related pre-aromatic N-heterocycles could be designed as thermodynamically good catalytic hydrogen reductants. If 0 < ΔGH2R(YH2) < 3.5 kcal mol−1, the related pre-aromatic N-heterocycles could be designed as thermodynamically better catalytic hydrogen reductants. If −3.5 kcal mol−1 < ΔGH2R(YH2) < 0 kcal mol−1, the related pre-aromatic N-heterocycles could be designed as the thermodynamically best catalytic hydrogen reductants.
In summary, this work focused on two new thermodynamic parameters, namely ΔGH2R(YH2) and ΔGH2A(Y), which is an important supplement to our previous work to offer precise insights into the chemical hydrogen storage of pre-aromatic N-heterocycles, and hydrogenation reactions.
This journal is © The Royal Society of Chemistry 2024 |