scholarly journals Isomerization and Fragmentation Reactions on the [C2SH4] Potential Energy Surface. The Metastable Thione-S-Methylide Isomer

Author(s):  
Zoi Salta ◽  
Marc E. Segovia ◽  
Aline Katz ◽  
Nicola Tasinato ◽  
Vincenzo Barone ◽  
...  

Thione S-methylide (TSM), the parent species of the thiocarbonyl ylide family, is a 1,3-dipolar, planar species on the [C2SH4] potential energy surface (PES), which has not shared the richness of studies dedicated to its isomers, the cyclic thiirane (THI), and the keto-enol pair vinyl thiol (VTH)/thioacetaldehyde (THA). While the conrotatory ring closure reaction toward THI was studied in the ‘90s, no complete analysis of the PES is available in the literature. In the present paper, we report a computational study of the reaction scheme linking all species on that PES. We employ several levels of calculation, ranging from density functional theory (DFT), through CCSD(T) based composite schemes, to CASSCF/CASPT2 multi-reference procedures, to find the best description of TSM, its isomers, and the transition states (TSs) ruling their interconversion. Fragmentation of TSM, THA and THI were investigated and compared to the available experimental information. We found that the B2PLYP-D3 functional, contrary to M06-2XD3 or B97X-D, describes well the geometry of both TSM and the transition state connecting it to THI. The reverse barrier, from THI to TSM, amounts to 52.2 kcal mol-1 (to be compared to 17.6 kcal mol-1 for the direct one), thus explaining why, in general, thiocarbonyl ylides cannot be prepared from thiiranes. Conversion of THI to VTH implies also a large barrier, explaining why the reaction has been observed only at high temperatures. The fragmentation of THI to S(3P) or S(1D) and ethylene was also explored, together with the decomposition to H2S plus acetylene. Open species, both in triplet and singlet states, were identified as intermediates in the fragmentations, and their energies were found to be lower than the transition state for the isomerization of THI to VTH, thus explaining the preference for fragmentation over isomerization at relatively low temperatures.

2020 ◽  
Author(s):  
Zoi Salta ◽  
Marc E. Segovia ◽  
Aline Katz ◽  
Nicola Tasinato ◽  
Vincenzo Barone ◽  
...  

Thione S-methylide (TSM), the parent species of the thiocarbonyl ylide family, is a 1,3-dipolar, planar species on the [C2SH4] potential energy surface (PES), which has not shared the richness of studies dedicated to its isomers, the cyclic thiirane (THI), and the keto-enol pair vinyl thiol (VTH)/thioacetaldehyde (THA). While the conrotatory ring closure reaction toward THI was studied in the ‘90s, no complete analysis of the PES is available in the literature. In the present paper, we report a computational study of the reaction scheme linking all species on that PES. We employ several levels of calculation, ranging from density functional theory (DFT), through CCSD(T) based composite schemes, to CASSCF/CASPT2 multi-reference procedures, to find the best description of TSM, its isomers, and the transition states (TSs) ruling their interconversion. Fragmentation of TSM, THA and THI were investigated and compared to the available experimental information. We found that the B2PLYP-D3 functional, contrary to M06-2XD3 or B97X-D, describes well the geometry of both TSM and the transition state connecting it to THI. The reverse barrier, from THI to TSM, amounts to 52.2 kcal mol-1 (to be compared to 17.6 kcal mol-1 for the direct one), thus explaining why, in general, thiocarbonyl ylides cannot be prepared from thiiranes. Conversion of THI to VTH implies also a large barrier, explaining why the reaction has been observed only at high temperatures. The fragmentation of THI to S(3P) or S(1D) and ethylene was also explored, together with the decomposition to H2S plus acetylene. Open species, both in triplet and singlet states, were identified as intermediates in the fragmentations, and their energies were found to be lower than the transition state for the isomerization of THI to VTH, thus explaining the preference for fragmentation over isomerization at relatively low temperatures.


2017 ◽  
Vol 16 (05) ◽  
pp. 1750046 ◽  
Author(s):  
Vijay M. Miriyala ◽  
Priya Bhasi ◽  
Zanele P. Nhlabatsi ◽  
Sanyasi Sitha

