Heat of hydrogenation of 1,5-dehydroquadricyclane. A computational and experimental study of a highly pyramidalized alkene

The radical anion of the highly pyramidalized alkene 1,5-dehydroquadricyclane (1) was generated in the gas phase from the Squires reaction of 1,5-bis(trimethylsilyl)quadricyclane with F-/F2. The electron binding energy and proton affinity of 1*- were determined by bracketing experiments to be 0.6 +/- 0.1 eV and 386 +/- 5 kcal/mol, respectively. These values are in good agreement with values predicted by density functional theory (B3LYP/6-31+G*) and ab initio (CASPT2/6-31+G*) calculations. The experimental heat of hydrogenation of 1, obtained from a thermochemical cycle, was found to be 91 +/- 9 kcal/mol. This value of deltaH(H2) leads to values of 67 +/- 9 kcal/mol for the olefin strain energy (OSE) of 1, 172 +/- 9 kcal/mol for its heat of formation, and 23 +/- 9 kcal/mol for its pi bond dissociation enthalpy. Since the retro-Diels-Alder reaction of neutral 1 is computed to be highly exothermic, the finding that 1*- apparently does not undergo a retro-Diels-Alder reaction is of particular interest. The B3LYP/6-31+G* optimized geometry of 1 suggests that the bonding in this alkene is partially delocalized, presumably because the highly pyramidalized double bond in 1 interacts with the distal cyclopropane bonds in a manner that eventually leads to a retro-Diels-Alder reaction. The good agreement of the B3LYP and (2/2)CASPT2 values for the heat of hydrogenation and OSE of 1 with the experimentally derived values provides indirect evidence for the correctness of the B3LYP prediction that the equilibrium geometry of 1 lies part way along the reaction coordinate to the transition structure for the retro-Diels-Alder reaction.     

Ab initio and experimental studies on the hetero-Diels-Alder and cheletropic additions of sulfur dioxide to (E)-1-methoxybutadiene: a mechanism involving three molecules of SO(2)

Kinetics on the cheletropic addition of sulfur dioxide to (E)-1-methoxybutadiene (1) to give the corresponding sulfolene 2 (2-methoxy-2,5-dihydrothiophene-1,1-dioxide) gave the rate law d[2]/dt = k[1][SO(2)](x)() with x = 2.6 +/- 0.2 at 198 K. Under these conditions, no sultine 3 [(2RS,6RS)-6-methoxy-3,6-dihydro-1,2-oxathiin-2-oxide] resulting from a hetero-Diels-Alder addition was observed, and the cheletropic elimination 2 --> 1 + SO(2) did not occur. Ab initio and DFT quantum calculations confirmed that the cheletropic addition 1 + SO(2) --> 2 follows two parallel mechanisms, one involving two molecules of SO(2) and the transition structure with DeltaG(++) = 18.2 +/- 0.2 kcal/mol at 198 K (exptl); 22.5-22.7 kcal/mol [B3LYP/6-31G(d,p)], the other one involving three molecules of SO(2) with DeltaG(++) = 18.9 +/- 0.1 kcal/mol at 198 K (exptl); 19.7 kcal/mol [B3LYP/6-31G(d,p)]. The mechanism involving only one molecule of SO(2) in the transition structure requires a higher activation energy, DeltaG(++) = 25.2 kcal/mol [B3LYP/6-31G(d,p)]. Comparison of the geometries and energetics of the structures involved into the 1 + SO(2) --> 2, 3 and 1 + 2SO(2) --> 2, 3 + SO(2) reactions obtained by ab initio and DFT methods suggest that the latter calculation techniques can be used to study the cycloadditions of sulfur dioxide. The calculations predict that the hetero-Diels-Alder addition 1 + SO(2) --> 3 also prefers a mechanism in which three molecules of SO(2) are involved in the cycloaddition transition structure. At 198 K and in SO(2) solutions, the entropy cost (TDeltaS(++)) is overcompensated by the specific solvation by SO(2) in the transition structures of both the cheletropic and hetero-Diels-Alder reactions of (E)-1-methoxybutadiene with SO(2).     

Theoretical proton affinity and fluoride affinity of nerve agent VX

