A quantum chemical investigation of pyrolysis reactions of glyoxylic acid ethylester

Two possible reaction paths for the pyrolysis of the ethylester of glyoxylic acid have been studied by ab initio molecular orbital calculations. The basis sets 3-21G and 6-31G * have been used, and electron correlation has been included by Møller–Plesset calculations up to fourth order. Our calculations indicate that the reaction leading to acid and ethylene through a 6-membered ring transition state is favored relative to a process involving a formyl hydrogen transfer via a 5-membered ring to the alkyl unit leading to ethane, CO, and CO₂. The predicted activation energies for these two reactions obtained at the highest level of calculation, MP 4(SDTQ )/6–31G *, are 50.4 and 71.7 kcal/mol, respectively. The transition states have RHF wave functions that are stable relative to UHF solutions using the 3–21G basis. The geometry of the transition states and IRC following indicate that both reactions are strongly asynchronous: The C? O bond rupture is virtually completed before hydrogen transfer occurs. For comparative purposes, analogous calculations have been performed for the ethylester of formic acid, where it is confirmed that a 6-membered ring transition state is preferred relative to a 4-membered one by around 42 kcal/mol at the highest level of calculation.     

Ground and excited states of the monomer and dimer of certain carboxylic acids

The ground-state properties of the monomer and the dimer of formic acid, acetic acid, and benzoic acid have been investigated using Hartree-Fock (HF) and density functional theory (DFT) methods using the 6-311++G(d,p) basis set. Some of the low-lying excited states have been studied using the time-dependent density functional theory (TDDFT) with LDA and B3LYP functionals and also employing complete-active-space-self-consistent-field (CASSCF) and multireference configuration interaction (MRCI) methodologies. DFT calculations predict the ground-state geometries in quantitative agreement with the available experimental results. The computed binding energies for the three carboxylic acid dimers are also in accord with the known thermodynamic data. The TDDFT predicted wavelengths corresponding to the lowest energy n-pi* transition in formic acid (214 nm) and acetic acid (214 nm) and the pi-pi* transition in benzoic acid (255 nm) are comparable to the experimentally observed absorption maxima. In addition, TDDFT calculations predict qualitatively correctly the blue shift (4-5 nm) in the excitation energy for the pi-pi* transition in going from the monomer to the dimer of formic acid and acetic acid and the red shift (approximately 19 nm) in pi-pi* transition in going from benzoic acid monomer to dimer. This also indicates that the electronic interaction arising from the hydrogen bonds between the monomers is marginal in all three carboxylic acids investigated.     

Theoretical insights into the mechanism of acetylcholinesterase-catalyzed acylation of acetylcholine

Acylation of acetylcholine (ACh) catalyzed by acetylcholinesterase (AChE) has been studied using high-level theoretical calculations on a model system that mimics the reaction center of the enzyme, and compared with uncatalyzed acylation reaction. The geometries of all the intermediates and transition states, activation energies, and solvent effects have been calculated. The calculations predict simultaneous formation of two short-strong hydrogen bonds (SSHB) in the rate-determining transition state structures [the first SSHB involves the hydrogen atom of Ser-200 (H(s)) and another involves the hydrogen atom of His-440 (H(h))]. In the intermediate states, the H-bond corresponding to H(h) involves SSHB, whereas the one corresponding to H(s) does not.     

Synergistic effect of Zr-incorporated framework and subsequent deposition of PEHA towards efficient and reusable HPW/PEHA/ZrSBA-15 composites

Heteropolyacid-based catalysts with molecular level distributions of surface functionalities were synthesized by anchoring phosphotungstic acid (HPW) onto the pentaethylenehexamine (PEHA) layers stabilized by platelet ZrSBA-15 with short mesochannels. Namely, the presence of zirconium species in pore walls could significantly improve the surface acidity to attract and stabilize PEHA species, leading to a stronger interaction between them and a superiority to make HPW species highly stable. Standard characterizations of XRD, SEM, TEM, FT-IR, TG, N₂ adsorption–desorption and N elemental analyses were preserved to demonstrate the synergistic effect. The resulting composite were investigated for the oxidation of benzaldehyde to benzoic acid in the H₂O₂-mediated environment, the catalyst 25 % HPW/PEHA/ZrSBA-15 screened high yield for benzoic acid (~99.7 %). Besides, various reaction parameters such as reaction temperature, reaction time, the amount of hydrogen peroxide and catalyst reusability were also investigated for the oxidation of benzaldehyde to benzoic acid with a target to confirm the validity of the catalyst.     

