MNDO calculations on diazirines

The MNDO molecular orbital method gives good agreement between experimental geometries and, when available, heats of formation of diazirines. MNDO and ab initio calculations have been performed on a novel intermediate, methylenediazirine, proposed in a reaction of 3-chloro-3-methyl-diazirine. The relation of this intermediate to its diazo-isomer, its dimer and the products of the reaction are discussed on the basis of MNDO calculations.     

5,5'-双(4-甲基苯)-2,2-联恶唑的晶体结构及构象

本文报道5,5'-双(4-甲基苯基)-2,2'-联恶唑(POOP)的晶体结构,并辅以分子力学法(MM-X程序)计算,探讨了POOP分子的基态构象.     

Theoretical study on the transition state of oxaphosphetane formation between ethylidenetriphenylphosphorane and acetaldehyde

The MNDO MO calculations were carried out for the formation of oxaphosphetanes from substituted ylides and aldehydes (Equations 2–4). It was found that the reactions proceed through a nearly planar cyclic four-membered ring transition state (TS) in all cases studied and that phenyl substitution on phosphorus has little effect on the TS geometry. This finding is not in line with the TS models that have been proposed to explain the experimentally observed cis-selectivity.     

Polycyclic aromatic hydrocarbons with five-membered rings: Modeling of experimental X-ray and neutron-diffraction structures

Several of the readily available theoretical programs are evaluated as tools for modeling the structures of polycyclic aromatic hydrocarbons with five-membered rings (CPAHs). The experimentally determined bond lengths and angles are compared to calculated values. Experimental bond lengths are also compared to Pauling and Huckel molecular orbital (HMO) bond orders. Previously published experimental X-ray and neutron-diffraction structures of acenaphthene, acenaphthylene, fluoranthene, cyclopent[o,p,q,r]benz[c]phenanthrene, and corannulene are modeled by the programs MMX, AM1, MNDO, and PM3, and previously reported STO-3G and 6-31G * data are also evaluated. In general, the error differences between the experimental and calculated results for all of the semiempirical programs were small. However, PM3 performed slightly better than AM1 and MMX, while MNDO generated structures which exhibited the largest deviation from experiment. Although the standard deviations for all programs are shown to be of comparable magnitude, a particular bond length or bond angle in any given theoretical calculation can exhibit significant error from the experimental data. The scatter in the bond order data computed from Huckel molecular orbital theory and valence bond theory is contrary to results obtained with alternant systems. It appears that these approaches are less successful at modeling accurately the nonalternant hydrocarbon systems described in this paper.     

Ab initio studies on infinite linear hydrogen fluoride chains

