Ab initio and molecular mechanics (MM2 and MM3) calculations of positively charged conjugated nitrogen-containing compounds

Ab initio calculations have been carried out on s-trans-N-vinylmethyleneammonium, pyridinium, and related compounds to obtain rotational barriers, structures, and vibrational frequencies. The restricted Hartree-Fock (RHF) level of theory with 6-31G** basis set was used for these calculations. In addition, the MM2(91) and MM3(94) force fields have been parameterized to calculate these positively charged nitrogen-containing compounds. A bond order term was incorporated in the force field to reproduce accurately the rotational barriers of s-trans-N-vinylmethyleneammonium and related compounds. Molecular mechanics geometries and vibrational frequencies compare well with those calculated by ab initio methods. © 1996 by John Wiley & Sons, Inc.     

Molecular mechanics studies (MM4) of sulfides and mercaptans

The MM4 force field has been extended to the title class of compounds. The vibrational spectra, structures, conformational equilibria, and heats of formation have been studied for 47 conformers of 29 compounds. In general, the properties may be calculated with accuracy that is competitive with that for hydrocarbons. The structures are better fit than previously because of the inclusion of a torsion–bend interaction term, which has its origin in the lone pair (Bohlmann) effect. Available experimental data do not suffice to yield detailed torsional potentials, or geometries as a function of torsion angle, and these quantities were determined by ab initio calculations at the MP2/6-31G* level. The rms error in the calculated frequencies of seven representative structures (with a total of 64 experimental and 96 ab initio frequencies) is 25 cm⁻¹. The heats of formation for 23 compounds have a weighted rms error of 0.36 kcal/mol. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1827–1847, 1997     

Molecular mechanics (MM3) parameterization for oxocarbenium ions

The physical properties of a diverse group of 12 oxocarbenium ions have been studied with ab initio calculations at the MP2/6‐31+G* level of theory. Based on theoretically derived properties such as molecular equilibrium geometry, dipole moment, and vibrational frequencies, a molecular mechanics (MM3) force field has been developed with the assistance of the programs TORSMART and MPMSR, components of our artificial parameter development and refinement method. The MM3 force field is now able to reproduce bond lengths, bond angles, moments of inertia, dipole moments, torsional energy profiles, and vibrational frequencies of oxocarbenium ions, which will allow further studies of glycoside hydrolysis and their rates of reaction. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 329–339, 2000     

Molecular mechanics calculations (MM3) on sulfones

The structures of several sulfones, including dimethyl sulfone, methyl ethyl sulfone, methyl vinyl sulfone, and diphenyl sulfone, have been fit with the MM3 force field to existing experimental data from electron diffraction and microwave spectroscopy. The vibrational spectra have also been fit for six of these compounds. The torsional parameters for the aliphatic sulfones were fit to ab initio 6-31G data. Heats of formation were also fit. © 1993 John Wiley & Sons, Inc.     

Molecular mechanics (MM3) calculations on azoxy compounds

The MM3 force field has been extended to include azoxy compounds and also the related amine oxides, both aliphatic and aromatic. The structures of nine molecules were all well fit. The heats of formation for the aliphatic compounds were also well fit, and the vibrational spectra of eight compounds were also fit to the accuracy expected for such calculations. Because many of the experimental data needed to derive the force field were either lacking or were inadequate, ab initio calculations on structures, optimized at the MP2/6-31G* level, were used as needed. © 1994 by John Wiley & Sons, Inc.     

Ab initio calculations and molecular mechanics (MM3) force field development for ammonium and protonated aliphatic amines

The geometries and vibrational frequencies of 11 training molecules containing the ammonium ion moiety were calculated at the MP2/6-31+G* level of theory. Various torsional energy profiles were also calculated using this basis set. From those ab initio calculations, a molecular mechanics (MM3) force field was developed using our Parameter Analysis and Refinement Toolkit System (PARTS). Using this set of parameters, the MM3 force field was found to well reproduce the molecular geometries and vibrational spectra for the all training molecules. CPU time was reduced from days to seconds. The availability of this new force field dramatically increases the feasibility of the computer-assisted drug design involving ammonium and protonated amino groups. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1371–1391, 1997     

