Ab initio molecular orbital calculations of 3,4-dihydro-1,2-dioxin, 3,6-dihydro-1,2-dioxin, 4H-1,3-dioxin (1,3-diox-4-ene), and 2,3-dihydro-1,4-dioxin (1,4-dioxene)

Optimized geometries and total energies for 3,4-dihydro-1,2-dioxin ( 1 ), 3,6-dihydro-1,2-dioxin ( 2 ), 4H-1,3-dioxin (1,3-diox-4-ene, 3 ), and 2,3-dihydro-1,4-dioxin (1,4-dioxene, 4 ) were calculated using ab initio 3-21G, 6-31G*, and MP2/6-31G*//6-31G* methods. The half-chair conformers of 1 (C₁), 2 (C₂), 3 (C₁), and 4 (C₂) are more stable than their respective planar structures [ 1 (C_s), 2 (C_2v), 3 (C_s), and 4 (C_2v)]. Among the four isomers 1 – 4 , the half-chair conformer of 3 is the most stable. It is 53.1, 54.6, and 3.4 kcal mol⁻¹ more stable than 1 , 2 , and 4 , respectively. The largest energy difference (19.0 kcal mol⁻¹) is observed between the half-chair and planar conformers of 2 . The boat conformers of 2 and 4 are less stable than their respective half-chair conformers, but are more stable than their planar structures. Hyperconjugative orbital interactions (anomeric effects) contribute to the greater stability of 3 (n_O(3) →σ*_C(2)—O(1), n_O(3)→σ*,n _O(3)→σ*) and of 4 (n_O(1)→ σ*). The ab initio calculated structural features of the half-chair conformations of the dihydrodioxins 1 – 4 are compared with the half-chair conformations of cyclohexene and the chair conformations of cyclohexane, oxacyclohexane (tetrahydropyran), 1,2-dioxacyclohexane (1,2-dioxane), 1,3-dioxacyclohexane (1,3-dioxane), and 1,4-dioxacyclohexane (1,4-dioxane) © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1392–1406, 1997     

Ab initio molecular orbital calculations for 3,6-dihydro-1,2-dithiin and 3,6-dihydro-1,2-dioxin

Optimized geometries and total energies for the conformers of 3,6-dihydro-1,2-dithiin ( 2 ) and 3,6-dihydro-1,2-dioxin ( 3 ) were calculated at several ab initio MO levels: RHF/3-21G(*), RHF/6-31G*, MP2/6-31G*, and MP2/6-31G*/ /RHF/3-21G(*). For the dioxin, in addition to the above levels the corresponding nonextended basis sets ab initio methods were also carried out. The dithiin results are compared with those of simple disulfanes, HSSH and (CH₃)₂S₂, whose optimized geometries agree closely with the observed structures, which is the gauche (C₂ symmetry). For the disulfanes, the gauche geometries from RHF/3-21G(*) are in good agreement with the observed structure while the RHF/3-21G results best fit the dioxin. Pertinent structural data at the RHF/3-21G(*) for the half-chair (C₂) dithiin are: bond lengths, ? SS? , ? CS? , ? CC?, and ? C?C? , 2.050, 1.817, 1.515, and 1.317 Å, respectively; bond angles, CSS, ?CCS, and C?CS, 98.0, 114.2, and 127.8°, respectively; CSSC dihedral angle of 63.2°; and twist angle of 36.5°. The total energy for half-chair dithiin at MP2/6-31G*//RHF/3-21G(*) is less than the planar (C_2v) and the half-boat (C_s) structures by 69.67 and 29.05 kJ/mol, respectively. The calculated structural data (vs. observed) at RHF/3-21G for the half-chair dioxin are: bond lengths, ? OO? , ? CO? , ? CC?, and C?C, 1.464 (1.463), 1.454, 1.509, and 1.313 Å (1.338 Å), respectively; bond angles, COO, ?CCO, and C?CO, 105.0, 109.8 (110.3), and 120.7° (119.9°), respectively; COOC dihedral angle of 79.7° (80 ± 2°); and twist angle of 39.0 (38.3°). The total energy for half-chair dioxin at MP2/6-31G//RHF/3-21G is less than the planar and the half-boat structures by 70.35 and 42.85 kJ/mol, respectively. The total energies calculated at the extended basis sets (*) ab initio levels for the C₂ symmetry dioxin are much lower than those of the nonextended basis sets. © John Wiley & Sons, Inc.     

