Next Article in Journal
Modeling and Optimization of Ellagic Acid from Chebulae Fructus Using Response Surface Methodology Coupled with Artificial Neural Network
Next Article in Special Issue
Kekulé Counts, Clar Numbers, and ZZ Polynomials for All Isomers of (5,6)-Fullerenes C52–C70
Previous Article in Journal
Castanea sativa Mill. By-Products: Investigation of Potential Anti-Inflammatory Effects in Human Intestinal Epithelial Cells
Previous Article in Special Issue
Triel Bonds between BH3/C5H4BX and M(MDA)2 (X = H, CN, F, CH3, NH2; M = Ni, Pd, Pt, MDA = Enolated Malondialdehyde) and Group 10 Transition Metal Electron Donors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Investigation of the Pyridinium-Inspired Catalytic Dehydration of Heptafluoro-Iso-Butyramide for the Synthesis of Environmentally Friendly Insulating Gas Heptafluoro-Iso-Butyronitrile

College of Chemistry and Molecular Science, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(16), 3952; https://doi.org/10.3390/molecules29163952
Submission received: 1 August 2024 / Revised: 16 August 2024 / Accepted: 20 August 2024 / Published: 21 August 2024

Abstract

:
Heptafluoro-iso-butyronitrile (i-C3F7CN) represents a feasible eco-friendly replacement gas for the most potent greenhouse gas sulfur hexafluoride in various high-voltage power transmission equipment. The reaction mechanisms for the in situ synthesis of i-C3F7CN from heptafluoro-iso-butyramide [i-C3F7C(O)NH2] in the presence of trifluoroacetic anhydride (TFAA) and pyridine (Py) in dimethylformamide solution have been studied within density functional theory with M06-2X exchange–correlation functional with the 6-311++G(d,p) basis set and the high-level ab initio complete basis set quadratic CBS-QB3 method. It is revealed that the unimolecular dehydration of i-C3F7C(O)NH2 can be catalyzed efficiently by TFAA in terms of both kinetic and thermodynamic aspects, producing i-C3F7CN and trifluoroacetic acid (TFA). Furthermore, Py is capable of reducing the energy barrier of the rate-determining step through hydrogen abstraction to form pyridinium hydrogen. The synergic effect of the TFAA/Py co-catalyst plays a pivotal role in the production of i-C3F7CN as the Gibbs free energy barrier can be lowered by more than 40 kcal/mol with the ratio of TFAA:2Py, in accordance with the experimental observation. The present theoretical work provides new insights into the rational design on the novel catalysts for large-scale synthesis of the perfluorinated nitriles.

Graphical Abstract

1. Introduction

Sulfur hexafluoride (SF6) has been widely used for decades in a variety of high-voltage power transmission equipment such as gas-insulated lines and switchgears due to its excellent dielectric and interruption performance [1,2]. However, SF6 is the most potent greenhouse gas with a global warming potential (GWP) of 24,300. Various restrictions on SF6 emissions have been implemented since the Kyoto Protocol [3]. To mitigate the greenhouse effects of SF6, two strategies have been considered: blending SF6 with buffer gases such as N2 and searching for new alternative gases [4,5]. Since the cumulative greenhouse gas effect remains as significant for the former, considerable attention has been paid to the identification of an eco-friendly replacement gas for SF6.
Recently, a novel gaseous compound heptafluoro-iso-butyronitrile (i-C3F7CN, abbr. C4) has been developed to replace SF6 in view of its excellent dielectric performance [6]. The dielectric strength of C4 is about twice of SF6 and its GWP is less than 10% of SF6. After mixing with a large amount of CO2, for instance, the 5%C4 + 95%CO2 gaseous admixture exhibits the same dielectric strength as SF6 and the liquefaction temperature is lowered to about −30 °C. Therefore, C4 appears to be the most promising alternative to SF6 as implemented in various high-voltage apparatus and thus the large-scale synthesis of C4 becomes critical for practical use [7,8].
Nitrile compounds can be obtained through ammoxidation, halobenzene cyanidation, diazonium salts, aldoxime dehydration, and so on [9,10,11,12]. However, under the mild experimental conditions, the dehydration of amides outperforms most other approaches with a high yield of nitriles [13]. C4 has been synthesized successfully via the dehydration of heptafluoro-iso-butyramide, i-C3F7C(O)NH2, as prepared by the reaction of i-C3F7COOH with NH3, in the presence of trifluoroacetic anhydride (TFAA) and pyridine (Py) [14]. In fact, no C4 could be produced by the unimolecular dehydration of i-C3F7C(O)NH2 in the dimethylformamide (DMF) solution. Adding Py to the solution does not make any difference. Even with the catalysis of TFAA, the yield of C4 is only 15% [14]. Interestingly, if both TFAA and Py are presented during dehydration, the yield of C4 can increase abruptly up to 86% with the TFAA to Py ratio of 1:2. Although the TFAA/Py-catalyzed dehydration of amides has been used extensively in the preparation of various nitriles, for instance, the aromatic and alkyl nitriles and their perfluorinated compounds [15], the detailed mechanism from amide to nitrile is still elusive especially in view of the synergic catalysis effect between TFAA and Py.
In contrast to the copper-catalyzed dehydration of the primary amides to nitriles [16], the metal-free catalysis of amide to nitrile involves the significant mechanisms as illustrated in Scheme 1. Herein, the TFAA/Py-catalyzed dehydration reactions of i-C3F7C(O)NH2 have been investigated systematically using density functional theory (DFT) and high-level ab initio methods. For the sake of comparison, five comparative reactions were examined, namely, i-C3F7C(O)NH2, i-C3F7C(O)NH2 + 2Py, i-C3F7C(O)NH2 + TFAA, i-C3F7C(O)NH2 + TFAA + 1Py, and i-C3F7C(O)NH2 + TFAA + 2Py, in order to clarify the respective role of TFAA and Py catalysis and the synergic effect, as observed experimentally [14] The present theoretical work does not only provide crucial insights and understanding on the experimental observations but also might stimulate the rational design on the industrial-scale preparation of C4 in the future.

2. Results and Discussion

2.1. The i-C3F7C(O)NH2→C4 + H2O Reaction

