Density Functional Theory Study of Competitive Reaction Pathways of Ti<sup>+</sup> with Fluorinated Acetone in the Gas Phase
Density Functional Theory Study of Competitive Reaction Pathways of Ti+ with Fluorinated Acetone in the Gas Phase
Journal of the Korean Chemical Society. 2012. Feb, 56(1): 14-19
Copyright © 2012, The Korean Chemical Society
  • Received : November 14, 2011
  • Accepted : November 25, 2011
  • Published : February 20, 2012
Export by style
Cited by
About the Authors
Kiryong, Hong
Tae Kyu, Kim

We investigate the doublet and quartet potential energy surfaces associated with the gas-phase reaction between Ti + and CF 3 COCH 3 for two plausible reaction pathways, TiF 2 + and TiO + formation pathways by using the density functional theory (DFT) method. The molecular structures of intermediates and transition states involved in these reaction pathways are optimized at the DFT level by using the PBE0 functional. All transition states are identified by using the intrinsic reaction coordinate (IRC) method, and the resulting reaction coordinates describe how Ti + activates CF 3 COCH 3 and produces TiF 2 + and TiO + as products. On the basis of presented results, we propose the most favorable reaction pathway in the reaction between Ti + and CF 3 COCH 3 .
Recent studies of the chemical reactions between transition metal ions and various organic molecules in the gas phase have yielded important insights, not only into the catalytic activity of transition metal ions, but into the reaction mechanisms and the structures of complexes that are present in many important organometallic reactions. 1 - 4 For this purpose, there have been many theoretical studies aiming at explaining reaction mechanisms of the first-row transition-metal ions with simple organic molecules. 5 - 10 Among these studies, the reactions of metal cations with acetone (ACT) have been extensively investigated and reaction pathways appear to be varied depending on the metal ion involved. 8 - 10 Due to the strong binding ability of the early transition metal ions such as Ti + and Sc + , the reaction of Ti + with acetone gives the reaction products of TiO + and CH 2 CHCH 3 . Recently it has been shown that density functional theory (DFT) has been used to study the oxidation pathway of acetone with Ti + and the resulting reaction coordinates describe how Ti + activates the C=O bond of CH 3 COCH 3 and yields TiO + and CH 2 CHCH 3 as products. 8 This theoretical study can explain experimental findings and give a deeper understanding to the reaction mechanism. 8
When H atoms in CH 3 COCH 3 are substituted by highly electronegative F atoms, it is expected that reaction pathways are altered significantly due to the strong binding affinity of Ti + ion to the F atoms. Therefore the reaction of Ti + with CF 3 COCH 3 (1,1,1-trifluoroacetone, TFA) provides a good model system to examine the influence of fluorine substituents on the chemical reactivity of Ti + with acetone. Our recent experimental results 11 show that both TiF 2 + and TiO + ions can be produced in the reaction of Ti + and TFA. This finding implies that the fluorine atom substitution can significantly alter the reaction pathways and related energetics in the reaction between Ti + and TFA. Here we have investigated reaction pathways and relevant potential energy surfaces (PES) of reactions between Ti + and TFA in the gas phase using the DFT method. Based on calculated results, we propose the most favorable reaction pathway in the reaction between Ti + and TFA.
The molecular geometries and harmonic vibrational frequencies were calculated by using DFT. The PBE0 (Hybrid, 25% of the exact exchange energy) functional 12 , 13 was used for the reaction of Ti + with TFA. The structures of all the reactants, products, intermediates, and transition states involved in the title reaction were fully optimized at the PBE0/6-311++G(d,p) levels of theory. Considering the possible spin-orbit coupling interactions that allow an intersystem crossing from the quartet to the doublet in the course of Ti + with TFA, both the doublet and quartet PESs for reaction of Ti + with TFA were searched. 8 - 10 The calculation details have been described in our previous work. 8 , 10 Once we found the association complexes between Ti + and TFA, we verified transition state to connect a desired pair of minima by tracing the intrinsic reaction coordinate (IRC). 14 Frequency calculations were performed at the same level to identify the stationary points and to estimate zero-point energy (ZPE) corrections that are applied to all reported energies. Natural bond orbital (NBO) analysis was also carried out for some important stationary pointsto gain electron configurations and bonding properties. All DFT calculations were performed using the GAUSSIAN 03 package. 15
Based on our experimental results for ion-molecule reaction between Ti + and TFA, here we estimated and performed electronic structure calculations for two plausible reaction pathways as followings:
PPT Slide
Lager Image
Unlike the Ti + + ACT reaction, in which only TiO + ions are generated by insertion of Ti + ions into the C=O bond of ACT molecules, 10 both TiF 2 + and TiO + ions are observed in the reaction of Ti + and TFA. 11 The optimized geometries of the stationary points in both doublet and quartet states for TiO + and TiF 2 + formation reaction pathways are depicted in . 1 and 2 , respectively. . 3 (A) and 3 (B) show the relevant PES along with the TiO + and TiF 2 + formation reaction pathways, respectively. In all cases, the superscript denotes the spin multiplicity. In typical ionmolecule experimental condition, 10 it is highly unlikely that the radical form of isomer can be generated. Therefore we described the most stable species as shown below. In this paper, we used the notations of IM n and TS nm for the intermediate n and the transition state between IM n and IM m , respectively. First, the excited state of Ti + ( 2 F) was calculated to lie 15.9 kcal·mol -1 above the ground state of Ti + ( 4 F) at the PBE0/6-311++G(d,p) levels of theory, which was in well agreement with the experimental value. This demonstrates the accuracy of our calculated results. 8
To form a stable association complex, initially a Ti + ion attacks the electron-rich oxygen of the TFA. 2 IM1 and 4 IM1 complexes are calculated to be more stable by 30.5 kcal/mol ( 4 IM1) and 45.2 kcal/mol ( 2 IM1), respectively, than the isolated reactants. The complex 4 IM 1 has a nearly linear binding of Ti + -O-C (170.2 o ), remarkably different from the bent geometry of the Ni + -CH 3 COCH 3 association complex (∠Ni + -O-C=138.8 o ), which produces Ni + CO+C 2 H 6 fragments. 9 Once the association complexes formed, the next step carries the association complex to the fluoridecontaining species IM2 via transition state TS12. As shown in . 1 , this process is strictly analogous to the fluorine atom-migration for HMg + -OCHCH 2 from Mg + -OCHCH 3 , in which the transition state has a five-membered structure. 16 It could be assumed that intersystem crossing (ISC) takes place between the ground spin state ( 4 F) and the doublet spin state ( 2 F) reaction pathways during process between TS12 and IM2. Similarly, intersystem crossing between the doublet and quartet spin states also has been observed in the gas phase ion-molecule reactions of Ti + with CH 4 , CH 3 OCH 3 , and CH 3 CHO. 5 , 7 , 17 After this intersystem crossing, two alternative pathways can be initiated along the reaction coordinate.
The PES for TiO + formation pathway is shown in . 3 (A). Following the TiO + formation pathway, the insertion of Ti + into the C1=O bond of TFA to produce the IM3 via 2 TS23 and 4 TS23 with a process barrier of 38.3 and 66.0 kcal·mol -1 , respectively. The bond distance of C1-O in TFA is substantially elongated in TS23 compared to those in IM2. The reaction barrier between 2 IM2 and 2 TS23 constituted the largest barrier on the doublet PES of TiO + formation pathway, suggesting that the Ti + insertion into C1=O bond is the rate-determining step of this pathway. After IM3, an electrostatic complex (IM4) between TiO + and CF 2 CFCH 3 , the direct precursor of products in this pathway, is produced by the intra-rotation of Ti + -F bond. The structure of IM4 is well characterized by five-membered ring structure. The dissociation of the (CF 2 CHCH 3 )-TiO + bond of IM4 gives rise to TiO + + CF 2 CFCH 3 and the TiO + pathway is computed to be exothermic by 38.5 kcal·mol -1 for TiO + ( 2 Δ) +CF 2 CFCH 3 and 46.5 kcal·mol -1 endothermic for producing TiO + ( 4 Δ) + CF 2 CFCH 3 .
Another reaction pathway is TiF 2 + formation pathway to account for final products of TiF 2 + + FCOCHCH 2 . The PES for this pathway is displayed in . 3 (B). Like the TiO + formation pathway, this pathway also starts from 4 IM1, passes through ISC and produces 2 IM2. Then, the C1-O-Ti + scissor vibration in IM2 transfers the second F atom to the metal center to form complex IM5 through 2 TS25 or 4 TS25. The five-membered doublet transition state 2 TS25 ( Erel = -82.0 kcal·mol -1 ) is more stable than 4 TS25. It is interesting to note that both 4 IM5 and 2 IM5 indeed have an almost identical stability ( 2 IM5 is 0.2 kcal·mol -1 stable than 4 IM5) as well as same structure whose Ti + -O bond is largely weakened ( RCO =1.957-1.960Å): thus effectively isolating the TiF 2 + from the association complex. The forward intermediate from IM5 is IM6, which involves a three-membered ring 2 TS56 and 4 TS56 with a barrier of 17.8 and 50.9 kcal·mol -1 , respectively. An interaction takes place between O and both C1 and C2 in the TS56. The energy of 2 IM6 is 10.3 kcal·mol -1 lower than that of the entrance channel. The H atom then migrates into C1 from C3 atoms, forming electrostatic complex, IM7 between TiF 2 + and FCOCHCH 2 , the direct precursor of products in the TiF 2 + formation pathway (TiF 2 + + FCOCHCH 2 ). Although the reaction barrier between IM7 ( Erel = -155.8 kcal·mol -1 ) and products represents the deepest barrier along the doublet surface of the TiF 2 + formation pathway, this effect can be attributed to the strong stabilization of both the constituent entities ( 2 TiF 2 + and FCOCHCH 2 ). Therefore, the actual rate determining step in this pathway exists between IM5 and TS56, which gives the reaction barrier of 17.8 kcal·mol -1 . Fragmentation of the (FCOCHCH 2 )-TiF 2 + bond of IM7 gives rise to the products of this pathway and the overall reaction energy of the 4 Ti + +CF 3 COCH 3 2 TiF 2 + + FCOCHCH 2 is calculated to be -100.2 kcal·mol -1 .
PPT Slide
Lager Image
Optimized geometrical parameters (bond lengths in Å and angles in degrees) for the stationary points found in both doublet and quartet electronic states for TiO+ formation pathway at the PBE0/6-311++G(d,p) level.
PPT Slide
Lager Image
Optimized geometrical parameters (bond lengths in Å and angles in degrees) for the stationary points found in both doublet and quartet electronic states for TiF2+ formation pathway at the PBE0/6-311++G(d,p) level.
PPT Slide
Lager Image
Potential energy diagrams along TiO+ (A) and TiF2+ (B) formation pathways.
In summary, a comparison between the values of the relative energies of the highest TSs for both reaction pathways shows that 2 TS23 (-38.3 kcal·mol -1 ) has a larger value than 2 TS56 (-17.8 kcal·mol -1 ), suggesting that TiF 2 + ions are primary reaction products. However, TiF 2 + formation pathway is rather complex than TiO + formation channel in terms of the number of involved reaction steps. Therefore two reaction pathways can be considered as competitive pathways, which is well in agreement with experimental findings. These reaction pathways account for how Ti + activates CF 3 COCH 3 in the gas phase and subsequently produces the competitive reaction products in detail.
This work was supported by for two years by Pusan National University Research Grant.
Roithová J. , Schröder D. 2010 Chem. Rev. 110 1170 -    DOI : 10.1021/cr900183p
Garvey J. F. , Peifer W. R. , Coolbaugh M. T. 1991 Acc. Chem. Res. 24 48 -    DOI : 10.1021/ar00002a004
Koo Y.-M. , Kim T. K. , Jung D. W. , Jung K.-W. 2006 J. Phys. Chem. A 248 810 -
Paul D. , Hong K. , Kim T. K. , Jung K.-W. 2010 Bull. Korean Chem. Soc. 31 271 -    DOI : 10.5012/bkcs.2010.31.02.271
Wang Y.-C. , Liu Z.-Y. , Geng Z.-Y. , Yang X.-Y. 2006 Chem. Phys. Lett. 427 271 -    DOI : 10.1016/j.cplett.2006.06.102
Zhao L. , Guo W. , Zhang R. , Wu S. , Lu X. 2006 Chem. Phys. Chem. 7 1345 -    DOI : 10.1002/cphc.200500720
Ding N. , Zhang S. , Chen X. 2008 Chem. Phys. Lett. 459 33 -    DOI : 10.1016/j.cplett.2008.05.021
Kim J. , Kim T. K. , Ihee H. 2009 J. Phys. Chem. A 113 11382 -    DOI : 10.1021/jp905158j
Chen X. , Guo W. , Zhao L. , Fu Q. 2006 Chem. Phys. Lett. 432 27 -    DOI : 10.1016/j.cplett.2006.10.031
Koo Y.-M. , Hong K. , Kim T. K. , Jung K.-W. 2010 Bull. Korean Chem. Soc. 31 953 -    DOI : 10.5012/bkcs.2010.31.04.953
Paul D. , Hong K. , Jung K.W. , Kim T. K. Chem. Phys. Lett. submitted
Adamo C. , Barone V. 1999 J. Chem. Phys. 110 6158 -    DOI : 10.1063/1.478522
Ernzerhof M. , Scuseria G. E. 1999 J. Chem. Phys. 110 5029 -    DOI : 10.1063/1.478401
Gonzalez C. , Schlegel H. B. 1990 J. Phys. Chem. 94 5523 -    DOI : 10.1021/j100377a021
Frisch M. J. 2004 Gaussian, Revision E.01. Gaussian, Inc. Wallingford, CT
Liu H.-C , Wang C.-S. , Guo W. , Wu Y.-D. , Yang S. 2002 J. Am. Chem. Soc. 124 3794 -    DOI : 10.1021/ja0122546
Sunderlin L. S. , Armentrout P. B. 1988 J. Phys. Chem. 92 1209 -    DOI : 10.1021/j100316a040