Calculated structures and relative stabilities of furoxan,some 1,2-dinitrosoethylenes and other isomers

The structures and relative stabilities of furoxan and some of its isomers, e.g., the 1,2-dinitrosoethylenes, have been determined by means of ab initio Hartee–Fock and Møller–Plesset calculations. Geometries were optimized at the HF/3-21G, HF/6-31G* and MP2/6-31G* levels, and subsequently used for computing MP2/6-31G*, MP3/6-31G*, and MP4/6-31G* energies. The results are markedly affected by the inclusion of electronic correlation, which renders three of the isomers unstable. It also emphasizes the importance of a zwitterionic contribution to the structure of furoxan, which promotes ring-opening through a cis 1,2-dinitrosoethylene intermediate/transition state that has an MP4/6-31G*//MP2/6-31G* energy that is 31.6 kcal/mol above furoxan.     

Structure,conformation and NMR studies on 1,2-dioxane and halogen substituted 1,2-dioxane molecules

Molecular structure and conformational stability of chair and twist conformers of 1,2-dioxane and halogen substituted compounds of the 1,2-dioxane have been studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 1,2-dioxane, 3,6-difluoro, 3,6-dichloro, 3,3,6,6-tetrafluoro and 3,3,6,6-tetrachloro 1,2-dioxane compounds were optimized at HF, MP2, B3LYP and B3PW91 levels of theory by implementing 6-31G* basis set. To study the effect of polar medium, self-consistent reaction field theory is used to optimize the conformers at B3LYP/6-31G* level of theory. The geometrical parameters of chair and twist conformers have been discussed in the light of interaction between lone pair electrons present in the oxygen and substituted halogen atoms. The relative stability of the conformers have been studied using relative energy, maximum hardness principle and thermodynamical quantities. The 13C-NMR chemical shift study for carbon atoms in the title compounds are calculated and the results have been discussed.     

Rearrangement of azirine intermediates to nitriles: Theoretical study of cleavage of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to cyanoketene–formaldimine complex

A detailed investigation of the reaction path for the thermal rearrangement of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to yield a cyanoketene–formaldimine complex is carried out at the MP2/6-31G* and B3LYP/6-31G* levels of theory. The ring opening of the five-membered pyrrolinone ring and the formation of the nitrile group takes place in a concerted manner, presenting a significant strain energy release and allowing for an electronic stabilization by coarctate conjugation of the transition structure (TS). These two factors make possible a moderate energy barrier. Although the structural features B3LYP/6-31G* theoretical levels, it is found that the MP2 energy barrier (28.8) CCSD(T)/6-31G*//MP2/6-31G* value (17.1 kcal/mol). The complex electronic rearrangement can be rationalized using the theory of coarctate transition structures developed by Herges as the evolution of an azirine structure without referring to a hypothetical vinyl nitrene intermediate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 912–922, 1998     

Synthesis, structure, and electrochemical properties of biferrocenes annulated with 1,2-dithiin and 1,2-dithiin 1,1-dioxides

2,2″-Bis(N,N-dimethylaminosulfonyl)-1,1″-biferrocene (6), a precursor of biferrocenes annulated with 1,2-dithiin and 1,2-dithiin 1,1-dioxides, was prepared by a sequence of selective ortho-lithiation and dimerization reaction from N,N-dimethylaminosulfonylferrocene. New biferrocenes annulated with 1,2-dithiin (1) and 1,2-dithiin 1,1-dioxides (2) and (3) were successfully synthesized in satisfactory yields by the reaction of compound 6 with lithium aluminum hydride followed by treatment with chlorotrimethylsilane. The electrochemical properties of the biferrocenes (1)-(3) were furnished by voltammetric studies.     

Predictive Abilities of Scaled Quantum Mechanical Molecular Force Fields: Application to 1,3-Butadiene

