PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
 
Comput Theor Chem. 2017 August 15; 1114: 140–145.
PMCID: PMC5521852

Theoretical study of Ni+ assisted C—C and C—H bond activations of propionaldehyde in the gas phase

Graphical abstract

An external file that holds a picture, illustration, etc.
Object name is fx1.jpg
Keywords: Density functional theory (DFT), Decomposition reaction, Charge decomposition analysis (CDA), Bonding analysis

Abstract

The reactions of Ni+ with propionaldehyde in the gas phase have been systematically investigated using density functional theory at the B3LYP/def2-TZVP level. The decomposition reaction mechanism has been identified. Our calculations indicated that Ni+ can assist decomposition of propionaldehyde to form Ni+CO and C2H6 through two types of reaction channel: C—C bond activation and C—H bond activation. In addition, charge decomposition analysis (CDA) was carried out to obtain a deeper understanding for orbital interaction of the initial complex. The bonding properties of the species involved were discussed by means of diverse analysis methods including electron localization function (ELF) and atoms in molecules (AIM).

1. Introduction

The transition metals play an important role in the catalysis, owing to their unique electronic configuration [1]. The last few decades, the reaction of transition metals with simple organic molecules in the gas phase is one of the important topics [2], [3], [4], [5], [6], [7], [8], [9], [10]. Not only because of the tremendous practical importance to the petroleum industry but also due to the fundamental importance of σ-bonds. A significant number of theoretical and experimental research have been extensively investigated at understanding the mechanisms C—C and C—H bond activations of small hydrocarbons by transition metal atoms, ions and oxide ions [11], [12], [13], [14], [15]. The experiment research has provides abundant information, such as intrinsic binding properties and thermodynamical data. The theoretical results not only rationalize the experimental findings but also provide molecular level details of catalytic processes.

Aldehydes and ketones have become promising candidates for characterizing C—C and C—H bond activations [16]. Over the last several years, the reactions of acetaldehyde and acetone with all kinds of transition metal ions have attracted significant attention, focusing mainly on the investigation of chemical bonds activation by experimental techniques and theoretical calculations. These researches indicated that the reaction product always consists of a stable, neutral organic molecule (methane, ethane, etc.) and the corresponding ionic fragment. At low energies, Cr+, Cu+ and Mn+ are unreactive with acetone, but Fe+, Co+, Ni+ react with acetone to produce C2H6 and M+CO [17], [18], [19], [20]. On the contrary, transition metal ions with high oxygen affinities react with acetone to form MO+ ions and C3H6, e.g. Sc+, V+, Ti+, Gd+ and Pr+ [21], [22]. The reactions of transition metal ions with acetaldehyde found that decarbonylation of acetaldehyde by Fe+, Co+, Cr+, Ni+ is a dominant process at low reaction energies. The reaction products are CH4 and M+CO [23], [24], [25], [26], [27], [28]. The 4Fe+, Co+, Cr+ mediated systems decarbonylation of CH3CHO only through C—C activation, but for 6Fe+, Ni+ decarbonylation reaction is take place via C—C activation and C—H activation [25].

Recently, D. J. Bellert and co-workers carried out experimental research using the single photon initiated decomposition rearrangement reaction (SPIDRR) technique and a simple theoretical research for Ni+ assisted decomposition of propionaldehyde [29]. On the basis of the experimental findings, the researchers proposed a possible reaction mechanism for the propionaldehyde decomposition. They guessed the reaction progressing along two competitive pathways. Ni+ insert into either the CH3CH2C(O)-H or CH3CH2-CHO σ-bond in the propionaldehyde, and the two different products were yielded.

Ni++C3H6ONi+(C3H6O)Ni+CO+C2H6
(1)

However, the details of the reaction pathways and the geometrical parameters of the intermediates, transition states involved have not been fully discussed. Hence theoretical study in the reaction is still needed. In this paper, we reported a theoretical study on the reaction of C—C and C—H bond activations of propionaldehyde mediated by Ni+ using density functional theory (DFT). The molecular structures and energetics properties of all the stationary points are optimized and calculated. Charge decomposition analysis (CDA) has been carried out to understand how fragment orbitals are mixed to form complex orbital. In addition, we also studied the bonding evolution during the reaction pathways by means of electron localization function (ELF) and atoms in molecules (AIM).

