Structural assignment of Diels-Alder adducts: an experimental and theoretical approach

A detailed NMR analysis with total assignment of (1)H and (13)C NMR data for the endo and the exo adducts, obtained by Diels-Alder reaction between 2-cyclohexenone and cyclopentadiene, is described. The unequivocal assignment of the endo and exo structures was performed by (1)H and (13)C NMR. These assignments were supported by theoretical chemical shift calculations at GIAO/HF level using 6-311 + g(2d, p) from optimized structures at the B3LYP/6-31g(d) level.     

A computational study of the Diels-Alder reaction of ethyl-S-lactyl acrylate and cyclopentadiene. Origins of stereoselectivity

The eight diastereoisomeric transition structures of the Diels-Alder addition of ethyl-S-lactyl acrylate and cyclopentadiene have been investigated in the gas phase and in solution by HF, MP2, and density functional theory (B3LYP and B3PW91) methods with the 6-31G(d,p) basis set. At all levels of theory used, the s-cis transition structures are more stable than the s-trans ones. The contribution of the s-trans transition structures increases in solution and, although still small, has to be taken in consideration for correct prediction of stereoselectivity. Diastereofacial selectivity is interpreted in terms of electrostatic (weak hydrogen bonding) C=O...H(C) interactions between the carbonyl group(s) of the dienophile and cyclopentadiene in the energetically favored transition structures. Endo/exo reaction selectivity is attributed to positive orbital interactions between the diene and the acrylate carbonyl oxygen in the endo s-cis transition structures. Ab initio methods reproduce well the experimentally observed trends in both endo/exo and diastereofacial selectivity. Density functional calculations in the gas phase correctly reproduce the observed trends in diastereofacial selectivity but single-point MP2 calculations are necessary to reproduce the experimental trend in endo/exo selectivity.     

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF,MP2, DFT,and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

Ab Initio Study of Endo/Exo and Diastereofacial Selectivities in Diels-Alder Reactions between Chiral Butenolides and Cyclopentadiene

Diels-Alder reactions of cyclopentadiene with crotonolactone and beta-angelica lactone have been studied using ab initio methods at the MP3/6-31G//RHF/3-21G and MP2/6-31G//RHF/6-31G levels. The transition states corresponding to the formation of different stereoisomers and diastereoisomers have been located. The results obtained show that the correct endo/exo selectivity is only obtained when polarization functions are included in the basis set. However, syn/antiselectivity is correctly described at all levels of calculation.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

Stereoselectivity in Diels-Alder reactions of diene-substituted N-alkoxycarbonyl-1,2-dihydropyridines

Diene substituent effects on the regiochemical and stereochemical outcomes of uncatalyzed Diels-Alder reactions of N-alkoxycarbonyl-1,2-dihydropyridines with both styrene and methyl vinyl ketone (MVK) were studied. Alkyl substitution on the diene in all cases examined resulted in a kinetic preference for 7-endo isomers (7-phenyl 51-96% exo and 7-acetyl 54-96% exo). For both dienophiles, the highest stereoselectivities (>or=89% endo) were observed with 5-methyl or 6-methyl substituents in the dihydropyridine. Theoretical calculations of the energies of gas phase endo and exo transition states at the RHF/3-21G(*) predict that total entropy, DeltaStotal, considerations favor endo cycloadducts for both dienophiles with DHP, while total energy considerations, DeltaEo, favor endo cycloadducts for styrene and exo cycloadducts for MVK. At this level, favored endo-phenyl isomers are correctly predicted for styrene reactions, but the calculation of 7-acetyl exo or endo isomer dominance is diene-substituent-dependent for MVK reactions. The preference for endo addition of MVK to the parent, 5-methyl, and 6-methyl-DHPs was successfully predicted by calculations at the B3LYP/6-31G* theory level.     

Radical ring expansion reactions of methylenecyclopropane derivatives: a theoretical study