Molecular mechanics (MM3) parameterization of hydroxylamine and methyl derivatives

Based on results of electron diffraction, gas phase infrared spectroscopy (IR), and MP2/6-31 + G* ab initio calculations, a set of molecular mechanics (MM3) parameters was developed for molecules containing the N(sp³)—O(sp³) moiety. Using this set of parameters, MM3 is able to reproduce structures (bond lengths and bond angles) and vibrational spectra satisfactorily. © 1994 by John Wiley & Sons, Inc.     

Molecular mechanics (MM3) calculations on alkyl radicals

The MM3 force field has been extended to cover alkyl radicals. Structures, conformational energies, vibrational spectra, and heats of formation have been well fit, mostly to ab initio data. © 1994 by John Wiley & Sons, Inc.     

Unusual conformational-determining interactions in oxymethylpyridines: An ab initio study and an improved method for refining molecular mechanics parameters

The conformational preferences of oxymethylpyridines have been investigated by ab initio calculations and compared to similar calculations for oxymethylbenzene. The C? O bond in the pyridine compounds was found to prefer eclipsing with a C? C bond in the ring, in agreement with previous observations but in disaccord with tentative MM2 calculations. The effect was most pronounced in the 2-substituted pyridine. The benzene compound, on the other hand, showed good agreement between the energies from MM2, MM3, and ab initio calculations. The conformational preferences are discussed in terms of stereoelectronic interactions. New MM2 and MM3 parameters were determined from ab initio calculations on nonstationary points on the energy hypersurface. The parameterization method is discussed. © 1995 by John Wiley & Sons, Inc.     

Molecular mechanics (MM4) calculations on carbonyl compounds part I: aldehydes

Aliphatic aldehydes have been studied with the aid of the MM4 force field. The structures, moments of inertia, vibrational spectra, conformational energies, barriers to internal rotation, and dipole moments have been examined for six compounds (nine conformations). MM4 parameters have been developed to fit the indicated quantities to the wide variety of experimental data. Ab initio (MP2) and density functional theory (B3LYP) calculations have been used to augment and/or replace experimental data, as appropriate. Because more, and to some extent, better, data have become available since MM3 was developed, it was anticipated that the overall accuracy of the information calculated with MM4 would be better than with MM3. The best single measure of the overall accuracy of a force field is the accuracy to which the moments of inertia of a set of compounds (from microwave spectroscopy) can be reproduced. For all of the 20 moments (seven conformations) experimentally known for the aldehyde compounds, the MM4 rms error is 0.30%, while with MM3, the most accurate force field presently available, the rms error over the same set is 1.01%. The calculation of the vibrational spectra was also improved overall. For the four aldehydes that were fully analyzed (over a total of 78 frequencies), the rms errors with MM4 and MM3 are 18 and 38 cm^?1, respectively. These improvements came from several sources, but the major ones were separate parameters involving the carbonyl carbon for formaldehyde, the alkyl aldehydes and the ketones, and new crossterms featured in the MM4 force field that are not present in the MM3 version. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1396–1425, 2001     

Rotational barriers of disilane,hexafluorodisilane, and hexamethyldisilane: Ab initio,density functional,and molecular mechanics (MM3) studies