The Gibbs free energy profiles of the dehydration of i-C3F7C(O)NH2 to form C4 and H2O are shown in Scheme 2. The optimized geometries for the stationary points are depicted in Figure 1. Evidently, the formation of C4 involves a stepwise H migration mechanism to eliminate the molecular water from i-C3F7C(O)NH2. One of the hydrogen atoms of the amino group is shifted to the carbonyl oxygen atom via a four-member-ring transition state TS1. The breaking N–H bond is elongated to 1.340 Å and the forming O–H bond is 1.317 Å, indicating a late barrier for TS1. The barrier height for H migration is calculated to be 50.9 kcal/mol and the formation of IM1 is endothermic by 15.7 kcal/mol.
Before the second H migration, a conformational change has to be made for IM1 to form IM3, in which the HNCO and HOCN moieties are in cis and trans geometries, respectively. Two energetic routes are presented via TS2/TS3 or TS4/TS5. One is that the terminal H atom of the amino group is rotated towards the O atom via TS2 to form IM2, followed by the internal rotation of the hydroxyl group around the C–O bond via TS3. The rotating N–H bond is shortened by 0.023 Å, and the NCO angle increases from 112.0° to 179.0° in TS2, tending to be a linear HNC geometry while the NCOH dihedral angle is 26.2°. The formation of IM2 is endothermic by 17.9 kcal/mol. As shown in Figure 1, the rotating O–H bond in TS3 remains unchanged but the NCOH dihedral angle changes to −93.7° in TS3. It is evident that the H–N rotation occurs more hardly than the H–O rotation since the energy barrier for TS3 is roughly 12 kcal/mol lower than that for TS2. Indeed, the other IM1→IM3 process takes place firstly via the rotation of the H–O bond to form IM4 by surmounting a barrier of only 7.5 kcal/mol (TS4), which is nearly identical to that for TS3. However, the subsequent H–N rotation involves a significant barrier TS5, which is 4.6 kcal/mol higher in energy than TS2. Therefore, the IM1→IM3 transformation should prefer the first mechanism.
The formation of IM3 is endothermic by 18.7 kcal/mol, as indicative of the unstable characteristics of imine relative to amine. Starting from IM3, the nitrile could be formed via the 1,3-H migration from N to O and simultaneously the molecular H2O is eliminated via the C–O bond breaking. The corresponding transition state is TS6, as shown in Figure 1, with a four-member-ring moiety. The breaking N–H bond is elongated to 1.458 Å and the forming H–O is shortened to 1.197 Å. The C–O bond is stretched to 1.566 Å. Meanwhile, the C–C–N bond angle increases considerably from 118.0 degrees in IM3 to 138.7 degrees in TS6, leading to the formation of the linear C–C≡N triple bond. The barrier height for TS6 was calculated to be as high as 63.1 kcal/mol, which is 12.2 kcal/mol higher than that for TS1. As a result, the conversion of RNH2 to C4 should be determined predominantly by the second H migration process as the rate-determining step. In addition, in the thermodynamic point of view, the RNH2→C4 reaction is endothermic by 20.3 kcal/mol, indicating that the unimolecular dehydration of RNH2 cannot occur spontaneously under ambient experimental condition, even though a weakly bound complex (IM5) might exist between RCN and H2O. In fact, no product has been observed experimentally for RNH2 in a DMF solution at 0 °C [14].

2.2. The i-C3F7C(O)NH2 + 2Py → C4 + H2O + 2Py Reaction

The free energy profile for the dehydration reaction of i-C3F7C(O)NH2 in the presence of two Py molecules in DMF solution is predicted in Scheme 3. The optimized geometries for the energetically preferable reaction paths are shown in Figure 2.
Rather than a co-solvent or a spectator, it is evident that both Py molecules take part in the dehydration process. A pre-reactive complex exists between the NH2 group of i-C3F7C(O)NH2 and two Py molecules due to the N–H...N hydrogen bonds. The binding energy of IM6 was calculated to be −12.2 kcal/mol but the relative Gibbs free energy was 6.6 kcal/mol above the reactants. Subsequently, the Py-assistant H migration transits via TS7 to form IM7. As can be seen in Figure 2, the shifting proton is essentially bonded to the N site of Py, namely, forming pyridinium, as indicated by the newly formed N–H bond of 1.036 Å in length. The breaking N–H and the forming O–H bonds are elongated significantly to 1.999 Å and 1.958 Å, respectively. In the formation of IM7, the proton is transferred back to carbonyl O to form the O–H bond, and thus the Py molecule is regenerated from pyridinium. The Py...HO hydrogen bond appears to be fairly strong in view of the N...H distance of as short as 1.552 Å, which is much shorter than that in the Py...HN hydrogen bond. As a result, IM7 corresponds to a Py-solvated imine. Such a pyridinium-inspired proton transfer path is energetically more preferable since the barrier height for TS7 is 25.2 kcal/mol, which is lowered by 25.7 kcal/mol with respect to that for TS1 in the unassisted dehydration. Although the formation of IM7 is endothermic by 18.2 kcal/mol, the energy of IM7 is only 11.6 kcal/mol higher than that of IM6, and the stabilization due to hydrogen bonding accounts for about 4 kcal/mol.
For the sake of H2O elimination, IM7 should transform to IM8 by rotating the remaining N–H bond via TS8. The reacting H–N=C moiety exhibits a nearly linear geometry. The Py molecule remains neutral as the N...HN hydrogen bond is 1.784 Å. The barrier height for TS8 is 18.8 kcal/mol, which is about 1 kcal/mol lower than that for TS2, implying that the Py molecules do not exert any influence on the internal rotation of N–H. The formation of IM8 is almost thermal neutral. It is interesting to note that the cis–cis H–N=C–O–H conformation of IM8 is considerably distorted due to the steric hindrance of two Py molecules. Therefore, the secondary proton transfer occurs straightforwardly starting from IM8 without the conformational conversion. The corresponding transition state is found to be TS9, which is a pyridinium-inspired path. The H atom is attached to the N site of Py by breaking the N–H bond to form pyridinium, which seems still to interact with the O site, as indicated by the H...O bond length of 1.792 Å. The barrier height for TS9 was calculated to be 35.2 kcal/mol with respect to IM8. Rather than the cleavage of the molecular H2O, the isomer IM9 appears to be a complex between pyridinium and RC(=N)–OH species, of which the energy is only about 2 kcal/mol below TS9. Finally, the H atom is kicked out of the pyridinium moiety to form the 2Py coordinated H2O via TS10, as the breaking C–O bond is elongated by 0.114 Å. The barrier height for such a proton-transfer path is negligible in view of the instability of IM9.
As shown in Scheme 3, a detour to C4 from IM7 also exists. The reaction takes place via the successive internal rotations, i.e., IM7→TS11→IM10→TS12→IM11, followed by the proton migration and the H2O elimination, i.e., IM11→TS13→IM12→TS14→C4. In terms of the energetic routes, the barrier height for TS13 is 52.4 kcal/mol, which is 17.2 kcal/mol higher than that for TS9. Therefore, such a detour to C4 is negligible and the production of C4 should be dominated by TS9.
In summary, two Py molecules act as the catalyst to facilitate the unimolecular reaction of i-C3F7C(O)NH2 through the pyridinium mechanism. The basic N site of the Py molecule is not only the acceptor for hydrogen bonds to NH2 and C=O groups but also the relay for the proton migration. The presence of pyridinium leads to a significant reduction in the barrier for the rate-determining steps from 63.1 kcal/mol to 35.2 kcal/mol. However, such an activation energy is still too high to be important for the production of C4. Indeed, the yield of C4 in the amine/2Py system is negligible.

