<Emphasis Type="Italic">Ab Initio</Emphasis> conformational analysis of tetrahydro-1,3-oxazine

Nonempirical RHF/STO-3G and 6-31G(d) studies on conformational behavior of tetrahydro-1,3- oxazine showed that interconversion between the axial (global minimum) and equatorial chair conformers can follow five independent pathways. The potential energy surface contains seven minima corresponding to chair and twist conformers, as well as seven transition states having sofa, half-chair, and symmetrical and unsym- metrical boat conformations. Additional potential barriers to interconversion, resulting from pyramidal inversion of the nitrogen atom, were also revealed.     

Conformational analysis of 1,3-oxathiane

Computer simulation of pathways of conformational isomerization of 1,3-oxathiane molecule carried out with the help of HF/6-31G(d), MP2.6-31G(d)/HF/6-31G(d), and PNE/3z quantum-chemical approximations showed that interconversion between the degenerate in energy chair conformers proceeds through seven independent pathways: directly and via six flexible forms. Potential energy surface contains eight minimum points including chair conformers and enantiomeric pairs of twist forms, and also five transition states, among them different modification of semi-chair, symmetric and unmmetrical boat. Molecular dynamics methods show that flexible forms at room temperature convert into one another and into the chair conformers.     

Stereochemistry of Seven-Membered Heterocycles: XLI. Stereoselective Synthesis of 3-R-1,5-Dihydro-2,4-benzothiepine 2-Oxides

The oxidation of 3-R-1,5-dihydro-2,4-benzothiepines with m-chloroperbenzoic acid occurs in a highly diastereoselective fashion. The resulting trans-sulfoxides at -60°C in CDCl3 exist as an equilibrium mixture of the chair and boat forms with the substituents in the equatorial position. The fraction of the boat conformation increases in the series R = Ph, Me, t-Bu.     

NBO and NICS analysis of the allylic rearrangements (the Cope and 3-aza-Cope rearrangements) of hexa-1,5-diene and N-vinylprop-2-en-1-amine: A DFT study

In this work, ab initio density functional theory (DFT) calculations have been performed on the 3,3-sigmatropic rearrangements of hexa-1,5-diene (Cope) and N-vinylprop-2-en-1-amine (3-aza-Cope) in the gas phase. The barrier heights and heats of reactions calculated at the B3LYP/6-311G** level of theory were in good agreement with experimental data. Transition states optimized with B3LYP/6-311G** theory were used for calculating the nucleus independent chemical shift (NICS) and, a natural bond orbital (NBO) analysis was also performed at the same level of theory. Our results indicate that the aromaticities of the transition states are controlled by the out-of-plane component and that the chair-like transition state of the Cope rearrangement exhibits the strongest aromatic character. Analysis of donor-acceptor (bonding and anti-bonding) interactions of σ_3–4 → π*_1–2 suggests that the TS structure in the hexa-1,5-diene reaction (the Cope rearrangement) has more aromatic character than the N-vinylprop-2-en-1-amine reaction (the 3-aza-Cope rearrangement). The NBO results show that in the hexa-1,5-diene and N-vinylprop-2-en-1-amine rearrangements, activation energies are controlled by σ_3–4 → π*_1–2 and σ_3–4 → π*_1–2 resonance energies.      

PO → Si intramolecular coordination in the derivatives of 1,4‐phosphasilacyclohexane 1‐oxides

The chair and boat conformers for a series of derivatives of 1,4‐phosphasilacyclohexane 1‐oxides have been calculated at the B3LYP/6‐311+G** level of theory in the gas phase and taking into account the effect of solvent polarity using the IEF‐PCM model. The stability of the boat conformers containing pentacoordinate silicon due to formation of the P?O→Si intramolecular coordination bond depends on the environment of the phosphorus atom and polarity of the solvent, and the strength of the transannular bond depends also on the nature of the substituents at the silicon atom. The highly polar boat conformers are strongly stabilized in the DMSO solution. NBO analysis showed the importance of the σ(C? Si) → σ*(H₃C? N) hyperconjugative interaction in the two H₃C? N? C? Si fragments of the ring favoring the formation of the pentacoordinate silicon atom. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Anomeric effect and theoretical vibrational spectra compared to the experiment in 2-chloro-1,3,2-dioxaphosphorinane-2-oxide, -sulfide,and -selenide. I. Anomeric effect and structures