We investigated structures, vibrational frequencies, and rotational barriers of disilane (Si₂H₆), hexafluorodisilane (Si₂F₆), and hexamethyldisilane (Si₂Me₆) by using ab initio molecular orbital and density functional theories. We employed four different levels of theories (i.e., HF/6–31G*, MP2/6–31G*, BLYP/6–31G*, and B3LYP/6–31G*) to optimize the structures and to calculate the vibrational frequencies (except for Si₂Me₆ at MP2/6–31G*). MP2/6–31G* calculations reproduce experimental bond lengths well, while BLYP/6–31G* calculations largely overestimate some bond lengths. Vibrational frequencies from density functional theories (BLYP/6–31G* and B3LYP/6–31G*) were in reasonably good agreement with experimental values without employing additional correction factors. We calculated the ΔG^‡(298 K) values of the internal rotation by correcting zero-point vibration energies, thermal vibration energies, and entropies. We performed CISD/6–31G*//MP2/6–31G* calculations and found the ΔG^‡(298 K) values for the internal rotation of Si₂H₆, Si₂F₆, and Si₂Me₆ to be 1.36, 2.06, and 2.69 kcal/mol, respectively. The performance of this level was verified by using G2 and G2(MP2) methods in Si₂H₆. According to our theoretical results, the ΔG^‡(298 K) values were marginally greater than the ΔE^‡(0 K) values in Si₂F₆ and Si₂Me₆ due to the contribution of the entropy. In Si₂H₆ the ΔE^‡(0 K) and ΔG^‡(298 K) values were coincidently similar due to a cancellation of two opposing contributions between zero-point and thermal vibrational energies, and entropies. Our calculated ΔG^‡(298 K) values were in good agreement with experimental values published recently. In addition, we also performed MM3 calculations on Si₂H₆ and Si₂Me₆. MM3 calculated rotational barriers and thermodynamic properties were compared with high level ab initio results. Based on this comparison, MM3 calculations reproduced high level ab initio results in rotational barriers and thermodynamic properties of Si₂H₆ derivatives including vibrational energies and entropies, although large errors exist in some vibrational frequencies. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1523–1533, 1997     

The valence isomers of (CH)8 and (SiH)8: An ab initio MO study

The 22 possible valence isomers of the (CH)₈ and (SiH)₈ systems have been studied by ab initio molecular orbital calculations at the MP2/6-31G*//6-31G* + ZPE level. Optimized geometries, relative energies, and, for some selected compounds, vibrational frequencies are reported. The systematic differences between the carbon and silicon compounds are analyzed. © 1994 by John Wiley & Sons, Inc.     

Force field calculations (MM3) on glyoxal,quinones, and related compounds

A general force field type of calculation has been devised in connection with MM3 to treat 1,2- and 1,4-diketones, both when they are not conjugated (as in derivatives of glyoxal) and when they are conjugated (as in derivatives of ortho- and para-benzoquinone). The molecular structures, moments of inertia, dipole moments, and vibrational spectra have been examined for about 15 compounds, some in several conformations. Ab initio calculations (6-31G*) have been used to determine quantities that have not been previously defined by experiment. In general, the force field permits the calculation of the structures with high accuracy, and the spectroscopic and conformational energy data with fair accuracy. © 1994 by John Wiley & Sons, Inc.     

Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide

Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree–Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by “breaking” bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. Results from tests of the method to find relative conformational energies and geometries of alanine dipeptides and alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1468–1494, 1999     

Conformational energies and rotational barriers in 3-methyl-1-butene and 1-butene: An ab initio and molecular mechanics study

The calculated result obtained with MM2(87) for the rotation of the isopropyl group in 3-methyl-1-butene is not in agreement with experimental data. In order to reparametrize the C_sp2-C_sp3-C_sp-C_sp3 torsional angle, 3-methyl-1-butene and 1-butene have been studied by molecular mechanics (MM2(87)) and ab initio (MP2/6-31G* and MP3/6-31G*) calculations. The reparametrization of the torsional angle gives calculated results from MM2(87) in agreement with experimental data and ab initio calculations for both 3-methyl-1-butene and 1-butene. The calculated barriers for the rotation of alkyl groups in alkylbenzenes are improved with these new parameters.     

Computation of the Electronic and Spectroscopic Properties of Carbohydrates Using Novel Density Functional and Vibrational Self-Consistent Field Methods

ABSTRACT