2. Computational details

The geometries and energies of reactants, intermediates, transition states and products involved in the reaction were fully optimized by using B3LYP hybrid density functional method in conjunction with def2-TZVP basis sets for all atoms [30], [31], [32]. This selected approach has been valuated to be of good validity for the Ni+-containing system in previous papers [33], [34], [35]. Vibrational frequencies were further calculated at the same level to determine that each optimized species was a minimum or a saddle point on the potential energy surface. Local minima on the potential energy surface (PES) have no negative eigenvalue, and saddle points have only one negative eigenvalue. Intrinsic reaction coordinate (IRC) calculations [36], [37] were performed at the B3LYP/def2-TZVP level of theory to confirm the correct connection of the transition state with the corresponding reactants and products. All the calculations were carried out using the Gaussian 09 software package [38].

In order to gain a deeper understanding of the chemical bonds evolution, the wavefunction files (.wfn) produced by Gaussian 09 were used as the inputs into Multiwfn to perform the electron localization function (ELF) and atoms in molecules (AIM) analysis. Charge decomposition analysis (CDA) was created, with the aim of gathering insights into the orbital interaction of the initial complex, as implemented in the Multiwfn code [39].

3. Results and analysis

The optimized geometries and structural parameters corresponding to all the intermediates and transition states involved along the two different reaction pathways are depicted in Fig. 1. The orbital interaction diagram of the initial complex can be directly plotted in Fig. 2. Besides, bonding analysis of all species involved in the reaction was performed. The analysis results are depicted in Fig. 3 and Table S1. The double potential energy surface and relevant energies of all species are shown in Fig. 4 and Table S2. It can be found that there are two different reaction pathways on the PES toward the reaction of Ni+ with CH3CH2CHO. A path is the C—H bond activation. The other path is activation of C—C bond. Ni+ insertion into the C—H, C—C bonds and the following rearrangement are key steps for this system. In the following sections, a detailed analysis of the reaction is described.

Fig. 1
Geometrical parameters of the intermediates and transition states involved on the double PES at the B3LYP/def2-TZVP level (bond lengths in angstrom and angles in degrees).
Fig. 2
Orbital interaction diagrams of the initial complex (molecular orbital energy in eV).
Fig. 3
ELF projection map of stationary points on the Ni+ + CH3CH2CHO reaction pathway at B3LYP/def2-TZVP level of theory.
Fig. 4
Potential energy profiles for the reactions of Ni+ + CH3CH2CHO calculated at the B3LYP/def2-TZVP level of theory.

3.1. Initial complexes

When Ni+ collides with CH3CH2CHO, the initial complex Ni+(C3H6O) with the Cs symmetry is formed. In order to gain a deeper understanding of how fragment orbitals are mixed to form complex orbital, charge decomposition analysis (CDA) was performed and the orbital interaction diagrams are plotted in Fig. 2. As shown in Fig. 2, the horizontal lines on the left, the right and the middle are corresponding to the CH3CH2CHO fragment orbitals, the Ni+ fragment orbitals and the I1a/I1b complex orbitals. The calculated results indicate that electrons have been transferred from CH3CH2CHO to Ni+ during the formation of the Ni+(C3H6O) complex. The net numbers of alpha and beta electrons of I1a are 3.9350 and 0.1254, respectively. So the net number of all electrons of I1a is 4.0604. The net numbers of alpha and beta electrons of I1b are 0.8895 and 0.2112, respectively. Thus, the net number of all electrons of I1b is 1.1001. From Fig. 2, it can be seen that the highest occupied molecular orbital (HOMO) of I1a has two orbitals of the same energy. These are both orbitals that are composed by 3dyz orbital of Ni, have a′ and a″ spatial symmetry, respectively. This orbital specifications stem from the Cs symmetry of the complex. the lowest unoccupied molecular orbital (LUMO) with a′ spatial symmetry is primarily composed by the 4s orbital of Ni. The highest occupied molecular orbital (HOMO) of I1b, which is primarily composed by Ni 3dz2 orbital, has a′ spatial symmetry. And the lowest unoccupied molecular orbital (LUMO) with a′ spatial symmetry is mainly composed of the empty 4s orbital of Ni.