The reaction of atomic oxygen with methanethiol. A theoretical study of the structures and the potential energy surface

Ab initio calculations at the Hartree-Fock, MP2 and MP4 levels were performed to find structures of the equilibrium and transition states and the reaction energies and energies of activation of several competing reaction pathways of O (3P)+CH3SH. A 6-31G* basis set was used in all calculations. The mechanism of hydrogen atom abstraction from the S-H group methanethiol was found to be very competitive with the oxygen atom addition to the sulfur atom.     

乙炔氢氯化反应铜系催化剂的反应机理

以电石乙炔法制备氯乙烯的非汞催化反应体系为研究对象,用量子化学密度泛函理论(DFT)研究了以Cu基催化剂为代表的非汞催化剂的反应机理.模拟了以石墨烯为载体的乙炔氢氯化反应在Cu基催化剂作用下的两条反应路径及其过渡态.将该反应机理应用到不同金属氯化物催化剂中,通过计算相应的活化自由能和反应速率常数讨论了不同金属氯化物的活性顺序,与实验结果相比较印证了反应机理的合理性.讨论了Cu基催化剂的失活原因以及在载体中掺杂氮原子或磷原子对反应活性的影响.为非贵金属无汞催化剂的研究提供了一定的理论指导.     

Comparison of AM1 and PM3 semiempirical to ab initio methods in the study of Diels-Alder reactions of butadiene and cyclopentadiene with cyanoethylenes

Assuming a concerted synchronous mechanism with one transition state of the Diels-Alder reactions, the structures of the transition states and the activation energies for the reactions of butadiene and cyclopentadiene with cyanoethylenes were calculated by AM1 and PM3 semiempirical methods. The structural parameters were compared with those obtained by high level Gaussian calculations, whereas the activation energies were compared both with the ab initio calculations and those obtained experimentally. The structural properties calculated with PM3 methods are in general in better agreement with the ab initio calculations. The low level ab initio calculations are in many cases worse than the semiempirical methods. All predicted activation energies with both semiempirical methods are up to 300% higher than the experimental values. The predicted reactivity is also opposite to the experimental data. Only the very high level Gaussian calculations are in good correlation with experimental results. The predicted selectivity of the reaction is also opposite to the experimental facts. Two explanations are offered for this discrepancy: AM1 and PM3 methods cannot handle the calculation of the concerted Diels-Alder transition states and are not recommended to be used for that purpose, or this Diels-Alder reaction is not concerted but is stepwise.     

Stereocartography: a computational mapping technique that can locate regions of maximum stereoinduction around chiral catalysts

A hypothesis concerning asymmetric induction by chiral catalysts is posited, tested, and found to be valid. The hypothesis states that chiral catalysts that are efficient at inducing asymmetry will have their region of maximum stereoinduction spatially congruent with the site of chemistry but inefficient catalysts will not. A simple mapping strategy (stereocartography) is used to assess where the region of maximum stereoinduction is located around a given catalyst. The protocol compares interaction energies between mirror image probes at each point in space around the catalyst being considered. The probes are models of the actual transition states of the reaction being catalyzed by a particular catalyst. The hypothesis was tested on three Diels-Alder reactions. Seventeen of the eighteen catalysts conform to the hypothesis. The idea of using this as a catalyst design tool is presented.     

Experimental and DFT studies of the conversion of ethanol and acetic acid on PtSn-based catalysts

