Solvent-mediated tautomerization of purine: single to quadruple proton transfer

We present calculations for the structures and the tautomerization reaction of purine and purine – (H₂O)_n (n=1–3) clusters. We find two pathways (via the carbene and the sp³-type intermediate) for the 9 ↔ 7 tautomerization of bare purine. The barrier heights for the 9 → 3 and 9 → 7 tautomerization of bare purine are calculated to be large (60–70 kcal/mol). Hydrogen bonding with the water molecule(s), however, dramatically lowers the 9 → 3 barrier by the concerted multiple proton transfer mechanism, favoring the formation of the conformer 3(H)- relative to the 7(H)-purine in the microsolvated environment, in contrast to the gas phase or the aqueous solution.     

A Theoretical Investigation on N7(H)→N9(H)Tautomerization of Isolated and Monohydrated 2,6-Dithiopurine Combined with IPCM Solvent Model

The tautomerization reaction mechanism has been reported between N7(H) and N9(H) of isolated and monohydrated 2,6‐dithiopurine using B3LYP/6‐311+G(d,p). The isodensity polarized continuum model (IPCM) in the self‐consistent reaction field (SCRF) method is employed to account for the solvent effect of water on the tautomerization reaction activation energies. The results show that the two pathways P(1) (via the carbene intermediate I1) and P(2) (via the sp³‐hybrid intermediate I2) are found in intramolecular proton transfer, and each pathway is composed by two primary steps. The calculated activation energy barriers of the rate‐determining steps in isolated 2,6‐dithiopurine N7(H)→N9(H) tautomerism are 308.2 and 220.0 kJ·mol^?1 in the two pathways, respectively. Interestingly, in one‐water molecule catalyst, it dramatically lowers the N7(H)→N9(H) energy barriers by the concerted double proton transfer mechanism in P(1), favoring the formation of 2,6‐dithiopurine N9(H). However, the single proton transfer mechanism assisted with out‐of‐plane water in the first step of P(2) increases the activation energy barrier from 220.0 to 232.3 kJ·mol^?1, while the second step is the out‐of‐plane concerted double proton transfer mechanism, indicating that they will be less preferable for proton transfer. Additionally, the results also show that all the pathways are put into the aqueous solution, and the activation energy barriers have no significant changes. Therefore, the long‐range electrostatic effect of bulk solvent has no significant impact on proton transfer reactions and the interaction with explicit water molecules will significantly influence proton transfer reactions.     

Excited-state proton transfer through water bridges and structure of hydrogen-bonded complexes in 1H-pyrrolo[3,2-h]quinoline: adiabatic time-dependent density functional theory study

Proton transfer reaction is studied for 1H-pyrrolo[3,2-h]quinoline-water complexes (PQ-(H(2)O)(n), n = 0-2) in the ground and the lowest excited singlet states at the density functional theory (DFT) level. Cyclic hydrogen-bonded complexes are considered, in which water molecules form a bridge connecting the proton donor (pyrrole NH group) and acceptor (quinoline nitrogen) atoms. To understand the effect of the structure and length of water bridges on the excited-state tautomerization in PQ, the potential energy profile of the lowest excited singlet state is calculated adiabatically by the time-dependent DFT (TDDFT) method. The S(0) --> S(1) excitation of PQ is accompanied by significant intramolecular transfer of electron density from the pyrrole ring to the quinoline fragment, so that the acidity of the N-H group and the basicity of the nitrogen atom of the quinoline moiety are increased. These excited-state acid-base changes introduce a driving force for the proton transfer reaction. The adiabatic TDDFT calculations demonstrate, however, that the phototautomerization requires a large activation energy in the isolated PQ molecule due to a high energy barrier separating the normal form and the tautomer. In the 1:1 cyclic PQ-H(2)O complex, the energy barrier is dramatically reduced, so that upon excitation of this complex the tautomerization can occur rapidly in one step as concerted asynchronous movements of the two protons assisted by the water molecule. In the PQ-(H(2)O)(2) solvate two water molecules form a cyclic bridge with sterically strained and unfavorable hydrogen bonds. As a result, some extra activation energy is needed for initiating the proton dislocation along the longer hydrogen-bond network. The full tautomerization in this complex is still possible; however, the cooperative proton transfer is found to be highly asynchronous. Large relaxation and reorganization of the hydrogen-bonded water bridge in PQ-(H(2)O)(2) are required during the proton translocation from the pyrrole NH group to the quinoline nitrogen; this may block the complete tautomerization in this type of solvate.     