3.2. Bonding analysis

In this section, the exploration of the bonding evolution of all species involved is based on the analysis of the electron localization function (ELF) and atoms in molecules (AIM). The ELF analysis shows a detailed picture of the bonding, which exhibits maxima at the most probable positions of local electron pairs and all special positions are surrounded by the basin in which there is an increased probability of finding an electron pair [40]. In the AIM analysis, the existence of the interaction positively correlated with the presence of the bond critical point (bcp) and the strength of the bond is indicated from the magnitude of the electron density at the bond critical point (bcp) [41].

The ELF projection figures of the stationary points of Ni+ + CH3CH2CHO are shown in Fig. 3. The AIM parameters [42], [43] of bond critical points (bcp) are listed in Table S1. As is shown in Table S1, Bond critical points (bcp) have been analyzed in terms of the values of the electron density ρ(r) and its Laplacian of electron density [nabla]2ρ(r) at the (3,−1) critical points, the Lagrangian kinetic energy G(r), the electron potential energy density V(r) and the total electron energy density E(r) or H(r), E(r) = G(r) + V(r) = H(r). The larger the electron density ρ at the bond critical points is, the stronger the interaction will be. The sign of Laplacian [nabla]2ρ(r) is a generally recognized criterion for discriminating bond types [44]. The positive [nabla]2ρ(r) implies the loose charge density at the critical point. Whereas the negative [nabla]2ρ(r) means the shared interaction [45]. The E(r) standard has been proven to be very applicable to characterize the nature of a chemical bond for heavy-atom systems [46], [47], [48]. Negative and positive signs suggest that the interactions are covalent character and closed-shell interactions, respectively [42], [49].

  • (1)
    I1a and I1b species. The ELF analysis of the complex reveals that the absence of a disynaptic valence basin between the Ni atom and O atom of CH3CH2CHO, uses it as proof that no covalent bond is formed, which confirms that the interaction can be considered as an electrostatic interaction. According to the AIM analysis, there is a (3,−1) critical point between the Ni and O atoms, but the corresponding charge density is really low (ρ(bcp) = 0.093 au for I1a. ρ(bcp) = 0.086 au for I1b). The [nabla]2ρ(bcp) value is small and positive ([nabla]2ρ(bcp) = 0.591 au for I1a. [nabla]2ρ(bcp) = 0.514 au for I1b). The energy density E(r) is negative but |E(r)| is quite small (E(r) = −0.018 au for I1a. E(r) = −0.017 au for I1b). All available evidence points to the fact that no covalent bond is formed between the Ni atom and O atom of CH3CH2CHO.
  • (2)
    I2a species. There is an increasing trend in the ρ(bcp) of the Ni—C1 and Ni—H1 bonds. The[nabla]2ρ(bcp) is decreased. E(r) is still negative (E(r) = −0.087 au for Ni—H1 bond. E(r) = −0.062 au for Ni—C1 bond), but |E(r)| is increased. The ELF analysis indicates that the formation of Ni—C1 and Ni—H1 bonds at this step, as evidenced by the strengthening of the V(Ni,C1) and (Ni,H1) basins, giving place to the disappeared of a disynaptic V(C1—H1) basin.
  • (3)
    I2b species. Our ELF analysis reveals that the C1—C2 bond is completely broken and Ni—C1, Ni—C2 bond formed at this process. This fact is evidenced by the lack of V(C1,C2) disynaptic valence basin and its replacement by a disynaptic V(Ni,C1) and V(Ni,C2) basins. This description is consistent with the results of AIM analysis, which shows the inexistence of the bond critical point between the C1 and C2 atoms. The ρ(bcp) and |E(r)| of the Ni—C1 and Ni—C2 bonds have a stable increase. The [nabla]2ρ(bcp) also has a marked decrease, although it is still positive.
  • (4)
    I3 species. There are two ways to form species I3. One is completely broken of the C1—C2 covalent bond of T2a. This is evidenced by the lack of a disynaptic V(C1,C2) basin in the ELF analysis. The fact is also confirmed by AIM analysis, which shows a bcp (3,−1) vanish between the C1 and C2 atoms. The other is H1-migration of T2b from C1 atom to Ni atom to form Ni—H1 bond, as evidenced that a bcp (3,−1) exists between the Ni and H1 atoms. The ρ(bcp) and [nabla]2ρ(bcp) are 0.135 au and −0.016 au, respectively.
  • (5)
    I4 species. At this stage, the precursor (I4) is generated with the transferring of H1 from the Ni atom to C2 atom. The I4 intermediate can directly dissociate into the Ni+CO and C2H6. The appearance of the V(C2,H1) basin as shown in the ELF analysis. According to the AIM analysis, the ρ(bcp) and E(r) of the C2—H1 bond is 0.254 au and −0.246 au, respectively, indicating the covalent nature of the C2—H1 bond.