2.3. The i-C3F7C(O)NH2 + TFAA→C4 + 2TFA Reaction

TFAA has been an efficient compound and widely used for dehydration in organic chemistry. The potential energy profiles for the reaction of i-C3F7C(O)NH2 with TFAA are illustrated in Scheme 4. The optimized geometries for some important stationary points are shown in Figure 3. Evidently, the reaction undergoes via multiple stepwise addition/elimination mechanisms. It is worth noting that the overall reaction becomes almost thermo-neutral and thus the formation of C4 is more thermodynamically feasible due to the catalysis of TFAA.
After a weak hydrogen-bond pre-reactive complex, the reaction of i-C3F7C(O)NH2 with TFAA takes place via the addition of the carbonyl C of TFAA to the carbonyl O of i-C3F7C(O)NH2. A total of three pathways have been found. The energetically most feasible path transits by surmounting a six-member-ring barrier TS15, leading to imine adduct IM14. The forming C–O bond is 1.497 Å and the C=O double bond in amine is stretched to 1.310 Å. Simultaneously, one of the H atoms of the NH2 group is shifted from the N atom to the adjacent carbonyl O atom of TFAA. The breaking N–H bond and the forming O–H bond are 1.208 Å and 1.268 Å, respectively. Meanwhile, the C=N double is formed in IM14. The barrier height for TS15 is 24.7 kcal/mol with respect to the initial reactants and the formation of IM14 is endothermic by 22.1 kcal/mol.
Instead of the intermediate route, the other two reaction routes proceed to form the esterified imine (IM15 in Scheme 4) by releasing one trifluoroacetic acid (TFA) via TS16 or TS17. The geometry of TS16 appears to be similar to that of TS15 but the reacting H atom is shifted to the central O atom of TFAA as accompanied by breaking the C–O bond (e.g., 2.196 Å). The energy difference between these two six-center barriers is 10.5 kcal/mol, excluding any significance of TS16. In contrast, the shifting H atom in TS17 is bonding to the carbonyl O atom of TFAA and the C–O addition occurs on the C atom of the other carbonyl group, resulting in a flexible eight-member-ring geometry. The central C–O bond is elongated to 1.899 Å to give the molecular TFA. The energy of TS17 is 3.3 kcal/mol higher than that of TS15.
The formation of IM15 and TFA is less endothermic than that of the adduct IM14. It should be noted that a feasible path was found for the decomposition of IM14 to form IM15. The H atom of the OH group is shifted to the terminal carbonyl O accompanied by the simultaneous C–O bond cleavage, leading to IM15 and TFA via TS18. The breaking O–H bond is as long as 1.314 Å and the new O–H bond appears to be stretched only slightly (e.g., 1.110 Å). Meanwhile, the breaking C–O bond is elongated by about 25%. Therefore, TS18 is a product-like late barrier, in accordance with the exothermicity of the IM14→IM15 + TFA reaction. The barrier height for TS18 was calculated to be 12.3 kcal/mol. Besides TS18, a detour to IM15 for IM14 involves TS19, which is a four-center barrier for H migration and C–O bond fission. However, the energy of TS19 is about 27 kcal/mol higher than that of TS18, and thus the significance of this detour can be ruled out.
As shown in Figure 3, IM14 is an adduct with intramolecular hydrogen bond occurring between the N atom and the newly formed H–O bond, representing the most stable conformer. The remaining H atom on N site could be rotated to the cis orientation via TS20 or TS21, forming the less stable conformer IM16 or IM17, where the hydrogen bond changes to the N–H...O geometry. In view of the nearly linear C=N–H geometry, the barriers for the internal rotations are around 20 kcal/mol, even higher than that for the direct decomposition of IM14 via TS18.
It was found that both IM16 and IM17 can decompose to C4 and dimeric TFA through three-body rupture. The barriers TS22 and TS23 involve the eight- and six-member-ring reacting moieties, respectively, in which the proton transfer from N to O occurs along with two C–O bond fissions. The corresponding barrier heights are 26.2 kcal/mol and 31.9 kcal/mol, respectively. Although the formation of C4 and TFA dimer is exothermic with respect to the IM14 adduct, the significant barriers preclude any importance of such a mechanism for the production of C4.
Alternatively, the secondary decomposition of IM15 seems to be more favorable in energy, as demonstrated in Scheme 4. Firstly, the H atom on the N site is rotated to form the less stable conformation IM18 via TS24. Secondly, the H atom migrates to the O site, and the C–O bond cleavage takes place simultaneously via TS25, forming C4 after TFA elimination. The barrier height for TS25 is only 26.5 kcal/mol. Thus, it is conceivable that the production of C4 can be promoted significantly in the presence of TFAA, together with the improved thermodynamic equilibrium. In fact, the 15% yield of C4 has been observed experimentally in the reaction of i-C3F7C(O)NH2 with TFAA [14].

2.4. The i-C3F7C(O)NH2 + TFAA + 1Py→C4 + TFA-Py + TFA Reaction