Ab initio study of gas phase and water-assisted tautomerization of maleimide and formamide

Maleimide serves as an important starting material in the synthesis of drugs and enzyme inhibitors. In the present paper, knowing the importance of tautomerization in maleimide for its drug action, potential energy surface of maleimide is studied and its tautomerization has been discussed and compared with tautomerization of formamide. Gas phase tautomerization of maleimide requires large amount of energy (23·21 kcal/mol) in comparison to formamide (15·05 kcal/mol) at HF/6-31+G* level. Thus making the proton transfer reaction a difficult process in gas phase. Water molecule lowers the energy barrier of tautomerization thus facilitating the tautomerization of maleimide to 5-hydroxy-pyrrol-2-one. Water assisted tautomerization of maleimide requires 19·60 kcal/mol energy at HF/6-31+G* and 17·63 kcal/mol energy at B3LYP/6-31+G* level, a decrease of 3·61 and 5·96 kcal/mol over gas phase tautomerization. Whereas, tautomerization of formamide requires 14·16 and 12·84 kcal/mol energy, a decrease of 0·89 and 2·01 kcal/mol energy over gas phase tautomerization at HF/6-31+G* and B3LYP/6-31+G* level, respectively. Water-assisted tautomerization in maleimide and formamide showed that difference in energy barrier reduces to 2·83 kcal/mol from 10·41 kcal/mol (in gas phase) at B3LYP level, which resulted that maleimide readily undergoes tautomerization in water molecule.     

Reaction paths of tautomerization between hydroxypyridines and pyridones

Tautomerization paths of 2(and 4)-hydroxypyridine (called here HP) to 2(and 4)-pyridone (called here PY) with water molecules were investigated by the use of density functional theory calculations. Potential energies were compared for a number of water molecules. The 2-HP molecule was found to be isomerized most readily and concertedly to the 2-PY one via proton relays with two water molecules. The reaction pattern is invariant even when outer water molecules are added. The 4-HP(H(2)O)(n) --> 4-PY(H(2)O)(n) reaction model did not give small activation energies. However, a reaction of (4-HP)(2)(H(2)O)(2) --> (4-PY)(2)(H(2)O)(2) was found to occur readily through a transient ion-pair intermediate. The conversion processes of (2-PY)(2) to the tautomerization reacting system were discussed. The hydrogen-bond directionality regulates the tautomerization paths.     

Density functional study on the reaction mechanism of proton transfer in 2-pyridone: effect of hydration and self-association

The proton-transfer mechanism in the isolated, mono, dehydrated forms and dimers of 2-pyridone and the effect of hydration or self-assistance on the transition state structures corresponding to proton transfer from the keto form to the enol form have been investigated using B3LYP and BH-LYP hybrid density functional methods at the 6-311++G (2d, 2p) basis set level. The barrier heights for both H2O-assisted and self-assisted reactions are significantly lower than that of the bare tautomerization reaction from 2-pyridone to 2-hydroxypyridine, implying the importance of the superior catalytic effect of H2O and (H2O)2 and the important role of 2-pyridone itself for the intramolecular proton transfer. Long-range solvent effects have also been taken into account by using the continuum model (Onsager model and polarizable continuum model (PCM)) of water. The tautomerization energies and the potential energy barriers are increased both for the water-assisted and for the self-assisted reaction because of the bulk solvent, which imply that the tautomerization of PY becomes less favorable in the polar solvent.     

Kinetic, DFT and TD-DFT studies on the mechanism of stabilization of pyramidal H3PO3 at the [Mo3M'S4(H2O)10]4+ clusters (M' = Pd, Ni)

The kinetics of reaction between the [Mo(3)M'S(4)(H(2)O)(10)](4+) clusters (M' = Pd, Ni) and H(3)PO(3) has been studied in 4.0 M Hpts/Lipts (pts(-) = p-toluenesulfonate). For both complexes there is an initial kinetic step with small absorbance changes that corresponds to substitution of the water coordinated to Pd by a molecule of tetrahedral H(3)PO(3). For the Pd complex, tautomerization of H(3)PO(3) occurs in a slower kinetic step with much larger absorbance changes; it leads to formation of [Mo(3)Pd(pyr-H(3)PO(3))S(4)(H(2)O)(9)](4+) in which H(3)PO(3) adopts a pyramidal structure, but the process is not as favored as for H(3)PO(2). The kinetics of this second step is independent of the concentration of H(3)PO(3) but dependent on the concentration of Hpts on the supporting electrolyte. For the Ni complex, the second step is severely hindered and its kinetics could not be studied. DFT calculations indicate that tautomerization of H(3)PO(3) is expected to be less favoured than that of H(3)PO(2), both processes being less favored at the Ni cluster than at its Pd analogue. With regard to the tautomerization mechanism, the calculations indicate that the mechanism previously proposed for H(3)PO(2) can be the same for H(3)PO(3), although the initial H-shift can also occur through a protonation-deprotonation sequence with participation of external protons instead of a second molecule of the phosphorus acid. TD-DFT studies have been also carried out to understand the similarity between the spectra of the starting complex and the reaction intermediate formed in the first kinetic step as well as the large spectral changes associated to the tautomerization process. Although it contains contributions from several transitions, the intense band observed for clusters containing coordinated pyr-H(3)PO(3) involves essentially ligand-to-metal charge-transfer (LMCT) transitions.     