Using computational calculations, we have revisited the potential energy surface (PES) of the reaction between OH and SO2, which is believed as the rate-limiting step in the atmospheric formation of H2SO4. In this work, we report for the first time the presence of a pre-reaction hydrogen-bonded complex between OH and SO2 in the reaction PES. Based on this finding, it has been shown that the reaction can be considered as a two-step process in which the first step is the formation of the pre-reaction complex and the second step is the transformation of this complex to the product. It was observed that due to the presence of this pre-reaction complex as a potential well in the reaction PES, the barrier height got increased by around two-fold for the second step. Based on this observation, it has been proposed that the kinetics of the reaction is going to be affected. Also based on the analysis of the geometries of this pre-reaction complex and the transition state, it has been argued that the step involving the transformation of this pre-reaction complex to the product via the transition state is going to be the slowest step as this transformation involves large structural changes of the stationary points involved.


2020 ◽  
Author(s):  
Shi Jun Ang ◽  
Wujie Wang ◽  
Daniel Schwalbe-Koda ◽  
Simon Axelrod ◽  
Rafael Gomez-Bombarelli

<div>Modeling dynamical effects in chemical reactions, such as post-transition state bifurcation, requires <i>ab initio</i> molecular dynamics simulations due to the breakdown of simpler static models like transition state theory. However, these simulations tend to be restricted to lower-accuracy electronic structure methods and scarce sampling because of their high computational cost. Here, we report the use of statistical learning to accelerate reactive molecular dynamics simulations by combining high-throughput ab initio calculations, graph-convolution interatomic potentials and active learning. This pipeline was demonstrated on an ambimodal trispericyclic reaction involving 8,8-dicyanoheptafulvene and 6,6-dimethylfulvene. With a dataset size of approximately</div><div>31,000 M062X/def2-SVP quantum mechanical calculations, the computational cost of exploring the reactive potential energy surface was reduced by an order of magnitude. Thousands of virtually costless picosecond-long reactive trajectories suggest that post-transition state bifurcation plays a minor role for the reaction in vacuum. Furthermore, a transfer-learning strategy effectively upgraded the potential energy surface to higher</div><div>levels of theory ((SMD-)M06-2X/def2-TZVPD in vacuum and three other solvents, as well as the more accurate DLPNO-DSD-PBEP86 D3BJ/def2-TZVPD) using about 10% additional calculations for each surface. Since the larger basis set and the dynamic correlation capture intramolecular non-covalent interactions more accurately, they uncover longer lifetimes for the charge-separated intermediate on the more accurate potential energy surfaces. The character of the intermediate switches from entropic to thermodynamic upon including implicit solvation effects, with lifetimes increasing with solvent polarity. Analysis of 2,000 reactive trajectories on the chloroform PES shows a qualitative agreement with the experimentally-reported periselectivity for this reaction. This overall approach is broadly applicable and opens a door to the study of dynamical effects in larger, previously-intractable reactive systems.</div>


2007 ◽  
Vol 06 (03) ◽  
pp. 549-562
Author(s):  
ABRAHAM F. JALBOUT

The transition states for the H 2 NO decomposition and rearrangements mechanisms have been explored by the CBS-Q method or by density functional theory. Six transition states were located on the potential energy surface, which were explored with the Quadratic Complete Basis Set (CBS-Q) and Becke's one-parameter density functional hybrid methods. Interesting deviations between the CBS-Q results and the B1LYP density functional theory lead us to believe that further study into this system is necessary. In the efforts to further assess the stabilities of the transition states, bond order calculations were performed to measure the strength of the bonds in the transition state.


2017 ◽  
Vol 70 (4) ◽  
pp. 362 ◽  
Author(s):  
Young J. Hong ◽  
Dean J. Tantillo

An energetically viable (on the basis of results from density functional theory computations) pathway to the diterpene variediene is described. Only one of the three secondary carbocations along this pathway is predicted to be a minimum on the potential energy surface.


2019 ◽  
Vol 21 (21) ◽  
pp. 11395-11404
Author(s):  
Abel Carreras ◽  
Luca Fuligni ◽  
Pere Alemany ◽  
Miquel Llunell ◽  
Josep Maria Bofill ◽  
...  

We present a computational study of a reduced potential energy surface (PES) to describe enantiomerization and internal rotation in three triptycyl-n-helicene molecules, centering the discussion on the issue of a proper reaction coordinate choice.


Sign in / Sign up

Export Citation Format

Share Document