Reaction kinetics studies were conducted for the conversions of ethanol and acetic acid over silica-supported Pt and Pt/Sn catalysts at temperatures from 500 to 600 K. Addition of Sn to Pt catalysts inhibits the decomposition of ethanol to CO, CH4, and C2H6, such that PtSn-based catalysts are active for dehydrogenation of ethanol to acetaldehyde. Furthermore, PtSn-based catalysts are selective for the conversion of acetic acid to ethanol, acetaldehyde, and ethyl acetate, whereas Pt catalysts lead mainly to decomposition products such as CH4 and CO. These results are interpreted using density functional theory (DFT) calculations for various adsorbed species and transition states on Pt(111) and Pt3Sn(111) surfaces. The Pt3Sn alloy slab was selected for DFT studies because results from in situ (119)Sn M?ssbauer spectroscopy and CO adsorption microcalorimetry of silica-supported Pt/Sn catalysts indicate that Pt-Sn alloy is the major phase present. Accordingly, results from DFT calculations show that transition-state energies for C-O and C-C bond cleavage in ethanol-derived species increase by 25-60 kJ/mol on Pt3Sn(111) compared to Pt(111), whereas energies of transition states for dehydrogenation reactions increase by only 5-10 kJ/mol. Results from DFT calculations show that transition-state energies for CH3CO-OH bond cleavage increase by only 12 kJ/mol on Pt3Sn(111) compared to Pt(111). The suppression of C-C bond cleavage in ethanol and acetic acid upon addition of Sn to Pt is also confirmed by microcalorimetric and infrared spectroscopic measurements at 300 K of the interactions of ethanol and acetic acid with Pt and PtSn on a silica support that had been silylated to remove silanol groups.     

Density Functional Theory Study of Reaction Equilibria in Signal Amplification by Reversible Exchange

An in-depth theoretical analysis of key chemical equilibria in Signal Amplification by Reversible Exchange (SABRE) is provided, employing density functional theory calculations to characterize the likely reaction network. For all reactions in the network, the potential energy surface is probed to identify minimum energy pathways. Energy barriers and transition states are calculated, and harmonic transition state theory is applied to calculate exchange rates that approximate experimental values. The reaction network energy surface can be modulated by chemical potentials that account for the dependence on concentration, temperature, and partial pressure of molecular constituents (hydrogen, methanol, pyridine) supplied to the experiment under equilibrium conditions. We show that, under typical experimental conditions, the Gibbs free energies of the two key states involved in pyridine-hydrogen exchange at the common Ir-IMes catalyst system in methanol are essentially the same, i. e., nearly optimal for SABRE. We also show that a methanol-containing intermediate is plausible as a transient species in the process.     

Hydrogen spillover mechanism on a Pd-doped Mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports-for the first time by ab initio simulation-the proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.     

Ring-closing olefin metathesis on ruthenium carbene complexes: model DFT study of stereochemistry

Ring-closing metathesis (RCM) is the key step in a recently reported synthesis of salicylihalamide and related model compounds. Experimentally, the stereochemistry of the resulting cycloolefin (cis/trans) depends strongly on the substituents that are present in the diene substrate. To gain insight into the factors that govern the observed stereochemistry, density functional theory (DFT) calculations have been carried out for a simplified dichloro(2-propylidene)(imidazole-2-ylidene)ruthenium catalyst I, as well as for the real catalyst II with two mesityl substituents on the imidazole ring. Four model substrates are considered, which are closely related to the systems studied experimentally, and in each case, two pathways A and B are possible since the RCM reaction can be initiated by coordination of either of the two diene double bonds to the metal center. The first metathesis yields a carbene intermediate, which can then undergo a second metathesis by ring closure, metallacycle formation, and metallacycle cleavage to give the final cycloolefin complex. According to the DFT calculations, the stereochemistry is always determined in the second metathesis reaction, but the rate-determining step may be different for different catalysts, substrates, and pathways. The ancillary N-heterocyclic carbene ligand lies in the Ru-Cl-Cl plane in the simplified catalyst I, but is perpendicular to it in the real catalyst II, and this affects the relative energies of the relevant intermediates and transition states. Likewise, the introduction of methyl substituents in the diene substrates influences these relative energies appreciably. Good agreement with the experimentally observed stereochemistry is only found when using the real catalyst II and the largest model substrates in the DFT calculations.     

Transition state docking: a probe for noncovalent catalysis in biological systems. Application to antibody-catalyzed ester hydrolysis

