Mechanisms of glycerol dehydration

Dehydration of neutral and protonated glycerol was investigated using quantum mechanical calculations (CBS-QB3). Calculations on neutral glycerol show that there is a high barrier for simple 1,2-dehydration, E(a)=70.9 kcal mol(-1), which is lowered to 65.2 kcal mol(-1) for pericyclic 1,3-dehydration. In contrast, the barriers for dehydration of protonated glycerol are much lower. Dehydration mechanisms involving hydride transfer, pinacol rearrangement, or substitution reactions have barriers between 20 and 25 kcal mol(-1). Loss of water from glycerol via substitution results in either oxirane or oxetane intermediates, which can inter-convert over a low barrier. Subsequent decomposition of these intermediates proceeds via either a second dehydration step or loss of formaldehyde. The computed mechanisms for decomposition of protonated glycerol are supported by the gas-phase fragmentation of protonated glycerol observed using a triple--quadrupole mass spectrometer.     

Computational study on kinetics and mechanisms of unimolecular decomposition of succinic acid and its anhydride

The mechanisms and kinetics of unimolecular decomposition of succinic acid and its anhydride have been studied at the G2M(CC2) and microcanonical RRKM levels of theory. It was shown that the ZsgsZ conformer of succinic acid, with the Z-acid form and the gauche conformation around the central C-C bond, is its most stable conformer, whereas the lowest energy conformer with the E-acid form, ECGsZ, is only 3.1 kcal/mol higher in energy than the ZsgsZ. Three primary decomposition channels of succinic acid producing H2O + succinic anhydride with a barrier of 51.0 kcal/mol, H2O + OCC2H3COOH with a barrier of 75.7 kcal/mol and CO2 + C2H5COOH with a barrier of 71.9 kcal/mol were predicted. The dehydration process starting from the ECGCZ-conformer is found to be dominant, whereas the decarboxylation reaction starting from the ZsgsZ-conformer is only slightly less favorable. It was shown that the decomposition of succinic anhydride occurs via a concerted fragmentation mechanism (with a 69.6 kcal/mol barrier), leading to formation of CO + CO2 + C2H4 products. On the basis of the calculated potential energy surfaces of these reactions, the rate constants for unimolecular decomposition of succinic acid and its anhydride were predicted. In addition, the predicted rate constants for the unimolecular decomposition of C2H5COOH by decarboxylation (giving C2H6 + CO2) and dehydration (giving H3CCHCO + H2O) are in good agreement with available experimental data.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

Computational study on the kinetics and mechanisms for the unimolecular decomposition of formic and oxalic acids

The kinetics and mechanisms for the unimolecular decomposition reactions of formic acid and oxalic acid have been studied computationally by the high-level G2M(CC1) method and microcanonical RRKM theory. There are two reaction pathways in the decomposition of formic acid: The dehydration process starting from the Z conformer is found to be the dominant, whereas the decarboxylation reaction starting from the E conformer is less competitive. The predicted rate constants for the dehydration and decarboxylation reactions are in good agreement with the experimental data. The calculated CO/CO2 ratio, 13.6-13.9 between 1300 and 2000 K, is in close agreement with the ratio of 10 measured experimentally by Hsu et al. (In The 19th International Symposium on Combustion; The Combustion Institute: Pittsburgh, PA, 1983; p 89). For oxalic acid, its isomer with two intramolecular hydrogen bonds is the most stable structure, similar to earlier reports. Two primary decomposition channels of oxalic acid producing CO2+HOCOH with barriers of 33-36 kcal/mol and CO2+CO+H2O with a barrier of 39 kcal/mol were found. At high temperatures, the latter process becomes more competitive. The rate constant predicted for the formation of CO2 and HOCOH (the precursor of HCOOH) agrees well with available experimental data. The mechanism for the isomerization of HOCOH to HCOOH is also discussed.     

Mechanism of HCS + O2 reaction: hydrogen- or oxygen-transfer?

In spite of the potential importance of the HCS radical in both combustion and interstellar processes, its chemical reactivity has not been tackled previously. In the present paper, the oxidation reaction of the HCS radical is theoretically investigated for the first time at the CCSD(T)/6-311++G(3df,2p)//BH&HLYP/6-311++G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels. It is shown that the most feasible pathway is the O2 addition to the HCS radical forming the intermediate SC(H)OO which can undergo a subsequent O-extrusion leading to SC(H)O + 3O. This features an indirect O-transfer mechanism with the overall barrier of 4.4 and 3.5 kcal mol(-1), respectively, at the two levels. However, formation of the H-transfer product CS + HO2 is kinetically much less feasible, i.e., the direct mechanism has barriers of 14.3 and 8.7 kcal mol(-1), whereas the indirect mechanism has barriers of 12.6 and 10.7 kcal mol(-1), respectively. This result is in sharp contrast to the analogous HCO + O2 reaction, where the direct (with a barrier of 2.98 kcal mol(-1)) and indirect (2.26 kcal mol(-1)) H-transfer processes are highly competitive over the indirect O-transfer process (the least endothermicity is 19.9 kcal mol(-1)). The possible explanations and implications of the present results are provided.     