In consideration of the positive role of the pyridinium mechanism to the production of C4, it is worthwhile investigating the potential significance of the Py molecule as a co-catalyst to the TFAA-mediated dehydration of i-C3F7C(O)NH2. The reaction of i-C3F7C(O)NH2 with the 1:1 ratio of TFAA to Py has been calculated first. Depending on the bonding modes of pyridinium toward the H atoms of the NH2 group, two mechanisms, denoted as E- and Z-pyridinium, have been studied. For the sake of simplicity, only the rate-determining mechanism in the i-C3F7C(O)NH2 + TFAA reaction, as outlined in Scheme 4, has been considered to evaluate the pyridinium-inspired routes.

2.4.1. The E-Pyridinium Mechanism

The potential energy profiles and the structures of the key intermediates and transition states for the reaction of E-i-C3F7C(O)NH2 + TFAA + 1Py are illustrated in Scheme 5. Geometric structures for the key species are shown in Figure 4.
When the Py molecule approaches the H atom of NH2, a weakly bound complex IM19 is formed via an N...H–N hydrogen bond of 1.833 Å, which is similar to the IM13 in the i-C3F7C(O)NH2 + TFAA reaction in geometry but the free energy increases to 7.7 kcal/mol. The hydrolysis undergoes the simultaneous C–O addition and proton migration via TS26. Evidently, the breaking N–H bond in TS26 is longer than that in the analogous TS15 (Figure 3), resulting in the shorter forming H–O bond. Meanwhile, the forming C–O bond is 1.462 Å, which is 0.035 Å longer than that in TS15. Therefore, TS26 turns to be a later barrier due to the participation of Py. In fact, the formation of IM20 becomes slightly more endothermic, as shown by the extended H–bond (Figure 4), and the barrier height for TS26 increases by roughly 1 kcal/mol as well. The subsequent decomposition of IM20 occurs via TS27. Since Py only acts as a spectator, the reacting geometries of TS27 are similar to those in TS18 but the barrier height increases to 14.8 kcal/mol for the former with respect to 12.3 kcal/mol for the latter. In addition, the release of TFA to form IM21 becomes even more endothermic by about 5 kcal/mol.
Interestingly, it was found that the further decomposition of IM21 takes place via a new mechanism, namely, the 3-body cleavage barrier TS28. As can be seen in Figure 4, the C–O bond is broken directly, accompanied by the proton transfer from NH to the basic N atom of Py, leading to C4, pyridinium, and CF3CO2 anion. In comparison with the intermolecular H migration in TS25, such a 3-body fission path involves a barrier height of only 7.1 kcal/mol, which is considerably lower than that for TS25. As a result, one Py molecule is capable of reverse tuning the rate-determining step for the hydrolysis of i-C3F7C(O)NH2 with TFAA. Although the relative free energy of TS27 increases slightly to 39.2 kcal/mol, the barrier height for the new rate-determining step, 14.8 kcal/mol, has been significantly reduced due to pyridinium. Note that the 3-body decomposition can occur via a conformer of IM21, namely, IM22, which is formed by the Py-assistant N–H rotation via TS30. The breaking C–O bond (1.638 Å) in TS30 is shorter than that in TS28 (1.761 Å). The reactant-like barrier induces a further 1 kcal/mol reduction in energy.
In fact, a detour to form pyridinium exists, as shown in Scheme 5. The addition of the carbonyl C of TFAA to the carbonyl O of i-C3F7C(O)NH2 goes via TS31 along with the simultaneous N–H bond breaking to form the pyridinium. The newly formed C–O bond in the adduct IM23 is 1.504 Å and the pyridinium is coordinated to the N site via an H–bond of 1.703 Å. The subsequent decomposition of IM23 occurs via TS32, which is a three-body fission process with a barrier of only 1.7 kcal/mol, leading to the pyridinium-coordinated product IM24 and CF3COO anion. The energy of TS31 is 6.7 kcal/mol lower than that of TS26, indicating that such a detour to pyridinium should be a dominant path. However, the further decomposition of IM24 to form C4 surmounts the significant barrier, i.e., TS33, which corresponds to a six-centered intramolecular TFA elimination process. The barrier height is as high as 28.9 kcal/mol, excluding any importance for such a detour mechanism in the production of C4.
It should be noted that the formation of C4 can be exothermic after the recombination of pyridinium with CF3COO anion. The overall reaction-free energy is calculated to be −2.8 kcal/mol. Therefore, the pyridinium-inspired TFAA-catalyzed hydrolysis of i-C3F7C(O)NH2 becomes not only kinetically feasible but also thermodynamically spontaneous under ambient conditions.

2.4.2. The Z-Pyridinium Mechanism

While the Py molecule approaches to the H atom of NH2 from the Z conformational direction, the hydrolysis of i-C3F7C(O)NH2 takes place via a completely different mechanism from that for the E-conformations, as illustrated in Scheme 6 and Figure 5.
Starting from the pre-reactive complex, the addition of TFAA molecule to the carbonyl O atom goes via TS34 with the forming C–O bond as long as 2.014 Å. Simultaneously, the H atom is transferred to Py to form pyridinium, which involves two hydrogen bonds towards the N site and the O site in the adduct IM26. In comparison with TS34, the H atom does not migrate from N to O. The other addition path occurs via TS36, which appears to be more product-like than the analogous E-conformer in view of the forming C–O bond the breaking N–H bond. As a result, the energy of TS36 increases to 22.7 kcal/mol. In contrast, the energy of TS34 drops to 19.1 kcal/mol, which is 6.5 kcal/mol lower than that of the corresponding E-conformer. Moreover, the adduct has been stabilized as well due to the presence of two H-bonds of pyridinium. Both adducts can decompose readily. The pyridinium is shifted to the terminal CF3COO group via TS35 or TS37, accompanied by the C–O bond cleavage. The relevant barriers are only marginal to produce IM15 and the ionized CF3COO-pyridinium complex. Evidently, the Z-directed pyridinium is able to catalyze the first H migration as an efficient proton carrier.
The decomposition of IM15 remains the same as that in the i-C3F7C(O)NH2 + TFAA reaction. The significant barrier TS25 leading to C4 makes it the rate-determining step. Therefore, the Z-directed pyridinium mechanism would not change the overall production of C4 even though the first-step H migration has been accelerated.

2.5. The i-C3F7C(O)NH2 + TFAA + 2Py→C4 + 2TFA-Py