Proton affinity and fluoride affinity of nerve agent VX at all of its possible sites were calculated at the RI-MP2/cc-pVTZ//B3LYP/6-31G* and RI-MP2/aug-cc-pVTZ//B3LYP/6-31+G* levels, respectively. The protonation leads to various unique structures, with H(+) attached to oxygen, nitrogen, and sulfur atoms; among which the nitrogen site possesses the highest proton affinity of -ΔE ～ 251 kcal/mol, suggesting that this is likely to be the major product. In addition some H(2), CH(4) dissociation as well as destruction channels have been found, among which the CH(4) + [Et-O-P(═O)(Me)-S-(CH(2))(2)-N(+)(iPr)═CHMe] product and the destruction product forming Et-O-P(═O)(Me)-SMe + CH(2)═N(+)(iPr)(2) are only 9 kcal/mol less stable than the most stable N-protonated product. For fluoridization, the S-P destruction channel to give Et-O-P(═O)(Me)(F) + [S-(CH(2))(2)-N-(iPr)(2)](-) is energetically the most favorable, with a fluoride affinity of -ΔE ～ 44 kcal. Various F(-) ion-molecule complexes are also found, with the one having F(-) interacting with two hydrogen atoms in different alkyl groups to be only 9 kcal/mol higher than the above destruction product. These results suggest VX behaves quite differently from surrogate systems.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Theoretical study of three predominant tautomers of 2-oxo-6-methylpurine and their two transition state structures

The three predominantly stable tautomers of 2-oxo-6-methylpurine were studied in gas phase and aqueous solution by means of quantum mechanical calculations. Two transition state structures connecting these three tautomeric forms on the free energy surface were determined. The activation free energy for the intramolecular proton transfer in gas phase was calculated to be considerably smaller than the bond energy of either N-H or O-H: 59.01 and 30.37 kcal/mol for N9→N3 and N1→O2, respectively, obtained at the QCISD(T)/6-31G+(d)//MP2(full)/6-31G(d) level of theory.     

Tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime: ab initio and NMR study

The tautomeric and conformational equilibrium of 2-nitrosophenol and 9,10-phenanthrenequinonemonooxime was studied by ab initio methods. The geometry optimizations of the structures investigated were done without any geometrical restrictions at HF/6-31G** and MP2/6-31G** levels of theory. The transition structures for tautomeric and rotameric conversions were located. To correct for electron correlation, single-point calculations were carried out up to MP4/6-311G*//MP2/6-31G* level of theory.

Ab initio calculations for 2-nitrosophenol in agreement with the available experimental data define the nitroso form as more stable. It was found that the influence of the correlation energy on the relative stabilities is smaller for the rotamers of the nitroso tautomer but substantially (4–6 kcal/mol) for the oxime forms. It was found that the barrier height of tautomerization reaction is 10.24 kcal/mol.

The structure of the 9,10-phenanthrenequinonemonooxime was studied by solid and liquid state NMR spectroscopy. Ab initio calculations in agreement with our experimental data predict that the compound exists as oxime tautomer and the syn-oxime is most stable. It was found that the solvent influence on the relative stabilities of both isomers: syn- and anti-oxime. While in chloroform solution the syn-oxime is preferred but in DMSO anti-oxime is more stable in energy.

At the MP4/6-311G*//MP2/6-31G**+ZPE level of theory the barrier of tautomerization was predicted to be 10.96 kcal/mol and the rotational barrier around the single C–O bond in the syn-oxime was found to be 7.57 kcal/mol. The rotation is facile and this explains the absence of nitroso tautomers in solution.     

Ab Initio Calculations of Structure and Energetics of the Hypermetalated Hydrogen Oxide Molecules: Li2NaOH and LiNa2OH

