Force field for computation of conformational energies,structures, and vibrational frequencies of aromatic polyesters

A CFF93¹ all-atom force field for aromatic polyesters based on ab initio calculations is reported. The force field parameters are derived by fitting to quantum mechanical data which include total energies, first and second derivatives of the total energies, and electrostatic potentials. The valence parameters and the ab initio electrostatic potential (ESP) derived charges are then scaled to correct the systematic errors originating from the truncation of the basis functions and the neglect of electron correlation in the HF/6-31G* calculations. Based on the force field, molecular mechanics calculations are performed for homologues of poly(p-hydroxybenzoic acid) (PHBA) and poly(ethylene terephthalate) (PET). The force field results are compared with available experimental data and the ab initio results. © 1994 by John Wiley & Sons, Inc.     

On the use of conformationally dependent geometry trends from ab initio dipeptide studies to refine potentials for the empirical force field CHARMM

Previous 4-21G ab initio geometry optimizations of various conformations of the model dipeptides (N-acetyl N'methyl amides) of glycine (GLY) and the alanine (ALA) have been used to help refine the empirical force constants and equilibrium geometry in the CHARMM force field for peptides. Conformationally dependent geometry trends from ab initio calculations and positions of energy minima on the ab initio energy surfaces have been used as guides in the parameter refinement, leading to modifications in the bond stretch, angle bending, and some torsional parameters. Preliminary results obtained with these refined empirical parameters are presented for the protein Crambin. Results for the cyclic (Ala-Pro-DPhe)₂ are compared with those from other calculations. It seems that the dihedral angle fit achieved by the new parameters is significantly improved compared with results from force fields whose derivation does not include ab initio geometry trends.     

A charge-conserving approximation method for ab initio calculations

A new method is presented for approximate ab initio calculations in quantum chemistry. It is called CCAM (charge conserving approximation method). The calculation method does not include the use of empirical parameters. We use Slater type orbitals as basis set, replacing STO's by STO-2G functions to evaluate three- and four-center integrals and making the STO-2G two-orbital charge distributions have the same total charge as STO. The results are presented for test calculations on five molecules. In view of these results, CCAM is better than ab initio calculations over STO-6G in the results on total energies, kinetic energies and occupied orbital energies. In atomic populations, dipole moments and unoccupied orbital energies, CCAM is also satisfactory. We estimate that CCAM would be as fast as ab initio calculations over STO-2G in evaluating molecular integrals.     

Interaction between polymers: Dispersion energy of two infinite linear chains

Using the newly developed discrete variant of the k · p perturbation theory for linear chains an expression is derived for the dispersion interaction between two infinite chains valid for both insulating and metallic polymers. It is shown that the formally divergent dispersion energy is also finite for metallic systems, but it may be considerably enhanced in polymers with partly filled bands as compared with insulting ones. The accuracy of the proposed k · p expansion is investigated on model chains through comparisons with direct ab initio Hartree–Fock crystal orbital calculations. Applying second order perturbation theory an explicit formula in terms of band parameters is given for the numerical calculation of the dispersion energy.     

QuickFF: A program for a quick and easy derivation of force fields for metal‐organic frameworks from ab initio input

QuickFF is a software package to derive accurate force fields for isolated and complex molecular systems in a quick and easy manner. Apart from its general applicability, the program has been designed to generate force fields for metal‐organic frameworks in an automated fashion. The force field parameters for the covalent interaction are derived from ab initio data. The mathematical expression of the covalent energy is kept simple to ensure robustness and to avoid fitting deficiencies as much as possible. The user needs to produce an equilibrium structure and a Hessian matrix for one or more building units. Afterward, a force field is generated for the system using a three‐step method implemented in QuickFF. The first two steps of the methodology are designed to minimize correlations among the force field parameters. In the last step, the parameters are refined by imposing the force field parameters to reproduce the ab initio Hessian matrix in Cartesian coordinate space as accurate as possible. The method is applied on a set of 1000 organic molecules to show the easiness of the software protocol. To illustrate its application to metal‐organic frameworks (MOFs), QuickFF is used to determine force fields for MIL‐53(Al) and MOF‐5. For both materials, accurate force fields were already generated in literature but they requested a lot of manual interventions. QuickFF is a tool that can easily be used by anyone with a basic knowledge of performing ab initio calculations. As a result, accurate force fields are generated with minimal effort. © 2015 Wiley Periodicals, Inc.     

Heuristic intermolecular potential function for the methane–water interaction based on ab initio quantum-mechanical calculations

An analytical potential function for the pairwise interaction of methane and water is reported. The function is representative of 225 ab initio quantum-mechanical calculations of the intermolecular interaction using 6–31G self-consistent-field molecular-orbital theory. The statistical parameters of the curve fitting are given and isoenergy contour maps of the interaction energy are presented and discussed.     

Modelisation of the Electric Field Linear Response of an Infinite Aggregate of All-Trans Hexatriene Chains by Electrostatic Interaction Model

A classical electrostatic polarization scheme using the additive distribution procedure has been applied to determine the longitudinal polarizability of an all-trans hexatriene molecule in an infinite stretched fiber. The parameters have been derived from ab initio CPHF/6-31G calculations and the electrostatic scheme has been validated via comparison with ab initio results on small clusters. Upon packing the polarizability of all-trans hexatriene increases by 7%. This small increase results from the balance between the enhancement of the polarizability due to collinear packing and the reduction associated with lateral packing.     

Accurate molecular electrostatic potentials based on modified PRDDO/M wave functions. I. Electrostatic potential derived atomic charges

A new approach for the calculation of electrostatic potential derived atomic charges is presented. Based on molecular orbital calculations in the PRDDO/M approximation, the new parametrized electrostatic potential (PESP) method is parametrized against ab initio MP2/6-31G** calculations. For a data set of 820 atoms in 145 molecules containing H, C, N. O, F, P, S, Cl, and Br (including hypervalent species), the PESP method achieves a mean absolute error of 0.037 e⁻ with a correlation coefficient of 0.990. Unlike other approximate approaches, no scaling factor is required to improve the agreement between PESP charges and the underlying ab initio results. PESP calculations are an order of magnitude faster than the simplest ab initio calculation (STO-3G) on large molecules while achieving a level of accuracy that rivals much more elaborate ab initio methods. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18: 955–969, 1997     

Conformational preferences for hydroxyl groups in substituted tetrahydropyrans

Quantum mechanical (ab initio and semiempirical) and force field calculations are reported for representative torsion potentials in several tetrahydropyran derivatives. The overall agreement between the various methods is quite good except that the AMBER torsion profiles are sensitive to the choice of atomic point charges. Using electrostatic potential (ESP) derived atomic point charges determined with the STO-3G basis set we find that AMBER is able to match the best quantum mechanical results quite well. However, when the point charges are derived using the 6-31G* basis set we find that scaling the intramolecular electrostatic nonbond interactions is necessary. AM1 does not work very well for these compounds when compared to the ab initio methods and, therefore, should only be used in cases when ab initio calculations would be prohibitive. Based upon our results we feel that any force field that makes use of 6-31G* ESP derived atomic point charges will need to scale intramolecular interactions. Implications of scaling intramolecular interactions to the development of force fields based on 6-31G* ESP derived atomic point charges are discussed. © 1992 by John Wiley & Sons, Inc.     

Force field parameters for carbohydrates

A new set of force field parameters for carbohydrates is reported. The parameter set is based on the CHARMM22 force field of Karplus and co-workers. The parameterization is based on newly performed high-level ab initio calculations [MP2/6-311 + G (2d, 2p)/ /6-31G**] of fragment molecules. A good agreement of the modified force field and ab initio data is achieved, which is demonstrated with a variety of molecules. © 1996 by John Wiley & Sons, Inc.     

Bond fukui indices: Comparison of frozen molecular orbital and finite differences through mulliken populations

Bond Fukui functions and matrices are introduced for ab initio levels of theory using a Mulliken atoms in molecules model. It is shown how these indices may be obtained from first‐order density matrix derivatives without need for going to second‐order density matrices as in a previous work. The importance of taking into account the nonorthogonality of the basis in ab initio calculations is shown, contrasting the present results with previous work based on Hückel theory. It is shown how the extension of Fukui functions to Fukui matrices allows getting more insight into the nature of bond Fukui functions. All presently introduced indices respect the necessary normalization conditions and include the classical single atom condensed Fukui functions. © 2013 Wiley Periodicals, 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.     

Derivation of class II force fields. I. Methodology and quantum force field for the alkyl functional group and alkane molecules

A new method for deriving force fields for molecular simulations has been developed. It is based on the derivation and parameterization of analytic representations of the ab initio potential energy surfaces. The general method is presented here and used to derive a quantum mechanical force field (QMFF) for alkanes. It is based on sampling the energy surfaces of 16 representative alkane species. For hydrocarbons, this force field contains 66 force constants and reference values. These were fit to 128,376 quantum mechanical energies and energy derivatives describing the energy surface. The detailed form of the analytic force field expression and the values of all resulting parameters are given. A series of computations is then performed to test the ability of this force field to reproduce the features of the ab initio energy surface in terms of energies as well as the first and second derivatives of the energies with respect to molecular deformations. The fit is shown to be good, with rms energy deviations of less than 7% for all molecules. Also, although only two atom types are employed, the force field accounts for the properties of both highly strained species, such as cyclopropane and methylcyclopropanes, as well as unstrained systems. The information contained in the quantum energy surface indicates that it is significantly anharmonic and that important intramolecular coupling interactions exist between internals. The representation of the nature of these interactions, not present in diagonal, quadratic force fields (Class I force fields), is shown to be important in accounting accurately for molecular energy surfaces. The Class II force field derived from the quantum energy surface is characterized by accounting for these important intramolecular forces. The importance of each 4.2 to 18.2%. This fourfold increase in the second derivative error dramatically demonstrates the importance of bond anharmonicity in the ab initio potential energy surface. The Class II force field derived from the quantum energy surface is characterized by accounting for these important intramolecular forces. The importance of each of the interaction terms of the potential energy function has also been assessed. Bond anharmonicity, angle anharmonicity, and bond/angle, bond/torsion, and angle/angle/ torsion cross-term interactions result in the most significant overall improvement in distorted structure energies and energy derivatives. The implications of each energy term for the development of advanced force fields is discussed. Finally, it is shown that the techniques introduced here for exploring the quantum energy surface can be used to determine the extent of transferability and range of validity of the force field. The latter is of crucial importance in meeting the objective of deriving a force field for use in molecular mechanics and dynamics calculations of a wide range of molecules often containing functional groups in novel environments. © 1994 by John Wiley & Sons, Inc.     

A theoretical study of tautomerism: 2- and 4-oxopyrimidine and some of their derivatives

Here we present a detailed study of the theoretical approach to the tautomerization processes considering 2- and 4-oxo-pyrimidine and several of their derivatives. We look into molecular relaxation, suitability of semiempirical methods, and the effect of basis sets size in pseudopotential ab initio calculations. We found that no semiempirical method is suited for studying the relative stability of the tautomers, even ab initio calculations with minimal basis sets. On the other hand MNDO appears to be very useful for molecular relaxation and can be used, as well as Ab initio calculations, for reasonable estimates of relative tautomerization. We also considered the quality of prediction of other parameters, in particular the ionization potentials. We considered the correlation of experimental and theoretical values as a means to adjust the theoretical results to obtain more reliable predictions.     

Molecular Structure of ortho-Fluoronitrobenzene Studied by Gas Electron Diffraction and Ab Initio MO Calculations

The molecular structure of ortho-fluoronitrobenzene (o-FNB) has been investigated by gas-phase electron diffraction and ab initio MO calculations. The geometrical parameters and force fields of o-FNB were calculated by ab initio and DFT methods. The obtained force fields were used to calculate vibrational amplitudes required as input parameters in an electron diffraction analysis. Within the experimental error limits, the geometrical parameters obtained from the gas-phase electron diffraction analysis are mostly in agreement with the results obtained from the ab initio calculations. The main results are: the molecular geometry of o-FNB is nonplanar with a dihedral angle about C–N of 38(3)°. The r _g (C–F) bond is shortened to 1.307(13) Å in comparison with r _g (C–F) = 1.356(4) Å in C₆H₅F.     

A comparative study of the structure and bonding of HOO,HOS, HSO,and HSS radicals by CNDO/2 and INDO methods

Equilibrium geometries, force constants, barriers to linearity, charge distributions, dipole moments, and electron spin density of HOO, HOS, HSO, and HSS radicals are calculated by CNDO/2 and INDO methods using respectively the original and some recently introduced scheme of parametrization. Three sets of calculations, namely, CNDO/2(sp), CNDO/2(spd), and INDO, are performed, and the results are compared with the ab initio and experimental values, wherever available. A good agreement is obtained for geometry in the case of CNDO/2 (sp) and INDO calculations. The performance of CNDO/2 (spd) calculations in this regard is quite unreliable. The stretching force constants are considerably overestimated by all the methods, while the bending force constants are in reasonable agreement with the ab initio values. With respect to dipole moments, the CNDO/2 values are in better agreement with the ab initio results than the INDO values. In all the cases, the dipole moment vector directions are in complete disagreement with the ab initio predictions.     

Ab initio SCF calculations on low-energy conformers of cyclohexaglycine

Results from ab initio self-consistent field (SCF) calculations with a 3-21G and a double-zeta-plus polarization (DZP) basis set on four low-energy conformations of cyclohexaglycine are reported. In agreement with results from semiempirical and molecular mechanics force field calculations, the lowest-energy conformation found at the DZP level is a conformation forming six C₇ turns. However, the energy difference to the β-turn conformers is significantly smaller at the ab initio DZP level than calculated by the other methods. In contrast to the results obtained with some of the other methods, the present ab initio calculations show that both the double-type-I β turn and the double-type-II β-turn conformer of cyclohexaglycine are stable low-energy structures. © 1995 by John Wiley & Sons, Inc.