A density-functional study of the mechanism for the diastereoselective epoxidation of chiral allylic alcohols by the titanium peroxy complexes

The epoxidation of three stereolabeled methyl-substituted chiral allylic alcohols with (1,2)A and/or (1,3)A allylic strain, namely 3-methylbut-3-en-2-ol (1a), pent-3-en-2-ol (1b), and 3-methylpent-3-en-2-ol (1c), have been studied by the density-functional theory method, B3LYP/6-31+G(d,p). For each substrate we calculated the two prereaction complexes with Ti(OH)(4)/MeOOH (the oxidant model for Ti(O-i-Pr)(4)/t-BuOOH), their threo and erythro transition states for oxygen transfer, and the corresponding product complexes. For substrate 1a, the erythro transition state is 0.91 kcal/mol of lower energy than the threo one; for substrates 1b and 1c, the threo compared to the erythro transition states are by 1.05 and 0.21 kcal/mol more favorable, respectively. The threo/erythro product ratios have been estimated from the computed free energies for the competing threo and erythro transition states 3a-c in CH(2)Cl(2) solution to be 12:88 (1a), 92:8 (1b), and 77:23 (1c), which are in good accordance with the experimental values 22:78 (1a), 91:9 (1b), and 83:17 (1c). The diastereoselectivity of this diastereoselective oxyfunctionalization is rationalized in terms of the competition between (1,3)A and (1,2)A strain and the electronic advantage for the spiro transition state. In addition, solvent effects are also play a role for the diastereoselectivity at the same time.     

Epoxidation of acyclic chiral allylic alcohols with peroxy acids: spiro or planar butterfly transition structures? A computational DFT answer

The mechanism of the epoxidation of two chiral allylic alcohols, i.e., 3-methyl-3-buten-2-ol and (Z)-3-penten-2-ol, with peroxyformic acid has been investigated by locating 20 transition structures with the B3LYP/6-31G* method and by evaluating their electronic energy also at the B3LYP/6-311+G**@B3LYP/6-31G* theory level. Relative stability of TSs, as far as electronic energy is concerned, is basis set dependent; moreover, it also depends on entropy and solvent effects. Free enthalpies, calculated by using electronic energy at the higher theory level and with inclusion of solvent effects, indicates that syn, exo TSs, where the olefinic OH group hydrogen bonds the peroxy oxygens of the peroxy acid, outweigh syn, endo TSs, where the peroxy acid carbonyl oxygen is involved in hydrogen bonding. In the former TSs the peroxy acid moiety maintains its planar geometry while in the latter ones a strong out-of-plane distortion of peroxy acid is observed. This distortion makes it viable an unprecedented 1,2-H shift, as a possible alternative to the 1,4-H shift, for the peroxy acid hydrogen. In fact, for one syn, endo TS IRC analysis demonstrated that the 1,2-H shift mechanism is actually operative. The geometry of all TSs substantially conforms to a spiro (i.e., with the peroxy acid plane almost perpendicular to the C=C bond axis) butterfly orientation of the reactants while no TS resembles, even loosely, the planar butterfly structure. Theoretical threo/erythro epoxide ratios are in fair accord with experimental data. Calculations indicate that threo epoxides derive mostly from TSs in which the olefinic OH assumes an outside conformation while erythro epoxides originate from TSs with the OH group in an inside position. Computational findings do not support the qualitative TS models recently proposed for these reactions.     

氮原子与卤代甲烷反应的直接氢抽提过程研究

利用B3LYP理论研究了N(~4S)+CH_3X(X = H, F, Cl)反应体系的直接氢抽提过 程,分别得到了各反应物、产物和过渡态的优化构型和谐振频率。同时应用了6- 31G(d), 6-311+G(d,p)和6-311+ + G(2d,2p)基组,考察其大小对反应体系中各物 种构型及能量的影响。理论计算表明,随着基组的增加,反应势垒逐渐降低,反应 吸热减少。对比取代甲烷的情形,结果表明反应过程中卤素原子具有典型的诱导效 应,降低了抽提势垒。     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

DFT Studies on the Isomerization of Butene Double Bond Catalyzed by 1-Butyl-3-methyl-imidazolium in Ionic Liquid

1 INTRODUCTION Butene and its isomers are important petroleum raw materials. Isomerization reaction of butene plays a key role in the course of C4 alkylation and its reaction mechanism has captured the attention of chemists all along[1, 2]. As a green so…     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

The role of trimethylsilyl group in highly stereoselective eposidation of allylic alcohols