Hybrid quantum mechanics/molecular mechanics-based molecular dynamics simulation of acid-catalyzed dehydration of polyols in liquid water

We use the conversion of protonated glycerol to acrolein for a case study of the mechanism of acid-catalyzed dehydration of polyols in aqueous environments. We employ hybrid Quamtum Mechanics/Molecular Mechanics Molecular Dynamics (QM/MM MD) simulations with biased sampling and perform free energy calculations for the elementary steps of the reaction. We investigate the effects of solvent dynamics and in particular the role of quantum mechanical water in the dehydration mechanism. We present results supporting a mechanism that proceeds via water-mediated proton transfers and thus through an enol intermediate. We find that the first dehydration may take place by two, low-energy pathways requiring, respectively, 20.9 and 18.8 kcal/mol of activation free energy. The second dehydration requires 19.9 kcal/mol of activation free energy while for the overall reaction we compute a free energy change of -8 kcal/mol.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Structures and reaction mechanisms of glycerol dehydration over H-ZSM-5 zeolite: a density functional theory study

The initial stage of glycerol conversion over H-ZSM-5 zeolite has been investigated using density functional theory (DFT) calculations on an embedded cluster model consisting of 128 tetrahedrally coordinated atoms. It is found that glycerol dehydration to acrolein and acetol proceeds favourably via a stepwise mechanism. The formation of an alkoxide species upon the first dehydration requires the highest activation energy (42.5 kcal mol(-1)) and can be considered as the rate determining step of the reaction. The intrinsic activation energies for the first dehydration are virtually the same for both acrolein and acetol formation, respectively, suggesting the competitive removal of the primary and secondary OH groups. A high selectivity to acrolein at moderate temperatures can be attributed to the selective activation of the stronger adsorption mode of glycerol through the secondary OH group and the kinetically favoured subsequent consecutive steps. In addition, the less reactive nature of acrolein relative to acetol precludes it from being converted to other products upon conversion to glycerol. In accordance with typical endothermic reactions, the forward rate constant for glycerol dehydration significantly increases with increasing reaction temperature.     

Computational study of the aminolysis of 2-benzoxazolinone

Three possible mechanisms (zwitterionic, neutral stepwise, and neutral concerted) of the ring-opening reaction of 2-benzoxazolinone (BO) upon aminolysis with methylamine were studied at the B3LYP/6-31G* level. In the gas phase, the neutral concerted mechanism is shown to be most favorable, which proceeds via a rate-determining barrier of 28-29 kcal/mol. The transition state, CTS, associated with this barrier is a four-centered one, where 1,2-addition of the N[bond]H of methylamine to the C[bond]O of BO ring occurs. The rate-determining barrier of the neutral stepwise pathway is found to be ca. 42 kcal/mol. The inclusion of solvent effects by a polarizable continuum model (PCM) does not change the conclusions based on the gas-phase study; the barrier at CTS is reduced to 20, 20, and 22 kcal/mol in water, ethanol, and acetonitrile, respectively.     

Hydrogen atom mediated Stone-Wales rearrangement of pyracyclene: a model for annealing in fullerene formation

We have investigated the Stone-Wales (SW) rearrangement of pyracyclene (C(14)H(12)) using quantum mechanical molecular modeling. Of particular interest in this study is the effect of an added hydrogen atom on the barriers to SW rearrangement. Hydrogen atoms are found in high abundance during combustion, and their effect upon isomerization of aromatic compounds to more stable species may play an important role in the combustion synthesis of fullerenes. We have calculated the barriers for the SW rearrangement in pyracyclene using density functional theory B3LYP/6-31G(d) and B3LYP/6-311G(d,p). Two mechanisms have been investigated: (i) a mechanism with two identical transition states of C(1) symmetry and a cyclobutyl intermediate and (ii) a mechanism with one transition state containing an sp(3) carbon (J. Am. Chem. Soc. 2003, 125, 5572-5580; Nature 1993, 366, 665-667). We find that the barriers for these mechanisms are 120.0 kcal mol(-1) for the cyclobutyl mechanism and 130.1 kcal mol(-1) for the sp(3) mechanism. Adding a hydrogen atom to the internal bridge carbon atoms of pyracyclene reduces the barrier of the cyclobutyl mechanisms to 67.0 kcal mol(-1) and the sp(3) mechanism to 73.1 kcal mol(-1). The bonding of carbon atoms in pyracyclene is similar to those found in isomers of C(60), and the barriers are low enough so that these reactions can become significant during fullerene synthesis in flames. Adding hydrogen atoms to the external bridge atoms on pyracyclene produces a smaller reduction in the SW barrier and adding hydrogen atoms to nonbridge external carbon atoms results in no reduction of the barrier.     