A strategy for pinpointing favorable noncovalent interactions between transition states and active sites of biological catalysts is described. This strategy combines high-level quantum mechanical calculations of transition state geometries with an automated docking procedure using AutoDock. By applying this methodology to antibody-catalyzed hydrolyses of aryl esters (by the 48G7, CNJ206, and 17E8 families of antibodies), varying levels of catalysis are explained in terms of specific hydrogen bonding interactions between combining site residues and transition states. Although these families of antibodies were produced in separate experiments by different researchers using related but different haptens, the mechanism of transition state stabilization appears to be highly conserved. Despite being elicited in response to anionic phosphonate haptens, the best catalysts often utilize hydrogen bond acceptors to stabilize transition states. A mutant of antibody CNJ206, designed based on this observation and predicted to be a better catalyst, is proposed. In the case of antibody 48G7, affinity maturation is shown to produce a catalyst that is highly selective for one of two enantiomeric transition states from a nonselective germline precursor.     

Reactivity of alkanes on zeolites: a computational study of propane conversion reactions

In this work, quantum chemical methods were used to study propane conversion reactions on zeolites; these reactions included protolytic cracking, primary hydrogen exchange, secondary hydrogen exchange, and dehydrogenation reactions. The reactants, products, and transition-state structures were optimized at the B3LYP/6-31G level and the energies were calculated with CBS-QB3, a complete basis set composite energy method. The computed activation barriers were 62.1 and 62.6 kcal/mol for protolytic cracking through two different transition states, 30.4 kcal/mol for primary hydrogen exchange, 29.8 kcal/mol for secondary hydrogen exchange, and 76.7 kcal/mol for dehydrogenation reactions. The effects of basis set for the geometry optimization and zeolite acidity on the reaction barriers were also investigated. Adding extra polarization and diffuse functions for the geometry optimization did not affect the activation barriers obtained with the composite energy method. The largest difference in calculated activation barriers is within 1 kcal/mol. Reaction activation barriers do change as zeolite acidity changes, however. Linear relationships were found between activation barriers and zeolite deprotonation energies. Analytical expressions for each reaction were proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts, as long as the deprotonation energies are first acquired.     

Study of solid-state catalytic isotope exchange of hydrogen inl-hydroxyproline under the action of spillover tritium

Reaction of high-temperature solid-state catalytic isotope exchange (HSCIE) of hydrogen in L-hydroxyproline was studied byab initio quantum-chemical calculations. A one-center synchronous mechanism of isotope exchange between the amino acid and the H₃O⁺ model acidic center was considered. The structures of transition states of the reaction and the activation energies were determined. Relative reactivity of the C−H bonds in the hydroxyproline molecule under conditions of HSCIE was studied. The results obtained are in agreement with experimental data on the stereoselectivity and regioselectivity of the HSCIE reaction,viz., the lower the calculated activation energy of isotope exchange, the larger the portion of hydrogen substituted by tritium in a given position of the amino acid molecule. The enhancement of the reactivity under conditions of solid-state isotope exchange can be associated with additional interaction between the exchanging H atoms and the electron-donor O and N atoms of the amino acid molecule in transition state. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1056–1060, June, 1999.     

Mechanistic study of asymmetric Michael addition of malonates to enones catalyzed by a primary amino acid lithium salt

A mechanistic study was carried out for the asymmetric Michael addition reaction of malonates to enones catalyzed by a primary amino acid lithium salt to elucidate the origin of the asymmetric induction. A primary β-amino acid salt catalyst, O-TBDPS β-homoserine lithium salt, exhibited much higher enantioselectivity than that achieved with the corresponding catalysts derived from α- and γ-amino acids for this reaction. Detailed studies of the transition states with DFT calculations revealed that the lithium cation and carboxylate group of the β-amino acid salt catalyst have important roles in achieving high enantioselectivity in the Michael addition reaction of malonates to enones.     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

Role of rare earth cations in Y zeolite for hydrocarbon cracking

