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.     

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.     

Application of the computationally efficient self-consistent-charge density-functional tight-binding method to magnesium-containing molecules

The geometric properties, ionization potentials, heats of formation, incremental binding energies, and protonation energies for up to 75 magnesium-containing compounds have been studied using the self-consistent-charge density-functional tight-binding method (SCC-DFTB), the complete-basis set (CBS-QB3) method, traditional B3LYP density-functional theory, and a number of modern semiempirical methods such as Austin Model 1 (AM1), modified neglect of diatomic overlap without and with inclusion of d functions (MNDO, MNDO/d), and the Parametric Method 3 (PM3) and its modification (PM5). The test set contains some widely varying chemical motifs including ionic or covalent, closed-shell or radical compounds, and many biologically relevant complexes. Geometric data are compared to experiment, if available, and otherwise to previous high-level ab initio calculations or the present B3LYP results. SCC-DFTB is found to predict bond lengths to high accuracy, with the root-mean-square (RMS) error being less than half that found for the other semiempirical methods. However, SCC-DFTB performs very poorly for absolute heats of formation, giving an RMS error of 29 kcal mol(-1), but for this property B3LYP and the other semiempirical methods also yield poor but useful results with errors of 12-22 kcal mol(-1). Nevertheless, SCC-DFTB does provide useful results for biologically relevant chemical-process energies such as protonation energies (RMS error 10 kcal mol(-1), with the range 6-19 kcal mol(-1) found for the other semiempirical methods) and ligation energies (RMS error 9 kcal mol(-1), less than the errors of 12-23 kcal mol(-1) found for the other semiempirical methods). SCC-DFTB is shown to provide a computationally expedient means of calculating properties of magnesium compounds, providing results with at most double the inaccuracy of the high-quality but dramatically more-expensive B3LYP method.     

The remarkably invariant interaction energies of lithium first-row compounds with water and with ammonia

Solvation energies of lithium first-row compounds LiX (X ? H, Li, BeH, BH₂, CH₃, NH₂, OH, F) and of the lithium cation with the model solvents, water and ammonia, have been calculated ab inito (MP2/6-31 + G^*//6-31G^* with zero-point vibrational energy corrections at 3-21G//3-21G). The solvation energies are found to be remarkably constant: ?18.0 ± 1.2 and ?21.5 ± 1.3 kcal/mol for the hydrates and ammonia solvates, respectively. This independence on the nature of X is due largely to the ionic character of the LiX compounds (dipole moments 4.7–6.6 debye). The unexpectedly high solvation energies of the lithium molecule (?14.3 and ?17.8 kcal/mol, respectively) are due to the polarizability of Li₂. At the same level, the lithium cation has interaction energies with H₂O and NH₃ of ?34.1 and ?39.7 kcal/mol, respectively. For the hydrates of LiOH and LiF cyclic structures with hydrogen bonds and somewhat increased solvation energies also are described.     

MNDO method for calculations of magnesium clusters

A new set of parameters for the magnesium atom has been developed within the MNDO method. In contrast to previously published parameters, the new parameters correctly describe molecules with different chemical natures: magnesium halides, organomagnesium compounds and the recently found small magnesium clusters Mg _n (n=2-8). The average errors in the calculated heats of formation and bond lengths of magnesium compounds, including clusters are: 10.7 kcal/mol and 0.167 Å, respectively.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1384–1388, August, 1994.     

Optimization of parameters for semiempirical methods II. Applications

MNDO/AM1-type parameters for twelve elements have been optimized using a newly developed method for optimizing parameters for semiempirical methods. With the new method, MNDO-PM3, the average difference between the predicted heats of formation and experimental values for 657 compounds is 7.8 kcal/mol, and for 106 hypervalent compounds, 13.6 kcal/mol. For MNDO the equivalent differences are 13.9 and 75.8 kcal/mol, while those for AM1, in which MNDO parameters are used for aluminum, phosphorus, and sulfur, are 12.7 and 83.1 kcal/mol, respectively. Average errors for ionization potentials, bond angles, and dipole moments are intermediate between those for MNDO and AM1, while errors in bond lengths are slightly reduced.     

Theoretical studies of [n]paracyclophanes and their valence isomers. I. Geometries,strain energies,and enthalpies of the inter-conversions of [n]paracyclophanes and their Dewar benzee isomers

Four semiempirical methods (AM1, MNDO, PM3, and MINDO/3) are used to calculate the deformation angles of [n]paracyclophanes and their Dewar benzene isomers for n = 3… 10. The results obtained by all these methods are in good agreement with data from X-ray studies. We have determined the strain energies that, in both series of compounds, are due to two components: (1) the strain energy of deformation of the cycle (aromatic or Dewar Benzene skeletons) and (2) the strain energy of the oligomethylene chain. In [6]paracyclophane, the strain energy [SE_ring(MNDO) ≈? 32.9 kcal/mol] almost compensates the resonance energy (E_resonance ≈ 36 kcal/mol) so that its chemical properties are closer to alkenes than to benzenic compounds. To better reproduce the enthalpy of the valence isomerization [n]Dewar bezene → [n]paracyclophane, which is poorly calculated with these methods, a correction is proposed and the reaction enthalpy of [6]paracyclophane is estimated to be about ΔH_r ≈ 15 ± 15 kcal/mol. It is found that MNDO and MINDO/3 need the smallest corrections, but MNDO leads to better geometries than MINDO/3. In conclusion, MNDO seems to be the best technique for further studies of these compounds. © 1992 by John Wiley & Sons, Inc.     

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.     

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.     

Extension of MNDO to d orbitals: Parameters and results for the halogens

A recently proposed extension of the MNDO formalism to d orbitals has been parameterized for the halogens CI, Br, and I. Extensive test calculations indicate slight consistent improvements for normalvalent molecules and dramatic improvements for hypervalent molecules, in comparison with established MNDO -type methods without d orbitals. The mean absolute errors in calculated heats of formation are 3.9 kcal/mol for 155 normalvalent compounds and 2.8 kcal/mol for 23 hypervalent compounds. The predicted structures of the hypervalent molecules are qualitatively correct, with a mean absolute error of 2° in 19 bond angles.     

Nonlocal density functional calculation of gas phase heats of formation

A general procedure is presented for computing the gas phase heats of formation of a wide variety of organic compounds. ΔE for the formation of the molecule from its elements at 0 K is obtained from density functional calculations (Gaussian 92/DFT) for optimized geometries. This result is converted to ΔH at 298 K by assuming ideal behavior and adding the translational, rotational, and vibrational energies. Additive correction terms corresponding to the various coordination states of the carbons, nitrogens, and oxygens were developed using a database of 54 compounds. The experimental ΔH values of these compounds are then reproduced with an average absolute error of 3 kcal/mol and a standard deviation of 4 kcal/mol. For a group of 10 test cases that were not part of the database, the average absolute error is 3.5 kcal/mol and the standard deviation is 4.1 kcal/mol. © 1995 by John Wiley & Sons, Inc.     

Improving description of hydrogen bonds at the semiempirical level: water–water interactions as test case

Hydrogen bonding is not well described by available semiempirical theories. This is an important restriction because hydrogen bonds represent a key feature in many chemical and biochemical processes, besides being responsible for the singular properties of water. In this study, we describe a possible solution to this problem. The basic idea is to replace the nonphysical gaussian correction functions (GCF) appearing in the core–core repulsion terms of most MNDO‐based semiempirical methods by a simple function exhibiting the correct physical behavior in the whole range of intermolecular separation distances. The parameterized interaction function (PIF) is the sum of atom‐pair contributions, each one having five adjustable parameters. In this work, the approach is used to study water–water interactions. The parameters are optimized to reproduce a reference ab initio intermolecular energy surface for the water–water dimer obtained at the MP2/aug‐cc‐pVQZ level. OO, OH, and HH parameters are reported for the PM3 method. The results of PM3‐PIF calculations remarkably improve qualitatively and quantitatively those obtained at the standard PM3 level, both for water–dimer properties and for water clusters up to the hexamer. For example, the root‐mean‐square deviation of the PM3‐PIF interaction energies, with respect to ab initio values obtained using 700 points of the water dimer surface, is only 0.47 kcal/mol. This value is much smaller than that obtained using the standard PM3 method (4.2 kcal/mol). The PM3‐PIF water dimer energy minimum (−5.0 kcal/mol) is also much closer to ab initio data (−5.0±0.01 kcal/mol) than PM3 (−3.50 kcal/mol). The method is therefore promising for the development of new semiempirical approaches as well as for application of combined quantum mechanics and molecular mechanics techniques to investigate chemical processes in water. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 572–581, 2000     

Evaluation of MNDO calculated proton affinities

A large set of charged species arising mainly from protonation or deprotonation of hydrocarbons, alcohols, ethers, carboxylic acids, amines, imines, and nitriles has been studied by means of the semiempirical self-consistent-field (SCF ) molecular orbital (MO ) MNDO method. From the calculated heats of formation of such charged species and those of neutral molecules, MNDO -estimated proton affinities have been obtained and the results compared with experimental gas-phase proton affinities. If the small size anions and acetylides, for which the method predicts heats of formation too large, are ruled out, the mean absolute error in calculated proton affinities is ca. 7 kcal/mol for hydrocarbons (22 acid-base pairs) and ca. 8 kcal/mol for oxygen-containing compounds (25 acid-base pairs). For nitrogen-containing molecules it is necessary to discard, in addition, the values corresponding to the protonation of alkylamines and imines in order to achieve a reasonable mean absolute error of 7–8 kcal/mol.     

Density functional study of chlorine–oxygen compounds related to the ClO self-reaction

Geometrical parameters, vibrational frequencies, relative stabilities, and dissociation energies of the three stable Cl₂O₂ isomers and the OClO and ClOO radicals were investigated by density functional theory (DFT). The present analysis shows that DFT using hybrid functionals is capable of describing these systems to at least the same degree of accuracy as ab initio methods. The average absolute bond-length deviation of ClClO₂, ClOOCl, and ClO₂ from experimental results is 0.024/0.027 Å, with a maximum deviation for the dichlorine peroxide O(SINGLE BOND)O bond equal to 0.072/0.063 Å, for the B3PW91 and B3LYP functionals, respectively. The average absolute bond-angle deviation for the hybrid functionals is 0.8°. Harmonic vibrational frequencies calculated with DFT give for all Cl(SINGLE BOND)O compounds good agreement with experiments. The dissociation energies of ClOOCl, OClO, and ClOO were found to be in good agreement with experiments, the average error being less than 1.2 kcal/mol. The two isomers chloryl chloride (ClClO₂) and dichlorine peroxide (ClOOCl) were found to be approximately 9 kcal/mol more stable than the chlorine chlorite (ClOClO) isomer. The ClOO isomer is predicted to be 3.0 kcal/mol more stable than OClO, in accordance with the experimental value of 4 kcal/mol. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 203–217, 1998     

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.     

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.     

Ground states of molecules. 53.MNDO calculations for molecules containing chlorine

MNDO calculations of heats of formation, dipole moments, ionization potentials, and structures are reported for a wide range of compounds containing chlorine in its characteristic valence state (Cl^I) and one or more of the elements H, B, Be, C, N, O, and F. The calculated errors in the heats of formation and the dipole moments are not significantly greater than those previously reported for compounds containing no chlorine. First vertical ionization potentials were on average 0.95 eV too high. The ordering of higher cationic states was found to be correct, even for species such as Cl₂O, Cl₂, and HOCl, where ab initio–Koopmans' theorem calculations predict the incorrect ordering. The calculated energies and geometries of compounds such as CIF₃ are qualitatively incorrect, probably because of the lack of 3d atomic orbitals in the orbital basis set.     

Improved prediction of heats of formation of energetic materials using quantum mechanical calculations

We present simple atom and group-equivalent methods that will convert quantum mechanical energies of molecules to gas phase heats of formation of CHNO systems. In addition, we predict heats of sublimation and vaporization derived from information obtained from the quantum-mechanically calculated electrostatic potential of each isolated molecule. The heats of sublimation and vaporization are combined with the aforementioned gas phase heats of formation to produce completely predicted condensed phase heats of formation. These semiempirical computational methods, calibrated using experimental information, were applied to a series of CHNO molecules for which no experimental information was used in the development of the methods. These methods improve upon an earlier effort of Rice et al. [Rice, B. M.; Pai, S. V.; Hare, J. Combust. Flame 1999, 118, 445] through the use of a larger basis set and the application of group equivalents. The root-mean-square deviation (rms) from experiment for the predicted group-equivalent gas phase heats of formation is 3.2 kcal/mol with a maximum deviation of 6.5 kcal/mol. The rms and maximum deviation of the predicted liquid heats of formation are 3.2 and 7.4 kcal/mol, respectively. Finally, the rms and maximum deviation of predicted solid heats of formation are 5.6 and 12.2 kcal/mol, respectively, an improvement in the rms of approximately 40% compared to the earlier Rice et al. predictions using atom equivalents and a smaller basis set (B3LYP/6-31G*).     

The heat of formation of C2F4

The CCSD(T) atomization energies are extrapolated to the complete basis set limit, and are corrected for zero-point energy, spin–orbit, core-valence, and scalar relativistic effects. Our best heats of formation at 298 K for CF₄ and C₂F₄ are −223.1±1.1 and −160.5±1.5 kcal/mol, respectively. The CF₄ value is in excellent agreement with experiment (−223.04±0.18 kcal/mol), while the C₂F₄ result suggests that the experimental value (−157.6±0.6 kcal/mol) has a larger error than believed. Our value for C₂F₄ also shows that the G3 value has the expected error of ±2 kcal/mol.     

Theoretical studies on the heats of formation and the interactions among the difluoroamino groups in polydifluoroaminocubanes

The heats of formation (HOFs) were calculated for a series of polydifluoroaminocubanes by using density functional theory (DFT), Hartree-Fock, and MP2 method with 6-31G basis set as well as semiempirical methods. The cubane skeleton was not broken in the process of designing isodesmic reactions; i.e., the cubane skeleton was chosen for a reference compound. The contribution of difluoroamino group to the heat of formation deviates from group additivity. The semiempirical MO (MNDO, AM1, and PM3) methods did not produce accurate and reliable results for the HOFs of the title compounds. The relationship between HOFs and molecular structures was discussed. It was found that the HOFs decreased dramatically initially and then gradually with each difluoroamino group attached to the cubane skeleton. The distance between difluoroamino groups influences the values of HOFs. The interacting energies of polydifluoroaminocubanes are in the range 14-20 kJ/mol. The interaction of neighbor difluoroamino groups discords with the group additivity. The average interaction energy between the nearest-neighbor NF(2) group in the most stable conformer of octadifluoroaminocubane is 13.94 kJ/mol at the B3LYP/6-31G level. The NF(2) group can rotate freely around the C-N bond. The relative stability of the title compounds was accessed on the basis of the calculated HOFs, the energy gaps between the frontier orbitals, and the bond order of C-NF(2). These results provide basic information for the molecular design of novel high energetic density materials.