Mechanisms of the Knoevenagel hetero Diels-Alder sequence in multicomponent reactions to dihydropyrans: experimental and theoretical investigations into the role of water

The role of water in a multicomponent domino reaction (MCR) involving styrene, 2,4-pentanedione, and formaldehyde was studied. Whereas anhydrous conditions produced no reaction, the MCR successfully proceeded in the presence of water, affording the targeted dihydropyran derivatives with good yield. The mechanism of this MCR (Knoevenagel hetero Diels-Alder sequence) was studied with and without explicit water molecules using the SMD continuum solvation model in combination with the B3LYP density functional and the 6-311++G** basis set to compute the water and acetone (aprotic organic solvent) solution Gibbs free energies. In the Knoevenagel step, we found that water acted as a proton relay to favor the formation of more flexible six-membered ring transition state structures both in concerted (direct H(2)O elimination) and stepwise (keto-enol tautomerization and dehydration) pathways. The inclusion of a water molecule in our model resulted in a significant decrease (-8.5 kcal mol(-1)ΔG(water)(?)) of the direct water elimination activation barrier. Owing to the presence of water, all chemical steps involved in the MCR mechanism had activation free energies barriers lower than 39 kcal mol(-1) at 25 °C in aqueous solvent (<21 kcal mol(-1) ZPE corrected electronic energies barriers). Consequently, the MCR proceeded without the assistance of any catalyst.     

Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine. A theoretical study

The influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine has been studied using density functional theory (DFT) methods. For the neutral system, two different pathways (with retention and inversion of configuration at the C1' anomeric carbon) have been found, both of them consisting of two steps and involving the formation of a dihydrofurane-like intermediate. The Gibbs free energy barrier for the first step is very high in both cases (53 and 46 kcal/mol for the process with inversion and with retention, respectively). However, the N7-protonated system shows a very different mechanism which consists of two steps. The first one leads to the formation of an oxacarbenium ion intermediate, with a Gibbs free energy barrier of 27 kcal/mol, and the second one corresponds to the nucleophilic attack of the water molecule to the oxacarbenium ion and takes place with a barrier of 1.3 kcal/mol. Thus, these results agree with a stepwise SN1 mechanism (DN*AN), with a discrete intermediate formed between the leaving group and the nucleophile approach, and show that N7 protonation strongly catalyzes the hydrolysis of the N-glycosidic bond, making the guanine a better leaving group. Finally, kinetic isotope effects have been calculated for the protonated system, and the results obtained are in very good agreement with experimental data for analogous systems.     

Barriers to rotation in methyl formate by dynamic NMR spectroscopy and barriers to 1,3 oxygen-to-oxygen migration in methyl formate and trifluoromethyl formate by ab initio calculations

Free-energy barriers of 9.85 and 11.91 +/- 0.15 kcal/mol at -70.8 degrees C were found by dynamic NMR spectroscopy for the E-to-Z and Z-to-E conversions, respectively, of methyl formate (1) enriched in 13C to 99% for the carbonyl carbon [methyl formate 13C (2)]. These barriers are higher than the literature values reported for -53 degrees C. The free-energy barrier to 1,3 oxygen-to-oxygen migration of the methyl group in methyl formate was determined by ab initio calculations at several levels. The value of 58.7 kcal/mol obtained at the MP2/6-311+G (df,pd) level was compared to a literature barrier for this process (MINDO/3) and to barriers for related compounds. A free-energy barrier of 63.0 kcal/mol for the oxygen - to - oxygen migration of the CF3 group in trifluoromethyl formate (3) was calculated at the MP2/6-31+G level.     

Insights into the mechanism of methionine oxidation catalyzed by metal (Cu2+, Zn2+, and Fe3+)—Amyloid beta (Aβ) peptide complexes: A computational study

In this DFT study, a mechanism of the oxidation of methionine (Met) amino acid residue catalyzed by the metal (Cu²⁺, Zn²⁺, and Fe³⁺) bound amyloid beta (Aβ) peptide has been proposed. Based on experimental information, two different mechanisms: (1) stepwise and (2) concerted mechanisms for this important process have been investigated. The B3LYP calculations suggest that in the stepwise mechanism, the two separate pathways leading to the same sulfoxide product [Met(O)] go through prohibitively high barriers of 27.3 and 35.1 kcal/mol, therefore it is ruled out. In the concerted mechanism, the Cu²⁺‐Aβ complex has been found to be the most efficient catalyst with the computed barrier of 14.3 kcal/mol. The substitutions of Cu²⁺ by Zn²⁺ and Fe³⁺ increase barriers to 19.6 and 16.9 kcal/mol, respectively and make the reaction thermodynamically less favorable. It was also found that, in comparison with the cysteine (Cys) residue, Met is more susceptible toward oxidation. Its substitution with Cys slightly increased the barrier to 15.8 kcal/mol for the Cu²⁺‐Aβ complex. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

Mechanistic pathways of the hydroxyl radical reactions of quinoline. 2. Computational analysis of hydroxyl radical attack at C atoms

Density functional theory (DFT) calculations are employed to compare the mechanism of the *OH attacks at all carbon atoms in quinoline. The computational analysis of the energy surface for the reaction of *OH with quinoline reveals that the formation of OH adducts proceeds through exothermic formation of pi-complexes/H-bonded complexes. The gas-phase reactions have activation energies ranging from <1.3 kcal/mol for the attack at positions C3 through C8 to 8.6 kcal/mol for the attack at the C2 position. Solvation, as described by the CPCM cavity model, lowers these activation barriers so that the attack at all carbon atoms except C2 is effectively barrierless. The *OH attack at C2 in solution is significantly different than at all other quinoline positions because it involves the only transition structure with energy higher than that of the starting materials and with an energetic barrier of 5.1 kcal/mol. The specific solvation approach also corroborates this finding because the attack at C2 was shown to have an energy barrier of 2.3 kcal/mol compared to the barrierless attack at C5. These results are in agreement with our recent experimental studies but differ from literature reports on the degradation of quinoline using the photo-Fenton reaction.     

On the mechanisms of oxidation of organic sulfides by H2O2 in aqueous solutions

The mechanism of oxidation of organic sulfides in aqueous solutions by hydrogen peroxide was investigated via ab initio calculations. Specifically, two reactions, hydrogen transfer of hydrogen peroxide to form water oxide and the oxidation of dimethyl sulfide (DMS) by hydrogen peroxide to form dimethyl sulfoxide, were studied as models of these processes in general. Solvent effects are included both via including explicitly water molecules and via the polarizable continuum model. The former was found to have a much more significant effect than the latter. When explicit water molecules are included, a mechanism different from those proposed in the literature was found. Specific interactions including hydrogen bonding with 2-3 water molecules can provide enough stabilization for the charge separation of the activation complex. The energy barrier of the oxidation of DMS by hydrogen peroxide was estimated to be 12.7 kcal/mol, within the experimental range of the oxidation of analogous compounds (10-20 kcal/mol). The major reaction coordinates of the reaction are the breaking of the O-O bond of H2O2 and the formation of the S-O bond, the transfer of hydrogen to the distal oxygen of hydrogen peroxide occurring after the system has passed the transition state. Reaction barriers of the hydrogen transfer of H2O2 are an average of 10 kcal/mol or higher than the reaction barriers of the oxidation of DMS. Therefore, a two-step oxidation mechanism in which, first, the transfer of a hydrogen atom occurs to form water oxide and, second, the transfer of oxygen to the substrate occurs is unlikely to be correct. Our proposed oxidation mechanism does not suggest a pH dependence of oxidation rate within a moderate range around neutral pH (i.e., under conditions in which hydronium and hydroxide ions do not participate directly in the reaction), and it agrees with experimental observations over moderate pH values. Also, without including a protonated solvent molecule, it has activation energies that correspond to measured activation energies.     

HZSM-5催化甲苯和甲醇烷基化反应机理的密度泛函理论研究

对二甲苯(PX)是重要的有机化工原料,主要用于生产对苯二甲酸(PTA)和对苯二甲酸二甲酯(DMT), PTA和 DMT可经缩聚生产化纤、合成树脂和塑料等聚酯产品. PX主要通过甲苯歧化、二甲苯异构化或甲苯与 C9芳烃烷基转移等方式生产.由于三种二甲苯和乙苯的沸点接近,需要经过吸附分离或深冷分离才能得到高纯度的 PX,传统工艺物料循环量大,设备庞大,操作费用高.而通过甲苯和甲醇烷基化反应直接高选择性生成 PX,可大大降低成本,具有非常高的经济效益和研究价值.自1970年代以来,国内外众多科研院所对甲苯和甲醇烷基化催化剂进行了广泛研究,但催化剂选择性和稳定性仍需进一步提高.为了加深对甲苯和甲醇烷基化反应的认识,指导催化剂开发,有必要对甲苯和甲醇烷基化生成二甲苯的反应机理进行深入研究.当前甲苯和甲醇烷基化机理研究主要存在以下问题：(1)计算得到的能量多为电子能,而非自由能;(2)所采用的模型多为团簇模型,使用 ONIOM方法,对长程作用力描述不充分;(3)认为甲苯只有一种吸附状态;(4)没有考虑偕烷基化反应.本文采用周期性模型,通过密度泛函理论研究了 HZSM-5分子筛上甲苯和甲醇烷基化反应机理,通过计算熵得到了反应自由能,并考虑了偕烷基化反应.由于甲基的存在,在甲苯的吸附态中,甲基会伸向孔道的不同方向,因此我们认为甲苯有多种吸附态,而不同的吸附态会生成不同的二甲苯.结果表明,甲苯可以在对位、间位、邻位和偕位上通过协同机理或分步机理发生烷基化反应.在协同机理中,甲苯在对位、间位、邻位和偕位发生烷基化反应的自由能垒分别为167,138,139和183 kJ/mol.在分步机理中,甲醇脱水生成甲氧基的自由能垒为145 kJ/mol,是决速步骤;而甲苯和甲氧基对位、间位、邻位和偕位烷基化的自由能垒分别为127,105,106和114 kJ/mol.两种机理中 PX的生成能垒均比 MX和 OX高,与文献报道的结果不同.文献均认为, PX的生成能垒最低.一方面这可能是由于所采用模型的不同,本文采用周期性模型,能更充分考虑长程作用力的影响;另一方面可能是由于对甲苯吸附态的不同处理,我们认为甲苯有多种吸附态,不同的吸附态会生成不同的二甲苯,而文献均只考虑了一种甲苯吸附态.但是,在实验中, PX选择性最高.这可能是由于：(1) PX在 HZSM-5孔道的扩散速率比 MX和 OX高2–3个数量级;(2)甲苯和甲醇烷基化生成的 MX和OX迅速发生异构化反应生成 PX,异构化反应速率高于甲苯烷基化速率.两种机理中, C8H11+都是重要的中间物种,它可以反馈一个质子给分子筛骨架,生成二甲苯;也可以脱烷基生成甲烷和乙烯等气相产物.研究发现,甲烷的生成是由于 C8H11+物种中的一个 H质子从苯环上的碳原子转移到甲基上的碳原子造成的,计算得到的对位、间位和邻位 C8H11+生成甲烷的能垒分别为136,132和134 kJ/mol.由于十元环孔道的限制, HZSM-5孔道中很难通过甲苯歧化反应生成苯;偕烷基化生成的碳正离子有可能脱烷基生成乙烯和乙烷等产物,进而生成苯.碳正离子脱烷基反应生成了大量气相产物,造成反应液收降低.碳正离子脱烷基反应与甲醇制烯烃过程的烃池机理相一致,因此甲苯和甲醇烷基化反应也遵循烃池机理.     

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.     

Exploring a reaction mechanism for acetato ligand replacement in paddlewheel tetrakisacetatodirhodium (II,II) complex by ammonia: computational density functional theory study

This study focuses on the first step of interaction between DNA and the paddle-wheel dirhodium complex. The ammonia molecule was used to model the oligonucleotide sequence. The reaction was considered in neutral and acidic conditions, in gas phase, and in solvent, using the COSMO model. Molecular structures of the complexes were optimized in both models at the B3PW91/6-31G(d) level. The B3LYP functional and aug-cc-pvdz basis set were employed for single-point energy determination and electron distribution analyses. It was shown that in neutral solution the replacement of axial aqua ligand is mildly exoergic. The reaction is characterized by a relatively low activation barrier (10-12 kcal/mol), and, according to Eyring transition state theory, it proceeds very quickly. The breaking of the Rh-O(ac) bond in neutral solution is mildly endoergic (less than 1 kcal/mol) with an activation barrier of about 21 kcal/mol. However, this process can occur much more spontaneously (ΔG of -14 kcal/mol) when the dirhodium complex is protonated at the acetyl oxygen in remote position.