Reaction kinetics data were collected for isobutane conversion over a series of ultra stable Y (USY) zeolite catalysts with and without rare earth cations and subjected to various extents of dealumination by steaming. We conducted these reaction studies at low temperatures (523-573 K) using isobutane feed streams containing known levels of isobutylene (100-400 ppm) so that the kinetics were controlled by bimolecular hydride transfer and oligomerization/beta-scission processes with little or no participation of monomolecular initiation reactions. These experimental conditions led to stable catalyst performance with the main products of isobutane conversion being propane, n-butane, and isopentane, with smaller amounts of propylene, trans-2-butene, and cis-2-butene. The rates of formation of these products per Br?nsted acid site (as counted by pyridine adsorption) depended exponentially on Br?nsted acid site density, regardless of whether the catalyst contained rare earth cations. Kinetic modeling showed an exponential dependence of hydride transfer and oligomerization/ beta-scission reaction rates on Br?nsted acid site density which translated into composite activation energies for these reactions having a linear relationship with site density. Based on results in the literature from theoretical calculations, we suggest that increasing Br?nsted acid site density in zeolite Y leads to larger zeolite elasticity, increased stabilization of cationic transition states, and lower composite activation barriers for hydride transfer and beta-scission steps. The role of rare earth cations, therefore, is to ensure the retention of high Br?nsted acid site density under hydrothermal conditions, such as in fluid catalytic cracking (FCC) regenerators, where steam would dealuminate the Y zeolite framework and reduce this site density. It is for this reason that hydride transfer reaction rates are high in the presence of rare earth cations and lead to higher yields of less olefinic gasoline during FCC.     

The isomerization of methoxy radical: intramolecular hydrogen atom transfer mediated through acid catalysis

The catalytic ability of water, formic acid, and sulfuric acid to facilitate the isomerization of the CH(3)O radical to CH(2)OH has been studied. It is shown that the activation energies for isomerization are 30.2, 25.7, 4.2, and 2.3 kcal mol(-1), respectively, when the reaction is carried out in isolation and with water, formic acid, or sulfuric acid as a catalyst. The formation of a doubly hydrogen bonded transition state is central to lowering the activation energy and facilitating the intramolecular hydrogen atom transfer that is required for isomerization. The changes in the rate constant for the CH(3)O-to-CH(2)OH isomerization with acid catalysis have also been calculated at 298 K. The largest enhancement in the rate, by over 12 orders of magnitude, is with sulfuric acid. The results of the present study demonstrate the feasibility of acid catalysis of a gas-phase radical isomerization reaction that would otherwise be forbidden.     

Activation Energy and Transition State Determination of the Olefin Insertion Process of Metallocene Catalysts Using a Semiempirical Molecular Orbital Calculation

The transition state of the olefin insertion process of metallocene catalysts can be determined by adopting the semiempirical PM3 model. In computational chemistry, the computational methods most employed are the ab initio method and density functional theory, which are very time consuming. The semiempirical molecular orbital method requires much less computational resources than the above methods. However, the accuracy and reliability of the semiempirical molecular orbital method remains to be determined. The PM3 model is the most recently developed the semiempirical molecular orbital method and can also be applied to transition metal calculations. This study is intended to investigate the reliability of computational results determined using semiempirical PM3 model on metallocene catalysts through comparison with published results on the density functional theory (DFT). The saddle point finding procedure is adopted to find the transition state of the ethylene insertion process of metallocene catalysts. Results on the geometry and energy trends of the ethylene insertion process of metallocene catalysts determined using the PM3 model are in good agreement with the DFT results. In addition, the saddle point of the potential energy surface of ethylene insertion is verified in accordance with the eigenvalue of the vibrational frequency spectrum. Correct eigenvalues indicate that the correct saddle point of the potential energy surface of ethylene insertion has been successfully located. Hence, the eigenvalue of the vibrational frequency spectrum is a valuable reference in terms of saddle point justification. Computational results and vibrational frequency spectrum analysis demonstrate that the PM3 model can be used to locate the correct saddle point of the potential energy surface. The results obtained using the PM3 model confirm that the eigenvalue of the transition state lies nearly on the vibrational frequency spectrum. The eigenvalues are also analyzed, providing a valuable reference for further studies of the transition state of olefin insertion of metallocene catalysts. The activation energies for the olefin insertion reaction are also studied for evaluation of the catalyst.