The epoxidation of allylic alcohols carrying Me₃Si group on the double bonds with VO(acac)2-tBuOOH or MCPBA provides erythro or threo epoxy alcohols, respectively with high stereoselectivity.     

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

Computational study of the aminolysis of anhydrides: effect of the catalysis to the reaction of succinic anhydride with methylamine in gas phase and nonpolar solution

Different possible pathways of the aminolysis reaction of succinic anhydride were investigated by applying high level electronic structure theory, examining the general base catalysis by amine and the general acid catalysis by acetic acid, and studying the effect of solvent. The density functional theory at the B3LYP/6-31G(d) and B3LYP/6-311++G(d,p) levels was employed to investigate the reaction pathways for the aminolysis reaction between succinic anhydride and methylamine. The single point ab initio calculations were based on the second-order M?ller-Plesset perturbation theory (MP2) with 6-31G(d) and 6-311++G(d,p) basis sets and CCSD(T)/6-31G(d) level calculations for geometries optimized at the B3LYP/6-311++G(d,p) level of theory. A detailed analysis of the atomic movements during the process of concerted aminolysis was further obtained by intrinsic reaction coordinate calculations. Solvent effects were assessed by the polarized continuum model method. The results show that the concerted mechanism of noncatalyzed aminolysis has distinctly lower activation energy compared with the addition/elimination stepwise mechanism. In the case of the process catalyzed by a second methylamine molecule, asynchronous proton transfer takes place, while the transition vectors of the acid-catalyzed transition states correspond to the simultaneous motion of protons. The most favorable pathway of the reaction was found through the bifunctional acid catalyzed stepwise mechanism that involves formation of eight-membered rings in the transition state structures. The difference between the activation barriers for the two mechanisms averages 2 kcal/mol at various levels of theory.     

Solvent effects on the excited state properties of 2-aminopurine—a theoretical study by the ONIOM and Supramolecular method

In this paper, the micro-solvated effects on the lowest-energy vertical transition state and adiabatic excited states of 2-aminopurine (2Ap) were studied by Supramolecular method (B3LYP/6-31++G(d)) and ONIOM (B3LYP/6-31++G(d):PM3) method. The results are as follows: (1) In 2Ap molecule surrounded by six water molecules the pyramidalization of amino group in 2Ap almost disappears, and the hex-atomic ring is obviously buckled. The adiabatic lowest-energy valence excitation of gaseous 2Ap also causes the disappearance of amino pyramidalization. (2) The energy for lowest-energy singlet π→π^* vertical transition in water is predicted as 3.99 and 4.29 eV by Supramolecular and ONIOM method, respectively. Both values are in good agreement with the reported experimental result, 4.11 eV. The energy for the second lowest-energy n→π^* transition, 4.72 eV, by the Supramolecular method is obviously deviated from the reported experimental value 4.46 eV. The corresponding value given by the two-layer ONIOM method, 4.43 eV, is in good agreement with the experimental value. (3) The optimized energy of the fluorescent emission state (S₁ state) are 3.61 and 3.87 eV by Supramolecular and ONIOM methods, respectively. The calculated oscillator strengths, in both gas and water clusters, were compared with reported experimental and theoretical results. These results indicated that both Supramolecular and ONIOM methods, combined with TD DFT B3LYP/6311++G(d), can provide good results of calculating excited state and spectra properties of 2Ap in condensed phase. This fact encouraged us to extend our study to 2Ap-T base pair and its solvated model so as to obtain the spectra properties of 2Ap in real DNA environment.     

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.     

Vinylogous Mannich reactions: some theoretical studies on the origins of diastereoselectivity

[reaction: see text] Ab initio calculations at the RHF/3-21G level were used to investigate the limiting transition states in the addition of 2-methoxyfuran to a pyrrolinium ion. Four stationary points were found on the potential energy surface with relative energies of threo Diels-Alder, 0.0 kcal/mol; erythro open, 0.9 kcal/mol; erythro Diels-Alder, 1.3 kcal/mol; and threo open, 1.8 kcal/mol.     

Relaxed energetic maps of κ‐carrabiose: A DFT study