The positions of some IR bands of the s-trans-1,3-butadiene-h ₆ and -1,1,2-d ₃ isotopomers in the gas phase have been measured using a Brucker IFS 120 HR spectrometer with a resolution of 2 cm^–1. The structural parameters of the s-trans- and s-gauche-1,3-butadiene conformers were optimized completely at the MP2/6-31G* theoretical level and their MP2/6-31G*//MP2/6-31G* quantum mechanical force fields (QMFFs) were calculated. Using only the experimental vibrational frequencies of s-trans-1,3-butadiene-h ₆ the QMFF of the s-trans conformer was corrected by Pulay's scaling method (eight scale factors were involved). The scaled QMFF was used to calculate the mean vibrational amplitudes and the Coriolis coupling constants of s-trans-1,3-butadiene-h ₆ and the vibrational frequencies of 12 of its deuterated isotopomers. The set of scale factors obtained for correction of the s-trans QMFF was transferred to the QMFF of the s-gauche conformer. Its theoretical vibrational spectrum and those of some deuterated and C¹³ isotopomers were calculated. The ability of this scaling approach (transferring of scale factors) to predict the vibrational frequencies of rotational conformers and their isotopomers, as well as other molecular characteristics, and to permit detection of perturbations of the experimental bands are discussed.     

Ab initio MO calculations on 1,2-dithietes and valence isomers

The valence isomerization of the 1,2-dithiete parent compound to the open-chain dithial was studied by CASSCF multiconfiguration methods including the CASPT2 perturbational treatment. The isomerization energy remains small at the highest level of theory. In agreement with Jonas and Frenking, the cyclic structure is only then preferred over the acyclic ones if f-functions on the sulfur atoms are considered. If they are included, the 1,2-dithiete is more stable by 3.8 kcal/mol and the barrier amounts to 24.9 kcal/mol at the CASPT2(8,8)/6-31G(2df)//MP2/6-31G* level of theory. According to MP2/6-31G* geometry optimizations, substitution of H by NH₂ and CH₃ reduces the stability of the 1,2-dithiete ring structures relative to the open-chain dithiocarbonyl structures, whereas the inverse holds for acceptor substitution by CN and CF₃. A higher stability of benzodithiete relative to ortho-dithiobenzoquinone is predicted at all employed levels of theory. This is in good agreement with conclusions drawn from experimental results. Whereas the experimental microwave geometry of 1,2-dithiete is well reproduced theoretically, the experimental electron diffraction geometry of 3,4-bis-(trifluoromethyl)-1,2-dithiete differs from the calculated one. © 1996 John Wiley & Sons, Inc.     