In consideration of the important role of Py to the i-C3F7C(O)NH2 + TFAA reaction, two Py molecules were introduced towards the naked NH2 group. The hydrolysis reaction routes are depicted in Scheme 7 and the corresponding transition states and adducts are illustrated in Figure 6.
Note that the energetic routes in Scheme 7 might be viewed to be either E-Py-mediated Z- i-C3F7C(O)NH2 + TFAA + Py reaction or Z-Py-mediated E- i-C3F7C(O)NH2 + TFAA + Py reaction. It is evident that the additional Py in the E-direction leads to the slight increase in the barrier heights for the first-step H migration. The lowest barrier is denoted as TS38. The addition between the carbonyl C site of TFAA to the carbonyl O site of i-C3F7C(O)NH2 goes along with the proton transfer from N to Py. The reacting geometrical parameters between TS38 and TS34 do show significant differences. The forming C–O bond is only 1.623 Å in the former, which is about 0.4 Å shorter than that in the latter. More significantly, the pyridinium has not been formed yet in view of the breaking N–H bond (1.238 Å) and the forming H–N bond (1.334 Å) in TS38 in comparison with the equilibrium geometry of pyridinium. The other Py molecule is H–bonded to the NH2 group as a spectator. Due to the steric hindrance of two Py rings, the energy of TS38 indeed increases by 2.7 kcal/mol.
The alternative path for the C–O addition occurs via TS40, in which only the E-directed Py molecule is reacting and the Z-directed Py acts as a spectator. It appears that TS40 is an early barrier in terms of the forming C–O bond (1.952 Å), in contrast to 1.875 Å in TS31 (Figure 4). However, the barrier height for TS40 is about 5 kcal/mol higher than for TS31 because of the π-stacking complexation of the two Py molecules. The subsequent decomposition of the adducts takes place rapidly in view of the negligible barrier heights for TS41 to kick the CF3COO anion out and TS39 to generate IM21. The latter channel is more exothermic than the former by nearly 8 kcal/mol, supporting that the formation of IM21 should be the major intermediate of hydrolysis.
It is noted that IM21 decomposes to form C4 according to the aforementioned mechanism in Scheme 5. The lowest free energy for the 3-body fission to C4, pyridinium, and anionic CF3COO is 7.1 kcal/mol. Therefore, the rate-determining step for the hydrolysis of i-C3F7C(O)NH2 + TFAA + 2Py should be the first-step C–O addition. In contrast, the pyridinium intermediate IM31 invokes a new barrier, i.e., TS42. The H atom is transferred from N to the terminal O in the form of pyridinium and the C–O bond is ruptured as the elongated distance of 1.810 Å, leading to C4, pyridinium, and the ionic CF3COO-pyridinium complex. The two Py rings appear to be well separated in TS42 and thus the barrier height is only 6.9 kcal/mol, which is nearly identical to that for the IM31 decomposition. It is conceivable that the two reaction pathways might be competitive in terms of the similar energetics.
For clarity, the hydrolysis reaction of i-C3F7C(O)NH2 to C4 as catalyzed by TFAA and 2Py is illustrated in Scheme 8. While the two Py molecules can be regenerated, TFAA has been hydrolyzed to two TFA molecules. It is worth noting that Py can form ionic pyridinium complex with TFA due to the strong H–bond interaction. As a result, the formation of C4+2TFA-Py can be exothermic by 8.2 kcal/mol with respect to the initial reactants, supporting that the hydrolysis of i-C3F7C(O)NH2 can proceed even at low temperatures. The calculated lowest Gibbs free energy barriers for the stepwise hydrolysis of i-C3F7C(O)NH2 to form C4 are summarized in Table 1. Depending on the ratios of the TFAA and Py molecules, the catalytic effect to hydrolysis appears to be in reasonable agreement with the experimental yields of C4. Therefore, the synergic effect of the TFAA/Py co-catalyst in the production of C4 by hydrolysis of i-C3F7C(O)NH2 has been revealed in the present theoretical calculations. More quantitative comparisons between theory and experiment can be made once the Kcat data are available for the title reactions.

3. Computational Methods

The geometries of reactants, intermediates (IM), transition states (TS), and products were fully optimized using the density functional M06-2X with a triple-ζ 6-311++G(d, p) basis set [17,18]. It has been demonstrated that the M06-2X functional performs well for main group thermochemistry and thermochemical kinetics [19,20]. Alternatively, the B3LYP-D3 and ωB97X-D were also employed for the purpose of comparison [21,22,23,24]. The polarizable continuum model using the integral equation formalism variant was employed in the self-consistent reaction field method [25,26,27]. Dimethylformamide was used as the solvent with the universal force-field atomic radii in accordance with the experimental conditions. Vibrational frequency was calculated at the same level to characterize the nature of the stationary points. The minimum has all real frequencies and TS has only one imaginary vibrational mode corresponding to the reaction coordinates (see Figure S1). In addition, the intrinsic reaction coordinates were calculated to confirm that the transition state links the designated reactants and products.
On the basis of the M06-2X/6-311++G(d,p) optimized geometrical parameters, the complete basis set model chemistry with quadratic mode, i.e., CBS-QB3, was employed to validate the DFT-energetics. The accuracy of CBS-QB3 approaches to 0.87 kcal/mol in terms of mean absolute error for the G2 test set, which is applicable to the calculations of transition states for chemical reactions [28]. Moreover, the T1 diagnostic values of the couple-cluster wave functions were examined to estimate the multi-configurational effect (Tables S1 and S2). It is evident that the T1 diagnostic values for all the species of concern are less than 0.02, showing a negligible multi-configurational effect [29]. For the sake of computational efficiency, the CBS-QB3 energetics were obtained using the double-layer (DL) scheme [30]. Because the i-C3F7 group in i-C3F7C(O)NH2 and the CF3 groups in TFAA are only spectators during the production of C4, they were treated as the low layers in the DL-CBS-QB3 calculations. The reactive C(O)NH2 and C(O)OC(O) moieties were considered to be the high layers. It has been demonstrated that the DL-CBS-QB3 energetics are in excellent agreement with the full-size CBS-QB3 data with the mean absolute deviation of 1.2 kcal/mol as calibrated by the C4 + OH reaction [30]. The computational efficiency could be improved by roughly two orders of magnitudes. The correlation between the DL-CBS-QB3 energetics and the M06-2X/6-311++G(d,p) data is demonstrated in Figure 7 on the basis of 29 key stationary points (Tables S1 and S2). Good agreement has been achieved in the range from −20 to 90 kcal/mol. The mean absolute deviation for all the energetics is only 1.0 kcal/mol. Moreover, the B3LYP-D3 and ωB97X-D optimized energetics are in good agreement with the M06-2X data (Tables S1 and S2). Therefore, the M06-2X/6-311++G(d,p) geometrical parameters and energies will be used in the discussion unless stated otherwise.
In all of the potential energy profiles, the solvation-corrected relative Gibbs free energies were presented. Since the ideal gas-phase model inevitably overestimates the entropic contributions to the free energies, the calculated free energies for two-to-one (or one-to-two) transformations at 298.15 K have been corrected by a scaling factor of 0.5 to the gas phase entropic contributions. All the quantum chemical calculations were performed using the Gaussian 16 suite of programs [31].