The B3LYP density function was used with the 6‐31G(d) basis set to perform relaxed energetic contour maps of the charged form of κ‐carrabiose in the gas phase and for the neutral form first in the gas phase and then by simulating the presence of water as solvent using the Onsager model. Only one starting conformation has been considered to perform all the calculations. Rigid energetic maps have been then constructed either by addition of diffuse or polarization functions to the basis set obtaining in that way 6‐31+G(d)//6‐31G(d), 6‐31+G(d,p)//6‐31G(d), and 6‐311++G(d,p)//6‐31G(d) energetic maps that have been carefully examined. The obtained structures corresponding to the lower energy conformers have been then fully optimized using different basis sets with the B3LYP method, a reversion in term of energy has been observed for the two first minima in the case of the charged disaccharide in the gas phase, this was attributed to the large grid of 30° that could lead to the exclusion of an intermediate value corresponding to the real minimum of energy. We thus suggest that after establishing potential energy maps it is essential to proceed to full optimizations of the lower energy conformers. Calculations using the more accurate correlated method MP2 with the 6‐31G(d) basis set have also been performed for conformers of the two disaccharides in the gas phase. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Diketoacid HIV-1 integrase inhibitors: An ab initio study

The stable tautomeric forms of two representative arene-substituted diketoacid HIV-1 integrase inhibitors, 5-ClTEP and L-731,988, were investigated by B3LYP with 6-31G*, 6-31G(d,p), and 6-31+G(d,p) basis sets. Optimization with MP2/6-31G* was also performed for 5-ClTEP. The solvation effect was considered using a conductor-like screening model. With the density functional theory method, the trans diketo conformations are more stable than the cis conformers. The difference is 14 kJ mol(-1) for 5-ClTEP and 33 kJ mol(-1) for L-731,988. Two trans diketo structures were obtained. The difference between these two trans diketo structures is less than 4 kJ mol(-1) calculated at the B3LYP/6-311+G(3df,2p) level. Two enol forms prevail over the diketo tautomers and are calculated to have the same free energy. Because there is no barrier observed between these two enol forms, they can interchange easily such that a delocalized transition state is suggested to be the observed form. Contradictory to the results of the MP2 method that predicts a preference for the trans diketo forms, the B3LYP method predicts a preference for the enol tautomers, which is in agreement with the experimental results.     

Structure, stability and vibrational spectrum of the hydrogen-bonded complex between HNO3 and H2O. Ab initio and DFT studies

The structure, stability and vibrational spectrum of the binary complex between HONO2 and H2O have been investigated using ab initio calculations at SCF and MP2 levels with different basis sets and B3LYP/6-31G(d,p) calculations. Full geometry optimization was made for the complex studied. It was established that the hydrogen-bonded H2O...HONO2 complex has a planar structure. The corrected values of the dissociation energy at the SCF and MP2 levels and B3LYP calculations are indicative of relatively strong OH...O hydrogen-bonded interaction. The changes in the vibrational characteristics (vibrational frequencies and infrared intensities) arising from the hydrogen bonding between HONO2 and H2O have been estimated by using the ab initio calculations at SCF and MP2 levels and B3LYP/6-31G(d,p) calculations. It was established that the most sensitive to the complexation is the stretching O-H vibration from HONO2. In agreement with the experiment, its vibrational frequency in the complex is shifted to lower wavenumbers. The predicted frequency shift with the B3LYP/6-31G(d,p) calculations (-439 cm(-1)) is in the best agreement with the experimentally measured (-498 cm(-1)). The intensity of this vibration increases dramatically upon hydrogen bonding. The ab initio calculations at the SCF level predict an increase up to five times; at the MP2 level up to 10 times and the B3LYP/6-31G(d,p) predicted increase is up to 17 times. The good agreement between the predicted values of the frequency shifts and those experimentally observed show that the structure of the hydrogen-bonded complex H2O...HONO2 is reliable.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Absolute configurations, predominant conformations, and tautomeric structures of enantiomeric tert-butylphenylphosphinothioic acid.

Vibrational absorption and circular dichroism spectra of dextrorotatory, levorotatory, and racemic mixture of tert-butylphenylphosphinothioic acid have been measured in CCl(4) solutions in the 2000-900 cm(-1) region. The conformations for both tautomeric structures of (S)-tert-butylphenylphosphinothioic acid are investigated using the B3LYP functional with the 6-31G* basis set. For the most stable conformation, the absorption and VCD spectra are predicted ab initio using the B3LYP functional with 6-31G*, 6-311G(2d, 2p), 6-31+G, and 6-311G(3df, 3pd) basis sets. A different functional, B3PW91, was also used with the 6-31G* basis set. The predicted spectra are compared to the experimental spectra. The comparison indicates that (-)-tert-butylphenylphosphinothioic acid is of the (S)-configuration and exists in only one tautomeric structure and one conformation in CCl(4) solution.