A novel density functional method is presented for the calculation of electronic and thermodynamical properties of oligosaccharides. This method, termed K2-BVWN, offers two advantages; it scales as N ³, where N is the number of basis functions, and there are only two adjustable parameters. The current density functional method is tested in terms of reproducing high level gas phase ab initio calculations in eleven low energy conformers of D-glucopyranose including exo-anomeric and different hydroxymethyl orientations (G ^?, G ⁺, and T). The K2-BVWN method is also tested in terms of reproducing the spectroscopic features of D-glucopyranose and D-mannopyranose (α/β) as compared with both a vibrational self-consistent field calculation (VSCF) as well as experimental infrared spectroscopy. The VSCF calculations offer the advantage that it is possible to include higher order mode coupling and anharmonic effects directly into the calculation of the vibrational frequencies. In general, the K2-BVWN method reproduces the ab initio energetic trends of the different conformers of D-glucose. While the absolute energies are not the same between the ab initio and the K2-BVWN method, both methods do predict a preference for the α-anomer in the gas phase (0.4 kcal/mol ab initio, 0.0 – 0.5 kcal/mol K2-BVWN). The K2-BVWN method was able to reproduce the experimental and VSCF calculated spectrum of both D-glucopyranose and D-mannopyranose in the frequency range between 1500 – 800 cm^?1. Because the current density functional method is both relatively quick and accurate, it represents a significant advancement in the development of oligosaccharide force fields.     

Chlorosilanes: Development and application of MM2 force field parameters

The geometries, relative conformational energies, and dipole moments of mono and polychlorosilanes have been calculated using ab initio molecular orbital (MO) theory. Calculations at the HF/3–21G(*) level, with the exception of dipole moments, give reasonable agreement with experimental data. A new MM2 force field for chlorosilanes, which includes terms for bond length shortening and bond angle compression due to the attachment of electronegative Cl atoms, has been developed on the basis of experimental and ab initio results. The new force field is generally successful in predicting structural parameters, but is unable to reproduce the dipole moments of several model systems. While dipole moment predictions are not the authors' main interest, this failure defines a shortcoming in the MM2 method. The new parameters have been applied to problems in the prediction of stereochemistries of cyclic systems, and compared with experimental results where data are available.     

Infrared and Raman Spectra,Conformational Stability,Ab Initio Calculations,and Vibrational Assignments for Cyclopropyldifluorosilane

The infrared (3200–30 cm^–1) spectra of gaseous and solid Cyclopropyldifluorosilane, c-C₃H₅SiF₂H, and the Raman spectra (3200–20 cm^–1) of the liquid with quantitative depolarization values and the solid have been recorded. Both the syn (cis) and skew (gauche) conformers have been identified in the fluid phases, but only the syn conformer remains in the solid. Variable temperature (–55 to –100°C) studies of the infrared spectra of the sample dissolved in liquid xenon have been carried out. From these data, the enthalpy difference has been determined to be 73 ± 10 cm^–1 (209 ± 29 cal mol^–1), with the syn conformer being the more stable rotamer, which is at variance with the predictions from ab initio calculations. A complete vibrational assignment is proposed for both conformers based on infared band contours, relative intensities, depolarization values, and group frequencies. The vibrational assignments are supported by normal coordinate calculations utilizing the force constants from ab initio MP2/6-31G* calculations. Utilizing the frequencies of the silicon–hydrogen sketch, the rm Si—H bond distances of 1.474 and 1.472 Å have been obtained for the syn and skew conformers, respectively. Complete equilibrium geometries have been determined for both rotamers by ab initio calculations employing the 6-31G* and 6-311 +G** basis sets at levels of restricted Hartree-Fock (RHF) and/or Moller–Plesset (MP) to second order. The potential energy terms for the conformer interconversion have been obtained from the MP2/6-31G* calculation. The results are discussed and compared to those obtained for some similar molecules.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

Optimization of the Lennard-Jones parameters for a combined ab initio quantum mechanical and molecular mechanical potential using the 3-21G basis set

A combined ab initio quantum mechanical and molecular mechanical (AI-QM/MM) potential for use in molecular modeling and simulation has been described. In this article, we summarize a procedure for deriving the empirical parameters embedded in a combined QM/MM model and suggest a set of Lennard-Jones parameters for the combined ab initio 3-21G and MM OPLS-TIP3P (AI-3/MM) potential. Interaction energies and geometrical parameters predicted with the AI-3/MM model for over 80 hydrogen-bonded complexes of organic compounds with water were found to be in good accord with ab initio 6-31G(d) results. We anticipate that the AI-3/MM potential should be reasonable for use in condensed phase simulations. © 1996 by John Wiley & Sons, Inc.