Can the ebselen derivatives catalyze the isomerization of peroxynitrite to nitrate?

The reaction of ebselen and its derivatives (1-7) with peroxynitrite anion (ONOO(-); PN) has been studied in gas phase and in aqueous, dichloromethane, benzene, and cyclohexane solutions using B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) and PCM-B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) approaches, respectively. It was shown that the reaction of 2 (R=H) with PN proceeds via 2 + PN --> 2-PN --> 2-TS1 (O-O activation) --> 2-O(NO(2)(-)()) --> 2-SeO + NO(2)(-) pathway with a rate-determining barrier of 25.3 (14.8) kcal/mol at the NO(2)(-) dissociation step (numbers presented without parentheses are enthalpies, and those in parentheses are Gibbs free energies). The NO(3)(-) formation process, starting from the complex 2-O(NO(2)(-)()), requires by (7.9) kcal/mol more energy than the NO(2)(-) dissociation process and is unlikely to compete with the latter. Thus, in the gas phase, the peroxynitrite --> nitrate isomerization catalyzed by complex 2 is unlikely to occur. It is shown that the NO(3)(-) formation process is slightly more favorably than the NO(2)(-) dissociation process for complex 4, with a strongest electron-withdrawing ligand R=CF(3). Therefore, complex 4 (as well as complex 6 with R=OH) is predicted to be a good catalyst for peroxynitrite <--> nitrite isomerization in the gas phase. Solvent effects (a) change the rate-determining step of the reaction 2 + PN from NO(2)(-) dissociation in the gas phase to O-O activation, which occurs with barriers of (13.9), (8.4), (8.4), and (8.2) kcal/mol in water, dichloromethane, benzene, and cyclohexane, respectively, and (b) significantly reduce the NO(2)(-) dissociation energy, while only slightly destabilizing the NO(3)(-) formation barrier, and make the peroxynitrite <--> nitrate isomerization process practically impossible, even for complex 4.     

Theoretical Study on the Tautomerization Reaction of Formamide to Formamidic Acid Encapsulated in a Nano-sized Molecular Container

The hybrid calculations with ONIOM(B3LYP/6-31G*:AM1) method were carried out on the tautomerization reaction of formamide to formamidic acid in the microcontainer-encapsulated state. The free-state tautomerization process was also investigated with B3LYP/6-31G**//B3LYP/6-31G* method for the purpose of comparison. Bare tautomerization, H2O-assisted(single-H2O or multiple-H2O) and self-assisted mechanisms were all taken into consideration for the encapsulated state. The results show that only bare tautomerization and single-H2O catalysis mechanisms are possible to the encapsulated for-mamide tautomerization owing to the container's size confinement effect. Geometrical changes in the complexed container and guest molecules are discussed to deeply understand the complex's structural properties. The bare tautomerization barrier in the encapsulated state increases by 23.826 kJ/mol, ac-counting for 12% of the corresponding total energy barrier in the free state, and the increased values for the single-H2O catalysis process are 12.958 kJ/mol, accounting for 16% of the corresponding total energy barrier, respectively. This finding suggests that the encapsulation can make the tautomerization process slightly difficult.     

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.     

Theoretical exploration of hydrogen loss from Al3H9