4. Conclusions

The mechanisms for the dehydration of i-C3F7C(O)NH2 to produce C4 have been theoretically clarified. It is confirmed that the formation of C4 involves a stepwise H migration mechanism via the imide intermediate. Without catalysis, the reaction is dominated by the second step with an activation barrier of 63.1 kcal/mol. Pyridine is capable of reducing the energy barrier for the rate-determining step by acting as the acceptor for hydrogen bonds to NH2 and C=O groups and the proton relay to form the pyridinium species, leading to a significant reduction in the barrier of about 28 kcal/mol. The production of C4 cannot be kinetically feasible unless TFAA is presented. The reaction of i-C3F7C(O)NH2 with TFAA goes via multiple stepwise addition/elimination mechanisms to form C4 and two TFA molecules and the overall reaction is almost thermo-neutral. Moreover, the dehydration of i-C3F7C(O)NH2 could be promoted considerably through the synergic effect of the TFAA/Py co-catalyst. The pyridine molecules can act differently in view of the E- and Z-directed reaction paths. In comparison with the uncatalyzed reaction, the pyridinium-inspired TFAA-catalyzed hydrolysis of i-C3F7C(O)NH2 appears to be dominated by the formation of imide via a Gibbs free energy barrier of only 21.8 kcal/mol and the formation of the C4 and ionic Py-TFA complex is exothermic by about 8 kcal/mol, as indicative of the kinetically feasible and thermodynamically spontaneous mechanisms under ambient conditions. The calculated activation barriers for the rate-determining steps support the experimental yields of C4 with the different ratios of pyridine to TFAA. The present theoretical work provides new insights into the rational design on the novel catalysts for large-scale synthesis of the perfluorinated nitriles.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29163952/s1, Figure S1: Imaginary frequency (in cm-1) and the corresponding atomic displacements for transition states of concern calculated at the M06-2X/6-311++G(d,p) level of theory; Table S1: Relative energies (in kcal/mol) for the species involved in the unimolecular dehydration reaction of i-C3F7C(O)NH2 calculated at the M06-2X/6-311++G(d, p), B3LYP-D3/6-311++G(d, p), ωB97X-D/6-311++G(d, p) and DL-CBS-QB3 levels of theory on the basis of the M06-2X/6-311++G(d,p) optimized geometries. The T1 diagnostic values are calculated at the CCSD(T)/6-31+G(d) level; Table S2: Relative energies (in kcal/mol) for the species involved in the unimolecular dehydration reaction of i-C3F7C(O)NH2 with TFAA calculated at the M06-2X/6-311++G(d, p), B3LYP-D3/6-311++G(d, p), ωB97X-D/6-311++G(d, p) and DL-CBS-QB3 levels of theory on the basis of the M06-2X/6-311++G(d,p) optimized geometries. The T1 diagnostic values are calculated at the CCSD(T)/6-31+G(d) level; Table S3: M06-2X/6-311++G(d,p) optimized Cartesian coordinates, zero point energies (ZPE), and various electronic total energies for the species involved in unimolecular dehydration reaction of i-C3F7C(O)NH2. Table S4: M06-2X/6-311++G(d, p) optimized Cartesian coordinates, zero point energies (ZPE), and electronic total energies for the species involved in the reaction of i-C3F7C(O)NH2+2Py; Table S5: M06-2X/6-311++G(d, p) optimized Cartesian coordinates, zero point energies (ZPE), and various electronic total energies for the species involved in the reaction of i-C3F7C(O)NH2+TFAA; Table S6: M06-2X/6-311++G(d, p) optimized Cartesian coordinates, zero point energies (ZPE), and electronic total energies for the species involved in the reaction of E-i-C3F7C(O)NH2+TFAA+Py; Table S7: M06-2X/6-311++G(d, p) optimized Cartesian coordinates, zero point energies (ZPE), and electronic total energies for the species involved in the reaction of Z-i-C3F7C(O)NH2+TFAA+Py; Table S8: M06-2X/6-311++G(d, p) optimized Cartesian coordinates, zero point energies (ZPE), and electronic total energies for the species involved in the reaction of i-C3F7C(O)NH2+TFAA+2Py.

Author Contributions