2-Chloro-1,3,2-dioxaphosphorinane-2-oxide, -sulfide, and -selenide are studied at the DFT/B3LYP level and several ab initio methods using a 6-311G** basis set. Our energy optimizations by all these methods show that for oxide DFT and ab initio methods are not much different, while for the sulfide and the selenide the DFT relative energies are higher by about a kcal/mol as compared to those of MP2, MP3, MP4(SDTQ)//MP2, and CCSD(T)//MP2 (//MP2 indicates that a single-point calculation based on the MP2 optimized geometry is performed). However, regardless of rather large relative energies, that does not change the fact that in all three cases the conformational equilibrium mixture contains more than 95% of the lowest, chair-equatorial conformer (this indicates that the P=X bond is in equatorial position). This one and the next higher conformer (chair-axial) are confirmed to be real conformers (energy minima) in all cases. The energetically much higher twist and boat forms are probably just stationary states and local maxima because in many cases, geometry optimizations do not converge to them. Only for MP2 and the selenide do all optimizations converge to the desired stationary state. The relative energies could all be explained in terms of anomeric effects and ring strains. The decreasing covalent character of the P=X bond, with X changing from O to S and to Se, shows itself in the increasing bond lengths and the decreasing strength of anomeric effects.     

All-trans- and t,T,t,C,t,T,t-deca-1,3,5,7,9-Pentaenes: Ab Initio Structures,Vibrational Analyses,and Some Regularities in the Series of Related Molecules

Total geometry optimization and calculation of the force constants for all-transand t,T,t,C,t,T,tdeca-1,3,5,7,9-pentaene were carried out at the ab initio, HF/6-31G level. The HF/6-31G//HF/ 6-31G force fields were modified using empirical scale factors transferred from trans-buta-1,3-diene augmented by an additional scale factor for the central formal carbon-carbon double bond coordinates (determined previously for all-trans-hexa-1,3,5-triene). The total number of scale factors was seven. The vibrational problems for both decapentaenes were solved using the respective scaled HF/6-31G//HF/6-31G force field. Infrared intensities and Raman activities were calculated from the unscaled HF/6-31G//HF/6-31G force fields. Complete assignment of all the fundamental vibrational frequencies is given. Geometrical parameters, vibrational frequencies and force constants are compared with the corresponding values of buta-1,3-diene, hexa-1,3,5-triene and octa-1,3,5,7-tetraene. Regularities in the properties of this molecular series are discussed. Special attention is given to the possibility of using the vibrational spectra for detection of distortions from the regular trans structure of these oligoenes.     

Stereochemistry of seven-membered heterocycles: XLV. Highly diastereoselective addition of dimethyl 1,2,4,5-tetrazine-3,6-dicarboxylate to 2-substituted 1,3-dithiacyclohept-5-enes

Conformationally heterogeneous 2-substituted 1,3-dithiacyclohept-5-enes (R = Ph, Me, t-Bu), which exist in solution as chair and boat conformers, react with dimethyl 1,2,4,5-tetrazine-3,6-dicarboxylate with high exo-diastereoselectivity: only the chair conformer is involved. The steric structure of 4-methyl-3,5-dithia-9,10-diazabicyclo[5.4.0]undeca-7,10-diene was determined by X-ray analysis. Its crystal packing and supramolecular structure were also analyzed.     

Effect of basic set quality and electron correlation on the scale factors of a harmonic force field

The completely optimized structure and harmonic force field of s-trans-buta-1,3-diene are reported at the MP2/6-31G and MP2/6-31G* levels of computation. Sets of empirical scale factors for the calculated force fields are derived and compared with the corresponding values computed at the RHF/4-31G and RHF/6-31G levels. Changes in the scale factors for this series of force fields are discussed. The vibrational frequencies are also reported for thirteen isotopomers of s-trans-buta-1,3-diene using the MP2/6-31G* force field. Some characteristics of the gauche and cis forms of buta-1,3-diene are also given.     

DFT法研究离子液中EMIM＋催化丁烯双键异构反应机理