3.3. Reaction mechanism

All geometrical parameters of the intermediates and transition states involved on the doublet state are shown in Fig. 1. The relevant energies of all species are listed in Table S2. As shown in Fig. 4. There are two different pathways. Six minima (noted as I1a, I1b, I2a, I2b, I3, I4, respectively) and five first-order saddle points (noted as T1a, T1b, T2a, T2b, and T3, respectively) have been located along the reaction coordinate. ‘intermediates’ are specified by an ‘I’, while ‘T’ is the abbreviation for ‘transition states’. Moreover, ‘a’ denotes the C—H bond activation path, while ‘b’ stands for the Ni+ insertion into the C—C bond.

An encounter complex Ni+(C3H6O) is initially formed by the linkage of Ni+ with the O atom of propionaldehyde. We found that the interaction of Ni+ and CH3CH2CHO leads to two different encounter complexes, suggesting that the metal ion interacts with the two different lone pair orbitals of the oxygen atom. Once the complex is formed. The C1—H1 and C1—C2 bonds can be activated by the metal ion. Structurally, the two encounter complexes have overall Cs symmetry, and the symmetry plane is defined as Ni+—O—C1—C2. The Ni+—O—C1 angle of I1a and I1b is 137.8 and 119.4, respectively. The equilibrium Ni+—O distance severally is 1.888 Å and 1.935 Å. Upon binding to Ni+, the most significant change in propionaldehyde is the lengthening of C1—O bond. Because of oxygen polarizing charge toward Ni+, the C1—O bond is weakened. Energetically, the two complexes lie at below the ground state reactants. The stability of two different structures is nearly the same (226.1 kJ/mol for I1a and 243.1 kJ/mol for I1b), indicating coexistence of them in the gas phase.

The species H3CCH2CO—Ni+—H (I2a) may be formed by the metal ion insertion into the C1—H1 bond of propionaldehyde. This complex I2a has C1 symmetry and lies at −96.02 kJ/mol below the ground state reactants. The bond lengths of Ni+—C1 and Ni+—H1 are calculated as 1.888 Å and 1.476 Å, respectively. And a C1—Ni+—H1 angle is 89.11. The intermediate I2a and the encounter complex I1a are connected by T1a, as confirmed by IRC calculations. The system needs to overcome the activation energy barrier of 132.6 kJ/mol. The imaginary frequency is 410.4i cm−1 corresponding primarily to the H1 atom migrates from the C1 atom to the Ni atom to form the Ni—H1 bond. The Ni—H1 distance is 1.495 Å. It can be found that the geometry of T1a takes a large similar to I2a. Thus it is indeed the late transition state on the PES.

