Dataset Information

Differences in the Abilities to Mechanically Eliminate Activation Energies for Unimolecular and Bimolecular Reactions.

ABSTRACT: Mechanochemistry, i.e. the application of forces, F, at the molecular level, has attracted significant interest as a means of controlling chemical reactions. The present study uses quantum chemical calculations to explore the abilities to mechanically eliminate activation energies, ?E(‡), for unimolecular and bimolecular reactions. The results demonstrate that ?E(‡) can be eliminated for unimolecular reactions by applying sufficiently large F along directions that move the reactant and/or transition state (TS) structures parallel to the zero-F reaction coordinate, S0. In contrast, eliminating ?E(‡) for bimolecular reactions requires the reactant to undergo a force-induced shift parallel to S0 irrespective of changes in the TS. Meeting this requirement depends upon the coupling between F and S0 in the reactant. The insights regarding the differences in eliminating ?E(‡) for unimolecular and bimolecular reactions, and the requirements for eliminating ?E(‡), may be useful in practical efforts to control reactions mechanochemically.

SUBMITTER: Kochhar GS

PROVIDER: S-EPMC4789786 | biostudies-literature | 2016 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Differences in the Abilities to Mechanically Eliminate Activation Energies for Unimolecular and Bimolecular Reactions.

Kochhar Gurpaul S GS Mosey Nicholas J NJ

Scientific reports 20160314

Mechanochemistry, i.e. the application of forces, F, at the molecular level, has attracted significant interest as a means of controlling chemical reactions. The present study uses quantum chemical calculations to explore the abilities to mechanically eliminate activation energies, ΔE(‡), for unimolecular and bimolecular reactions. The results demonstrate that ΔE(‡) can be eliminated for unimolecular reactions by applying sufficiently large F along directions that move the reactant and/or transi ...[more]

PMID: 26972114

Similar Datasets

Project description:A detailed computational study of the atmospheric reaction of the simplest Criegee intermediate CH2OO with methane has been performed using the density functional theory (DFT) method and high-level calculations. Solvation models were utilized to address the effect of water molecules on prominent reaction steps and their associated energies. The structures of all proposed mechanisms were optimized using B3LYP functional with several basis sets: 6-31G(d), 6-31G (2df,p), 6-311++G(3df,3pd) and at M06-2X/6-31G(d) and APFD/6-31G(d) levels of theory. Furthermore, all structures were optimized at the B3LYP/6-311++G(3df,3pd) level of theory. The intrinsic reaction coordinate (IRC) analysis was performed for characterizing the transition states on the potential energy surfaces. Fifteen different mechanistic pathways were studied for the reaction of Criegee intermediate with methane. Both thermodynamic functions (ΔH and ΔG), and activation parameters (activation energies Ea, enthalpies of activation ΔHǂ, and Gibbs energies of activation ΔGǂ) were calculated for all pathways investigated. The individual mechanisms for pathways A1, A2, B1, and B2, comprise two key steps: (i) the formation of ethyl hydroperoxide (EHP) accompanying with the hydrogen transfer from the alkanes to the terminal oxygen atom of CIs, and (ii) a following unimolecular dissociation of EHP. Pathways from C1 → H1 involve the bimolecular reaction of EHP with different atmospheric species. The photochemical reaction of methane with EHP (pathway E1) was found to be the most plausible reaction mechanism, exhibiting an overall activation energy of 7 kJ mol-1, which was estimated in vacuum at the B3LYP/6-311++G(3df,3pd) level of theory. All of the reactions were found to be strongly exothermic, expect the case of the sulfur dioxide-involved pathway that is predicted to be endothermic. The solvent effect plays an important role in the reaction of EHP with ammonia (pathway F1). Compared with the gas phase reaction, the overall activation energy for the solution phase reaction is decreased by 162 and 140 kJ mol-1 according to calculations done with the SMD and PCM solvation models, respectively.

