Exploring The Sequence of Electron Density Along The Chemical Reactions Between Carbonyl Oxides And Ammonia/Water Using Bond Evolution Theory

Full-dimensional automated potential energy surface development and detailed dynamics for the CH2OO + NH3 reaction

With the help of the ROBOSURFER program package, a global full-dimensional potential energy surface (PES) for the reaction of the Criegee intermediate, CH2OO, with the NH3 molecule is developed iteratively using different ab initio methods and the monomial symmetrization fitting approach. The final permutationally-invariant analytical PES is constructed based on 23447 geometries and the corresponding ManyHF-based CCSD(T)-F12b/cc-pVTZ-F12 energies. The accuracy of the PES is confirmed by the excellent agreement of its stationary-point properties and one-dimensional potential energy curves compared with the corresponding ab initio data. The reaction probabilities and integral cross sections are calculated for the ground-state and several vibrationally excited-state reactions by quasi-classical trajectory simulations. Remarkable is that the maximum impact parameter b where reactivity vanishes is almost independent of collision energy ranging from 1 to 40 kcal mol-1, and the reaction probability increases with increasing collision energy for this negative-barrier reaction. At the same time, a slight mode-specificity effect is observed. In addition, the deuterium effect is investigated and the sudden vector projection is discussed.

Exploring the Mechanism of the Intramolecular Diels-Alder Reaction of (2E,4Z,6Z)-2(allyloxy)cycloocta-2,4,6-trien-1-one Using Bonding Evolution Theory

In the present work, the bond breaking/forming events along the intramolecular Diels-Alder (IMDA) reaction of (2E,4Z,6Z)-2(allyloxy)cycloocta-2,4,6-trien-1-one have been revealed within bonding evolution theory (BET) at the density functional theory level, using the M05-2X functional with the cc-pVTZ basis set. Prior to achieving this task, the energy profiles and stationary points at the potential energy surface (PES) have been characterized. The analysis of the results finds that this rearrangement can proceed along three alternative reaction pathways (a-c). Paths a and b involve two steps, while path c is a one-step process. The first step in path b is kinetically favored, and leads to the formation of an intermediate step, Int-b. Further evolution from Int-b leads mainly to 3-b1. However, 2 is the thermodynamically preferred product and is obtained at high temperatures, in agreement with the experimental observations. Regarding the BET analysis along path b, the breaking/forming process is described by four structural stability domains (SSDs) during the first step, which can be summarized as follows: (1) the breaking of the C-O bond with the transfer of its population to the lone pair (V(O)), (2) the reorganization of the electron density with the creation of two V(C) basins, and (3) the formation of a new C-C single bond via the merger of the two previous V(C) basins. Finally, the conversion of Int-b (via TS2-b1) occurs via the reorganization of the electron density during the first stage (the creation of different pseudoradical centers on the carbon atoms as a result of the depopulation of the C-C double bond involved in the formation of new single bonds), while the last stage corresponds to the non-concerted formation of the two new C-C bonds via the disappearance of the population of the four pseudoradical centers formed in the previous stage. On the other hand, along path a, the first step displays three SSDs, associated with the depopulation of the V(C2,C3) and V(C6,C7) basins, the appearance of the new monosynaptic basins V(C2) and V(C7), and finally the merging of these new monosynaptic basins through the creation of the C2-C7 single bond. The second step is described by a series of five SSDs, that account for the reorganization of the electron density within Int-a via the creation of four pseudoradical centers on the C12, C13, C3 and C6 carbon atoms. The last two SSDs deal with the formation of two C-C bonds via the merging of the monosynaptic basins formed in the previous domains.

Decoding the reaction mechanism of the cyclocondensation of ethyl acetate2-oxo-2-(4-oxo-4H-pyrido [1.2-a] pyrimidin-3-yl) polyazaheterocycle and ethylenediamine using bond evolution theory

We investigated the flow of electron density along the cyclocondensation reaction between ethyl acetate 2-oxo-2-(4-oxo-4H-pyrido[1.2-a]pyrimidin-3-yl) polyazaheterocycle (1) and ethylenediamine (2) at the ωB97XD/6-311++G(d,p)computational method within of bond evolution theory (BET). The exploration of potential energy surface shows that this reaction has three channels (1-3) with the formation of product 3 via channel-2 (the most favorable one) as the main product and this is in good agreement with experimental observations. The BET analysis allows identifying unambiguously the main chemical events happening along channel-2. The mechanism along first step (TS2-a) is described by a series of four structural stability domains (SSDs), while five SSDs for the last two steps (TS2-b and TS2-c). The first and third steps can be summarized as follows, the formation of N1-C6 bond (SSD-II), then, the restoration of the nitrogen N1 lone pair (SSD-III), and finally, the formation of the last O1-H1 bond (SSD-IV). For the second step, the formation of hydroxide ion is noted, as a result of the disappearance of V(C6,O7) basin and the transformation of C6-N1 single bond into double one (SSD-IV). Finally, the appearance of V(O7,H2) basin lead to the elimination of water molecule within the last domain is observed.

Unveiling the [3+2] cycloaddition between difluoromethyl diazomethane and 3-ylideneoxindole from the perspective of molecular electron density theory

Evolution of some key ELF basins along the IRC of the most favorable ortho/endo reaction path.

Topological unraveling of the [3+2] cycloaddition (32CA) reaction between N-methylphenylnitrone and styrene catalyzed by the chromium tricarbonyl complex using electron localization function and catastrophe theory

We have investigated the reaction mechanisms of [3+2] cycloaddition (32CA) between N-methylphenylnitrone and styrene catalyzed by the chromium tricarbonyl complex at the MPWB1K/6-311G(d,p) level of approximation.