The complex I2a can be rearranged to form species I3 after the system surpasses the second transition state T2a, which is labeled by imaginary frequency 221.3i cm−1. The energy of T2a amounts to −16.30 kJ/mol. And corresponds to the highest barrier along the path, i.e., it is a rate-determining transition state in path a. The Ni+—H1, Ni+—C1, Ni+—C2 distances of species I3 are 1.492 Å, 1.834 Å, 1.997 Å, respectively, indicating that the Ni+ center interacts with all of these atoms.

Along the dissociation coordinates, the complex I3 can be rearranged to form a direct precursor (I4) of the dissociation products via transition state T3, which has an imaginary frequency of 498.4i cm−1. The species T3 lies at −66.13 kJ/mol below the reactants’ asymptote. Nonreactive dissociation of the Ni+(CO)-C2H6 would give rise to the decomposition products: Ni+CO and C2H6. The overall decomposition process of Ni+ + CH3CH2CHO is calculated to be exothermic by 149.5 kJ/mol. Our results indicated that the energy of the whole reaction pathway on this PES is lower than the reactant’s asymptote.

Next, we turn to the alternative C—C activation path for the formation of species I3. This branch involves the initial complex I1b. Subsequent Ni+ inserts into the C1—C2 bond to form species H3CH2C-Ni+-CHO (I2b). In H3CH2C-Ni+-CHO (I2b) structure, the Ni+—C1, Ni+—C2 bond distances are 1.832 Å, 1.899 Å, and a C1—Ni+—C2 angle is 107.3. Energetically, species I2b is calculated to be −74.47 kJ/mol below the separated reactants, or 21.55 kJ/mol higher than the corresponding C—H insertion minimum I2a. The encounter complex I1b and species I2b are connected by transition state T1b. The energy of the transition state is −51.83 kJ/mol. Note that T1b is rate-determining transition states in path b.

Intermediate I2b can also be converted to the species I3 through transition state T2b. The hydrogen atom is shifted from C1 atom to the transition metal atom to form species I3a. Energetically, the transition state T2b is calculated to be −55.04 kJ/mol lower than the corresponding transition state along the C—H bond activation pathway (T2a).

4. Conclusions

In this work, a specific theoretical calculation has been performed to further understand the gas phase reaction mechanism details of Ni+ with propionaldehyde using density functional theory. All geometries along the potential energy surfaces were fully optimized at B3lYP/def2-TZVP level. The orbital interaction of the initial complex has been performed by CDA. Bonding evolution was analyzed on the basis of ELF and AIM calculations. Those conclusions can be summarized as follows:

  • (1)
    We found that the decomposition reaction of Ni+ + CH3CH2CHO has two types of reaction channels: C—C activation and C—H activation. This reaction contains four elementary steps of encounter complexation, C—C and C—H bond activations, H-shift and nonreactive dissociation. The second step in the C—H bond activation process has to surmount a highest barrier, which constitutes the rate-determining for this path. The rate-determining step of the C—C bond activation path is the C—C bond break, which has an energy barrier of 44.18 kJ/mol. The formation of stable precursor (H6C2-Ni+(CO)) is the reaction with the higher exothermicity, −281.9 kJ/mol, and the C-C activation path is the most favorable pathway.
  • (2)
    The ELF and AIM analysis give important insights into the chemical bonds evolution of the reactions. The results of the analysis indicated that the initial complex, Ni+(C3H6O), is formed by electrostatic interaction, because no bond is formed between Ni atom and O atom of CH3CH2CHO.

Acknowledgments

We are grateful to the financial support from the National Natural Science Foundation of China-China (Grant No. 21263023) and support from the Supercomputing Center of Gansu Province.

Footnotes

Appendix ASupplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.comptc.2017.05.030.

Appendix A. Supplementary material

Supplementary data 1:

.

References