A Computational Study of Conformers of 2-Thiaoxacyclohexane (1,2-Oxathiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(2d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), and 6-311G(d,p) basis sets and the hybrid density functionals B3LYP, B3P86, and B3PW91 have been used to calculate the optimized geometries and relative energies of the chair, half-chair, sofa, twist, and boat structures of 2-thiaoxacyclohexane (1,2-oxathiane). The values of the energy difference (E, kcal/mol) between the chair and 3,6-twist structures of 1,2-oxathiane were 4.92 (HF), 4.73 (MP2), and 4.66 (DFT). The HF chair–twist energy difference (G _c–t ^o) for 1,2-oxathane was 5.16 kcal/mol. Intrinsic reaction coordinate (IRC) calculations connected a transition state (TS-A) between the chair conformation and the less stable 2,5-twist form and connected two transition states (TS-B, TS-C) between the chair conformation and the more stable 3,6-twist conformer. The DFT energy differences between the chair and TS-A, TS-B, and TS-C were 11.4, 10.8, and 12.6 kcal/mol, respectively. Hyperconjugative stereoelectronic interactions were observed in the chair (n _o and ) and 3,6-twist (n _S and n _O ) conformers of 1,2-oxathiane. The chair conformation of 1,2-oxthiane is 9.6 and 10.0 kcal/mol, respectively, less stable than the chair conformations of 3-thiaoxacyclohexane (1,3-oxathiane) and 4-thiaoxacyclohexane (1,4-oxathiane, thioxane).     

Theoretical predictions of the structure, gas-phase acidity and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure, acidity (free energy of deprotonation, ΔG^o), and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid (3,4-dithiohydroxy-3-cyclobutene-1,2-diselenone, H₂C₄Se₂S₂) are reported. The global minimum found on the potential energy surface of 1,2-diseleno-3,4-dithiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanediselenone, and cyclobutenedithiol. The computed aromatic stabilization energy (ASE) by homodesmotic reaction (Eq 1) is −20.1 kcal/mol (MP2(fu)/6-311+G** //RHF/6-311+G**) and −14.9 kcal/mol (B3LYP//6-311+G**//B3LYP/6-311+G**). The aromaticity of 1,2-diseleno-3,4-dithiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −17.91 (CSGT(IGAIM)-RHF/6-311+G**//RHF/6-311+G**) and −31.01 (CSGT(IGAIM)-B3LYP/6-311+G**//B3LYP/6-311+G**). Thus, 1,2-diseleno-3,4-dithiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The calculated theoretical gas-phase acidity is ΔG^o _1(298K)=302.7 kcal/mol and ΔG^o _2(298K)=388.4 kcal/mol. Hence, 1,2-diseleno-3,4-dithiosquaric acid is a stronger acid than squaric acid(3,4-dihydroxy-3-cyclobutene-1,2-dione, H₂C₄O₄). Received: 11 April 2000 / Accepted: 7 July 2000 / Published online: 27 September 2000     

Proton affinities of molecules containing nitrogen and oxygen: comparing ab initio and semiempirical results to experiments

We report a comparison of theoretical and experimental proton affinities at nitrogen and oxygen sites within a series of small molecules. The calculated proton affinities are determined using the semiempirical methods AM 1, MNDO , and PM 3; the ab initio Hartree–Fock method at the following basis levels: 3-21G //3-21G , 3-21+G //3-21G , 6-31G *//6-31G *, and 6-31+G (d, p)//6-31G *; and Møller–Plesset perturbation calculations: MP 2/6-31G *//6-31G *, MP 3/6-31G *//6-31G *, MP 2/6-31G +(d, p)//6-31G *, MP 3/6-31G +(d, p)//6-31G *, and MP 4(SDTQ )/6-31G +G (d, p)//6-31G *. The semiempirical methods have more nonsystematic scatter from the experimental values, compared to even the minimal 3-21G level ab initio calculations. The thermodynamically corrected 6-31G *//6-31G * proton affinities provide acceptable results compared to experiment, and we see no significant improvement over 6-31G *//6-31G * in the proton affinities with any of the higher-level calculations. © 1992 John Wiley & Sons, Inc.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Conformational analysis of carbonic anhydrase inhibitors using ab initio molecular orbital methods. 1. Rotational isomerism in methane sulfonamide anion, CH3-SO2-NH−

Ab initio calculations have been performed on methane sulfonamide anion. Geometries have been optimized using Hartree-Fock basis sets up to 6-31+G*, and single-point calculations employing those Hartree-Fock geometries have been performed at levels up to MP2/6-311++G**. In addition, geometry optimizations for the 0°, 90°, 150°, and 180° conformers have been carried out at the MP2/6-31G*, MP2/6-31+G*, and MP2/6-311++G** levels. Vibrational frequencies have been calculated using the HF/4-31G*, MP2/6-31G*, and MP2/6-31+G* geometries. All calculations at or above the 4-31G* level agree that H—N—S—C ˜90° is the global minimum. The H—N—S—C = 180° conformer is clearly higher in energy although the relative energy of this conformer varies from 0.36 to 1.03 kcal/mol for the post-HF calculations depending on basis set. The H—N—S—C = 180° conformer appears to be a very shallow local minimum. However, the potential energy surface is quite flat in this region, and the highest-level calculations, including MP2 optimizations and vibrational frequency analysis, are ambiguous on this point. The conformer with an H—N—S—C torsion of 0° is a transition state with a relative energy ˜8 kcal/mol. Received: 3 December 1996 / Accepted: 2 January 1997     

Bifunctionalized allenes. Part XXIV. Competitive electrophilic cyclization of 5-(dimethoxyphosphoryl)-alka-3,4-dienoates leading to 2,5-dihydro-1, 2-oxaphospholes and 5,6-dihydro-2H-pyranes

Abstract

We report herein a study on the competitive electrophilic cyclization of 5-(dimethoxyphosphoryl)-alka-3,4-dienoates involving 5-endo-trig and 6-endo-trig mode cyclizations. Reaction with electrophiles produces mixtures of the 2-(2-oxo-2,5-dihydro-1,2-oxaphosphol-5-yl)-alkanoates and (6-oxo-5,6-dihydro-2H-pyran-2-yl)-phosphonates by competitive electrophilic cyclization due to the participation of the neighboring phosphonate and carboxylate groups.     

Ab initio study of the conformational equilibrium of ethylene glycol

Summary The conformational equilibrium of ethylene glycol (CH₂OHCH₂OH) has been examined by performing geometry optimizations at the 6-31G*, MP2/6-31G* and 6-31G** levels. Final energies have been calculated at the MP3 level with the optimized geometries. The two most stable conformers are atGg andgGg but it is verified that the inclusion of electronic correlations reduces their energy difference of 0.6 kcal/mol at the HF level to less than 0.2 kcal/mol. The possible coexistence of the two most stable conformers is in agreement with some previous studies of Frei et al. For thetXg conformer a detailed analysis of the intramolecular potential as a function of rotation around the C-C bond is also reported.     

Experimental and theoretical studies on the thermal decomposition of heterocyclic nitrosimines

A series of substituted 2-nitrosiminobenzothiazolines (2) were synthesized by the nitrosation of the corresponding 2-iminobenzothiazolines (6). Thermal decomposition of 2a--f and of the seleno analogue 7 in methanol and of 3-methyl-2-nitrosobenzothiazoline (2a) in acetonitrile, 1,4-dioxane, and cyclohexane followed first-order kinetics. The activation parameters for thermal deazetization of 2a were measured in cyclohexane (Delta H(++) = 25.3 +/- 0.5 kcal/mol, Delta S(++) = 1.3 +/- 1.5 eu) and in methanol (Delta H(++) = 22.5 +/- 0.7 kcal/mol, Delta S(++) = -12.9 +/- 2.1 eu). These results indicate a unimolecular decomposition and are consistent with a proposed stepwise mechanism involving cyclization of the nitrosimine followed by loss of N(2). The ground-state conformations of the parent nitrosiminothiazoline (9a) and transition states for rotation around the exocyclic C==N bond, electrocyclic ring closure, and loss of N(2) were calculated using ab initio molecular orbital theory at the MP2/6-31G* level. The calculated gas-phase barrier height for the loss of N(2) from 9a (25.2 kcal/mol, MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE) compares favorably with the experimental barrier for 2a of 25.3 kcal/mol in cyclohexane. The potential energy surface is unusual; the rotational transition state 9a-rot-ts connects directly to the orthogonal transition state for ring-closure 9aTS. The decoupling of rotational and pseudopericyclic bond-forming transition states is contrasted with the single pericyclic transition state (15TS) for the electrocyclic ring-opening of oxetene (15) to acrolein (16). For comparison, the calculated homolytic strength of the N--NO bond is 40.0 kcal/mol (MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE).     

Conformational analysis of 1,2-difluoroethane, an ab-initio investigation

Contrary to expectation, the gauche conformer of 1,2-difluoroethane is more stable than the trans conformer in the gas phase. In order to understand the underlying causes of the “gauche effect”, a complete geometry relaxation was performed for the gauche and trans conformers of 1,2-difluoroethane with the 4-21G, 4-31G and 4-31G** basis sets. The 4-31G** optimized geometry of the gauche conformer compares well with the experimental values obtained from a number of electron-diffraction studies. A correction for the correlation energy, calculated by means of second-order Møller—Plesset perturbation theory with the 6-31G** basis set, proves to be essential to obtain a correct estimate of the energy difference between the gauche and trans conformers 1,2-difluoroethane.     

吡啶-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^*)。     

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.