Project description:We present here a greatly updated version of an earlier study on the conformational energies of protein-ligand complexes in the Protein Data Bank (PDB) [Nicklaus et al. Bioorg. Med. Chem. 1995, 3, 411-428], with the goal of improving on all possible aspects such as number and selection of ligand instances, energy calculations performed, and additional analyses conducted. Starting from about 357,000 ligand instances deposited in the 2008 version of the Ligand Expo database of the experimental 3D coordinates of all small-molecule instances in the PDB, we created a "high-quality" subset of ligand instances by various filtering steps including application of crystallographic quality criteria and structural unambiguousness. Submission of 640 Gaussian 03 jobs yielded a set of about 415 successfully concluded runs. We used a stepwise optimization of internal degrees of freedom at the DFT level of theory with the B3LYP/6-31G(d) basis set and a single-point energy calculation at B3LYP/6-311++G(3df,2p) after each round of (partial) optimization to separate energy changes due to bond length stretches vs bond angle changes vs torsion changes. Even for the most "conservative" choice of all the possible conformational energies-the energy difference between the conformation in which all internal degrees of freedom except torsions have been optimized and the fully optimized conformer-significant energy values were found. The range of 0 to ~25 kcal/mol was populated quite evenly and independently of the crystallographic resolution. A smaller number of "outliers" of yet higher energies were seen only at resolutions above 1.3 Å. The energies showed some correlation with molecular size and flexibility but not with crystallographic quality metrics such as the Cruickshank diffraction-component precision index (DPI) and R(free)-R, or with the ligand instance-specific metrics such as occupancy-weighted B-factor (OWAB), real-space R factor (RSR), and real-space correlation coefficient (RSCC). We repeated these calculations with the solvent model IEFPCM, which yielded energy differences that were generally somewhat lower than the corresponding vacuum results but did not produce a qualitatively different picture. Torsional sampling around the crystal conformation at the molecular mechanics level using the MMFF94s force field typically led to an increase in energy.

Project description:Many of the radical-molecule reactions are nonelementary reactions with negative activation energies, which usually proceed through two steps. They exist extensively in the atmospheric chemistry and hydrocarbon fuel combustion, so they are extensively studied both theoretically and experimentally. At the same time, various models, such as a two transition state model, a steady-state model, an equilibrium-state model, and a direct elementary dynamics model are proposed to get the kinetic parameters for the overall reaction. In this paper, a conversion temperature T C1 is defined as the temperature at which the standard molar Gibbs free energy change of the formation of the reaction complex is equal to zero, and it is found that when T ? T C1, the direct elementary dynamics model with an inclusion of the tunneling correction of the second step reaction is applicable to calculate the overall reaction rate constants for this kind of reaction system. The reaction class of hydroxyl radical addition to alkenes is chosen as the objects of this study, five reactions are chosen as the representative for the reaction class, and their single-point energies are calculated using the method of CCSD(T)/CBS, and it is shown that the highest conversion temperature for the five reactions is 139.89 K, far below the usual initial low-temperature (550 K) oxidation chemistry of hydrocarbon fuels; therefore, the steady-state approximation method is applicable. All geometry optimizations are performed at the BH&HLYP/6-311+G(d,p) level, and the result shows that the geometric parameters in the reaction centers are conserved; hence, the isodesmic reaction method is applicable to this reaction class. To validate the accuracy of this scheme, a comparison of electronic energy difference at the BH&HLYP/6-311+G(d,p) level and the corrected electronic energy difference with the electronic energy difference at the CCSD(T)/CBS level is performed for the five representative reactions, and it is shown that the maximum absolute deviation of electronic energy difference can be reduced from 2.54 kcal·mol-1 before correction to 0.58 kcal·mol-1 after correction, indicating that the isodesmic reaction method is applicable for the accurate calculation of the kinetic parameters for large-size molecular systems with a negative activation energy reaction. The overall rate constants for 44 reactions of the reaction class of hydroxyl radical addition to alkenes are calculated using the transition-state theory in combination with the isodesmic correction scheme, and high-pressure limit rate rules for the reaction class are developed. In addition, the thermodynamic parameter is calculated and the results indicate that our dynamics model is applicable for our studied reaction class. A chemical kinetic modeling and sensitivity analysis using the calculated kinetic data is performed for the combustion of ethene, and the results indicate the studied reaction is important for the low-to-medium temperature combustion modeling of ethene.

Dataset Information

Differences in the Abilities to Mechanically Eliminate Activation Energies for Unimolecular and Bimolecular Reactions.

Publications

Differences in the Abilities to Mechanically Eliminate Activation Energies for Unimolecular and Bimolecular Reactions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure