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

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

Theoretical studies on heats of formation for cubylnitrates using density functional theory B3LYP method and semiempirical MO methods

The heats of formation (HOF) have been calculated for all the 21 cubylnitrate compounds using the semiemprical molecular orbital (MO) methods (MINDO/3, MNDO, AM1, and PM3) and for 8 of 21 cubylnitrates containing 1–4 ? ONO₂ groups using the density functional theory (DFT) method at the B3LYP/6‐31G* level by means of designed isodesmic reactions. The cubane cage skeletons in cubylnitrate molecules have been kept in setting up isodesmic reactions to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs of the 8 cubylnitrates calculated using B3LYP/6‐31G* and two semiempirical MO (PM3 and AM1) methods, and the linear correlation coefficients of PM3 and AM1 methods are 0.9901 and 0.9826, respectively. Subsequently, the accurate HOFs at B3LYP/6‐31G* level of other 13 cubylnitrates containing 4–8 ? ONO₂ groups are obtained by systematically correcting their PM3‐calculated HOFs. Compared with noncaged nitrates, all the 21 cubylnitrates have high heats of formation implying that they may be very powerful energetic materials and have highly exploitable value. The relationship between the HOFs and the molecular structures of cubylnitrates has been discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

PDDG/PM3 and PDDG/MNDO: improved semiempirical methods

Two new semiempirical methods employing a Pairwise Distance Directed Gaussian modification have been developed: PDDG/PM3 and PDDG/MNDO; they are easily implemented in existing software, and yield heats of formation for compounds containing C, H, N, and O atoms with significantly improved accuracy over the standard NDDO schemes, PM5, PM3, AM1, and MNDO. The PDDG/PM3 results for heats of formation also show substantial improvement over density functional theory with large basis sets. The PDDG modifications consist of a single function, which is added to the existing pairwise core repulsion functions within PM3 and MNDO, a reparameterized semiempirical parameter set, and modified computation of the energy of formation of a gaseous atom. The PDDG addition introduces functional group information via pairwise atomic interactions using only atom-based parameters. For 622 diverse molecules containing C, H, N, and O atoms, mean absolute errors in calculated heats of formation are reduced from 4.4 to 3.2 kcal/mol and from 8.4 to 5.2 kcal/mol using the PDDG modified versions of PM3 and MNDO over the standard versions, respectively. Several specific problems are overcome, including the relative stability of hydrocarbon isomers, and energetics of small rings and molecules containing multiple heteroatoms. The internal consistency of PDDG energies is also significantly improved, enabling more reliable analysis of isomerization energies and trends across series of molecules; PDDG isomerization energies show significant improvement over B3LYP/6-31G* results. Comparison of heats of formation, ionization potentials, dipole moments, isomer, and conformer energetics, intermolecular interaction energies, activation energies, and molecular geometries from the PDDG techniques is made to experimental data and values from other semiempirical and ab initio methods.     

硝酸酯几何构型、生成热和电子结构的PM3 研究

用PM3 SCF-MO方法,通过能量梯度全优化计算,得到39个硝酸酯化合物的分子几何构型、生成热和电子结构。与实验结果和AM₁计算结果进行了比较和评估。     

Improved semiempirical heats of formation through the use of bond and group equivalents

Deficiencies in energetics obtained using the common semiempirical methods, AM1, PM3, and MNDO, may partly be traced to the use of pseudoatomic equivalents for conversion of molecular energies to heats of formation at 298 K. We present an alternative scheme based on the use of bond and group equivalents. Values for the 61 bond and group equivalents necessary for treatment of molecules containing the common organic elements, hydrogen, carbon, nitrogen, and oxygen have been derived. For a set of 583 neutral, closed-shell molecules mean absolute errors in AM1, PM3, and MNDO heats of formation are reduced from 6.6, 4.2, and 8.2 kcal/mol to 2.3, 2.2, and 3.0 kcal/mol, respectively. Several systematic problems are overcome in the present scheme including relative stabilities of branched hydrocarbons, energetics of conjugated systems, heats of formation of long chain hydrocarbons, and enthalpies of molecules containing multiple heteroatoms. Although the approach is restricted to molecules with well-defined functional groups, the equivalents are easy to incorporate and are chemically relevant. This revised procedure allows semiempirical methods to be used for far more reliable evaluations of heats of reactions. Estimates are made of the errors inherent in these semiempirical formalisms, arising from integral approximations and the neglect of explicit treatment of electron correlation effects, while excluding those from inadequate parameterization.     

Comparison of SCC-DFTB and NDDO-based semiempirical molecular orbital methods for organic molecules

Extensive testing of the SCC-DFTB method has been performed, permitting direct comparison to data available for NDDO-based semiempirical methods. For 34 diverse isomerizations of neutral molecules containing the elements C, H, N, and O, the mean absolute errors (MAE) for the enthalpy changes are 2.7, 3.2, 5.0, 5.1, and 7.2 kcal/mol from PDDG/PM3, B3LYP/6-31G(d), PM3, SCC-DFTB, and AM1, respectively. A more comprehensive test was then performed by computing heats of formation for 622 neutral, closed-shell H, C, N, and O-containing molecules; the MAE of 5.8 kcal/mol for SCC-DFTB is intermediate between AM1 (6.8 kcal/mol) and PM3 (4.4 kcal/mol) and significantly higher than for PDDG/PM3 (3.2 kcal/mol). Similarly, SCC-DFTB is found to be less accurate for heats of formation of ions and radicals; however, it is more accurate for conformational energetics and intermolecular interaction energies, though none of the methods perform well for hydrogen bonds with strengths under ca. 7 kcal/mol. SCC-DFTB and the NDDO methods all reproduce MP2/cc-pVTZ molecular geometries with average errors for bond lengths, bond angles, and dihedral angles of only ca. 0.01 A, 1.5 degrees , and 3 degrees . Testing was also carried out for sulfur containing molecules; SCC-DFTB currently yields much less accurate heats of formation in this case than the NDDO-based methods due to the over-stabilization of molecules containing an SO bond.     

Correlation of singlet-triplet gaps for aryl carbenes calculated by MINDO/3, MNDO,AM1, and PM3 with Hammett-type substituent constants

Heats of formation, atomic charges, and geometries of some 110 structures involving substituted singlet and triplet phenyl and 4,4-dimethyl-1,4-dihydronaphthalene carbenes and the corresponding diazomethanes were calculated by MINDO/3, MNDO, AM1, and PM3 semiempirical molecular orbital methods. The singlet-triplet gaps for AM1 and PM3 calculations for the para derivatives in both systems have been successfully correlated with Brown σ⁺ constants. Good correlations with σ⁺ were found for the charges on the carbenic centers of the singlets as well as with the energy barrier for rotation of the aryl group about the C-C single bond in substituted singlet phenylcarbenes. Comparisons of these results with experimental data indicate that AM1 and PM3 are much better than MNDO and MINDO/3 in predicting the intrinsic substituent effects in singlet carbenes.     

Extension of the PDDG/PM3 and PDDG/MNDO semiempirical molecular orbital methods to the halogens

The new semiempirical methods, PDDG/PM3 and PDDG/MNDO, have been parameterized for halogens. For comparison, the original MNDO and PM3 were also reoptimized for the halogens using the same training set; these modified methods are referred to as MNDO' and PM3'. For 442 halogen-containing molecules, the smallest mean absolute error (MAE) in heats of formation is obtained with PDDG/PM3 (5.6 kcal/mol), followed by PM3' (6.1 kcal/mol), PDDG/MNDO (6.6 kcal/mol), PM3 (8.1 kcal/mol), MNDO' (8.5 kcal/mol), AM1 (11.1 kcal/mol), and MNDO (14.0 kcal/mol). For normal-valent halogen-containing molecules, the PDDG methods also provide improved heats of formation over MNDO/d. Hypervalent compounds were not included in the training set and improvements over the standard NDDO methods with sp basis sets were not obtained. For small haloalkanes, the PDDG methods yield more accurate heats of formation than are obtained from density functional theory (DFT) with the B3LYP and B3PW91 functionals using large basis sets. PDDG/PM3 and PM3' also give improved binding energies over the standard NDDO methods for complexes involving halide anions, and they are competitive with B3LYP/6-311++G(d,p) results including thermal corrections. Among the semiempirical methods studied, PDDG/PM3 also generates the best agreement with high-level ab initio G2 and CCSD(T) intrinsic activation energies for S(N)2 reactions involving methyl halides and halide anions. Finally, the MAEs in ionization potentials, dipole moments, and molecular geometries show that the parameter sets for the PDDG and reoptimized NDDO methods reduce the MAEs in heats of formation without compromising the other important QM observables.     

Synthesis and conformational study of acridine derivatives related to 1,4-dihydropyridines

Novel acridine derivatives have been synthesized from dimedone and different aromatic aldehydes by following the classical Hantzsch's procedure. The particular substitution pattern of these compounds is responsible for the observed strong push-pull effect. Quantum chemical calculations were carried out on these molecules by using the AM1 method with complete geometry optimization. The calculated heats of formation reveal two equally favoured conformations. The parameter of planarity and the charge density calculations are in agreement with the ¹³C nmr spectroscopic data.     

Accurate parametrized variational calculations of the molecular electronic polarizability by NDDO‐based methods

The variational method for the calculation of the electronic polarizability of molecules within the NDDO‐based semiempirical MO methods MNDO, AM1, and PM3 was parametrized to improve its accuracy. A training set of 156 compounds was used to fit 34 parameters simultaneously for 12 elements using a simplex optimization. The resulting parameters were tested for a test set of 83 molecules and the calculated polarizabilities compared with the experimental data. For AM1, the RMS deviation between experimental and calculated polarizabilities was reduced from 2.99 (using the original variational treatment) to 0.70 ų for the test set and from 2.81 to 0.40 ų for the training set. MNDO and PM3 gave similar improvements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 17–31, 1999     

Quantitative Correlations for the Estimation of Thermodynamic and Molecular Properties of Pyridines and 2,2′-Bipyridyines

Standard heats of formation, entropies, ionization potentials, and molecular dipole moments of a series of pyridines have been calculated by MNDO, AM1 and PM3 methods. Linear relationship have been established which permit a priori estimation of thermodynamic and molecular characteristics of pyridines. Correlation have been found between the values of pK_a for 2,2′-bipyridines for aqueous solutions and their gas phase proton affinities.__________Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 3, 391–402, March, 2005.     

Evaluation of PM3, AM1, and MNDO for calculation of higher energy ionization potentials

Higher ionization energies were calculated with PM3, AM1, and MNDO for three series of molecules, representative small molecules, molecules containing heteroatoms, and sterically congested alkenes. Values from PM3, AM1, and MNDO were compared to experimental values. In most instances, the semiempirical calculations correctly predict the ordering of higher ionization energies. In the absence of steric hindrance, MNDO is the method of choice. Within groups of molecules, AM1 performs better on hydrocarbons, especially twisted hydrocarbons, than PM3. PM3 commonly gives sigma orbitals which are too high in energy compared to related pi orbitals. PM3 performed better than AM1 with molecules containing oxygen, but failed to give the correct geometry for hydrogen peroxide.     

NDDO semiempirical approximations coupled with Green's function technique—a reliable approach for calculating ionization potentials

The ionization potentials of different molecules have been calculated with the outer valence Green's function (OVGF) technique, coupled with semiempirical MNDO, AM1 and PM3 methods. It is found that the OVGF method gives significantly better agreement with the experimental data than do results obtained with semiempirical calculations using Koopman's theorem including a new SAM1 and MNDO/d methods. Of the three semiempirical methods tested (MNDO, AM1, PM3) the OVGF (AM1) method gives the best agreement with experiment.     

RM1: a reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I

Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.     

Correlation of the acidity of substituted phenols,anilines, and benzoic acids calculated by MNDO,AM1, and PM3 with Hammett-type substituent constants

Heats of formation and net atomic charges of some 120 structures involving substituted phenols, anilines, and benzoic acids and the corresponding anions were calculated by MNDO, AM1, and PM3 semiempirical methods. The gas phase acidities of substituted phenols and anilines and the net atomic charges on the anionic heteroatoms of the corresponding anions have been successfully correlated with σ^? constants. Moreover, good correlations with σ were found for the charges on the acidic hydrogens of substituted phenols and anilines. In contrast, the gas phase acidities of substituted benzoic acids and the charges on the anionic oxygens of the corresponding anions are better correlated with Taft σ° constants. Comparisons of these results with experimental data and ab initio theoretical calculations indicate that AM1 and PM3 methods are much better than MNDO in predicting the acidity of aromatic compounds.     

PM3 study of the proton affinities of 2-, 3-, and 4-monosubstituted pyridines in the gas phase