The Al(3)H(9) and Al(3)H(7) potential energy surfaces were explored using quantum chemistry calculations to investigate the H(2) loss mechanism from Al(3)H(9), which provide new insights into hydrogen production from bulk alane, [AlH(3)](x), a possible energy storage material. We present results of B3LYP/6-311++G(d,p) calculations for the various Al(3)H(9) and Al(3)H(7) optimized local minima and transition state structures along with some reaction pathways for their interconversion. We find the energy for Al(3)H(9) decomposition into Al(2)H(6) and AlH(3) is slightly lower than that for H(2) loss and Al(3)H(7) formation, but the calculations show that H(2) loss from Al(3)H(9) is a lower energy process than for losing hydrogen from either Al(2)H(6) or AlH(3). We found four transition state structures and reaction pathways for Al(3)H(9) → Al(3)H(7) + H(2), where the lowest energy activation barrier is around 25-73 kJ/mol greater than the experimental value for H(2) loss from bulk alane. Intrinsic reaction coordinate calculations show that the H(2) loss pathway involves considerable rearrangement of the H atom positions around a single Al center. Three of the pathways start with the formation of an AlH(3) moiety, which then enables a terminal H on the AlH(3) to get within 1.1 to 1.2 ? of a nearby bridging H atom. The bridging and terminal H atoms eventually combine to form H(2) and leave Al(3)H(9). One implication of these H(2) loss reaction pathways is that, since the H atoms in bulk alanes are all at bridging positions, if a similar H(2) loss mechanism were to apply to bulk alane, then H(2) loss would most likely occur on the bulk alane surface or at a defect site where there should be more terminal H atoms available for reaction with nearby bridging H atoms.     

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.     

Gas-phase reactions of organic radicals and diradicals with ions

Reactions of polyatomic organic radicals with gas phase ions have been studied at thermal energy using a flowing afterglow-selected ion flow tube (FA-SIFT) instrument. A supersonic pyrolysis nozzle produces allyl radical (CH2CHCH2) and ortho-benzyne diradical (o-C6H4) for reaction with ions. We have observed: [CH2CHCH2 + H3O+ --> C3H6+ + H2O], [CH2CHCH2 + HO- --> no ion products], [o-C6H4 + H3O+ --> C6H5+ + H2O], and [o-C6H4 + HO- --> C6H3- + H2O]. The proton transfer reactions with H3O+ occur at nearly every collision (kII approximately with 10(-9) cm3 s(-1)). The exothermic proton abstraction for o-C6H4 + HO- is unexpectedly slow (kII approximately with 10(-10) cm3 s(-1)). This has been rationalized by competing associative detachment: o-C6H4 + HO- --> C6H5O + e-. The allyl + HO- reaction proceeds presumably via similar detachment pathways.     

Photodissociation and photoisomerization of alpha-fluorotoluene and 4-fluorotoluene in a molecular beam

The photodissociation of jet-cooled alpha-fluorotoluene and 4-fluorotoluene at 193 and 248 nm was studied using vacuum ultraviolet (vuv) photoionization/multimass ion imaging techniques as well as electron impact ionization/photofragment translational spectroscopy. Four dissociation channels were observed for alpha-fluorotoluene at both 193 and 248 nm, including two major channels C6H5CH2F-->C6H5CH2 (or C7H7)+F and C6H5CH2F-->C6H5CH (or C7H6)+HF and two minor channels C6H5CH2F-->C6H5CHF+H and C6H5CH2F-->C6H5+CH2F. The vuv wavelength dependence of the C7H7 fragment photoionization spectra indicates that at least part of the F atom elimination channel results from the isomerization of alpha-fluorotoluene to a seven-membered ring prior to dissociation. Dissociation channels of 4-fluorotoluene at 193 nm include two major channels C6H4FCH3-->C6H4FCH2+H and C6H4FCH3-->C6H4F+CH3 and two minor channels C6H4FCH3-->C6H5CH2 (or C7H7)+F and C6H4FCH3-->C6H5CH (or C7H6)+HF. The dissociation rates for alpha-fluorotoluene at 193 and 248 nm are 3.3 x 10(7) and 5.6 x 10(5) s(-1), respectively. The dissociation rate for 4-fluorotoluene at 193 nm is 1.0 x 10(6) s(-1). An ab initio calculation demonstrates that the barrier height for isomerization from alpha-fluorotoluene to a seven-membered ring isomer is much lower than that from 4-fluorotoluene to a seven-membered ring isomer. The experimental observed differences of dissociation rates and relative branching ratios between alpha-fluorotoluene and 4-fluorotoluene may be explained by the differences in the six-membered ring to seven-membered ring isomerization barrier heights, F atom elimination threshold, and HF elimination threshold between alpha-fluorotoluene and 4-fluorotoluene.     

(H_2O)_n(n=0～2)催化HS和HOCl反应机理的理论研究