1. Ertl G. Wiley; Hoboken: 2009. Reactions at Solid Surfaces.
2. Eller K, Schwarz H. Organometallic chemistry in the gas phase. Chem. Rev. 1991;91:1121–1177.
3. Haynes Chris L, Fisher Ellen R, Armentrout PB. Probing the [CoC2H6]+ Potential energy surface: a detailed guided-ion beam study. J. Am. Chem. Soc. 1996;118:3269–3280.
4. Cho Han-Gook, Andrews Lester. Infrared spectra of metallacyclopropane, insertion, and dihydrido complex products in reactions of laser-ablated group 6 metal atoms with ethylene molecules. J. Phys. Chem. A. 2008;112:12071–12080. [PubMed]
5. Richard M Cox, P.B. Armentrout, Wibe A. de Jong, Activation of CH4 by Th+ as Studied by Guided Ion Beam Mass Spectrometry and Quantum Chemistry, Inorg. Chem. 54 (2015) 3584-3599. [PubMed]
6. Sievers MR, Jarvis LM, Armentrout PB. Transition-metal ethene bonds: thermochemistry of M+(C2H4)n (M = Ti-Cu, n = 1 and 2) complexes. J. Am. Chem. Soc. 1998;120:1891–1899.
7. Zhang Qiang, Bowers Michael T. activation of methane by MH+ (M  =  Fe Co, and Ni): a combined mass spectrometric and DFT Study. J. Phys. Chem. A. 2004;108:9755–9761.
8. Andrews Lester. Matrix preparation and spectroscopic and theoretical investigations of simple methylidene and methylidyne complexes of group 4–6 transition metals. Organometallics. 2006;25:4040.
9. Holthausen Max C., Koch Wolfram. A theoretical view on Co+-mediated C-C and C-H bond activations in ethane. J. Am. Chem. Soc. 1996;118:9932–9940.
10. Zhang Dongju, Liu Chengbu, Bi Siwei, Yuan Shiling. A Comprehensive theoretical study on the reactions of Sc+ with CnH2n+2 (n = 1–3): structure, mechanism, and potential-energy surface. Chem. Eur. J. 2003;9:484. [PubMed]
11. Zhou Shaodong, Li Jilai, Schlangen Maria, Schwarz Helmut. thermal activation of methane by [HfO]+ and [XHfO]+ (X = F, Cl, Br, I) and the origin of a remarkable ligand effect. Angew. Chem. Int. Ed. 2016;55:7685–7688. [PubMed]
12. Schwarz Helmut, González-Navarrete Patricio, Li Jilai, Schlangen Maria, Sun Xiaoyan, Weiske Thomas, Zhou Shaodong. unexpected mechanistic variants in the thermal gas-phase activation of methane. Organometallics. 2017;36:8–17.
13. Steinmetz Marc, Grimme Stefan. benchmark study of the performance of density functional theory for bond activations with (ni, pd)-based transition-metal catalysts. ChemistryOpen. 2013;2:115–124. [PubMed]
14. Jana Roithova′, Detlef Schro¨der, Selective activation of alkanes by gas-phase metal ions, Chem. Rev. 110 (2010) 1170-1211. [PubMed]
15. Kang H, Beauchamp JL. Gas-phase studies of alkene oxidation by transition-metal oxides. Ion-beam studies of CrO+ J. Am. Chem. Soc. 1986;108:5663. [PubMed]
16. Schroden Jonathan J, Teo Maurice, Floyd H, Davis Dynamics of CO elimination from reactions of yttrium atoms with formaldehyde, acetaldehyde, and acetone. J. Chem. Phys. 2002;117:9258.
17. Sunderlin LS, Armentrout PB. Reactions of manganese(1+) with isobutene, neopentane, acetone, cyclopropane and ethylene oxide, bond energies for MnCH2+, MnH, and MnCH3. J. Phys. Chem. 1990;94:3589.
18. Carpenter Catherine J., van Koppen Petra A.M., Bowers Michael T. Details of potential energy surfaces involving C-C bond activation: reactions of Fe+, Co+, and Ni+ with acetone. J. Am. Chem. Soc. 1995;117:10976.
19. Jason Dee S., Castleberry Vanessa A., Villarroel Otsmar J., Laboren Ivanna E., Frey Sarah E., Ashley Daniel, Bellert Darrin J. Rate-limiting step in the low-energy unimolecular decomposition reaction of Ni+ acetone into Ni+CO+ ethane. J. Phys. Chem. A. 2009;113:14074–14080. [PubMed]
20. Castleberry Vanessa A., Jason Dee S., Villarroel Otsmar J., Laboren Ivanna E., Frey Sarah E., Bellert Darrin J. The low-energy unimolecular reaction rate constants for the gas phase, Ni+-mediated dissociation of the C-C σ bond in acetone. J. Phys. Chem. A. 2009;113:10417–10424. [PubMed]
21. Chen Xiangfeng, Zang Hao, Yeung Hoi-Sze, Xiaoqing Lu., Dominic Chan T-W. Reaction pathways of Sc+(3D, 1D) and Fe+(6D, 4F) with acetone in the gas phase: metal ion oxidation and acetone deethanization. J. Mass Spectrom. 2012;47:1518–1525. [PubMed]
22. Schilling JB, Beauchamp JL. Hydrocarbon activation by gas-phase lanthanide cations: interaction of praseodymium (Pr+), europium (Eu+), and gadolinium (Gd+) with small alkanes, cycloalkanes, and alkenes. J. Am. Chem. Society. 1988;110:15–24.
23. Zhao Lianming, Zhang Rongrong, Guo Wenyue, Wu Shujuan, Lu Xiaoqing. Does the Co+-assisted decarbonylation of acetaldehyde occur via C-C or C-H activation ?, a theoretical investigation using density functional theory. Chem. Phys. Lett. 2005;414:28.
24. Zhao Lianming, Guo Wenyue, Zhang Rongrong, Wu Shujuan, Lu Xiaoqing. Theoretical Investigation of the decarbonylation of acetaldehyde by Fe+and Cr+ ChemPhysChem. 2006;7:1345. [PubMed]
25. Halle L.F., Crowe W.E., Armentrout P.B., Beauchamp J.L. Reactions of atomic cobalt ions with aldehydes and ketones. Observation of decarbonylation processes leading to formation of metal alkyls and metallacycles in the gas phase. Organometallics. 1984;3:1694.
26. Chen Xiangfeng, Guo Wenyue, Zhao Lianming, Fu Qingtao, Ma Yan. Reaction of acetaldehyde with Ni+: an extended theoretical study of the decarbonylation mechanism of acetaldehyde by first-row transition metal ions. J. Phys. Chem. A. 2007;111:3566–3570. [PubMed]
27. Sonnenfroh' D.M., Farrar J.M. Crossed-beam studies of energy and angular distributions of organometallic reactions: decarbonylation of acetaldehyde by iron (I) and chromium (I) J. Am. Chem. Soc. 1986;108:3521–3522.
28. Jason Dee S., Castleberry Vanessa A., Villarroel Otsmar J., Laboren Ivanna E., Bellert Darrin J. Low-energy reaction rate constants for the ni+-assisted decomposition of acetaldehyde: observation of C-H and C-C activation. J. Phys. Chem. A. 2010;114:1783–1789. [PubMed]
29. Mansell A., Theis Z., Gutierrez M.G., Nieto O., Faza, Silva Lopez C., Bellert D.J. Submerged barriers in the Ni+ assisted decomposition of propionaldehyde. J. Phys. Chem. A. 2016;120:2275–2284. [PubMed]
30. Becke Axel D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993;98:5648.
31. Becke A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A. 1988;38:3098–3100. [PubMed]
32. Lee Chengteh, Yang Weitao, Parr Robert G. Development of the colle-salvetti correlationenergy formula into a functional of the electron density. Phys. Rev. B. 1988;37:785–789. [PubMed]
33. Chen Xiangfeng, Guo Wenyue, Zhao Lianming, Fu Qingtao. Theoretical survey of the potential energy surface of Ni+ + acetone reaction. Chem. Phys. Lett. 2006;432:27–32.
34. Rodríguez-Santiago L., Noguera M., Sodupe M., Salpin J.Y., Tortajada J. Gas phase reactivity of Ni+ with urea. Mass spectrometry and theoretical studies. J. Phys. Chem. A. 2003;107:9865–9874.
35. Luna Alberto, Alcamí Manuel, Mó Otilia, Yáñez Manuel, Tortajada Jeanine. A theoretical study of the interaction between Ni+ and small oxygen- and nitrogen-containing bases. Int. J. Mass Spectrom. 2002;217:119.
36. Gonzalez Carlos, Schlegel H.Bernhard. An improved algorithm for reaction path following. J. Chem. Phys. 1989;90:2154.
37. Gonzalez Carlos, Schlege H.Bernhard. Reaction path following in mass-weighted internal coordinates. J. Phys. Chem. 1990;94:5523–5527.
38. M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G.A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H.P. Hratchian, A.F. Izmaylov, J. Bloino, G. Zheng, J.L. Sonnenberg, M.Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J.A. Montgomery Jr., J.E. Peralta, F. Ogliaro, M. Bearpark, J.J. Heyd, E. Brothers, K.N. Kudin, V.N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J.C. Burant, S.S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J.M. Millam, M. Klene, J.E. Knox, J.B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R.E. Stratmann, O. Yazyev, A.J. Austin, R. Cammi, C. Pomelli, J.W. Ochterski, R.L. Martin, K. Morokuma, V.G. Zakrzewski, G.A. Voth, P. Salvador, J.J. Dannenberg, S. Dapprich, A.D. Daniels, Ö. Farkas, J.B. Foresman, J.V. Ortiz, J. Cioslowski, D.J. Fox, Gaussian 09, Revision D.01, Gaussian Inc., Wallingford, CT, 2009.
39. Lu Tian, Chen Feiwu. Multiwfn: a multifunctional wavefunction analyzer. J. Comput. Chem. 2012;33:580. [PubMed]
40. Li Peng, Niu Wenxia, Tian Xiaofeng, Gao Tao, Wang Hongyan. Ab Initio molecular dynamics study of the reaction of U+ and U2+ with H2O in the gas phase: direct classical trajectory calculations. J. Phys. Chem. A. 2013;117:3761–3770. [PubMed]
41. Shakerzadeh Ehsan. A DFT study on the formaldehyde (H2CO and (H2CO)2) monitoring using pristine B12N12 nanocluster. Physica E. 2016;78:1–9.
42. Cremer Dieter, Kraka Elfi. Chemical bonds without bonding electron density-does the difference electron-density analysis suffice for a description of the chemical bond ? Angew. Chem. Int. Ed. Engl. 1984;23:627–628.
43. Di Santo Emanuela, del Carmen Maria, Michelini Nino Russo. Methane C-H bond activation by gas-phase Th+ and U+: reaction mechanisms and bonding analysis. Organometallics. 2009;28:3716–3726.
44. Stalke Dietmar. Meaningful structural descriptors from charge density. Chem. Eur. J. 2011;17:9264–9278. [PubMed]
45. Bader RFW. A Quantum Theory; Clarendon, Oxford: 1990. Atoms in Molecules.
46. K.J. de Almeida, T.C. Ramalho, J.L. Neto, R.T. Santiago, V.C. Felicíssimo, H.A. Duarte, Neto, Methane Dehydrogenation by Niobium Ions: a First-Principles Study of the Gas-Phase Catalytic Reactions, Organometallics. 32 (2013) 989-999.
47. Du Jiguang, Sun Xiyuan, Chen Jun, Zhang Li, Jiang Gang. An icosahedral Ta122+ cluster with spherical aromaticity. Dalton Trans. 2014;43:5574–5579. [PubMed]
48. Wang Xiaoli, Wang Yongcheng, Li Shuang, Zhang Yuwei, Ma Panpan. Theoretical study on the reaction mechanism of Ti with CH3CN in the Gas Phase. J. Phys. Chem. A. 2016;120:5457–5463. [PubMed]
49. Niu Wenxia, Zhang Hong, Li Peng, Gao Tao. Gas-Phase ammonia activation by Th, Th+, and Th2+: reaction mechanisms, bonding analysis, and rate constant calculations. Int. J. Quantum. Chem. 2015;115:6–18.