Ab initio crystal orbital calculations on linear infinite chains of hydrogen fluoride and MO calculations on HF and the linear dimer have been performed. Equilibrium geometries, force constants, band structures, densities of states and longitudinal phonon dispersions are presented, and compared with the available data. In agreement with experiment the most common features of hydrogen bonding, elongation of HX bond length (ΔR_HX and decrease in HF stretching force constants, are much more pronounced in the solid state than in the isolated dimer.     

共轭烯烃的分子几何构型、电子结构和生成热的DMM和AM1、PM3的比较研究

运用Ｄｅｌｆｔ分子力学（ＤＭＭ）力场和程序以及半经验分子轨道ＡＭ１和ＰＭ３方法计算研究了丁二烯、苯、甲苯、联苯、苯乙烯、富烯、、环辛四烯、［２，２］对环烷和菲等１０个共轭烯烃分子的几何构型、电子结构和生成热．ＤＭＭ计算的几何构型和生成热与实验结果相吻合，电荷分布结果与从头计算结果较接近．ＡＭ１和ＰＭ３计算的几何构型较好，但计算的生成热与实验结果偏差较大．ＰＭ３计算值比ＡＭ１的稍好．     

Quantum study of the ground state of a series of unsaturated boron–nitrogen compounds by the MNDO method. I. Geometries and stabilities

A MNDO study of molecular geometries, enthalpies, formation, atomization, bond separation, and hydrogenation of a series of unsaturated boron–nitrogen compounds, linear with two and three members and cyclic with three, four, five, and six members are presented. For all these molecules, MNDO calculations are in excellent agreement with available ab initio calculations or experimental data. The high rotational barrier in aminoborane H₂BNH₂, 30.6 kcal/mol, and the length of bond in iminoborane HBNH, 1.183 Å, imply strong double and triple BN bond character in these two molecules. In the odd membered heterocycles, examination of the molecular geometries and energies of equilibrium states shows that in all cases, the stability of the compounds grows with the number of boron atoms and decreases with that of the nitrogen atoms. Moreover, compared study of the two BN-fulvenes with their homologous hydrocarbon shows that only BN (BB)-fulvene has a polyenic structure similar to that of fulvene.     

Pyrrolizidine alkaloids necine bases: Ab initio,semiempirical, and molecular mechanics approaches to molecular properties

The structural stabilities of endo and exo conformations of retronecine and heliotridine molecules were analyzed using different ab initio, semiempirical, and molecular mechanics methods. All electron and pseudopotential ab initio calculations at the Hartree-Fock level of theory with 6-31G* and CEP-31G* basis sets provided structures in excellent agreement with available experimental results obtained from X-ray crystal structure and ¹H-NMR (nuclear magnetic resonance) studies in D₂O solutions. The exo conformations showed a greater stability for both molecules. The most significant difference between the calculations was found in the ring planarity of heliotridine, whose distortion was associated with the interaction between the O(11)H group and the C(1)-C(2) double bond as well as with a hydrogen bond between O(11)H and N(4). The discrepancy between pseudopotential and all-electron optimized geometries was reduced after inclusion of the innermost electrons of C(1), C(2), and N(4) in the core potential calculation. The MNDO, AM1, and PM3 semiempirical results showed poor agreement with experimental data. The five-membered rings were observed to be planar for AM1 and MNDO calculations. The PM3 calculations for exo-retronecine showed a greater stability than the endo conformer, in agreement with ab initio results. A good agreement was observed between MM3 and ab initio geometries, with small differences probably due to hydrogen bonds. While exo-retronecine was calculated to be more stable than the endo conformer, the MM3 calculations suggested that endo-heliotridine was slightly more stable than the exo form. © 1996 by John Wiley & Sons, Inc.     

A study of hydration effects on the conformational aspects of gaba mediators

The conformation of numerous chemical compounds is strongly influenced by solvents. Knowledge of their structure in solution is necessary, especially for a discussion of the biological and pharmacological activity of the molecules. The neurotransmitters and their agonists and antagonists are known to be flexible molecules that interact with quite distinct receptors. The conformational properties of several GABA (γ-aminobutyric acid) mediators have been studied by the MNDO technique. The optimized geometries of the molecules have been used to study the solvent effects on their conformational properties considering the supermolecule approach for their first hydration shell. A conjectural pharmacophoric pattern for several GABA inhibitors has already been suggested from the molecular electrostatic potentials (MEP ) of several molecules using a localized bond orbital technique. In the present work, MEP calculations have been carried out considering a solvent effect on the MNDO optimized geometries to investigate any deviation from the earlier results.     

MNDO Calculations of silicon-containing molecules

A comparison is made of MNDO and MINDO /3 calculations for saturated silicon-containing molecules, and with experimental values, for heats of formation, molecular geometries, charge distributions, and ionization potentials. Except for bond angles, it is found that with the published parameter values the MINDO /3 program gives more reliable results than MNDO . For unsaturated molecules, a comparison of bond lengths and stabilities of Si multiple bonds as given by the two programs and ab initio methods is made, and large discrepancies between predicted structures are pointed out. Some reasons for the dicrepancies are discussed.     

Comparison of the use of the MNDO and MINDO/3 methods with ab initio methods to study tautomerism in histamine, 2- and 4-methylhistamine

The semiempirical MNDO and MINDO/3 methods are used to study the various tautomeric forms of histamine, 2-methylhistamine, and 4-methylhistamine. Comparisons of the optimized structures and tautomerization energies are made with values obtained from ab initio Hartree-Fock calculations using the 3-21G and STO-3G basis sets. Based on these results and previous comparisons of STO-3G results with x-ray structures, the present results indicate that while there are some differences in the values of the structural parameters, the changes in structure upon tautomerization and/or protonation are very similar. Further analysis of the MNDO and MINDO/3 structures by means of their utilization in 3-21G and STO-3G calculations indicates that either of these semiempirical methods provides reliable values for the structural parameters. Both methods give good qualitative agreement with the ab initio calculations for the relative energies of the various tautomers in the three compounds. In these studies the MNDO method appears to give better quantitative agreement with the 3-21G and STO-3G results than the MINDO/3 method.     

Theoretical study of the mechanism of the Wittig reaction: Ab initio and MNDO-PM3 treatment of the reaction of unstabilized,semistabilized, and stabilized ylides with acetaldehyde

In this study, we describe the results of ab initio (HF and MP2) and MNDO-PM3 calculations on the model reactions of unstabilized (Me₃P=CH–CH₃), semi-stabilized (Me₃P=CH–C≡CH), and stabilized (Me₃P=CH–C≡N) ylides with acetaldehyde to form their respective Z and E olefins and trimethylphosphine oxide. These reactions occur in three stages: oxaphosphetane formation, oxaphosphetane pseudorotation, and oxaphosphetane decomposition. The calculated barriers for these processes vary considerably depending on the level of theory employed (ab initio vs. MNDO-PM3 or HF vs. MP2 at the ab initio level). However, self-consistent geometries of reactants, intermediates, transition states and products are obtained at all levels. Oxaphosphetane formation is best described as very asynchronous cycloaddition (borderline two-step mechanisms). The geometries of the transition states are near planar with respect to P, C, C, and O atoms. Analysis of the bond indices of these reactions shows that the C–C bonds are between 44% (unstabilized case) and 60% (stabilized case) formed, whereas the corresponding P–O bonds have not been formed to any significant degree. Oxaphosphetane decomposition can be described as a very asynchronous retrocycloaddition where P–C bond breakage runs ahead of C–O bond breakage. These results are compared with experimental findings for the Wittig reaction, and its relevance to the overall mechanism of the olefination is discussed. © 1997 John Wiley & Sons, Inc. Heteroatom Chem 8: 557–569, 1997     

硝酸酯分子几何构型的量子化学研究

运用MINDO / 3、MNDO 和AM1 三种半经验分子轨道(MO)方法, 通过SCF计算, 首次系统地获得了32个硝酸酯化合物分子的全优化几何构型。三种方法的计算结果与已报道的四个化合物(硝酸甲酯、吉纳、硝化甘油和太安)的实验结果相比, AM1法较好。所有硝酸酯的酯基(-ONO~2)具有近似不变的几何参数。直链烷基硝酸酯的键长和键角极为相近, 全部重原子均共平面。二元直链和四元硝酸酯具有对称的分子构型。     

Computational structural determination and energy landscape analysis of the hepatic carcinogen 2-(acetylamino)fluorene

 2-(Acetylamino)fluorene (AAF), a potent mutagen and a prototypical example of the mutagenic aromatic amines, forms covalent adducts to DNA after metabolic activation in the liver. A benchmark study of AAF is presented using a number of the most widely used molecular mechanics and semiempirical computational methods and models. The results are compared to higher-level quantum calculations and to experimentally obtained crystal structures. Hydrogen bonding between AAF molecules in the crystal phase complicates the direct comparison of gas-phase theoretical calculations with experiment, so Hartree–Fock (HF) and Becke–Perdew (BP) density functional theory (DFT) calculations are used as benchmarks for the semiempirical and molecular mechanics results. Systematic conformer searches and dihedral energy landscapes were carried out for AAF using the SYBYL and MMFF94 molecular mechanics force fields; the AM1, PM3 and MNDO semiempirical quantum mechanics methods; HF using the 3-21G*and 6-31G* basis sets; and DFT using the nonlocal BP functional and double numerical polarization basis sets. MMFF94, AM1, HF and DFT calculations all predict the same planar structures, whereas SYBYL, MNDO and PM3 all predict various nonplanar geometries. The AM1 energy landscape is in substantial agreement with HF and DFT predictions; MMFF94 is qualitatively similar to HF and DFT; and the MNDO, PM3 and SYBYL results are qualitatively different from the HF and DFT results and from each other. These results are attributed to deficiencies in MNDO, PM3 and SYBYL. The MNDO, PM3 and SYBYL models may be unreliable for compounds in which an amide group is immediately adjacent to an aromatic ring. Received: 26 May 2002 / Accepted: 12 December 2002 / Published online: 14 February 2003     

An empirical MM2 augmented force field for the cycloimmonium ylides

An original procedure approach taking into account the implementation of parameters determined, using calculations based on density functional theory, for the amidocyano-pyridinium methylide in the MM2 augmented harmonic potential function has been proposed. A good agreement between theoretical force field calculation and X-ray diffraction data has been observed. Thus an empirical force field for cycloimmonium ylides has been established. It provides good quality geometries for cycloimmonium ylide molecules by energy minimization. In this study we proposed a new MM2 augmented atom-type for the ylidic carbon atom. To our knowledge no attempt has been done in this way for such organic systems. Thus, we have shown that parameterization established by the DFT method is able to reproduce or to predict with good accuracy the structures of the cycloimmonium ylide compounds. This study also includes a full conformational analysis.     

A density functional treatment of organolithium compounds: Comparison to ab initio,semiempirical, and experimental results

Density functional calculations on several classes of organolithium compounds are described. The compounds studied include lithium bonds to carbon, oxygen, and nitrogen and are representative of most types of organolithium compounds that have appeared in the recent literature. The computational results are compared to those using MNDO, which has been shown to have some serious deficiencies in compounds involving carbon–lithium bonds, and to PM3 results, which offer some improvement over MNDO for many organolithium compounds. Most of the density functional calculations with a large basis set are in good agreement with available ab initio and experimental data. Calculated carbon–lithium bond lengths were slightly shorter than those calculated by other ab initio methods and were substantially longer than those calculated by MNDO, which is known to underestimate carbon–lithium bond lengths severely. Dimerization energies of methyllithium, calculated by DMol, were also in good agreement with those of other ab initio calculations. Lithium–nitrogen bonds in lithium amides were calculated to be slightly shorter by DMol than by MNDO, although the two methods were in qualitative agreement for this type of compound. © 1995 by John Wiley & Sons, Inc.     

Dynamic NMR study of the oxaphosphetane complexation with lithium during the Wittig reaction

Dynamic NMR spectroscopy at very low temperatures (148–182 K) reveal the dynamic behavior of the (2‐tri(3‐furyl)‐3‐methyl‐4‐cyclopropenyl‐oxaphosphetane) generated during a Wittig reaction between tri(3‐furyl)ethylphosphonium iodine and cyclopropylaldehyde. The possibility of formation of different adducts between Li⁺ ions and oxaphosphetane or betainic intermediates was checked calculating the formation enthalpies using the MNDO, AM1, and PM3 semiempirical MO methods. The observed species are interpreted as oxaphosphetane complexes with lithium ions present in solution. Quantum mechanical calculations confirm the spectroscopic results. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Bond length calculations for norbornane,bicyclo[2.2.2]octane,and related olefins using a consistent force field. Influence of stretch-bend crossterms

It is shown that a valence force field CCC stretch-bend crossterm which is consistent in sign and magnitude with vibrational spectroscopic calculations leads to bond stretchings in norbornane, norbornenes, and bicyclo[2.2.2]octenes. The bond length effects are in satisfactory agreement with most of the available experimental results and have not been obtained by earlier force field calculations.     

Evaluation of MINDO/3 calculated structures. II. Branching errors in alkanes and cycloalkanes

MINDO/3 calculations have been carried out for a series of branched chain alkanes in order to assess effects of branching on calculated geometries and heats of formation (ΔH_f). With vicinal branching, MINDO/3 calculates the central C? C bond to be too long. Bond angles are also found to be distorted. Errors in calculated heats of formation are large when geminal branching is present and significant with vicinal branching. Branching error corrections for ΔH_f have been derived and applied to a separate series of branched acyclic and cyclic compounds. For the test sample, application of the branching error corrections gave calculated structures of acyclic branched hydrocarbons with heats of formation having an average absolute error of 1.3 kcal/mole rather than 17.3 kcal/mole before correction. Cyclic branched hydrocarbons are shown to be less well corrected. Calculations of heats of reaction have also been carried out for some isomerization and cyclization reactions using the MINDO/3 and MNDO methods. It is clear from the comparisons that MNDO calculations give less severe errors for highly branched compounds but the errors are still substantial. For prediction of heats of reaction, the error-corrected calculations are shown to be superior to the “raw” calculations obtained by MINDO/3 or MNDO.