Heats of formation (ΔH_f) and proton affinities (PA) of 2-, 3-, and 4-monosubstituted pyridines in the gas phase are calculated using the AM1 and PM3 semiempirical methods. The following substitutents are considered: NO₂, CN, CF₃, CHO, F, Cl, COCH₃, H, CH₃, OCH₃, SCH₃, NH₂, and N(CH₃)₂. The results are compared with the experimental data. Both methods reproduce the ΔH_f with comparble accuracy; the rms deviations are 4.1 (AM1) and 4.5 kcal/mol (PM3) for the free bases and 9.5 (AM1) and 9.7 kcal/mol (PM3) for their conjugated acids. The PA are systematically underestimated by both methods, but AM1 appears to be clearly better than PM3 for reproducing the experimental values. The rms deviations for AM1 and PM3 are 5.1 and 9.6 kcal/mol, respectively. This is due to a cancellation of systematic errors in the calculated ΔH_f in the AM1 case and to a summation of the errors in the PM3 case. Both methods correctly reproduce conformations of the molecules under consideration.     

PM3-compatible zinc parameters optimized for metalloenzyme active sites

Recent studies have shown that semiempirical methods (e.g., PM3 and AM1) for zinc-containing compounds are unreliable for modeling structures containing zinc ions with ligand environments similar to those observed in zinc metalloenzymes. To correct these deficiencies a reparameterization of zinc at the PM3 level was undertaken. In this effort we included frequency corrected B3LYP/6-311G* zinc metalloenzyme ligand environments along with previously utilized experimental data. Average errors for the heats of formation have been reduced from 46.9 kcal/mol (PM3) to 14.2 kcal/mol for this new parameter set, termed ZnB for "Zinc, Biological." In addition, the new parameter sets predict geometries for the Bacillus fragilis active site model and other zinc metalloenzyme mimics that are qualitatively in agreement with high-level ab initio results, something existing parameter sets failed to do.     

Comparative semiempirical and DFT study of styrylnaphthalenes and styrylquinolines and their photocyclization products

Semiempirical molecular orbital (PM3, PM6, and RM1) and density functional theory (DFT) (B3LYP/6‐31G*) studies are carried out for 1‐ and 2‐styrylnaphthalenes and their aza‐derivatives—2‐ and 4‐styrylquinolines. Relative stabilities of three isomeric forms: E‐ and Z‐isomers and the closed‐ring dihydrocyclophotoproduct (derivative of dihydrophenanthrene) are calculated. Compared to PM3, PM6 and especially RM1 understate heats of formation; in some cases, PM6 and RM1 even place Z‐isomer in energy below E‐isomer. PM3 rather close to DFT predicts heats of isomerization reaction, whereas PM6 and especially RM1 underestimate these values. Semiempirical methods in comparison with DFT markedly underestimate heats of cyclization reaction; however, reproduce trends in relative stabilities of different isomers in dependence on the structure of styrylnaphthalenes and styrylquinolines. Qualitative correlation is found between calculated relative stabilities of the closed‐ring forms (heats of cyclization reaction) and experimental data: cyclized products with low heats of cyclization are observed in steady‐state photolysis and those with high heats of cyclization are not. In the latter case, the closed‐ring compounds, if formed in the excited state, due to thermal instability decompose rapidly with ring opening in the ground state that prevents their observation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Optimization and application of lithium parameters for PM3

Lithium parameters have been optimized for Stewart's standard PM3 method. The average deviation of the heats of formation calculated for 18 reference compounds is 6.2 kcal/mol from the experimental or high-level ab initio data; the average deviation with Li/MNDO is 18.9 kcal/mol. The average error in bond lengths is also reduced by a factor of two to three. Ionization potentials and dipole moments are reproduced with comparable accuracy than Li/MNDO. However, the mean deviation for the heats of formation of both methods increases when being applied to other systems, especially to small inorgnic molecules. The applicability of the new parameter set is demonstrated further for various compounds not included in the reference set, for the calculation of the activation barriers of several lithiation reactions, as well as for the estimation of oligomerization energies of methyl lithium (including the tetramer). Li/PM3 gives reliable results even for large dimeric complexes, like [{4-(CH₃CR)C₅H₄N}Li]₂, containing TMEDA or THF as coligands and reproduces the haptotropic interaction between Li⁺ and π-systems (e.g., in benzyl lithium) as well as the relative energies and structural features of compounds with “hypervalent” atoms (e.g., in lithiated sulfones). © John Wiley & Sons, Inc.     

Investigation of the properties of some furylpyridine and thienylpyridine derivatives in the gas and aqueous phases using semiempirical methods

The gas and aqueous phase geometries, heats of formation, entropies, dipole moments, proton affinities and methyl cation affinities of some furylpyridines and thienylpyridines have been studied to predict their basicity and pK_a values, using semiempirical AM1, PM3 and MNDO quantum-chemical calculations at SCF level in gas and aqueous phases, with full geometry optimisation.