We made ab initio electronic calculations of the structure and energetics of mixed hypermetalated hydrogen oxides, Li₂NaOH and LiNa₂OH. There exist five equilibrium geometries for each complex. In all levels of calculation the global minimum structure for Li₂NaOH has C_2v symmetry and a large distance between sodium and oxygen, 4.24 Å (MP2/6-31G*). The dissociation energies to all possible products were also calculated. Li₂NaOH → Na + Li₂OH δH = +25.33 kcal/mol (at MP4/6-311++G**//6-31G* + ZPE scaled by 0.9). All other dissociation processes are highly endothermic. Similar procedures were applied to LiNa₂OH. The global minimum structure for LiNa₂OH belongs to point group C_s. It is also endothermic to all possible dissociation paths. LiNa₂OH →Na + LiNaOH δH = +12.72 kcal/mol (at MP4/6-311++G*//6-31G* + ZPE scaled by 0.9). The nuclear repulsion energy is crucial in energetics of the structures. The distribution of electron density and bonding properties for these equilibrium structures were analyzed.     

Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

Force field parameters for sulfates and sulfamates based on ab initio calculations: Extensions of AMBER and CHARMm fields

Ab initio self-consistent field (SCF) Hartree-Fock calculations of sulfates R? O? SO₃(?1) (R = Me, Et, i-Pr) and sulfamates R? NHSO₃(?1) (R = H, Me, Et, i-Pr) were performed at the 4-31G(*S*N) //3-21G(*S*N) basis set levels, where asterisks indicate d functions on sulfur and nitrogen atoms. These standard levels were determined by comparing calculation results with several basis sets up to MP2/6-31G*//6-31G*. Several conformations per compound were studied to obtain molecular geometries, rotational barriers, and potential derived point charges. In methyl sulfate, the rotational barrier around the C? O bond is 1.6 kcal/mol at the MP2 level and 1.4 kcal/mol at the standard level. Its ground state has one of three HCOS torsion angles trans and one of three COSO torsion angles trans. Rotation over 60° around the single O? S bond in the sulfate group costs 2.5 kcal/mol at the MP2 and 2.1 kcal/mol at the standard level. For ethyl sulfate, the calculated rotational barrier in going from the ground state, which has its CCOS torsion angle trans, to the syn-periplanar conformation (CCOS torsion angle cis) is 4.8 kcal/mol. However, a much lower barrier of 0.7 kcal/mol leads to a secondary gauchelike conformation about 0.4 kcal/mol above the ground state, with the CCOS torsion angle at 87.6°. Again, one of the COSO torsion angles is trans in the ground state, and the rotational barrier for a 60° rotation of the sulfate group amounts to 1.8 kcal/mol. For methyl sulfamate, the rotational barriers are 2.5 kcal/mol around the C? N bond and 3.3 kcal/mol around the N? S bond. This is noteworthy because sulfamate itself has a calculated rotational barrier around the N? S bond of only 1.7 kcal/mol. These and other data were used to parameterize the well-known empirical force fields AMBER and CHARMm. When the new fields were tested by means of vibrational frequency calculations at the 6-31G*//6-31G* level for methyl sulfate, sulfamate, and methyl sulfamate ground states, the frequencies compared favorably with the AMBER and CHARMm calculated frequencies. The transferability of the force parameters to β-D -glucose-6-sulfate and isopropyl sulfate appears to be better than to isopropyl sulfamate. © 1995 by John Wiley & Sons, Inc.     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Ab initio studies on Al(+)(H(2)O)(n), HAlOH(+)(H(2)O)(n-1), and the size-dependent H(2) elimination reaction

We report computational studies on Al(+)(H(2)O)(n), and HAlOH(+)(H(2)O)(n-1), n = 6-14, by the density functional theory based ab initio molecular dynamics method, employing a planewave basis set with pseudopotentials, and also by conventional methods with Gaussian basis sets. The mechanism for the intracluster H(2) elimination reaction is explored. First, a new size-dependent insertion reaction for the transformation of Al(+)(H(2)O)(n), into HAlOH(+)(H(2)O)(n-1) is discovered for n > or = 8. This is because of the presence of a fairly stable six-water-ring structure in Al(+)(H(2)O)(n) with 12 members, including the Al(+). This structure promotes acidic dissociation and, for n > or = 8, leads to the insertion reaction. Gaussian based BPW91 and MP2 calculations with 6-31G* and 6-31G** basis sets confirmed the existence of such structures and located the transition structures for the insertion reaction. The calculated transition barrier is 10.0 kcal/mol for n = 9 and 7.1 kcal/mol for n = 8 at the MP2/6-31G** level, with zero-point energy corrections. Second, the experimentally observed size-dependent H(2) elimination reaction is related to the conformation of HAlOH(+)(H(2)O)(n-1), instead of Al(+)(H(2)O)(n). As n increases from 6 to 14, the structure of the HAlOH(+)(H(2)O)(n-1) cluster changes into a caged structure, with the Al-H bond buried inside, and protons produced in acidic dissociation could then travel through the H(2)O network to the vicinity of the Al-H bond and react with the hydride H to produce H(2). The structural transformation is completed at n = 13, coincident approximately with the onset of the H(2) elimination reaction. From constrained ab initio MD simulations, we estimated the free energy barrier for the H(2) elimination reaction to be 0.7 eV (16 kcal/mol) at n = 13, 1.5 eV (35 kcal/mol) at n = 12, and 4.5 eV (100 kcal/mol) at n = 8. The existence of transition structures for the H(2) elimination has also been verified by ab initio calculations at the MP2/6-31G** level. Finally, the switch-off of the H(2) elimination for n > 24 is explored and attributed to the diffusion of protons through enlarged hydrogen bonded H(2)O networks, which reduces the probability of finding a proton near the Al-H bond.     

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.     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

Molecule-induced alkane homolysis with dioxiranes.

The mechanisms of C-H and C-C bond activations with dimethyldioxirane (DMD) were studied experimentally and computationally at the B3LYP/6-311+G**//B3LYP/6-31G* density functional theory level for the propellanes 3,6-dehydrohomoadamantane (2) and 1,3-dehydroadamantane (3). The sigma(C-C) activation of 3 with DMD (Delta G(*) = 23.9 kcal mol(-1) and Delta G(r) = -5.4 kcal mol(-1)) is the first example of a molecule-induced homolytic C-C bond cleavage. The C-H bond hydroxylation observed for 2 is highly exergonic (Delta G(r) = -74.4 kcal mol(-1)) and follows a concerted pathway (Delta G(*) = 34.8 kcal mol(-1)), in contrast to its endergonic molecule-induced homolysis (Delta G(*) = 28.8 kcal mol(-1) and Delta G(r) = +9.2 kcal mol(-1)). The reactivities of 2 and 3 with CrO(2)Cl(2), which follow a molecule-induced homolytic activation mechanism, parallel the DMD results only for highly reactive 3, but differ considerably for more stable propellanes such as 4-phenyl-3,6-dehydrohomoadamantane (1) and 2.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Aspirin. An ab initio quantum-mechanical study of conformational preferences and of neighboring group interactions.

The potential energy surface of acetylsalicylic acid, aspirin, has been explored at the RHF/6-31G* and B3LYP/6-31G* levels, and single-point calculations were performed at levels up to B3LYP/6-311G**//B3LYP/6-31G*. All conformational isomers have been located, the thermochemical functions have been computed, and relative energies and free enthalpies were determined. The conformational space of aspirin is spanned by three internal coordinates, and these are the carboxylic acid C-O conformation (s-trans preferred by about 7 kcal/mol), the C-COOH conformation (Z preferred unless there are H-bonding opportunities), and the ester C-O conformation (s-trans preferred by about 4 kcal/mol). There are nine aspirin isomers since one of the conformers realizes hydrogen-bonding structure isomerism as well. Neighboring group interactions are discussed with reference to the intrinsic properties of benzoic acid and phenyl acetate. The intrinsic conformational preference energies for benzoic acid and phenyl acetate are not additive. The acid s-trans preference energies differ by as much as 9 kcal/mol depending on the Ph-COOH and ester conformations. Similarly, the E-preference energies about the Ph-COOH bond vary by as much as 6 kcal/mol depending on the ester conformation. The structural discussion suggests an overall ortho repulsion between the functional groups in all aspirin isomers including the intramolecularly hydrogen-bonded isomers. The isodesmic reaction between the most stable conformers of benzoic acid and phenyl acetate to form aspirin and benzene is found to be endothermic by 2.7 kcal/mol and provides compelling evidence for and a quantitative measure of ortho repulsion. The ortho repulsion of 2.7 kcal/mol is a lower limit, and the ortho repulsion can increase to as much as 6 kcal/mol in some aspirin isomers.     

Theoretical study of the reactions of the 1Σ+ ground state of MS+ (M = Sc,Y, and La) with oxygen‐transfer reagent MS+ + CO → ScO+ + CS in the gas phase

The reaction mechanisms of the ¹Σ⁺ ground state of MS⁺ (M = Sc, Y, and La) with oxygen‐transfer reagent MS⁺ + CO → MO⁺ + CS in the gas phase has been proposed and investigated by ab initio methods with the 6‐31G* basis set for nonmetal atoms and the effective core potentials of Lanl2dz for the metal atoms. A carbon migration from oxygen atom to sulfur atom via a four‐center transition state is involved on the reaction potential surface. The activation energies of the reactions are 34.0, 24.1, and 36.7 kcal/mol relative to their corresponding reactants and the reaction heats are 15.7, 18.6, and 18.0 kcal/mol (respectively, for M = Sc, Y, and La) at the MP4 (SDTQ)/6‐31G*//MP2/6‐31G* level plus zero‐point energy, which indicates that the cationic yttrium sulfide is more favorable for this type of reaction. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Multiple bond migration with participation of a protophilic agent

The pathways of migration of the double bond in the 1-methoxy-2-propene molecule with participation of a hydroxide ion were investigated by theab initio RHF/6-31+G^* and MP2/6-31+G^* methods. Stationary points corresponding to complexes between the molecule under study and the hydroxide ion and between the corresponding carbanion and a water molecule were found on the potential energy surfaces of the proton transfer reactions. As in the propene molecule, migration of the double bond in the 1-methoxy-2-propene molecule can occur in the gas phase by the mechanism of intramolecular proton transfer involving the proton of the hydroxide ion. Conformational isomerism of the initial molecule and reaction products was considered. The proposed mechanism of 1,3-hydrogen shift involving the proton-containing base suggests the formation of bothE- andZ-products with predominance of the latter irrespective of the nature of the substituent. In this case the direction of multiple bond migration will be completely determined by the energy difference between the initial reagents and final products.     

Proton and sodium cation affinities of harpagide: a computational study

The aim of this work was to estimate the proton and sodium cation affinities of harpagide (Har), an iridoid glycoside responsible for the antiinflammatory properties of the medicinal plant Harpagophytum. Monte Carlo conformational searches were performed at the semiempirical AM1 level to determine the most stable conformers for harpagide and its protonated and Na+-cationized forms. The 10 oxygen atoms of the molecule were considered as possible protonation and cationization sites. Geometry optimizations were then refined at the DFT B3LYP/6-31G level from the geometries of the most stable conformers found. Final energetics were obtained at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G level. The proton and sodium ion affinities of harpagide have been estimated at 223.5 and 66.0 kcal/mol, respectively. Since harpagide mainly provides HarNa+ ions in electrospray experiments, the DeltarG298 associated with the reaction of proton/sodium exchange between Har and methanol, MeOHNa+ + HarH+ --> MeOH2+ + HarNa+ (1), has been calculated; it has been estimated to be 1.9 kcal/mol. Complexing a methanol molecule to each reagent and product of reaction 1 makes the reaction become exothermic by 1.7 kcal/mol. These values are in the limit of the accuracy of the method and do not allow us to conclude definitely whether the reaction is endo- or exothermic, but, according to these very small values, the cation exchange reaction is expected to proceed easily in the final stages of the ion desolvation process.     

The reaction of triplet flavin with indole. A study of the cascade of reactive intermediates using density functional theory and time resolved infrared spectroscopy

As a model for riboflavin, lumiflavin was investigated using density functional theory methods (B3LYP/6-31G* and B3LYP/6-31+G**) with regard to the proposed cascade of intermediates formed after excitation to the triplet state, followed by electron-transfer, proton-transfer, and radical[bond]radical coupling reactions. The excited triplet state of the flavin is predicted to be 42 kcal/mol higher in energy than the singlet ground state, and the pi radical anion lies 45.1 kcal/mol lower in energy than the ground-state flavin and a free electron in the gas phase. The former value compares to a solution-phase triplet energy of 49.8 kcal/mol of riboflavin. For the radical anion, the thermodynamically favored position to accept a proton on the flavin ring system is at N(5). A natural population analysis also provided spin density information for the radicals and insight into the origin of the relative stabilities of the six different calculated hydroflavin radicals. The resulting 5H-LF* radical can then undergo radical[bond]radical coupling reactions, with the most thermodynamically stable adduct being formed at C(4'). Vibrational spectra were also calculated for the transient species. Experimental time-resolved infrared spectroscopic data obtained using riboflavin tetraacetate are in excellent agreement with the calculated spectra for the triplet flavin, the radical anion, and the most stable hydroflavin radical.