1
|
Liu X, Feng H, Li R, Zhang Q, Wu Y, Pang B. Mechanistic Insights into the Proton Transfer and Substitution Dynamics of N-Atom Center Reactions: A Study of CH 3O - with NH 2Cl. J Phys Chem A 2024. [PMID: 38502933 DOI: 10.1021/acs.jpca.3c08447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024]
Abstract
Bimolecular substitution reactions involving N as the central atom have continuously improved our understanding of substitution dynamics. This work used chemical dynamics simulations to investigate the dynamics of NH2Cl with N as the central atom and the multiatomic nucleophile CH3O- and compared these results with the F- + NH2Cl reaction. The most noteworthy difference is in the competition between proton transfer (PT) and the SN2 pathways. Our results demonstrate that, for the CH3O- + NH2Cl system, the PT pathway is considerably more favorable than the SN2 pathway. In contrast, no PT pathway was observed for the F- + NH2Cl system at room temperature. This can be attributed to the exothermic reaction of the PT pathway for the CH3O- + NH2Cl reaction and is coupled with a more stable transition state compared to the substitution pathway. Furthermore, the bulky nature of the CH3O- group impedes its participation in SN2 reactions, which enhances both the thermodynamic and the dynamic advantages of the PT reaction. Interestingly, the atomic mechanism reveals that the PT pathway is primarily governed by indirect mechanisms, similar to the SN2 pathway, with trajectories commonly trapped in the entrance channel being a prominent feature. These trajectories are often accompanied by prolonged and frequent proton exchange or proton abstraction processes. This current work provides insights into the dynamics of N-centered PT reactions, which are useful in gaining a comprehensive understanding of the dynamics behavior of similar reactions.
Collapse
Affiliation(s)
- Xu Liu
- College of Chemistry, Liaoning University, Shenyang 110036, P. R. China
| | - Huining Feng
- College of Chemistry, Liaoning University, Shenyang 110036, P. R. China
| | - Rui Li
- College of Chemistry, Liaoning University, Shenyang 110036, P. R. China
| | - Qiuju Zhang
- Key Laboratory of Advanced Fuel Cells and Electrolyzers Technology of Zhejiang Province, Ningbo Institute of Material Technology and Engineering, Chinese Academy of Sciences, Ningbo, Zhejiang 315201, China
| | - Yang Wu
- College of Chemistry, Liaoning University, Shenyang 110036, P. R. China
| | - Boxue Pang
- Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China
| |
Collapse
|
2
|
Giricz A, Czakó G, Papp D. Alternating Stereospecificity upon Central-Atom Change: Dynamics of the F - +PH 2 Cl S N 2 Reaction Compared to its C- and N-Centered Analogues. Chemistry 2023; 29:e202302113. [PMID: 37698297 DOI: 10.1002/chem.202302113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Indexed: 09/13/2023]
Abstract
Central-atom effects on bimolecular nucleophilic substitution (SN 2) reactions are well-known in chemistry, however, the atomic-level SN 2 dynamics at phosphorous (P) centers has never been studied. We investigate the dynamics of the F- +PH2 Cl reaction with the quasi-classical trajectory method on a novel full-dimensional analytical potential energy surface fitted on high-level ab initio data. Our computations reveal intermediate dynamics compared to the F- +CH3 Cl and the F- +NH2 Cl SN 2 reactions: phosphorus as central atom leads to a more indirect SN 2 reaction with extensive complex-formation with respect to the carbon-centered one, however, the title reaction is more direct than its N-centered pair. Stereospecificity, characteristic at C-center, does not appear here either, due to the submerged front-side-attack retention path and the repeated entrance-channel inversional motion, whereas the multi-inversion mechanism discovered at nitrogen center is also undermined by the deep Walden-well. At low collision energies, 6 % of the PH2 F products form with retained configuration, mostly through complex-mediated mechanisms, while this ratio reaches 24 % at the highest energy due to the increasing dominance of the direct front-side mechanism and the smaller chance for hitting the deep Walden-inversion minimum. Our results suggest pronounced central-atom effects in SN 2 reactions, which can fundamentally change their (stereo)dynamics.
Collapse
Affiliation(s)
- Anett Giricz
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| | - Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged, H-6720, Hungary
| |
Collapse
|
3
|
Zhen W, Zhao S, Fu G, Wang H, Sun J, Yang L, Zhang J. Effects of Methyl Substitution and Leaving Group on E2/S N2 Competition for Reactions of F - with RY (R = CH 3, C 2H 5, iC 3H 7, tC 4H 9; Y = Cl, I). Molecules 2023; 28:6269. [PMID: 37687098 PMCID: PMC10488877 DOI: 10.3390/molecules28176269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 08/22/2023] [Accepted: 08/24/2023] [Indexed: 09/10/2023] Open
Abstract
The competition between base-induced elimination (E2) and bimolecular nucleophilic substitution (SN2) is of significant importance in organic chemistry and is influenced by many factors. The electronic structure calculations for the gas-phase reactions of F- + RY (R = CH3, C2H5, iC3H7, tC4H9, and Y = Cl, I) are executed at the MP2 level with aug-cc-pVDZ or ECP/d basis set to investigate the α-methyl substitution effect. The variation in barrier height, reaction enthalpy, and competition of SN2/E2 as a function of methyl-substitution and leaving group ability has been emphasized. And the nature of these rules has been explored. As the degree of methyl substitution on α-carbon increases, the E2 channel becomes more competitive and dominant with R varying from C2H5, iC3H7, to tC4H9. Energy decomposition analysis offers new insights into the competition between E2 and SN2 processes, which suggests that the drop in interaction energy with an increasing degree of substitution cannot compensate for the rapid growth of preparation energy, leading to a rapid increase in the SN2 energy barrier. By altering the leaving group from Cl to I, the barriers of both SN2 and E2 monotonically decrease, and, with the increased number of substituents, they reduce more dramatically, which is attributed to the looser transition state structures with the stronger leaving group ability. Interestingly, ∆E0‡ exhibits a positive linear correlation with reaction enthalpy (∆H) and halogen electronegativity. With the added number of substituents, the differences in ∆E0‡ and ∆H between Y = Cl and I likewise exhibit good linearity.
Collapse
Affiliation(s)
| | | | | | | | | | - Li Yang
- State Key Laboratory of Urban Water Resource and Environment, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, China; (W.Z.); (S.Z.); (G.F.); (H.W.); (J.S.)
| | - Jiaxu Zhang
- State Key Laboratory of Urban Water Resource and Environment, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, China; (W.Z.); (S.Z.); (G.F.); (H.W.); (J.S.)
| |
Collapse
|
4
|
Li J, Tu Z, Xiang H, Li Y, Song H. Theoretical studies on the kinetics and dynamics of the BeH + + H 2O reaction: comparison with the experiment. Phys Chem Chem Phys 2023; 25:20997-21005. [PMID: 37503894 DOI: 10.1039/d3cp02322b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
The reaction of BeH+ with background gaseous H2O may play a role in qubit loss for quantum information processing with Be+ as trapped ions, and yet its reaction mechanism has not been well understood until now. In this work, a globally accurate, full-dimensional ground-state potential energy surface (PES) for the BeH+ + H2O reaction was constructed by fitting a total of 170 438 ab initio energy points at the level of RCCSD(T)-F12/aug-cc-pVTZ using the fundamental invariant-neural network method. The total root-mean-square error of the final PES was 0.178 kcal mol-1. For comparison, quasi-classical trajectory calculations were carried out on the PES at an experimental temperature of 150 K. The obtained thermal rate constant and product branching ratio of the BeD+ + H2O reaction agreed quite well with experimental results. In addition, the vibrational state distributions and energy disposals of the products were calculated and rationalized using the sudden vector projection model.
Collapse
Affiliation(s)
- Jiaqi Li
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China.
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Zhao Tu
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
- School of Chemical and Environmental Engineering, Hubei Minzu University, Enshi 445000, China
| | - Haipan Xiang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
- School of Physics and Electronics, Hunan University, Changsha 410082, China
| | - Yong Li
- College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China.
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| |
Collapse
|
5
|
Zhao S, Fu G, Zhen W, Yang L, Sun J, Zhang J. Reaction mechanism conversion induced by the contest of nucleophile and leaving group. Phys Chem Chem Phys 2022; 24:24146-24154. [PMID: 36168813 DOI: 10.1039/d2cp01987f] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Direct dynamic simulations have been employed to investigate the OH- + CH3Cl reaction with the chosen B3LYP/aug-cc-pVDZ method. The calculated rate coefficient for the bimolecular nucleophilic substitution reaction (SN2), 1.0 × 10-9 cm3 mol-1 s-1 at 300 K, agrees well with the experimental result of (1.3-1.6) × 10-9 cm3 mol-1 s-1. The simulations reveal that the majority of the SN2 reactions are temporarily trapped in the hydrogen-bonded complex at Ecoll = 0.89 kcal mol-1. Importantly, the influences of the leaving group and nucleophile have been discussed by comparisons of X- + CH3Y (X = F, OH; Y = Cl, I) reactions. For the X = F- reactions, the reaction probability of SN2 increases along the increased leaving group ability Cl < I, suggesting that the thermodynamic factor plays a key role. The indirect mechanisms were found to be dominant for both reactions. In contrast, for X = OH-, the fraction of SN2 drops with the enhanced leaving group ability. In particular, a dramatic transition occurs for the dominant atomic reaction mechanisms, i.e., from complex-mediated indirect to direct, implying an interesting contest between the leaving group and the nucleophile and the importance of the dynamic factors, i.e., the dipole moment, steric hindrance, and electronegativity.
Collapse
Affiliation(s)
- Siwei Zhao
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Gang Fu
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Wenqing Zhen
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Li Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China. .,State Key Lab of Advanced Welding and Joining, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Jianmin Sun
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China. .,State Key Lab of Urban Water Resource and Environment, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Jiaxu Zhang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China. .,State Key Lab of Advanced Welding and Joining, Harbin Institute of Technology, Harbin 150001, P. R. China
| |
Collapse
|
6
|
Li Y, Li C, Gao D, Wang D. Atomic-Level Mechanism, Solvent Effect, and Potential of the Mean Force of the F - + CH 3CH 2Cl S N2 Reaction in Aqueous Solution. J Phys Chem A 2022; 126:5527-5533. [PMID: 35947789 DOI: 10.1021/acs.jpca.2c02105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the bimolecular nucleophilic substitution (SN2) reaction of F- with CH3CH2Cl in aqueous solution using combined multilevel quantum mechanism (ML-QM) theories with molecular mechanics (MM). The synchronized, atomic-level structural and charge evolutions are analyzed along the reaction path. The potential mean force along the reaction path in water is calculated at high-accuracy CCSD(T)/aug-cc-pVTZ/MM level of theory with a free energy barrier of 16.8 kcal/mol and a free energy of reaction of -9.7 kcal/mol. The water solvent hinders the reactivity by raising its reaction barrier by 15.1 kcal/mol, of which 13.6 kcal/mol comes from solvent energy contribution and 1.5 kcal/mol comes from the polarization effect. This indicates that the water solvent plays an essential role on this reaction in aqueous solution. We also predict the potential mean force profile based on the gas-phase reaction path and the solvation free energies of the stationary points; the comparison between our calculated result at CCSD(T)/MM level shows an excellent agreement with the predicted one with the free energy barrier at 16.2 kcal/mol and the free energy of reaction at -8.3 kcal/mol.
Collapse
Affiliation(s)
- Yixuan Li
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Chen Li
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Delu Gao
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| |
Collapse
|
7
|
Unexpected steric hindrance failure in the gas phase F - + (CH 3) 3CI S N2 reaction. Nat Commun 2022; 13:4427. [PMID: 35907925 PMCID: PMC9338938 DOI: 10.1038/s41467-022-32191-6] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Accepted: 07/18/2022] [Indexed: 11/13/2022] Open
Abstract
Base-induced elimination (E2) and bimolecular nucleophilic substitution (SN2) reactions are of significant importance in physical organic chemistry. The textbook example of the retardation of SN2 reactivity by bulky alkyl substitution is widely accepted based on the static analysis of molecular structure and steric environment. However, the direct dynamical evidence of the steric hindrance of SN2 from experiment or theory remains rare. Here, we report an unprecedented full-dimensional (39-dimensional) machine learning-based potential energy surface for the 15-atom F− + (CH3)3CI reaction, facilitating the reliable and efficient reaction dynamics simulations that can reproduce well the experimental outcomes and examine associated atomic-molecular level mechanisms. Moreover, we found surprisingly high “intrinsic” reactivity of SN2 when the E2 pathway is completely blocked, indicating the reaction that intends to proceed via E2 transits to SN2 instead, due to a shared pre-reaction minimum. This finding indicates that the competing factor of E2 but not the steric hindrance determines the small reactivity of SN2 for the F− + (CH3)3CI reaction. Our study provides new insight into the dynamical origin that determines the intrinsic reactivity in gas-phase organic chemistry. Base-induced elimination (E2) and bimolecular nucleophilic substitution (SN2) are of significant importance in physical organic chemistry. Here, the authors show that the competing factor of E2 as opposed to steric hindrance determines the low reactivity of SN2 in the F− + (CH3)3CI reaction.
Collapse
|
8
|
Valverde D, Georg HC, Canuto S. Free-Energy Landscape of the S N2 Reaction CH 3Br + Cl - → CH 3Cl + Br - in Different Liquid Environments. J Phys Chem B 2022; 126:3685-3692. [PMID: 35543431 DOI: 10.1021/acs.jpcb.1c10282] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
This work describes in detail the reaction path of the well-known SN2 reaction CH3Br + Cl- → CH3Cl + Br-, whose reaction rate has a huge variation with the solvent in the gas phase and in protic and aprotic liquid environments. We employed the ASEC-FEG method to optimize for minima (reactants and products) and saddle points (transition states) in the in-solution free-energy hypersurface. The method takes atomistic details of the solvent into account. A polarizable continuum model (PCM) has also been employed for comparison. The most perceptive structural changes are noted in aqueous solution by using the ASEC-FEG approach. The activation energies in all solvents, estimated by means of free-energy perturbation calculations, are in good agreement with the experimental data. The total solute-solvent hydrogen bonds play an important role in the increased barrier height observed in water and are therefore crucial to explain the huge decrease in the kinetic constant. It is also found that the hydration shell around the ions breaks itself spontaneously to accommodate the molecule, thus forming minimum energy complexes.
Collapse
Affiliation(s)
- Danillo Valverde
- Instituto de Física, Universidade de São Paulo, Rua do Matão 1371 Cidade Universitária, CEP 05508-090 São Paulo, São Paulo, Brazil
| | - Herbert C Georg
- Instituto de Física, Universidade Federal de Goiás, Avenida Esperança, Campus Samambaia, CEP 74690-900 Goiânia, Goiás, Brazil
| | - Sylvio Canuto
- Instituto de Física, Universidade de São Paulo, Rua do Matão 1371 Cidade Universitária, CEP 05508-090 São Paulo, São Paulo, Brazil
| |
Collapse
|
9
|
Tachikawa H. Reaction mechanism of an intracluster S N2 reaction induced by electron capture. Phys Chem Chem Phys 2022; 24:3941-3950. [PMID: 35098286 DOI: 10.1039/d1cp04697g] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Bimolecular nucleophilic substitution (SN2) reactions have been widely investigated from both experimental and theoretical points of view because they represent one of the simplest organic reactions. Most studies on SN2 reactions have been focused on bimolecular collision. In contrast, information on intracluster SN2 reactions is limited. In this study, an intracluster SN2 reaction of NF3-CH3Cl triggered by electron attachment was investigated using a direct ab initio molecular dynamics (AIMD) method. In the structure of NF3-CH3Cl, the N-F bond in NF3 is oriented collinearly toward the carbon atom of CH3Cl. After electron capture by NF3-CH3Cl, the F- ion that is generated from the (NF3)- moiety collides with the carbon atom of CH3Cl. The intracluster SN2 reaction occurs as follows: (NF3-CH3Cl)- (electron capture state) → NF2-(F-)-CH3Cl (pre-reaction complex) → transition state (TS) → NF2-CH3F-Cl- (post-reaction complex) → NF2 + CH3F + Cl- (product state). The reaction energy is efficiently transferred to the translational mode of Cl-, and the Cl- ion with a high translational energy is then removed from the system. This energy is significantly larger than that of Cl- formed in the bimolecular SN2 reaction (F- + CH3Cl). The reaction mechanism is discussed based on the theoretical results.
Collapse
Affiliation(s)
- Hiroto Tachikawa
- Division of Applied Chemistry, Faculty of Engineering, Hokkaido University, Sapporo 060-8628, Japan.
| |
Collapse
|
10
|
Zhu Y, Li R, Song H. Kinetic and dynamic studies of the NH 2+ + H 2 reaction on a high-level ab initio potential energy surface. Phys Chem Chem Phys 2022; 24:25663-25672. [DOI: 10.1039/d2cp03859e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The dynamics and kinetics of the NH2+ + H2 reaction are investigated on a newly developed ab initio potential energy surface using the quasi-classical trajectory method.
Collapse
Affiliation(s)
- Yongfa Zhu
- School of Chemistry and Chemical Engineering, Hubei Polytechnic University, Huangshi 435003, China
| | - Rui Li
- School of Chemistry and Chemical Engineering, Hubei Polytechnic University, Huangshi 435003, China
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
| |
Collapse
|
11
|
Dékány AÁ, Kovács GZ, Czakó G. High-Level Systematic Ab Initio Comparison of Carbon- and Silicon-Centered S N2 Reactions. J Phys Chem A 2021; 125:9645-9657. [PMID: 34709818 PMCID: PMC8591615 DOI: 10.1021/acs.jpca.1c07574] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We characterize the
stationary points along the Walden inversion,
front-side attack, and double-inversion pathways of the X– + CH3Y and X– + SiH3Y [X,
Y = F, Cl, Br, I] SN2 reactions using chemically accurate
CCSD(T)-F12b/aug-cc-pVnZ [n = D,
T, Q] levels of theory. At the carbon center, Walden inversion dominates
and proceeds via prereaction (X–···H3CY) and postreaction (XCH3···Y–) ion-dipole wells separated by a usually submerged
transition state (X–H3C–Y)−, front-side attack occurs over high barriers, double inversion is
the lowest-energy retention pathway for X = F, and hydrogen- (F–···HCH2Y) and halogen-bonded
(X–···YCH3) complexes
exist in the entrance channel. At the silicon center, Walden inversion
proceeds through a single minimum (X–SiH3–Y)−, the front-side attack is competitive via a usually
submerged transition state separating pre- and postreaction minima
having X–Si–Y angles close to 90°, double inversion
occurs over positive, often high barriers, and hydrogen- and halogen-bonded
complexes are not found. In addition to the SN2 channels
(Y– + CH3X/SiH3X), we report
reaction enthalpies for proton abstraction (HX + CH2Y–/SiH2Y–), hydride substitution
(H– + CH2XY/SiH2XY), XH···Y– complex formation (XH···Y– + 1CH2/1SiH2), and halogen
abstraction (XY + CH3–/SiH3– and XY– + CH3/SiH3).
Collapse
Affiliation(s)
- Attila Á Dékány
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gyula Z Kovács
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
12
|
Tajti V, Győri T, Czakó G. Detailed quasiclassical dynamics of the F - + CH 3Br reaction on an ab initio analytical potential energy surface. J Chem Phys 2021; 155:124301. [PMID: 34598562 DOI: 10.1063/5.0065209] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Dynamics and mechanisms of the F- + CH3Br(v = 0) → Br- + CH3F (SN2 via Walden inversion, front-side attack, and double inversion), F- + inverted-CH3Br (induced inversion), HF + CH2Br- (proton abstraction), and FH⋯Br- + 1CH2 reactions are investigated using a high-level global ab initio potential energy surface, the quasiclassical trajectory method, as well as non-standard configuration- and mode-specific analysis techniques. A vector-projection method is used to identify inversion and retention trajectories; then, a transition-state-attack-angle-based approach unambiguously separates the front-side attack and the double-inversion retention pathways. The Walden-inversion SN2 channel becomes direct rebound dominated with increasing collision energy as indicated by backward scattering, initial back-side attack preference, and the redshifting of product internal energy peaks in accord with CF stretching populations. In the minor retention and induced-inversion pathways, almost the entire available energy transfers into product rotation-vibration, and retention mainly proceeds with indirect, slow double inversion following induced inversion with about 50% probability. Proton abstraction is dominated by direct stripping (evidenced by forward scattering) with CH3-side initial attack preference, providing mainly vibrationally ground state products with significant zero-point energy violation.
Collapse
Affiliation(s)
- Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Tibor Győri
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
13
|
Pan M, Xiang H, Li Y, Song H. Study on the kinetics and dynamics of the H 2 + NH 2- reaction on a high-level ab initio potential energy surface. Phys Chem Chem Phys 2021; 23:17848-17855. [PMID: 34612274 DOI: 10.1039/d1cp02423j] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gas-phase ion-molecule reactions play major roles in many fields of chemistry and physics. The reaction of an amino radical anion with a hydrogen molecule is one of the simplest proton transfer reactions involving anions. A globally accurate full-dimensional potential energy surface (PES) for the NH2- + H2 reaction is developed by the fundamental invariant-neural network method, resulting in a root mean square error of 0.116 kcal mol-1. Quasi-classical trajectory calculations are then carried out on the newly developed PES to give integral cross sections, differential cross sections and thermal rate coefficients. This reaction has two reaction channels, proton transfer and hydrogen exchange. The reactivity of the proton transfer channel is about one or two orders of magnitude stronger than that of the hydrogen exchange channel in the energy range studied. Vibrational excitation of H2 promotes the proton transfer reaction, while fundamental excitation of each vibrational mode of NH2- has a negligible effect. In addition, the theoretical rate coefficients of the proton transfer reaction on the PES show inverse temperature dependence from 150 to 750 K, in accordance with the available experimental results.
Collapse
Affiliation(s)
- Mengyi Pan
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | | | | | | |
Collapse
|
14
|
Alkorta I, Elguero J. The SN2 reaction and its relationship with the Walden inversion, the Finkelstein and Menshutkin reactions together with theoretical calculations for the Finkelstein reaction. Struct Chem 2021. [DOI: 10.1007/s11224-021-01805-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
AbstractThis communication gives an overview of the relationships between four reactions that although related were not always perceived as such: SN2, Walden, Finkelstein, and Menshutkin. Binary interactions (SN2 & Walden, SN2 & Menshutkin, SN2 & Finkelstein, Walden & Menshutkin, Walden & Finkelstein, Menshutkin & Finkelstein) were reported. Carbon, silicon, nitrogen, and phosphorus as central atoms and fluorides, chlorides, bromides, and iodides as lateral atoms were considered. Theoretical calculations provide Gibbs free energies that were analyzed with linear models to obtain the halide contributions. The M06-2x DFT computational method and the 6-311++G(d,p) basis set have been used for all atoms except for iodine where the effective core potential def2-TZVP basis set was used. Concerning the central atom pairs, carbon/silicon vs. nitrogen/phosphorus, we reported here for the first time that the effect of valence expansion was known for Si but not for P. Concerning the lateral halogen atoms, some empirical models including the interaction between F and I as entering and leaving groups explain the Gibbs free energies.
Collapse
|
15
|
Wester R. Fifty years of nucleophilic substitution in the gas phase. MASS SPECTROMETRY REVIEWS 2021; 41:627-644. [PMID: 34060119 PMCID: PMC9291629 DOI: 10.1002/mas.21705] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 05/05/2021] [Accepted: 05/05/2021] [Indexed: 06/12/2023]
Abstract
Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.
Collapse
Affiliation(s)
- Roland Wester
- Institut für Ionenphysik und Angewandte PhysikUniversität InnsbruckTechnikerstraße 256020 InnsbruckAustria
| |
Collapse
|
16
|
Czakó G, Győri T, Papp D, Tajti V, Tasi DA. First-Principles Reaction Dynamics beyond Six-Atom Systems. J Phys Chem A 2021; 125:2385-2393. [PMID: 33631071 PMCID: PMC8028310 DOI: 10.1021/acs.jpca.0c11531] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Revised: 02/05/2021] [Indexed: 11/29/2022]
Abstract
Moving beyond the six-atomic benchmark systems, we discuss the new age and future of first-principles reaction dynamics, which investigates complex, multichannel chemical reactions. We describe the methodology starting from the benchmark ab initio characterization of the stationary points, followed by full-dimensional potential energy surface (PES) developments and reaction dynamics computations. We highlight our composite ab initio approach providing benchmark stationary-point properties with subchemical accuracy, the Robosurfer program system enabling automatic PES development, and applications for the Cl + C2H6, F + C2H6, and OH- + CH3I post-six-atom reactions focusing on ab initio issues and their solutions as well as showing the excellent agreement between theory and experiment.
Collapse
Affiliation(s)
- Gábor Czakó
- MTA-SZTE Lendület
Computational Reaction Dynamics Research Group, Interdisciplinary
Excellence Centre and Department of Physical Chemistry and Materials
Science, Institute of Chemistry, University
of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Tibor Győri
- MTA-SZTE Lendület
Computational Reaction Dynamics Research Group, Interdisciplinary
Excellence Centre and Department of Physical Chemistry and Materials
Science, Institute of Chemistry, University
of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Dóra Papp
- MTA-SZTE Lendület
Computational Reaction Dynamics Research Group, Interdisciplinary
Excellence Centre and Department of Physical Chemistry and Materials
Science, Institute of Chemistry, University
of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Viktor Tajti
- MTA-SZTE Lendület
Computational Reaction Dynamics Research Group, Interdisciplinary
Excellence Centre and Department of Physical Chemistry and Materials
Science, Institute of Chemistry, University
of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Domonkos A. Tasi
- MTA-SZTE Lendület
Computational Reaction Dynamics Research Group, Interdisciplinary
Excellence Centre and Department of Physical Chemistry and Materials
Science, Institute of Chemistry, University
of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
17
|
Papp D, Czakó G. Facilitated inversion complicates the stereodynamics of an S N2 reaction at nitrogen center. Chem Sci 2021; 12:5410-5418. [PMID: 34168784 PMCID: PMC8179618 DOI: 10.1039/d1sc00490e] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Bimolecular nucleophilic substitution (SN2) reactions at carbon center are well known to proceed with the stereospecific Walden-inversion mechanism. Reaction dynamics simulations on a newly developed high-level ab initio analytical potential energy surface for the F− + NH2Cl nitrogen-centered SN2 and proton-transfer reactions reveal a hydrogen-bond-formation-induced multiple-inversion mechanism undermining the stereospecificity of the N-centered SN2 channel. Unlike the analogous F− + CH3Cl SN2 reaction, F− + NH2Cl → Cl− + NH2F is indirect, producing a significant amount of NH2F with retention, as well as inverted NH2Cl during the timescale within the unperturbed NH2Cl molecule gets inverted with only low probability, showing the important role of facilitated inversions via an FH…NHCl−-like transition state. Proton transfer leading to HF + NHCl− is more direct and becomes the dominant product channel at higher collision energies. Multiple-inversion, the analogue of the double-inversion pathway recently revealed for SN2@C, is the key mechanism in SN2 at N center undermining stereospecificity.![]()
Collapse
Affiliation(s)
- Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged Rerrich Béla tér 1 Szeged H-6720 Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged Rerrich Béla tér 1 Szeged H-6720 Hungary
| |
Collapse
|
18
|
Ji X, Zhao C, Xie J. Investigating the role of halogen-bonded complexes in microsolvated Y−(H2O)n + CH3I SN2 reactions. Phys Chem Chem Phys 2021; 23:6349-6360. [DOI: 10.1039/d0cp06299e] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A halogen-bonded complex pathway is computed for Y−(H2O)n + CH3I (Y = HO, F, Cl, Br, and I) ion–molecule nucleophilic substitution reactions and is compared with back-side and front-side attack pathways.
Collapse
Affiliation(s)
- Xiaoyan Ji
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- Beijing 100081
- China
| | - Chongyang Zhao
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- Beijing 100081
- China
| | - Jing Xie
- Key Laboratory of Cluster Science of Ministry of Education
- School of Chemistry and Chemical Engineering
- Beijing Institute of Technology
- Beijing 100081
- China
| |
Collapse
|
19
|
Tasi DA, Győri T, Czakó G. On the development of a gold-standard potential energy surface for the OH− + CH3I reaction. Phys Chem Chem Phys 2020; 22:3775-3778. [DOI: 10.1039/c9cp07007a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
We develop the first accurate full-dimensional ab initio PES for the OH− + CH3I SN2 and proton-transfer reactions treating the failure of CCSD(T) at certain geometries.
Collapse
Affiliation(s)
- Domonkos A. Tasi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Tibor Győri
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| |
Collapse
|
20
|
Czakó G, Győri T, Olasz B, Papp D, Szabó I, Tajti V, Tasi DA. Benchmark ab initio and dynamical characterization of the stationary points of reactive atom + alkane and SN2 potential energy surfaces. Phys Chem Chem Phys 2020; 22:4298-4312. [DOI: 10.1039/c9cp04944d] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review composite ab initio and dynamical methods and their applications to characterize stationary points of atom/ion + molecule reactions.
Collapse
Affiliation(s)
- Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Tibor Győri
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Balázs Olasz
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - István Szabó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| | - Domonkos A. Tasi
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
| |
Collapse
|
21
|
Li L, Fu B, Yang X, Zhang DH. A global ab initio potential energy surface and dynamics of the proton-transfer reaction: OH− + D2 → HOD + D−. Phys Chem Chem Phys 2020; 22:8203-8211. [DOI: 10.1039/d0cp00107d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction mechanisms of OH− + D2 → HOD + D− were first revealed by theory, based on an accurate full-dimensional PES.
Collapse
Affiliation(s)
- Lulu Li
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei 230026
- China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Xueming Yang
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei 230026
- China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
| | - Dong H. Zhang
- Department of Chemical Physics
- University of Science and Technology of China
- Hefei 230026
- China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
| |
Collapse
|
22
|
Niu X, Liu P, Wang D. Multilevel Quantum Mechanics and Molecular Mechanics Study of the Double-Inversion Mechanism at Nitrogen: F- + NH2Cl in Aqueous Solution. J Phys Chem A 2019; 124:141-147. [PMID: 31820988 DOI: 10.1021/acs.jpca.9b09689] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We employ a multilevel quantum mechanics and molecular mechanics method to investigate the double-inversion mechanism of the nucleophilic substitution reaction at the N center: the F- + NH2Cl reaction in aqueous solution. We find that the structures of the stationary points along the reaction path are quite different from the ones in the gas phase owing to the hydrogen-bond interactions between the solute and the surrounding water molecules. The atomic-level evolutions of the structures and charge transfer along the reaction path show that this double-inversion mechanism consists of an upside-down proton inversion process and a Walden-inversion process. The computed potential of mean force at the coupled-cluster singles and doubles with perturbative triples (CCSD(T))/molecular mechanics (MM) level of theory has the two-inversion barrier heights, and reaction free energy at 11.7, 29.6, and 12.6 kcal/mol, agreeing well with the predicted ones at 12.6, 32.5, and 12.2 kcal/mol obtained on the basis of the gas-phase reaction path and the solvation free energies of the stationary points.
Collapse
Affiliation(s)
- Xiao Niu
- College of Physics and Electronics , Shandong Normal University , Jinan , Shandong 250014 , China
| | - Peng Liu
- College of Physics and Electronics , Shandong Normal University , Jinan , Shandong 250014 , China
| | - Dunyou Wang
- College of Physics and Electronics , Shandong Normal University , Jinan , Shandong 250014 , China
| |
Collapse
|
23
|
Győri T, Czakó G. Automating the Development of High-Dimensional Reactive Potential Energy Surfaces with the robosurfer Program System. J Chem Theory Comput 2019; 16:51-66. [DOI: 10.1021/acs.jctc.9b01006] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tibor Győri
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
24
|
Tian L, Zhu Y, Song H, Yang M. Theoretical study of the F( 2P) + NH 3→ HF + NH 2 reaction on an accurate potential energy surface: dynamics and kinetics. Phys Chem Chem Phys 2019; 21:11385-11394. [PMID: 31111138 DOI: 10.1039/c9cp02113b] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The highly exothermic hydrogen abstraction reaction of the F atom with NH3 is investigated using the quasi-classical trajectory method on a newly developed potential energy surface (PES) for the ground electronic state. The full-dimensional PES is constructed by fitting 41 282 ab initio energy points at the level of UCCSD(T)-F12/aug-cc-pVTZ. The flexible fundamental invariant-neural network method is applied in the fitting, resulting in a total root mean square error of 0.13 kcal mol-1. On one hand, the calculated differential cross sections agree reasonably well with the experimental results and indicate that the reaction is dominated by the direct abstraction and stripping mechanisms while a considerable amount of reaction takes place by the indirect "yo-yo" mechanism. The product energy partition also reproduces well the experimental result, which can be understood according to the geometry change along the minimum energy path. On the other hand, the obtained vibrational state distribution of the product HF follows PνHF=2≈PνHF=1 > PνHF=0 > PνHF=3, less consistent with the scattered experimental results. In addition, the calculated thermal rate coefficients have practically no temperature dependence within the statistical errors.
Collapse
Affiliation(s)
- Li Tian
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China. and College of Physical Science and Technology, Huazhong Normal University, Wuhan 430079, China
| | - Yongfa Zhu
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China. and University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Minghui Yang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China.
| |
Collapse
|
25
|
Effect of solvent polarity on the potential energy surface in the SN2 reaction of F− + CH3Cl. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.01.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|
26
|
Pratihar S, Nicola Barbosa Muniz MC, Ma X, Borges I, Hase WL. Pronounced changes in atomistic mechanisms for the Cl - + CH 3I S N2 reaction with increasing collision energy. Phys Chem Chem Phys 2019; 21:2039-2045. [PMID: 30633280 DOI: 10.1039/c8cp06198j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
In a previous direct dynamics simulation of the Cl- + CH3I → ClCH3 + I- SN2 reaction, predominantly indirect and direct reaction was found at collision energies Erel of 0.20 and 0.39 eV, respectively. For the work presented here, these simulations were extended by studying the reaction dynamics from Erel of 0.15 to 0.40 eV in 0.05 eV intervals. A transition from a predominantly indirect to direct reaction is found for Erel of 0.27-0.28 eV, a finding consistent with experiment. The simulation results corroborate the understanding that in experiments indirect reaction is characterized by small product translational energies and isotropic scattering, while direct reaction has higher translational energies and anisotropic scattering. The traditional statistical theoretical model for the Cl- + CH3I SN2 reaction assumes the Cl--CH3I pre-reaction complex (A) is formed, followed by barrier crossing, and then formation of the ClCH3-I- post-reaction complex (B). This mechanism is seen in the dynamics, but the complete atomistic dynamics are much more complex. Atomistic SN2 mechanisms contain A and B, but other dynamical events consisting of barrier recrossing (br) and the roundabout (Ra), in which the CH3-moiety rotates around the heavy I-atom, are also observed. The two most important mechanisms are only formation of A and Ra + A. The simulation results are compared with simulations and experiments for Cl- + CH3Cl, Cl- + CH3Br, F- + CH3I, and OH- + CH3I.
Collapse
Affiliation(s)
- Subha Pratihar
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA.
| | | | | | | | | |
Collapse
|
27
|
Tasi DA, Fábián Z, Czakó G. Rethinking the X− + CH3Y [X = OH, SH, CN, NH2, PH2; Y = F, Cl, Br, I] SN2 reactions. Phys Chem Chem Phys 2019; 21:7924-7931. [DOI: 10.1039/c8cp07850e] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Explicitly-correlated ab initio computations reveal novel inversion and retention pathways for several SN2 reactions with different nucleophiles and leaving groups.
Collapse
Affiliation(s)
- Domonkos A. Tasi
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
- Hungary
| | - Zita Fábián
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
- Hungary
| | - Gábor Czakó
- Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science
- Institute of Chemistry
- University of Szeged
- Szeged H-6720
- Hungary
| |
Collapse
|
28
|
Caballero-García G, Mondragón-Solórzano G, Torres-Cadena R, Díaz-García M, Sandoval-Lira J, Barroso-Flores J. Calculation of VS,max and Its Use as a Descriptor for the Theoretical Calculation of p Ka Values for Carboxylic Acids. Molecules 2018; 24:molecules24010079. [PMID: 30587832 PMCID: PMC6337188 DOI: 10.3390/molecules24010079] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Revised: 12/21/2018] [Accepted: 12/25/2018] [Indexed: 01/30/2023] Open
Abstract
The theoretical calculation of pKa values for Brønsted acids is a challenging task that involves sophisticated and time-consuming methods. Therefore, heuristic approaches are efficient and appealing methodologies to approximate these values. Herein, we used the maximum surface electrostatic potential (VS,max) on the acidic hydrogen atoms of carboxylic acids to describe the H-bond interaction with water (the same descriptor that is used to characterize σ-bonded complexes) and correlate the results with experimental pKa values to obtain a predictive model for other carboxylic acids. We benchmarked six different methods, all including an implicit solvation model (water): Five density functionals and the Møller–Plesset second order perturbation theory in combination with six different basis sets for a total of thirty-six levels of theory. The ωB97X-D/cc-pVDZ level of theory stood out as the best one for consistently reproducing the reported pKa values, with a predictive power of 98% correlation in a test set of ten other carboxylic acids.
Collapse
Affiliation(s)
- Guillermo Caballero-García
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| | - Gustavo Mondragón-Solórzano
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| | - Raúl Torres-Cadena
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| | - Marco Díaz-García
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| | - Jacinto Sandoval-Lira
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| | - Joaquín Barroso-Flores
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM, Unidad San Cayetano, Carretera Toluca⁻Atlacomulco km 14.5, Personal de la UNAM, Toluca 50200, Mexico.
| |
Collapse
|
29
|
Ma X, Di Liberto G, Conte R, Hase WL, Ceotto M. A quantum mechanical insight into SN2 reactions: Semiclassical initial value representation calculations of vibrational features of the Cl−⋯CH3Cl pre-reaction complex with the VENUS suite of codes. J Chem Phys 2018; 149:164113. [DOI: 10.1063/1.5054399] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xinyou Ma
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Giovanni Di Liberto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - William L. Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
30
|
Olasz B, Czakó G. Mode-Specific Quasiclassical Dynamics of the F - + CH 3I S N2 and Proton-Transfer Reactions. J Phys Chem A 2018; 122:8143-8151. [PMID: 30230832 DOI: 10.1021/acs.jpca.8b08286] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Mode-specific quasiclassical trajectory computations are performed for the F- + CH3I( v k = 0, 1) SN2 and proton-transfer reactions at nine different collision energies in the range of 1.0-35.3 kcal/mol using a full-dimensional high-level ab initio analytical potential energy surface with ground-state and excited CI stretching ( v3), CH3 rocking ( v6), CH3 umbrella ( v2), CH3 deformation ( v5), CH symmetric stretching ( v1), and CH asymmetric stretching ( v4) initial vibrational modes. Millions of trajectories provide statistically definitive mode-specific cross sections, opacity functions, scattering angle distributions, and product internal energy distributions. The excitation functions reveal slight vibrational SN2 inversion inhibition/enhancement at low/high collision energies ( Ecoll), whereas large decaying-with- Ecoll vibrational enhancement effects for the SN2 retention (double inversion) and proton-transfer channels. The most efficient vibrational enhancement is found by exciting the CI stretching (high Ecoll) for SN2 inversion and the CH stretching modes (low Ecoll) for double inversion and proton transfer. Mode-specific effects do not show up in the scattering angle distributions and do blue-shift the hot/cold SN2/proton-transfer product internal energies.
Collapse
Affiliation(s)
- Balázs Olasz
- Department of Physical Chemistry and Materials Science, Institute of Chemistry , University of Szeged , Rerrich Béla tér 1 , Szeged H-6720 , Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry , University of Szeged , Rerrich Béla tér 1 , Szeged H-6720 , Hungary
| |
Collapse
|
31
|
Wang Y, Wang D. Quantum dynamics calculations reveal temperature independence of kinetic isotope effect of the OH + HBr/DBr reaction. J Chem Phys 2018; 149:034302. [DOI: 10.1063/1.5037542] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Yuping Wang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong, China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong, China
| |
Collapse
|
32
|
Tasi DA, Fábián Z, Czakó G. Benchmark ab Initio Characterization of the Inversion and Retention Pathways of the OH– + CH3Y [Y = F, Cl, Br, I] SN2 Reactions. J Phys Chem A 2018; 122:5773-5780. [DOI: 10.1021/acs.jpca.8b04218] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Domonkos A. Tasi
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Zita Fábián
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
33
|
Hamlin TA, Swart M, Bickelhaupt FM. Nucleophilic Substitution (S N 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent. Chemphyschem 2018; 19:1315-1330. [PMID: 29542853 PMCID: PMC6001448 DOI: 10.1002/cphc.201701363] [Citation(s) in RCA: 118] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Indexed: 11/12/2022]
Abstract
The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.
Collapse
Affiliation(s)
- Trevor A. Hamlin
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
| | - Marcel Swart
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Institut de Química Computacional I Catàlisi and Department de QuímicaUniversitat de Girona17003GironaSpain
- ICREAPg. Lluís Companys 2308010BarcelonaSpain
| | - F. Matthias Bickelhaupt
- Department of Theoretical Chemistry andAmsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081 HVAmsterdamThe Netherlands
- Institute of Molecules and Materials (IMM)Radboud UniversityHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| |
Collapse
|
34
|
Hajdu B, Czakó G. Benchmark ab Initio Characterization of the Complex Potential Energy Surfaces of the X - + NH 2Y [X, Y = F, Cl, Br, I] Reactions. J Phys Chem A 2018; 122:1886-1895. [PMID: 29360360 DOI: 10.1021/acs.jpca.7b11927] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report a comprehensive high-level explicitly correlated ab initio study on the X- + NH2Y [X,Y = F, Cl, Br, I] reactions characterizing the stationary points of the SN2 (Y- + NH2X) and proton-transfer (HX + NHY-) pathways as well as the reaction enthalpies of various endothermic additional product channels such as H- + NHXY, XY- + NH2, XY + NH2-, and XHY- + NH. Benchmark structures and harmonic vibrational frequencies are obtained at the CCSD(T)-F12b/aug-cc-pVTZ(-PP) level of theory, followed by CCSD(T)-F12b/aug-cc-pVnZ(-PP) [n = Q and 5] and core correlation energy computations. In the entrance and exit channels we find two equivalent hydrogen-bonded C1 minima, X-···HH'NY and X-···H'HNY connected by a Cs first-order saddle point, X-···H2NY, as well as a halogen-bonded front-side complex, X-···YNH2. SN2 reactions can proceed via back-side attack Walden inversion and front-side attack retention pathways characterized by first-order saddle points, submerged [X-NH2-Y]- and high-energy [H2NXY]-, respectively. Product-like stationary points below the HX + NHY- asymptotes are involved in the proton-transfer processes.
Collapse
Affiliation(s)
- Bálint Hajdu
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
35
|
Liu X, Zhao C, Yang L, Zhang J, Sun R. Indirect dynamics in S N2@N: insight into the influence of central atoms. Phys Chem Chem Phys 2018; 19:22691-22699. [PMID: 28816323 DOI: 10.1039/c7cp04199c] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Central atoms have a significant influence on the reaction kinetics and dynamics of nucleophilic substitution (SN2). Herein, atomistic dynamics of a prototype SN2@N reaction F- + NH2Cl is uncovered employing direct dynamics simulations that show strikingly distinct features from those determined for a SN2@C congener F- + CH3Cl. Indirect scattering is found to prevail, which proceeds predominantly through a hydrogen-bonded F--HNHCl complex in the reactant entrance channel. This unexpected finding of a pronounced contribution of indirect reaction dynamics, even at a high collision energy, is in strong contrast to a general evolution from indirect to direct dynamics with enhanced energy that characterizes SN2@C. This result suggests that the relative importance of different atomic-level mechanisms may depend essentially on the interaction potential of reactive encounters and the coupling between inter- and intramolecular modes of the pre-reaction complex. For F- + NH2Cl the proton transfer pathway is less competitive than SN2. A remarkable finding is that the more favorable energetics for NH2Cl proton transfer, as compared to that for CH3Cl, does not manifest itself in the reaction dynamics. The present work sheds light on the underlying reaction dynamics of SN2@N, which remain largely unclear compared to well-studied SN2@C.
Collapse
Affiliation(s)
- Xu Liu
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Chenyang Zhao
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Li Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Jiaxu Zhang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | - Rui Sun
- Department of Chemistry, University of Hawaii at Manoa, Honolulu, HI 96822-2275, USA
| |
Collapse
|
36
|
Liu P, Wang D, Xu Y. A new, double-inversion mechanism of the F - + CH 3Cl S N2 reaction in aqueous solution. Phys Chem Chem Phys 2018; 18:31895-31903. [PMID: 27844085 DOI: 10.1039/c6cp06195h] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Atomic-level, bimolecular nucleophilic substitution reaction mechanisms have been studied mostly in the gas phase, but the gas-phase results cannot be expected to reliably describe condensed-phase chemistry. As a novel, double-inversion mechanism has just been found for the F- + CH3Cl SN2 reaction in the gas phase [Nat. Commun., 2015, 6, 5972], here, using multi-level quantum mechanics methods combined with the molecular mechanics method, we discovered a new, double-inversion mechanism for this reaction in aqueous solution. However, the structures of the stationary points along the reaction path show significant differences from those in the gas phase due to the strong influence of solvent and solute interactions, especially due to the hydrogen bonds formed between the solute and the solvent. More importantly, the relationship between the two double-inversion transition states is not clear in the gas phase, but, here we revealed a novel intermediate complex serving as a "connecting link" between the two transition states of the abstraction-induced inversion and the Walden-inversion mechanisms. A detailed reaction path was constructed to show the atomic-level evolution of this novel double reaction mechanism in aqueous solution. The potentials of mean force were calculated and the obtained Walden-inversion barrier height agrees well with the available experimental value.
Collapse
Affiliation(s)
- Peng Liu
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Yulong Xu
- School of Science, Qilu University of Technology, Jinan, 250353, China
| |
Collapse
|
37
|
Carrascosa E, Meyer J, Michaelsen T, Stei M, Wester R. Conservation of direct dynamics in sterically hindered S N2/E2 reactions. Chem Sci 2017; 9:693-701. [PMID: 29629138 PMCID: PMC5869569 DOI: 10.1039/c7sc04415a] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2017] [Accepted: 11/09/2017] [Indexed: 11/21/2022] Open
Abstract
Nucleophilic substitution (SN2) and base-induced elimination (E2), two indispensable reactions in organic synthesis, are commonly assumed to proceed under stereospecific conditions. Understanding the way in which the reactants pre-orient in these reactions, that is its stereodynamics, is essential in order to achieve a detailed atomistic picture and control over such processes. Using crossed beam velocity map imaging, we study the effect of steric hindrance in reactions of Cl- and CN- with increasingly methylated alkyl iodides by monitoring the product ion energy and scattering angle. For both attacking anions the rebound mechanism, indicative of a direct SN2 pathway, is found to contribute to the reaction at high relative collision energies despite being increasingly hindered. An additional forward scattering mechanism, ascribed to a direct E2 reaction, also contributes at these energies. Inspection of the product energy distributions confirms the direct and fast character of both mechanisms as opposed to an indirect reaction mechanism which leads to statistical energy redistribution in the reaction complex. This work demonstrates that nonstatistical dynamics and energetics govern SN2 and E2 pathways even in sterically hindered exchange reaction systems.
Collapse
Affiliation(s)
- Eduardo Carrascosa
- Institut für Ionenphysik und Angewandte Physik , Universität Innsbruck , Technikerstraße 25 , 6020 Innsbruck , Austria . ; Tel: +43 512 507 52620
| | - Jennifer Meyer
- Institut für Ionenphysik und Angewandte Physik , Universität Innsbruck , Technikerstraße 25 , 6020 Innsbruck , Austria . ; Tel: +43 512 507 52620
| | - Tim Michaelsen
- Institut für Ionenphysik und Angewandte Physik , Universität Innsbruck , Technikerstraße 25 , 6020 Innsbruck , Austria . ; Tel: +43 512 507 52620
| | - Martin Stei
- Institut für Ionenphysik und Angewandte Physik , Universität Innsbruck , Technikerstraße 25 , 6020 Innsbruck , Austria . ; Tel: +43 512 507 52620
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik , Universität Innsbruck , Technikerstraße 25 , 6020 Innsbruck , Austria . ; Tel: +43 512 507 52620
| |
Collapse
|
38
|
Szabó I, Czakó G. Dynamics and Novel Mechanisms of S N2 Reactions on ab Initio Analytical Potential Energy Surfaces. J Phys Chem A 2017; 121:9005-9019. [PMID: 28985079 DOI: 10.1021/acs.jpca.7b08140] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We describe a novel theoretical approach to the bimolecular nucleophilic substitution (SN2) reactions that is based on analytical potential energy surfaces (PESs) obtained by fitting a few tens of thousands high-level ab initio energy points. These PESs allow computing millions of quasi-classical trajectories thereby providing unprecedented statistical accuracy for SN2 reactions, as well as performing high-dimensional quantum dynamics computations. We developed full-dimensional ab initio PESs for the F- + CH3Y [Y = F, Cl, I] systems, which describe the direct and indirect, complex-forming Walden-inversion, the frontside attack, and the new double-inversion pathways as well as the proton-transfer channels. Reaction dynamics simulations on the new PESs revealed (a) a novel double-inversion SN2 mechanism, (b) frontside complex formation,
Collapse
Affiliation(s)
- István Szabó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
39
|
Szabó I, Czakó G. Benchmark ab Initio Characterization of the Complex Potential Energy Surface of the Cl– + CH3I Reaction. J Phys Chem A 2017; 121:5748-5757. [DOI: 10.1021/acs.jpca.7b05503] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- István Szabó
- Department of Physical Chemistry
and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry
and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
40
|
Szabó I, Olasz B, Czakó G. Deciphering Front-Side Complex Formation in S N2 Reactions via Dynamics Mapping. J Phys Chem Lett 2017; 8:2917-2923. [PMID: 28598635 DOI: 10.1021/acs.jpclett.7b01253] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Due to their importance in organic chemistry, the atomistic understanding of bimolecular nucleophilic substitution (SN2) reactions shows exponentially growing interest. In this publication, the effect of front-side complex (FSC) formation is uncovered via quasi-classical trajectory computations combined with a novel analysis method called trajectory orthogonal projection (TOP). For both F- + CH3Y [Y = Cl,I] reactions, the lifetime distributions of the F-···YCH3 front-side complex revealed weakly trapped nucleophiles (F-). However, only the F- + CH3I reaction features strongly trapped nucleophiles in the front-side region of the prereaction well. Interestingly, both back-side and front-side attack show propensity to long-lived FSC formation. Spatial distributions of the nucleophile demonstrate more prominent FSC formation in case of the F- + CH3I reaction compared to F- + CH3Cl. The presence of front-side intermediates and the broad spatial distribution in the back-side region may explain the indirect nature of the F- + CH3I reaction.
Collapse
Affiliation(s)
- István Szabó
- Department of Chemistry, King's College London , London SE1 1DB, United Kingdom
| | - Balázs Olasz
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged , Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
41
|
Abstract
The competition between bimolecular nucleophilic substitution and base-induced elimination is of fundamental importance for the synthesis of pure samples in organic chemistry. Many factors that influence this competition have been identified over the years, but the underlying atomistic dynamics have remained difficult to observe. We present product velocity distributions for a series of reactive collisions of the type X− + RY with X and Y denoting the halogen atoms fluorine, chlorine and iodine. By increasing the size of the residue R from methyl to tert-butyl in several steps, we find that the dynamics drastically change from backward to dominant forward scattering of the leaving ion relative to the reactant RY velocity. This characteristic fingerprint is also confirmed by direct dynamics simulations for ethyl as residue and attributed to the dynamics of elimination reactions. This work opens the door to a detailed atomistic understanding of transformation reactions in even larger systems. The competition between chemical reactions critically affects our natural environment and the synthesis of new materials. Here, the authors present an approach to directly image distinct fingerprints of essential organic reactions and monitor their competition as a function of steric substitution.
Collapse
|
42
|
Affiliation(s)
- Jennifer Meyer
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, 6020 Innsbruck, Austria
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, 6020 Innsbruck, Austria
| |
Collapse
|
43
|
Yang L, Liu X, Zhang J, Xie J. Effects of microsolvation on a S N2 reaction: indirect atomistic dynamics and weakened suppression of reactivity. Phys Chem Chem Phys 2017; 19:9992-9999. [PMID: 28362011 DOI: 10.1039/c7cp00294g] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Systematic studies of microsolvation in the gas phase have enriched our knowledge of solvent effects. Here, the dynamics of a prototype SN2 reaction of a hydrated fluoride ion with methyl iodide is uncovered employing direct dynamics simulations that show strikingly distinct features from those determined for an unsolvated system. An indirect scattering is found to prevail, which occurs dominantly by forming hydrated F-(H2O)-HCH2I and F-(H2O)-CH3I pre-reaction complexes at low energies, but proceeds through their water-free counterparts at higher energies. This finding is in strong contrast to a general evolution from indirect to direct dynamics with enhancing energy for the unsolvated substitution reactions, and this discrepancy is understood by the substantial steric hindrance introduced by a water molecule. As established in experiments, solvation suppresses the reactivity, whereas we find that this depression is remarkably frustrated upon raising the energy given that collision-induced dehydration essentially diminishes the water block for reactive collisions. The present study sheds light on how solute-solvent interactions affect the underlying dynamics at a deeper atomic level, thereby promoting our understanding of the fundamental solvent effects on chemical reactions in solution.
Collapse
Affiliation(s)
- Li Yang
- MIIT Key Laboratory of Critical Materials Technology for New Energy Conversion and Storage, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China.
| | | | | | | |
Collapse
|
44
|
Tajti V, Czakó G. Benchmark ab Initio Characterization of the Complex Potential Energy Surface of the F– + CH3CH2Cl Reaction. J Phys Chem A 2017; 121:2847-2854. [DOI: 10.1021/acs.jpca.7b01572] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Viktor Tajti
- Department of Physical Chemistry
and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry
and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
45
|
Li Y, Wang Y, Wang D. Quantum Dynamics Study of the Potential Energy Minima Effect on Energy Efficiency for the F– + CH3Cl → FCH3 + Cl– Reaction. J Phys Chem A 2017; 121:2773-2779. [PMID: 28346779 DOI: 10.1021/acs.jpca.7b01547] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yida Li
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong China
| | - Yuping Wang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong China
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, Jinan 250014, Shandong China
| |
Collapse
|
46
|
Olasz B, Szabó I, Czakó G. High-level ab initio potential energy surface and dynamics of the F - + CH 3I S N2 and proton-transfer reactions. Chem Sci 2017; 8:3164-3170. [PMID: 28507692 PMCID: PMC5413972 DOI: 10.1039/c7sc00033b] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2017] [Accepted: 02/15/2017] [Indexed: 11/21/2022] Open
Abstract
Bimolecular nucleophilic substitution (SN2) and proton transfer are fundamental processes in chemistry and F- + CH3I is an important prototype of these reactions. Here we develop the first full-dimensional ab initio analytical potential energy surface (PES) for the F- + CH3I system using a permutationally invariant fit of high-level composite energies obtained with the combination of the explicitly-correlated CCSD(T)-F12b method, the aug-cc-pVTZ basis, core electron correlation effects, and a relativistic effective core potential for iodine. The PES accurately describes the SN2 channel producing I- + CH3F via Walden-inversion, front-side attack, and double-inversion pathways as well as the proton-transfer channel leading to HF + CH2I-. The relative energies of the stationary points on the PES agree well with the new explicitly-correlated all-electron CCSD(T)-F12b/QZ-quality benchmark values. Quasiclassical trajectory computations on the PES show that the proton transfer becomes significant at high collision energies and double-inversion as well as front-side attack trajectories can occur. The computed broad angular distributions and hot internal energy distributions indicate the dominance of indirect mechanisms at lower collision energies, which is confirmed by analyzing the integration time and leaving group velocity distributions. Comparison with available crossed-beam experiments shows usually good agreement.
Collapse
Affiliation(s)
- Balázs Olasz
- Department of Physical Chemistry and Materials Science , Institute of Chemistry , University of Szeged , Rerrich Béla tér 1 , Szeged H-6720 , Hungary .
| | - István Szabó
- Department of Physical Chemistry and Materials Science , Institute of Chemistry , University of Szeged , Rerrich Béla tér 1 , Szeged H-6720 , Hungary .
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science , Institute of Chemistry , University of Szeged , Rerrich Béla tér 1 , Szeged H-6720 , Hungary .
| |
Collapse
|
47
|
Piccini G, McCarty JJ, Valsson O, Parrinello M. Variational Flooding Study of a S N2 Reaction. J Phys Chem Lett 2017; 8:580-583. [PMID: 28071915 DOI: 10.1021/acs.jpclett.6b02852] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We have studied the reaction dynamics of a prototypical organic reaction using a variationally optimized truncated bias to accelerate transitions between educt and product reactant states. The asymmetric SN2 nucleophilic substitution reaction of fluoromethane and chloromethane CH3F + Cl- ⇌ CH3Cl + F- is considered, and many independent biased molecular dynamics simulations have been performed at 600, 900, and 1200 K, collecting several hundred transitions at each temperature. The transition times and relative rate constants have been obtained for both reaction directions. The activation energies extracted from an Arrhenius plot compare well with standard static calculations.
Collapse
Affiliation(s)
- GiovanniMaria Piccini
- Department of Chemistry and Applied Biosciences, ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltá di Informatica, Instituto di Scienze Computationali, Universitá della Svizzera italiana (USI) , Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| | - James J McCarty
- Department of Chemistry and Applied Biosciences, ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltá di Informatica, Instituto di Scienze Computationali, Universitá della Svizzera italiana (USI) , Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| | - Omar Valsson
- Department of Chemistry and Applied Biosciences, ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltá di Informatica, Instituto di Scienze Computationali, Universitá della Svizzera italiana (USI) , Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich , c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
- Facoltá di Informatica, Instituto di Scienze Computationali, Universitá della Svizzera italiana (USI) , Via Giuseppe Buffi 13, CH-6900 Lugano, Ticino, Switzerland
| |
Collapse
|
48
|
Carrascosa E, Meyer J, Wester R. Imaging the dynamics of ion–molecule reactions. Chem Soc Rev 2017; 46:7498-7516. [DOI: 10.1039/c7cs00623c] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A range of ion–molecule reactions have been studied in the last years using the crossed-beam ion imaging technique, from charge transfer and proton transfer to nucleophilic substitution and elimination.
Collapse
Affiliation(s)
- Eduardo Carrascosa
- Institut für Ionenphysik und Angewandte Physik
- Universität Innsbruck
- 6020 Innsbruck
- Austria
| | - Jennifer Meyer
- Institut für Ionenphysik und Angewandte Physik
- Universität Innsbruck
- 6020 Innsbruck
- Austria
| | - Roland Wester
- Institut für Ionenphysik und Angewandte Physik
- Universität Innsbruck
- 6020 Innsbruck
- Austria
| |
Collapse
|
49
|
Szabó I, Czakó G. Mode-specific multi-channel dynamics of the F - + CHD 2Cl reaction on a global ab initio potential energy surface. J Chem Phys 2016; 145:134303. [PMID: 27782409 DOI: 10.1063/1.4963664] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
We report a detailed quasiclassical trajectory study for the dynamics of the ground-state and CH/CD stretching-excited F- + CHD2Cl(vCH/CD = 0, 1) → Cl- + CHD2F, HF + CD2Cl-, and DF + CHDCl- SN2, proton-, and deuteron-abstraction reactions using a full-dimensional global ab initio analytical potential energy surface. The simulations show that (a) CHD2Cl(vCH/CD = 1), especially for vCH = 1, maintains its mode-specific excited character prior to interaction, (b) the SN2 reaction is vibrationally mode-specific,
Collapse
Affiliation(s)
- István Szabó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| | - Gábor Czakó
- Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Rerrich Béla tér 1, Szeged H-6720, Hungary
| |
Collapse
|
50
|
Proenza YG, de Souza MAF, Longo RL. Dynamical Bifurcation in Gas-Phase XH - + CH 3 Y S N 2 Reactions: The Role of Energy Flow and Redistribution in Avoiding the Minimum Energy Path. Chemistry 2016; 22:16220-16229. [PMID: 27651104 DOI: 10.1002/chem.201602976] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2016] [Indexed: 11/10/2022]
Abstract
The gas-phase reactions of XH- (X=O, S) + CH3 Y (Y=F, Cl, Br) span nearly the whole range of SN 2 pathways, and show an intrinsic reaction coordinate (IRC) (minimum energy path) with a deep well owing to the CH3 XH⋅⋅⋅Y- (or CH3 S- ⋅⋅⋅HF) hydrogen-bonded postreaction complex. MP2 quasiclassical-type direct dynamics starting at the [HX⋅⋅⋅CH3 ⋅⋅⋅Y]- transition-state (TS) structure reveal distinct mechanistic behaviors. Trajectories that yield the separated CH3 XH+Y- (or CH3 S- +HF) products directly are non-IRC, whereas those that sample the CH3 XH⋅⋅⋅Y- (or CH3 S- ⋅⋅⋅HF) complex are IRC. The IRCIRC/non-IRC ratios of 90:10, 40:60, 25:75, 2:98, 0:100, and 0:100 are obtained for (X, Y)=(S, F), (O, F), (S, Cl), (S, Br), (O, Cl), and (O, Br), respectively. The properties of the energy profiles after the TS cannot provide a rationalization of these results. Analysis of the energy flow in dynamics shows that the trajectories cross a dynamical bifurcation, and that the inability to follow the minimum energy path arises from long vibration periods of the X-C⋅⋅⋅Y bending mode. The partition of the available energy to the products into vibrational, rotational, and translational energies reveals that if the vibrational contribution is more than 80 %, non-IRC behavior dominates, unless the relative fraction of the rotational and translational components is similar, in which case a richer dynamical mechanism is shown, with an IRC/non-IRC ratio that correlates to this relative fraction.
Collapse
Affiliation(s)
- Yaicel G Proenza
- Departamento de Química Fundamental, CCEN, Universidade Federal de Pernambuco, 50.740-560, Recife, PE, Brazil
| | - Miguel A F de Souza
- Instituto de Química, CCET, Universidade Federal do Rio Grande do Norte, 59.072-970, Natal, RN, Brazil
| | - Ricardo L Longo
- Departamento de Química Fundamental, CCEN, Universidade Federal de Pernambuco, 50.740-560, Recife, PE, Brazil.
| |
Collapse
|