Criegee中间体气相反应热力学的G2理论计算

Criegee intermediate is believed to play an important role in the atmospheric chemistry. Because of its short life and the difficulty in experimental study, we carried out ah initio calculations on the thermochemistry of the Criegee involving reactions in this study. Thermochemistry data of reaction enthalpies and Gibbs free energies for four different stable structures of the Criegee intermediates (singlet CH2OO ①1 A1 in C2v, triplet CH2OO ②3B1 in C2v, singlet CH2OO ③1A' in Cs and triplet CH2OO ④ in C1 symmetry) involved in some of the gas-phase reactions were calculated at the standard Gaussian-2 [G2(MP2) and G2] and a modified G2, G2(fu1)[10],levels of theory. Relative energies among those Criegees and formic acid were compared. Chemical reactions include the formation of Criegees, re-arrangement from Criegee to formic acid, dissociations (producing CH2(3B1)+O2, CH2(1A1)+O2, CO2+H2, CO2+2H, CO+H2O, OH+HCO) and the reactions between Criegee and NO/H2O. Standard equilibrium constants for some reactions were investigated and may be obtained for all of the rest reactions involved in this study by the standard Gibbs free energies. It is shown that the formation of Criegee ①-④ by ethylene and ozone, the re-arrangement from any Criegee to formic acid, the dissociation in producing CO2+O2and CO+H2O and the reactions between any Criegee and NO/H2O are all favourable thermodynamicaly. The dissociation in forming CO2+2H and OH+HCO is less favourable. While the dissociation in forming carbene (either in 3B1 or 1A1 state) is not allowed by ΔrGm? values. Standard enthalpies of formation at 298 K for the four Criegees were predicted at the G2(ful) level of theory. Each value is the average value from ten of the above reactions and they are -4.3, 74.8,98.9 and 244.6 kJ mol-1 at the G2(ful) level for Criegee ① to Criegee ④, respectively. In addition, tile standard enthalpy of formation at 298 K for HOCH2OOH is further predicted to be -315.6 kJ mol-1 at the G2(MP2) level.     

Criegee中间体气相反应热力学的G2理论计算

Criegee中间体已被认为是大气中一些重要污染物的来源，是引发大量次级反应的关键物质，与光化学氧化剂、有机过氧化物的产生，以及与酸雨、大气光化学烟雾和森林衰退等许多环境问题有关[1-4]．但由于Criegee中间体寿命短，对其可能参与的一些大气化学反应进行实验研究有一定难     

量子化学的Gaussian－2理论及其应用与化合物的标准生成焓预测

本文较为全面地综述了Ｇａｕｓｓｉａｎ－１，Ｇａｕｓｓｉａｎ－２（简称Ｇ１，Ｇ２）理论以及简化的Ｇ２（ＭＰ２），Ｇ２（ＭＰ３）理论，将其主要结果进行了比较分析。关于Ｇ２理论的应用，除了较为详细地综述了几年来理论在重现实验数据、评价实验数据、预测实验数据及研究化学反应途径等方面的应用外，还结合我们近期研究结果的主要结论讨论了该理论在研究等电子－等自旋，价层等电－等旋，等旋及非等旋化学反应的能量计算中的应用情况，以及该理论在预测化合物的标准生成焓方面的应用情况。     

The Energetics of Large Deformations of a Single Polyimide Molecular Chain: DFT and MO Calculations

The large‐deformation energetics of a single molecular chain of the rod‐like polyimide PMDA‐PDA was investigated using DFT, ab initio MO and semi‐empirical MO methods. The force/displacement curves were calculated from tensile testing simulations along the axis of the molecular chain, allowing a discussion of the distribution and change of local strain of the molecular chain. The deformation behavior of a single PMDA‐PDA molecular chain under finite deformations as functions of bending angle and dihedral angle between PMDA and PDA groups are compared. It is found that the semi‐empirical MO calculations provide sufficient accuracy to express the energetics of large deformations except for compressive deformation.

     

Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data

Reactions involving Criegee intermediates (CIs, R₁R₂COO) are important in atmospheric ozonolysis models. In recent years, density functional theory (DFT) and CCSD(T)-based ab initio methods are increasingly being used for modeling reaction profiles involving CIs. We obtain highly accurate CCSDT(Q)/CBS reaction energies and barrier heights for ring-closing reactions involving atmospherically important CIs (R₁/R₂ = H, Me, OH, OMe, F, CN, cyclopropene, ethylene, acetaldehyde, and acrolein). We use this benchmark data to evaluate the performance of DFT, double-hybrid DFT (DHDFT), and ab initio methods for the kinetics and thermodynamics of these reactions. We find that reaction energies are more challenging for approximate theoretical procedures than barrier heights. Overall, taking both reaction energies and barrier heights into account, only one of the 58 considered DFT methods (the meta-GGA MN12-L) attains near chemical accuracy, with root-mean-square deviations (RMSDs) of 3.5 (barrier heights) and 4.7 (reaction energies) kJ mol⁻¹. Therefore, MN12-L is recommended for investigations where CCSD(T)-based methods are not computationally feasible. For reaction barrier heights performance does not strictly follow Jacob's Ladder, for example, DHDFT methods do not perform better than conventional DFT methods. Of the ab initio methods, the cost-effective CCSD(T)/CBS(MP2) approach gives the best performance for both reaction energies and barrier heights, with RMSDs of 1.7 and 1.4 kJ mol⁻¹, respectively. All the considered Gaussian-n methods show good performance with RMSDs below the threshold of chemical accuracy for both reaction energies and barrier heights, where G4(MP2) shows the best overall performance with RMSDs of 2.9 and 1.5 kJ mol⁻¹, respectively. © 2019 Wiley Periodicals, Inc.     

多面体碳烷的分子轨道方法研究

利用INDO自洽场方法和Edmiston-Ruedenberg定域化方法,计算了多面体碳烷C_2nH_2n(n=2,3,4,5和10)及其骨架C_2n,讨论了它们的电子结构、稳定性和化学键性质。     

含芯电子相关能修正的G2理论——G2(ful)

A general method in considering the core electronic correlation energies has been proposed and introduced into the standard Gaussian-2 (G2)[7] theory by small post-Hartree-Fock calculations. In this paper an additional MP2(FC)/6-31G(d) calculation over the G2 procedures is employed and examined in modification in modification to the flaw of Frozen-Core (FC) approximation of G2 vai eq.:
ΔE(full)= E[MP2(full)/6-31G(d)]-E[MP2(FC)/6-31G(d)]
where the MP2(full)/6-31G(d) energy has been obtained in the molecular geometry optimizations. This energy, ΔE(full), is directly added into the total G2 energy of a molecule in facilitating the effect of core electronic correlations for each molecule in chemical reactions. It has been shown that the over-all average absolute deviation for the 125 reaction energies of the G2 test set (test set 1) is slightly reduced from 5.09 to 5.01 kJ, mol(-1) while for the 55 D0 values, which have been used for the derivation of the A coefficient of the empirical High-Level...更多-Correction (HLC), it is also reduced from 4.99 [for both G2 and G2(COMPLETE)[8]]to 4.77 kJ• mol(-1). In addition, larger errors (greater than ±8.4 kJ•mol(-1) for the D0 energies are improved, especially for the largest error of the D0 of SO2 This error is reduced from 21.3 to 15.4 kJ. mol(-1), in which the experimental geometry would further reduce it by 7.1kJ.mol(-1)[8]. Another improvement is the absolute value of the A coefficient in HLC being reduced from 4.81 for G2 to 4.34 milli-hartrees which is believed to be useful in isolating the relationship between the HLC and the FC approximation. Modifications to the original G2 from this work is denoted as G2(fu 1) and thus the G2 (fu 1) total energy for a molecule is
E[G2(fu 1)]= E[G2]+Δ E(full)h
with a new ΔE[HLC] =-0.19α- 4.34nβ milli-hartree.     

An Ab Initio Molecular Orbital Theory and Density Functional Theory (DFT) Study of Conformers and Rotamers of 4-Substituted (Methyl, Hydroxymethyl, Methanoyl, Ethanoyl, Cyano, Fluoro, Chloro, Bromo, Acetoxy) Tetrahydro-2H-thiopyran-1,1-dioxides

The structures and energies of axial and equatorial conformers and rotamers of 4-substituted tetrahydro-2H-thiopyran-1,1-dioxides (tetrahydrothiopyran-1,1-dioxides, thiacyclohexane-1,1-dioxides, thiane-1,1-dioxides, and 1,1-dioxothianes; CH₃, CH₂OH, CHO, COCH₃, CN, F, Cl, Br, and OCOCH₃) were calculated using the hybrid density functionals B3LYP, B3P86, and B3PW91, as well as MP2 and the 6-31G(d), 6-31G(2d), 6-31G(3d), 6-31G(d,p), and 6-31+G(d) basis sets. MP2/6-31+G(d)/ /HF/6-31+G(d) [–G° = 1.73 kcal/mol], B3P86/6-31G(d) [–G° = 1.75 kcal/mol], and B3PW91/6-31G(d) [–G° = 1.85 kcal/mol] gave conformational free energy (G°) values at 180 K for 4-methyltetrahydro-2H-thiopyran-1,1-dioxide which were similar to the reported experimental values for methylcyclohexane (–G° = 1.80 kcal/mol), 4-methyltetrahydro-2H-thiopyran (–G° = 1.80 kcal/mol), and other 4-methyl-substituted heterocycles. All levels of theory showed that the conformational preferences of the 4-methanoyl (4-formyl), 4-ethanoyl (4-acetyl), and 4-cyano substituents were small. The HF calculations gave conformational free energy (G°) values for 4-chlorotetrahydro-2H-thiopyran-1,1dioxide which were closer to the experimental value than the MP2 and density functional methods. The best agreement with available experimental data for 4-bromotetrahydro-2H-thiopyran-1,1-dioxide was obtained from the HF/6-31G(2d), HF/6-31G(3d), and B3LYP/6-31G(2d) calculations, and, for 4-acetoxytetrahydro-2H-thiopyran-1,1-dioxide, from the HF/6–31G(3d) calculations. The conformational free energies (G°) and relative energies (E) of the conformers and rotamers have been compared with the correspondingly substituted cyclohexanes and tetrahydro-2H-thiopyrans and are discussed in terms of dipole–dipole (electrostatic) interactions and repulsive nonbonded interactions (steric) in the most stable axial and equatorial conformers. The axial S=O bond lengths are shorter than the equatorial S=O bond lengths and the C2–C3 bond lengths in the substituents with carbon-bonded to the ring are shorter than the C3–C4 and C4–C-5 bond lengths. In contrast, the C2–C3 bond lengths in the 4-halogen and 4-acetoxy substituents are longer than the C3–C4 and C4–C-5 bond lengths.     

Conformational Features of Monochloro-substituted Tetrahydropyrans According to Data from <Emphasis Type="Italic">ab initio</Emphasis> Calculations and <Superscript>35</Superscript>Cl NQR

Nonempirical quantum-chemical calculations of 2-, 3-, and 4-chloro-substituted tetrahydropyrans by the RHF/6-31G(d) and MP2/6-31G(d) methods showed that their conformational energy increases in the transition from the 4- to the 3-chloro-substituted compound, while in the transition to the 2-substituted compound the sign changes. In each of these isomers the axial C-Cl and C-H bonds, situated in the geminal position in relation to the oxygen atom, are longer, while the electron density at their Cl and H atoms is higher than in the corresponding equatorial bonds.__________Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 2, pp. 196–201, February, 2005.     

Theoretical studies on the structures and isomerization of methylene lithium‐chlorosilylenoid H2CSiLiCl

The methylene lithium‐chlorosilylenoid H₂C?SiLiCl was studied with ab initio calculations at the G2(MP2) level. Its four equilibrium structures, p‐complex, three‐membered ring, σ complex and silene, and three isomerization transition states were located. The calculations show that the nonplanar p‐complex structure is the lowest in energy among four equilibrium structures of H₂C?SiLiCl and should be experimentally detectable. The silene and σ complex structures with high energies are unstable and easy to isomerize to the most stable p‐complex structure via three‐membered ring one. Also, the geometric characteristics and bonding properties of various structures were analyzed and discussed. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Isomers and [2,3]-Sila-wittig Rearrangement of [（Allyloxy）silyl]lithium Silylenoid in the Gas Phase and THF Solvent

Theoretical calculations of the [2,3]-sila-wittig rearrangement of isomers of [（allyloxy）silyl]lithium （C3H5O）HzSiLi have been performed in the gas phase and THF solvent using the G3MP2B3 method. Seven isomers of silylenoid （C3H5O）H2SiLi, 1-7, are found. The [2,3]-silawittig rearrangement paths are followed using two isomers, 2 and 4, to yield the transition states as well as the products. In the transition state, the silicon center functions as a nucleophile and the aUyl as an electrophile. The interaction between the silicon and allylic sites leads to the formation of SiC（3） bond and the break of O-C（1） bond. Finally, the （allylsilyl）oxylithium （C3H5）H2SiOLi is obtained. The rearrangement paths are confirmed by the intrinsic reaction coordinate （IRC） calculations. The rearrangement mechanisms of reactions of 2 and 4 are similar, and the latter reaction is more favored in the gas phase and THF solvent. Also, the solvent effects are analyzed in this work.     

Isomers and [2,3]-Sila-wittig Rearrangement of [(Allyloxy)sily]lithium Silylenoid in the Gas Phase and THF Solvent

Theoretical calculations of the [2,3]-sila-wittig rearrangement of isomers of [(allyl-oxy)silyl]lithium (C3H50)H2SiLi have been performed in the gas phase and THF solvent using theG3MP2B3 method.Seven isomers of silylenoid (C3H5O)H2SiLi, 1～7, are found.The [2,3]-sila-wittig rearrangement paths are followed using two isomers, 2 and 4, to yield the transition states as well as the products.In the transition state, the silicon center functions as a nucleophile and the allyl as an electrophile.The interaction between the silicon and allylic sites leads to the formation of Si-C[3] bond and the break of O-C[1] bond.Finally, the (allylsilyl)oxylithium (C3H5)H2SiOLi is obtained.The rearrangement paths are confirmed by the intrinsic reaction coordinate (IRC)calculations.The rearrangement mechanisms of reactions of 2 and 4 are similar, and the latter reaction is more favored in the gas phase and THF solvent.Also, the solvent effects are analyzed in this work.     

4-Oxo-3,4-dihydroquinazolinyl- and Benzimidazolylacetonitriles in Annelation Reactions of a Haloquinoline Ring

The interaction of 4-oxo-3,4-dihydroquinazolinyl- and benzimidazolylacetonitriles with 2,6-dihalobenzaldehydes leads to 3-(2,6-dihalophenyl)-2-(4-oxo-3,4-dihydro-2-quinazolinyl)acrylonitriles and 2-(1H-benzo[d]imidazol-2-yl)-3-(2,6-dihalophenyl)acrylonitriles respectively. As a result of intramolecular cyclization of these nitriles 4-halo-12-oxo-12H-quino[2,1-b]quinazoline-6-carbonitriles and 4-halobenzo[4,5]imidazo[1,2-a]quinoline-6-carbonitriles respectively are formed.     

Accurate ab initio prediction of propagation rate coefficients in free-radical polymerization: Acrylonitrile and vinyl chloride

A systematic methodology for calculating accurate propagation rate coefficients in free-radical polymerization was designed and tested for vinyl chloride and acrylonitrile polymerization. For small to medium-sized polymer systems, theoretical reaction barriers are calculated using G3(MP2)-RAD. For larger systems, G3(MP2)-RAD barriers can be approximated (to within 1 kJ mol⁻¹) via an ONIOM-based approach in which the core is studied at G3(MP2)-RAD and the substituent effects are modeled with ROMP2/6-311+G(3df,2p). DFT methods (including BLYP, B3LYP, MPWB195, BB1K and MPWB1K) failed to reproduce the correct trends in the reaction barriers and enthalpies with molecular size, though KMLYP showed some promise as a low cost option for very large systems. Reaction rates are calculated via standard transition state theory in conjunction with the one-dimensional hindered rotor model. The harmonic oscillator approximation was shown to introduce an error of a factor of 2–3, and would be suitable for “order-of-magnitude” estimates. A systematic study of chain length effects indicated that rate coefficients had largely converged to their long chain limit at the dimer radical stage, and the inclusion of the primary substituent of the penultimate unit was sufficient for practical purposes. Solvent effects, as calculated using the COSMO model, were found to be relatively minor. The overall methodology reproduced the available experimental data for both of these monomers within a factor of 2.     

Discrete Metal Catalysts for Stereoselective Ring‐Opening Polymerization of Chiral Racemic β‐Lactones

Poly(β‐hydroxyalkanoate)s (PHAs) are a class of aliphatic polyesters that can be efficiently synthesized by ring‐opening polymerization (ROP) of β‐lactones. The case of chiral racemic β‐substituted β‐lactones is particularly appealing since these monomers open the way to original tacticities and materials different from those biotechnologically produced. In this overview, after briefly surveying general considerations associated to the ROP of β‐lactones and metal‐based catalysts used in stereoselective ROP of racemic β‐butyrolactone, special emphasis is given to discrete rare earth catalysts that have allowed the preparation of highly syndiotactic poly(3‐hydroxybutyrate)s. Recent developments – such as preparation of stereocontrolled PHAs with pendant structural groups via (co)polymerization of functional β‐substituted β‐lactones, and highly alternating copolymers obtained by ROP of mixtures of enantiomerically pure but different monomers – are also discussed.

     

Scale‐up of Microwave‐Assisted Polymerizations in Batch Mode: The Cationic Ring‐Opening Polymerization of 2‐Ethyl‐2‐oxazoline

Summary: The use of microwave heating in polymer science is a rapidly growing field of research leading to faster and cleaner polymerization procedures. However, the majority of the investigations are performed at small scales (≈1 mL), which is far away from potential commercial applications of microwave‐assisted polymerizations. In addition, it has been shown in organic chemistry that microwave‐assisted reaction protocols can be directly scaled without the need for process optimization. In this contribution, we have investigated the direct scaling of microwave‐assisted polymerization procedures under pressure conditions using the cationic ring‐opening polymerization of 2‐ethyl‐2‐oxazoline as the model system. This polymerization was performed at scales ranging from 4.0 mmol (1 mL) to 1.0 mol (250 mL) in different microwave synthesizers covering both monomode and multimode devices.

Scale‐up of microwave‐assisted polymerizations in batch mode: the cationic ring‐opening polymerization of 2 ethyl‐2‐oxazoline.     

A model of the chemical bond must be rooted in quantum mechanics, provide insight, and possess predictive power

In this response to the preceding paper by Bader, we show that the core arguments and statements presented in the latter are flawed. We argue that it is insufficient for a model of the chemical bond to be rooted in quantum mechanics. A good model must in addition provide insight and possess predictive power. Our molecular orbital (MO) model of the chemical bond, in particular, the associated energy-decomposition approach satisfies all these conditions. On the other hand, Atoms-in-Molecules (AIM) theory is only rooted in quantum mechanics as far as its mathematical framework is concerned. The physical status of its central concepts is not so clear. In particular, "bond paths" and "bond critical points" are once more confirmed not to be indicators of a stabilizing interaction. Moreover, AIM theory lacks any predictive power. We also address specific questions raised in the preceding paper. Finally, interpreting chemical bonding implies choosing a perspective on this phenomenon. That there are many perspectives is a matter of fact and this is in no way unphysical. What is unscientific is to claim uniqueness and truth for one of these choices, namely AIM, and to dismiss on this ground all other approaches.