利用密度泛函方法(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了离子液体中1-乙基-3-甲基咪唑阳离子(EMIM＋)催化丁烯双键异构反应过程中的反应物、产物以及过渡态的几何构型,分析了反应过程中键参数的变化.通过振动分析对平衡态和过渡态进行了验证,并得到了零点能.通过计算内禀反应坐标(IRC),确认了对应于过渡态的反应物和产物. 计算结果表明,EMIM＋催化丁烯双键异构可以基元反应的方式一步完成,1-丁烯异构化为2-丁烯的活化能约为192 kJ•mol－1, 逆反应活化能约为208 kJ•mol－1, 可在室温或高于室温条件下进行.     

Stereochemistry of seven-membered heterocycles: XLVII. Synthesis,NMR, XRD,and theoretical investigation of monosulfoxides of (<Emphasis Type="Italic">Z</Emphasis>)-4,7-dihydro-1λ<Superscript>4</Superscript>,3-dithiepine series

By oxidation of (Z)-4,7-dihydro-1,3-dithiepines with meta-chloroperbenzoic acid the corresponding 4,7-dihydro-1λ⁴,3-dithiepine 1-oxides were synthesized in ∼85% yield. A complete assignment of signals in the ¹H and ¹³C NMR spectra of the compounds obtained was established. According to XRD analysis the seven-membered heterocycles crystallized in conformations chair or boat. The calculations by DFT method in the B3LYP/6-31G(d,p) version showed the possibility of conformational equilibrium of these forms in the gas phase.     

Molecular modeling of ring-chain equilibria for the ring-opening cross-metathesis of cis,cis-1,5-dimethyl-cycloocta-1,5-diene with ethylene at T=298.15 K

The molecular modeling of ring-opening cross-metathesis of cis,cis-1,5-dimethyl-cycloocta-1,5-diene (1,5-DM-COD), cis,cis-1,6-dimethyl-cycloocta-1,5-diene (1,6-DM-COD) and cis,cis-cycloocta-1,5-diene (COD) with ethylene (ethenolysis) at T=298.15 K using the B3LYP/6-31G(d,p) level of theory reveals that ring-chain equilibrium constants are dependent on the nature of cyclic diene. The ring-chain equilibria for the ethenolysis of 1,5-DM-COD is completely shifted to the formation of monomeric 2-methyl-hexa-1,5-diene.     

Synthesis and conformational study of P‐heterocyclic androst‐5‐ene derivatives

The reactions of (20R)‐3β‐acetoxy‐21‐hydroxymethylpregn‐5‐en‐20‐ol ( 2 ) and (20R)‐3β‐acetoxypregn‐5‐ene‐20,21‐diol ( 11 ) with phenylphosphonic dichloride 3 and aryl dichlorophosphates 4–6 afforded novel types of P‐heterocyclic androst‐5‐ene derivatives 7–10 and 12 as epimeric pairs. The diastereomers were separated by column chromatography and were characterized by NMR spectroscopy. Estimation of the stereostructures of the corresponding epimers by B3LYP/631G(d) DFT ab initio calculations suggested that the six‐membered hetero ring in compounds 7b and 8a–10a adopts predominantly a chair conformation, with the P‐substituents in their preferred orientation. The cyclic phosphonate moiety in 7a or 8b–10b , however, seems to exist as an equilibrium mixture of chair–distorted‐boat or chair–chair forms. The theoretical calculations indicate that the conformational equilibrium is shifted toward the distorted‐boat conformer for 7a , with a pseudoequatorial P‐phenyl substituent, whereas for 8b–10b the chair conformer with an equatorial P‐phenoxy group predominates. © 2008 Wiley Periodicals, Inc. Heteroatom Chem 19:7–14, 2008; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20372     

Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths

The minimum energy path (MEP) of the reaction, CF₃CHFCF₃ + H → transition state (TS) → CF₃CFCF₃ + H₂, has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6‐31++G**, BH&HLYP/cc‐pVDZ, BMK/6‐31++G**, M05/6‐31+G**, M05‐2X/6‐31+G**, UMP2/6‐31++G**, PUMP2/6‐31++G**//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVDZ//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVTZ(spd,sp)//UMP2//6‐31++G**, RCCSD(T)/CBS//M05/6‐31+G**, and RCCSD(T)/CBS//UMP2/6‐31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero‐curvature, and small‐curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000–1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6‐31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (～1500 K), SCT corrections are significant at low temperatures (～300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections. © 2012 Wiley Periodicals, Inc.     

Ab initio study of possible intermediates in the no(4Π) catalyzed geometric isomerizations of olefins

The geometries of acyclic and three-membered ring (nitroxide) H₄C₂NO radicals in their ground ²Π electronic states have been optimized completely at ab initio UHF and ROHF theoretical levels with the STO-3G and the 6-31G** basis sets. The optimizations favour the cyclic nitroxide structure energetically. However ΔE(acyclic - cyclic) at the UHF and ROHF/6-31G** levels are only 3.2 and 1.9 kcal mol^-1, respectively. Incomplete MP2/6-311G** optimizations support these results. The zero-point energy computed at the ROHF/6-31G** level for the nitroxide radical is 2.5 kcal mol^-1 higher than that for the acyclic structure, thus reversing the relative energies by 0.6 kcal mol^-1. The energies of the two radical structures, relative to the sum of those for ethylene and NO, are very close to literature values of the activation energies for the thermal, NO catalyzed geometrical isomerizations of olefins. Thus cyclic nitroxide intermediates may play a role not only in the Hg 6(³P₁) photosensitized, but also in the thermal, NO-catalyzed geometric isomerizations of olefins. Paper dedicated to Professor Otto P. Strausz; presented in part at the 75th Canadian Chemical Congress and Exhibition, Edmonton, May 31 – June 4, 1992.     

In-depth quantum chemical investigation of electro-optical and charge-transport properties of trans-3-(3,4-dimethoxyphenyl)-2-(4-nitrophenyl)prop-2-enenitrile

The structural, electro-optical and charge-transport properties of compound trans-3-(3,4-dimethoxyphenyl)-2-(4-nitrophenyl)prop-2-enenitrile (DMNPN) were studied using quantum chemical methods. The neutral, cation and anion molecular geometries were optimized in the ground state using density functional theory (DFT) at the restricted and unrestricted B3LYP/6-31G** level of theory. The excited state geometries were optimized by applying time-dependent DFT at the TD-B3LYP/6-31G** level of theory. The absorption and fluorescence wavelengths were calculated at the TD-CAM-B3LYP/6-31G** and TD-LC-BLYP/6-31G** levels of theory. The distribution pattern of the charge densities on the highest occupied molecular orbitals (HOMOs) and lowest unoccupied molecular orbitals (LUMOs) are discussed. Intramolecular charge transfer was observed from the dimethoxyphenyl to (nitrophenyl)prop-2-enenitrile moieties. The detailed charge-transport behavior of the DMNPN molecule is investigated based on its ionization potential, electron affinity, hole and electron reorganization energies, hole and electron-transfer integrals, and hole and electron intrinsic mobilities. The total/partial densities of states and structure–property relationship are discussed in detail. The higher computed hole intrinsic mobility than electron intrinsic mobility reveals that DMNPN is an efficient hole-transport material.     

Estimation of the relative stability of chair forms of hexahydropyrimidine

Nonempirical approximations HF/6-31G(d), MP2/6-31G(d)//HF/6-31G(d), PBE/3z, and RI-MP2/λ2 were applied to the calculation of the relative stability and parameters of the conformational equilibrium of three chair forms of the unsubstituted hexahydropyrimidine with different orientation of the N-H bonds. It was established that in the framework of the used methods the most stable are the conformers of diaxial and axial-equatorial chairs. The occupancy of the diequatorial chair form at the room temperature does not exceed 1.7%.     

DFT法研究3-羟基丙烯醛的双键旋转异构反应机理

利用密度泛函理论(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了基态3-羟基丙烯醛分子在双键旋转异构反应过程中的平衡态以及过渡态的几何构型，分析了反应过程中键参数的变化，计算了该反应的内禀反应坐标(IRC)，发现在重排反应途径上存在一个四元环骨架的中间体.通过振动分析对平衡态和过渡态进行了确认，并得到了零点能.计算结果表明，基态3-羟基丙烯醛分子的双键旋转异构反应经过两步完成，第一步反应位垒稍高，第二步反应位垒较低，存在着发生重排反应的可能性.     

3-羟基丙烯醛单键旋转异构光反应机理的从头算研究

应用从头算分子轨道法分别在RHF/6-31G**和UHF/6-31G**水平上对3腔┤┓肿拥幕?S0)和三重激发态(T1)单键旋转异构反应机理进行了研究,优化出反应物和产物在S0态和T1态的4种平衡态和过渡态的几何构型,通过振动分析得到的虚频和计算的内禀反应坐标对过渡态进行了确认,并得到了零点能,根据基态和激发态反应途径分析了光反应机理。计算结果表明,基态和激发态的3腔┤┓肿拥ゼ旃狗从ξ焕萁系停嘏湃菀追⑸永砺凵涎橹ち斯夥从κ笛橹屑钢忠旃固宓拇嬖凇?     

An SCF and CI Study of the 1,3 Shift in the HX?CHY⇌XCH?YH Isoelectronic Series: X,YCH2, NH,and O

Ab initio molecular orbital calculations at SCF level with the 3-21G, 6-31G, and 6-31G** basis sets and CI level with the 6-31G basis set have been carried out for an isoelectronic series HX? CH?Y and X?CH? YH, where X, Y can be CH₂, NH, and O. Optimized structures (3-21G and 6-31G**) for both tautomers and the 1,3 hydrogen shift transition states are reported. The relative stabilities of the isomers and the barriers of the 1,3 shift are discussed in terms of proton affinities and bond orders. It is shown that both the relative stabilities of the tautomers and the relative barrier heights can be explained qualitatively using simple proton affinity arguments and that the barrier heights are quantitatively related to bond orders.