The evolution of the primary radicals from 1-(3-bromopropyl)-2-ethyl-3-methylenecyclopropane, 1-(3-bromopropyl)-1-trimethylsilyl-2-methylenecyclopropane, 1-(3-bromobutyl)-2-ethyl-3-methylenecyclopropane, and 1-(3-bromobutyl)-1-trimethylsilyl-2-methylenecyclopropane was theoretically studied at the ROMP2/6-311++G(d,p)//UB3LYP/6-31G(d,p) theory level taking into account the effect of solvent through a PCM-UAHF model. For the propyl-substituted radicals, the attack of the radical center on the double bond takes place most favorably in an exo fashion. The subsequent ring expansions yield the product corresponding to the rupture of the endo C-C bond as the most favorable one in accordance with the experimental results. In the case of 1-(3-bromobutyl)-2-ethyl-3-methylenecyclopropane, the Gibbs energy barriers for the endo and exo attacks are the same, and the subsequent reversible evolution yields the product corresponding to the rupture of the exo C-C bond as the most favorable one through thermodynamic control in agreement with experiment. Finally, for 1-(3-bromobutyl)-1-trimethylsilyl-2-methylenecyclopropane, our calculations predict that the endo attack is 0.8 kcal/mol more favorable than the exo one. In the subsequent reversible ring expansion, the product corresponding to the rupture of the endo C-C bond is kinetically the most favored one in reasonable agreement with the experimental observations.     

A new insight on hydrogen bond donor property of bridged B–H–B protons in the interaction of diborane with some hydrogen bond acceptor molecules

Structural Chemistry - Ab initio calculations at the MP2/6-311++G(d,p) computational level were used to analyze the interactions between a molecule of B2H6 with some small molecules which often...     

A theoretical study of the kinetics of the benzylperoxy radical isomerization

The rate constant of the benzylperoxy isomerization reaction has been computed using 54 different levels of theory and has been compared to the experimental value reported at 773 K. The aim of this methodology work is to demonstrate that standard theoretical methods are not adequate to obtain quantitative rate constants for the reaction under study. The use of the elaborated CASPT2 method is essential to estimate a quantitative rate constant. Geometry optimizations and vibrational frequency calculations are performed using three different methods (B3LYP, MPW1K, and MP2) and six different basis sets (6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), and cc-pVDZ). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the 6-31G(d,p) basis set, and with the CASPT2 level of theory with the ANO-L-VDZP basis set. Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. We recommend the use of the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory to compute the temperature dependence of the rate constant of the four-center isomerization of the benzylperoxy radical. It is given by the following relation: k(600-2000 K) (in s (-1)) = (1.29 x 10 (10)) T (0.79) exp[(-133.1 in kJ mol (-1))/ RT]. These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. This computational procedure can be extended to predict rate constants for other similar reactions where no available experimental data exist.     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Theoretical study of the regioselectivity of the cycloaddition reaction between cyclopentadiene and methyleneketene

Three possible reaction schemes for the cycloaddition reaction between methyleneketene and cyclopentadiene were studied by DFT (density functional theory) and ab initio calculations. All of these cycloaddition reactions are exothermic, concerted but nonsynchronous processes. The computed activation energies indicate that the norbornene product yielded from a 1,2-addition of methyleneketene with cyclopentadiene (reaction (1)) is the primary product. The performance of various computational methodologies, MP2, MP4, and CCSD(T), in conjunction with a wide array of basis sets, 6-31G(d), 6-311+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ, in obtaining reliable activation and reaction energies of the reactions under investigation has been critically analyzed.     

Ab initio investigation of the diels-alder reaction between 2H-phosphole and phosphaethene: A model for phosphole dimerization

All four possible Diels-Alder reactions between 2H-phosphole and phosphaethene were examined at various theoretical levels, including HF, MP4SDQ, CCSD(T), and CASSCF. MP2/6-31G^* geometry optimizations could not be employed since the potential energy surface is qualitatively incorrect at this level of theory, due to the inherent underestimation of the activation energies (ameliorated at higher-order MP or coupled-cluster levels). Solvent effects were examined employing the Onsager, polarized continuum, and isodensity and surface polarized continuum models. At MP4SDQ/6-31G^*//HF/6-31G^* these reactions are exothermic by 34–38 kcal mol⁻¹ and have very low activation energies, 5–7 kcal mol⁻¹. The P P/C C regioisomer products are lower in energy than the C P isomers and, within each pair, the exo isomer is lower in energy. At low computational levels the smallest activation energy is for the reaction leading to the C P endo product. Larger basis sets, electron correlation, and solvent favor the transition state leading to the experimentally observed P P/C P endo isomer. The dimerization of phosphole is, therefore, kinetically controlled. Based on geometric and electron density analysis, the reactions are concerted and synchronous. © 1997 by John Wiley & Sons, Inc.