Conceptualization, B.W.; methodology, J.X.; software, J.X.; validation, J.X.; investigation, B.W.; writing—original draft preparation, J.X.; writing—review and editing, B.W.; visualization, J.X.; supervision, B.W. and H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2021YFB2401400).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in this study are included in the article material; further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wooton, R.E.; Kegelman, M.R. Gases Superior to SF6 for Insulation and Interruption; Final Report; Electric Power Research Institute (EPRI)EL⁃2620: Washington, DC, USA, 1982; pp. 4–37. [Google Scholar]
  2. Boggs, S.A. Sulphur Hexafluoride-A Complex Dielectric. IEEE Electr. Insul. Mag. 1989, 5, 16–21. [Google Scholar] [CrossRef]
  3. Beroual, A.; Haddad, A. Recent Advances in the Quest for a New Insulation Gas with a Low Impact on the Environment to Replace Sulfur Hexafluoride (SF6) Gas in High-Voltage Power Network Applications. Energies 2017, 10, 1216. [Google Scholar] [CrossRef]
  4. Rokunohe, T.; Yagihashi, Y.; Endo, F.; Oomori, T. Fundamental Insulation Characteristics of Air; N2, CO2, N2/O2, and SF6/N2 Mixed Gases. Electr. Eng. Jpn. 2006, 155, 9–17. [Google Scholar] [CrossRef]
  5. Katagiri, H.; Kasuya, H.; Mizoguchi, H.; Yanabu, S. Investigation of the Performance of CF3I Gas as a Possible Substitute for SF6. IEEE Trans. Dielectr. Electr. Insul. 2008, 15, 1424–1429. [Google Scholar] [CrossRef]
  6. Kieffel, Y.; Biquez, F.; Ponchon, P.; Irwin, T. SF6 Alternative Development for High Voltage Switchgears. In Proceedings of the 2015 IEEE Power & Energy Society General Meeting, Denver, CO, USA, 26–30 July 2015; pp. 1–6. [Google Scholar]
  7. Kieffel, Y.; Irwin, T.; Ponchon, P.; Owens, J. Green Gas to Replace SF6 in Electrical Grid. IEEE Power Energy Mag. 2016, 14, 32–39. [Google Scholar]
  8. Rabie, M.; Franck, C.M. Assessment of Eco-Friendly Gases for Electrical Insulation to Replace the Most Potent Industrial Greenhouse Gas SF6. Environ. Sci. Technol. 2018, 52, 369–380. [Google Scholar] [CrossRef]
  9. Bulánek, R.; Novoveská, K.; Wichterlová, B. Oxidative Dehydrogenation and Ammoxidation of Ethane and Propane over Pentasil Ring Co-Zeolites. Appl. Catal. A-Gen. 2002, 235, 181–191. [Google Scholar] [CrossRef]
  10. Sundermeier, M.; Zapf, A.; Beller, M. Palladium-Catalyzed Cyanation of Aryl Halides: Recent Developments and Perspectives. Eur. J. Inorg. Chem. 2003, 2003, 3513–3526. [Google Scholar] [CrossRef]
  11. Hodgson, H.H. The Sandmeyer reaction. Chem. Rev. 1947, 40, 251–277. [Google Scholar]
  12. Wang, L.; Shen, C.; Wang, H.; Zhou, W.; Sun, F.; He, M.-Y.; Chen, Q. Selective Conversion of Aldehydes into Nitriles and Primary Amides in aQueous Media. J. Chem. Res. 2012, 36, 460–462. [Google Scholar] [CrossRef]
  13. Campagna, F.; Carotti, A.; Casini, G. A Convenient Synthesis of Nitriles from Primary Amides under Mild Conditions. Tetrahedron. Lett. 1977, 18, 1813–1815. [Google Scholar] [CrossRef]
  14. Gao, Z.; Wang, M.; Wang, S.; Wang, Y.; Peng, R.; Yu, P.; Luo, Y. Novel and Efficient Synthesis of Insulating Gas- Heptafluoroisobutyronitrile from Hexafluoropropylene. R. Soc. Open Sci. 2019, 6, 181751. [Google Scholar] [CrossRef] [PubMed]
  15. Sureshbabu, V.V.; Naik, S.A.; Nagendra, G. Synthesis of Boc-Amino Tetrazoles Derived from α-Amino Acids. Synth. Commun. 2009, 39, 395–406. [Google Scholar] [CrossRef]
  16. Liu, R.Y.; Bae, M.; Buchwald, S.L. Mechanistic Insight Facilitates Discovery of a Mild and Efficient Copper-Catalyzed Dehydration of Primary Amides to Nitriles Using Hydrosilanes. J. Am. Chem. Soc. 2018, 140, 1627–1631. [Google Scholar] [CrossRef]
  17. Frisch, M.J.; Pople, J.A.; Binkley, J.S. Self-consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets. J. Chem. Phys. 1984, 80, 3265–3269. [Google Scholar]
  18. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  19. Mardirossian, N.; Head-Gordon, M. Thirty Years of Density Functional Theory in Computational Chemistry: An Overview and Extensive Assessment of 200 Density Functionals. Mol. Phys. 2017, 115, 2315–2372. [Google Scholar] [CrossRef]
  20. Zhao, Y.; Truhlar, D.G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157–167. [Google Scholar] [CrossRef]
  21. Becke, A.D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  22. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed]
  23. Chai, J.-D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom–Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef] [PubMed]
  24. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  25. Cossi, M.; Barone, V.; Cammi, R.; Tomasi, J. Ab Initio Study of Solvated Molecules: A New Implementation of the Polarizable Continuum Model. Chem. Phys. Lett. 1996, 255, 327–335. [Google Scholar] [CrossRef]
  26. Barone, V.; Cossi, M.; Tomasi, J. A New Definition of Cavities for the Computation of Solvation Free Energies by the Polarizable Continuum Model. J. Chem. Phys. 1997, 107, 3210–3221. [Google Scholar] [CrossRef]
  27. Barone, V.; Cossi, M.; Tomasi, J. Geometry Optimization of Molecular Structures in Solution by the Polarizable Continuum Model. J. Comput. Chem. 1998, 19, 404–417. [Google Scholar] [CrossRef]
  28. Montgomery, J.A.; Frisch, M.J.; Ochterski, J.W.; Petersson, G.A. A Complete Basis Set Model Chemistry. VI. Use of Density Functional Geometries and Frequencies. J. Chem. Phys. 1999, 110, 2822–2827. [Google Scholar] [CrossRef]
  29. Lee, T.J.; Taylor, P.R. A Diagnostic for Determining the Quality of Single-Reference Electron Correlation Methods. Int. J. Quantum Chem. 2009, 36, 199–207. [Google Scholar] [CrossRef]
  30. Yu, X.; Hou, H.; Wang, B. Double-Layered Composite Methods Extrapolating to Complete Basis-Set Limit for the Systems Involving More than Ten Heavy Atoms: Application to the Reaction of Heptafluoroisobutyronitrile with Hydroxyl Radical. J. Phys. Chem. A 2017, 121, 9020–9032. [Google Scholar] [CrossRef]
  31. Frisch, M.J.; Truchs, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