在CCSD(T)/6-311+G(3df,2p)//M06-2X/6-311+G(3df,2p)水平上研究了(H_2O)n(n=0~2)催化HS和HOCl的反应机理.结果表明,HS与HOCl反应中HS夺取HOCl上的H原子形成产物H_2S和ClO.在无水催化时,该反应存在2种不同的路径(分别经过过渡态TS1和TS2,二者互为顺反结构),对应的能垒分别为100.28和100.91kJ/mol,到达产物(H_2S+ClO)需吸收18.99kJ/mol能量,反应不易发生;在单个水分子参与时,水分子可通过形成弱相互作用或者作为H原子转移桥梁影响反应机理,获得了4种水催化路径,能垒(间于53.97~92.39kJ/mol之间)均低于无水催化过程.同时发现,在反应到达产物前,水分子可以与产物形成中间体IM,IM相对能仅为0.46kJ/mol,有利于产物形成;有2个水分子参与反应时,找到了3条催化路径,最优反应路径过渡态TS7的能垒为45.05kJ/mol,低于无水催化过程,相比单个水分子最优路径能垒(53.97kJ/mol)并无显著降低.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

A computational investigation of the sulphuric acid‐catalysed 1,4‐hydrogen transfer in higher Criegee intermediates

Criegee intermediates (CIs) are formed during the ozonolysis of unsaturated hydrocarbons in the troposphere. The fate of CIs is of critical importance to tropospheric oxidation chemistry, particularly in the context of radical and secondary organic aerosol formation. Using the high‐level ab initio G4(MP2) method, we investigate the 1,4 hydrogen shift reaction in CIs formed from ozonolysis of two common biogenic hydrocarbons: isoprene and α‐pinene. We consider the uncatalysed reaction, as well as the reaction catalysed by a water molecule and by sulphuric acid. We show that sulphuric acid is a very effective catalyst, leading to a barrierless tautomerization relative to the free reactants and to very low reaction barrier heights relative to the reactant complexes. In particular, we obtain reaction barrier heights of Δ = 24.5 (isoprene CI) and 8.4 (α‐pinene CI) kJ mol⁻¹ relative to the reactant complexes. Given the reaction of OH radicals with SO₂ in the troposphere can ultimately yield sulphuric acid, these findings may have significant consequences for current atmospheric chemical models for regions of high sulphur concentrations.     

Dimers of and tautomerism between 2-pyrimidinethiol and 2(1H)-pyrimidinethione: a density functional theory (DFT) study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.     

General rules for predicting where a catalytic reaction should occur on metal surfaces: a density functional theory study of C-H and C-O bond breaking/making on flat,stepped, and kinked metal surfaces

To predict where a catalytic reaction should occur is a fundamental issue scientifically. Technologically, it is also important because it can facilitate the catalyst's design. However, to date, the understanding of this issue is rather limited. In this work, two types of reactions, CH(4) <--> CH(3) + H and CO <--> C + O on two transition metal surfaces, were chosen as model systems aiming to address in general where a catalytic reaction should occur. The dissociations of CH(4) --> CH(3) + H and CO --> C + O and their reverse reactions on flat, stepped, and kinked Rh and Pd surfaces were studied in detail. We find the following: First, for the CH(4) <--> Ch(3) + H reaction, the dissociation barrier is reduced by approximately 0.3 eV on steps and kinks as compared to that on flat surfaces. On the other hand, there is essentially no difference in barrier for the association reaction of CH(3) + H on the flat surfaces and the defects. Second, for the CO <--> C + O reaction, the dissociation barrier decreases dramatically (more than 0.8 eV on Rh and Pd) on steps and kinks as compared to that on flat surfaces. In contrast to the CH(3) + H reaction, the C + O association reaction also preferentially occurs on steps and kinks. We also present a detailed analysis of the reaction barriers in which each barrier is decomposed quantitatively into a local electronic effect and a geometrical effect. Our DFT calculations show that surface defects such as steps and kinks can largely facilitate bond breaking, while whether the surface defects could promote bond formation depends on the individual reaction as well as the particular metal. The physical origin of these trends is identified and discussed. On the basis of our results, we arrive at some simple rules with respect to where a reaction should occur: (i) defects such as steps are always favored for dissociation reactions as compared to flat surfaces; and (ii) the reaction site of the association reactions is largely related to the magnitude of the bonding competition effect, which is determined by the reactant and metal valency. Reactions with high valency reactants are more likely to occur on defects (more structure-sensitive), as compared to reactions with low valency reactants. Moreover, the reactions on late transition metals are more likely to proceed on defects than those on the early transition metals.     

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.