Scheme 1. Mechanistic routes for the dehydration of i-C3F7C(O)NH2 to produce i-C3F7CN as catalyzed by trifluoroacetic anhydride and pyridine.
Scheme 1. Mechanistic routes for the dehydration of i-C3F7C(O)NH2 to produce i-C3F7CN as catalyzed by trifluoroacetic anhydride and pyridine.
Molecules 29 03952 sch001
Scheme 2. Gibbs free energy profiles (in kcal/mol) for the unimolecular dehydration reaction of i-C3F7C(O)NH2 to form i-C3F7CN calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 2. Gibbs free energy profiles (in kcal/mol) for the unimolecular dehydration reaction of i-C3F7C(O)NH2 to form i-C3F7CN calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch002
Figure 1. Geometrical parameters for various species involved in the dehydration reaction of i-C3F7C(O)NH2 optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 1. Geometrical parameters for various species involved in the dehydration reaction of i-C3F7C(O)NH2 optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g001
Scheme 3. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2+2Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 3. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2+2Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch003
Figure 2. Geometrical parameters for the species involved in the energetically favorable paths for the i-C3F7C(O)NH2+2Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 2. Geometrical parameters for the species involved in the energetically favorable paths for the i-C3F7C(O)NH2+2Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g002
Scheme 4. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2 with TFAA in DMF solution calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 4. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2 with TFAA in DMF solution calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch004
Figure 3. Geometrical parameters for various species involved in the reaction of i-C3F7C(O)NH2 with TFAA optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 3. Geometrical parameters for various species involved in the reaction of i-C3F7C(O)NH2 with TFAA optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g003
Scheme 5. Gibbs free energy profiles (in kcal/mol) for the reaction of E-i-C3F7C(O)NH2+TFAA+1Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 5. Gibbs free energy profiles (in kcal/mol) for the reaction of E-i-C3F7C(O)NH2+TFAA+1Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch005
Figure 4. Optimized geometrical parameters for various species involved in the reaction of E-i-C3F7C(O)NH2+TFAA+1Py at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 4. Optimized geometrical parameters for various species involved in the reaction of E-i-C3F7C(O)NH2+TFAA+1Py at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g004
Scheme 6. Gibbs free energy profiles (in kcal/mol) for the reaction of Z-i-C3F7C(O)NH2+TFAA+1Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 6. Gibbs free energy profiles (in kcal/mol) for the reaction of Z-i-C3F7C(O)NH2+TFAA+1Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch006
Figure 5. Geometrical parameters for the species involved in the energetically favorable paths for the Z-i-C3F7C(O)NH2+TFAA+1Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 5. Geometrical parameters for the species involved in the energetically favorable paths for the Z-i-C3F7C(O)NH2+TFAA+1Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g005
Figure 6. Geometrical parameters for the species involved in the energetically favorable paths for the i-C3F7C(O)NH2+TFAA+2Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Figure 6. Geometrical parameters for the species involved in the energetically favorable paths for the i-C3F7C(O)NH2+TFAA+2Py reaction optimized at the M06-2X/6-311++G(d,p) level of theory. Bond distances are in Å and angles are in degrees.
Molecules 29 03952 g006
Scheme 7. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2+TFAA+2Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Scheme 7. Gibbs free energy profiles (in kcal/mol) for the reaction of i-C3F7C(O)NH2+TFAA+2Py calculated at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch007
Scheme 8. Pyridine regeneration mechanism in the i-C3F7C(O)NH2+TFAA+2Py reaction at the M06-2X/6-311++G(d,p) level of theory.
Scheme 8. Pyridine regeneration mechanism in the i-C3F7C(O)NH2+TFAA+2Py reaction at the M06-2X/6-311++G(d,p) level of theory.
Molecules 29 03952 sch008
Figure 7. Correlations between the DL-CBS-QB3 and the M06-2X/6-311++G(d,p) energies (in kcal/mol). Solid triangles: for the reaction of i-C3F7C(O)NH2+TFAA. Open squares: for the dehydration reaction of i-C3F7C(O)NH2.
Figure 7. Correlations between the DL-CBS-QB3 and the M06-2X/6-311++G(d,p) energies (in kcal/mol). Solid triangles: for the reaction of i-C3F7C(O)NH2+TFAA. Open squares: for the dehydration reaction of i-C3F7C(O)NH2.
Molecules 29 03952 g007
Table 1. Summary of reactant ratio, reaction yield, and the H migration step, for a series of reactions about i-C3F7C(O)NH2 by the M06-2X/6-311++G(d,p) method.
Table 1. Summary of reactant ratio, reaction yield, and the H migration step, for a series of reactions about i-C3F7C(O)NH2 by the M06-2X/6-311++G(d,p) method.
i-C3F7C(O)NH2: TFAA: PyFirst StepSecond StepExperimental Yield [14]
1:0:050.963.10%
1:0:225.235.20%
1:1:024.725.615%
1:1:125.6(E)/19.1(Z)7.1(E)/25.6(Z)75%
1:1:221.87.186%
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

Xiong, J.; Hou, H.; Wang, B. Theoretical Investigation of the Pyridinium-Inspired Catalytic Dehydration of Heptafluoro-Iso-Butyramide for the Synthesis of Environmentally Friendly Insulating Gas Heptafluoro-Iso-Butyronitrile. Molecules 2024, 29, 3952. https://doi.org/10.3390/molecules29163952

AMA Style

Xiong J, Hou H, Wang B. Theoretical Investigation of the Pyridinium-Inspired Catalytic Dehydration of Heptafluoro-Iso-Butyramide for the Synthesis of Environmentally Friendly Insulating Gas Heptafluoro-Iso-Butyronitrile. Molecules. 2024; 29(16):3952. https://doi.org/10.3390/molecules29163952

Chicago/Turabian Style

Xiong, Jiageng, Hua Hou, and Baoshan Wang. 2024. "Theoretical Investigation of the Pyridinium-Inspired Catalytic Dehydration of Heptafluoro-Iso-Butyramide for the Synthesis of Environmentally Friendly Insulating Gas Heptafluoro-Iso-Butyronitrile" Molecules 29, no. 16: 3952. https://doi.org/10.3390/molecules29163952

APA Style

Xiong, J., Hou, H., & Wang, B. (2024). Theoretical Investigation of the Pyridinium-Inspired Catalytic Dehydration of Heptafluoro-Iso-Butyramide for the Synthesis of Environmentally Friendly Insulating Gas Heptafluoro-Iso-Butyronitrile. Molecules, 29(16), 3952. https://doi.org/10.3390/molecules29163952

Article Metrics

Back to TopTop