Conformational analysis—131: A molecular mechanics treatment of phosphines

The molecular mechanics method (1973 force field) has been extended to include the general class of phosphaalkanes. Structural data on simple molecules were fit well, and equilibria between conformations can be calculated in general in agreement with experiment, in sofar as experimental data are known. Phosphacyclohexane and a few simple derivatives have been studied in some detail. It is concluded that the force field developed is adequate for the prediction of structures and energy differences, but it is not possible to calculate heats of formation, due to the lack of experimental data.     

Molecular dynamics effects on luminescence properties of oligothiophene derivatives: a molecular mechanics-response theory study based on the CHARMM force field and density functional theory

CHARMM force field parameter values for a class of oligothiophene derivatives have been derived with reference to density functional theory/B3LYP potential energy surfaces. The force field parametrization of these luminescent conjugated polyelectrolytes includes the electronic ground state as well as the strongly light absorbing first excited state. In conjunction with quantum chemical response theory calculations of transition state properties, a molecular dynamical model of the Stokes shift is obtained. The theoretical model is benchmarked against experimental data recorded at room temperature which refer to sodium salts of p-HTAA and p-FTAA with distilled water as a solvent. For p-HTAA the theoretically predicted Stokes shift of 112 nm is in good agreement with the experimental result of 124 nm, given the approximations about exciton localization that were introduced to obtain a force field for the excited state.     

分子力场进展

分子力学（简称ＭＭ）是近年来化学家常用的一种计算方法。与量子力学从头计算和半经验方法相比，用分子力学处理大分子可以大大节省计算时间，而且，在大多数情况下，用分子力学方法计算得到的分子几何构型参数与实验值之间的差值可在实验误差范围之内。所以，分子力学是研究生物化学体系的有效和可行的手段。分子力学的核心是分子力场。本文介绍了分子力场的量子力学背景、分子力场和光谱力场之间的关系。分子力场的一般形式、分力     

Derivation of class II force fields: V. Quantum force field for amides,peptides, and related compounds

As the field of biomolecular structure advances, there is an ever-growing need for accurate modeling of molecular energy surfaces to simulate and predict the properties of these important systems. To address this need, a second generation amide force field for use in simulations of small organics as well as proteins and peptides has been derived. The critical question of what accuracy can be expected from calculations in general, and with this class II force field in particular, is addressed for structural, dynamic, and energetic properties. The force field is derived from a recent methodology we have developed that involves the systematic use of quantum mechanical observables. Systematic ab initio calculations were carried out for numerous configurations of 17 amide and related compounds. Relative energies and first and second derivatives of the energy of 638 structures of these compounds resulted in 140,970 ab initio quantum mechanical observables. The class II peptide quantum mechanical force field (QMFF), containing 732 force constants and reference values, was parameterized against these observables. A major objective of this work is to help establish the role of anharmonicity and coupling in improving the accuracy of molecular force fields, as these terms have not yet become an agreed upon standard in the ever more extensive simulations being used to probe biomolecular properties. This has been addressed by deriving a class I harmonic diagonal force field (HDFF), which was fit to the same energy surface as the QMFF, thus providing an opportunity to quantify the effects of these coupling and anharmonic contributions. Both force field representations are assessed in terms of their ability to fit the observables. They have also been tested by calculating the properties of 11 stationary states of these amide molecules. Optimized structures, vibrational frequencies, and conformational energies obtained from the quantum calculations and from both the QMFF and the HDFF are compared. Several strained and derivatized compounds including urea, formylformamide, and butyrolactam are included in these tests to assess the range of applicability (transferability) of the force fields. It was found that the class II coupled anharmonic force field reproduced the structures, energies, and vibrational frequencies significantly more faithfully than the class I harmonic diagonal force field. An important measure, rms energy deviation, was found to be 1.06 kcal/mol with the class II force field, and 2.30 kcal/mol with the harmonic diagonal force field. These deviations represent the error in relative configurational energy differences for strained and distorted structures calculated with the force fields compared with quantum mechanics. This provides a measure of the accuracy that might be expected in applications where strain may be important such as calculating the energy of a system as it approaches a (rotational) barrier, in ligand binding to a protein, or effects of introducing substituents into a molecule that may induce strain. Similar results were found for structural properties. Protein dynamics is becoming of ever-increasing interest, and, to simulate dynamic properties accurately, the dynamic behavior of model compounds needs to be well accounted for. To this end, the ability of the class I and class II force fields to reproduce the vibrational frequencies obtained from the quantum energy surface was assessed. An rms deviation of 43 cm⁻¹ was achieved with the coupled anharmonic force field, as compared to 105 cm⁻¹ with the harmonic diagonal force field. Thus, the analysis presented here of the class II force field for the amide functional group demonstrates that the incorporation of anharmonicity and coupling terms in the force field significantly improves the accuracy and transferability with regard to the simulation of structural, energetic, and dynamic properties of amides. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 430–458, 1998     

Derivation of force field parameters for SnO2-H2O surface systems from plane-wave density functional theory calculations

Plane-wave density functional theory (DFT-PW) calculations were performed on bulk SnO2 (cassiterite) and the (100), (110), (001), and (101) surfaces with and without H2O present. A classical interatomic force field has been developed to describe bulk SnO2 and SnO2-H2O surface interactions. Periodic density functional theory calculations using the program VASP (Kresse et al., 1996) and molecular cluster calculations using Gaussian 03 (Frisch et al., 2003) were used to derive the parametrization of the force field. The program GULP (Gale, 1997) was used to optimize parameters to reproduce experimental and ab initio results. The experimental crystal structure and elastic constants of SnO2 are reproduced reasonably well with the force field. Furthermore, surface atom relaxations and structures of adsorbed H2O molecules agree well between the ab initio and force field predictions. H2O addition above that required to form a monolayer results in consistent structures between the DFT-PW and classical force field results as well.     

Evaluation of a new force field for describing the adsorption behavior of alkanes in various pure silica zeolites

The recently proposed united atom force field by Dubbeldam et al. (Phys. Rev. Lett. 2004, 93, 088302) for the adsorption of alkanes in MFI-type zeolites was extended to other zeolites in this work. Its applicability to FER-type zeolites was evaluated in detail, for which the Henry coefficients, the isosteric heat of adsorption, the adsorption isotherms, as well as the locations of alkanes in the FER-type zeolites were computed and compared to experimental values. The results show that the new force field works well for FER zeolites. Furthermore, its applicability to MWW-, MTW-, CFI-, LTA-, and STF-type zeolites was investigated, and we found that the experimental isotherms could be accurately predicted except for STF-type zeolites. This work shows that the new united atom force field proposed by Dubbeldam et al. is applicable to most pure silica zeolites.     

Design and application of a molecular mechanics force field for alkyl iodides including an electrostatic polarization model

We describe a force field extending Boyd and Kesner's molecule mechanics treatment of alkyl chlorides, fluorides and bromides to alkyl iodides: accordingly, as well as parameters for the potential energy of bond stretching and bending and for interactions between nonbonded atoms, the field includes parameters taking into account permanent and induced C-I dipoles, whose values have been estimated by fitting the model to experimental dipole moment data. The results of applying this force field to 22 organic iodides comprising several substitution series are reported, and the structural and energy trends in the various series are discussed. These results agree satisfactorily with available experimental data and with those obtained by A. Y. Meyer using a different kind of force field.     

Influence of cation Na/Ca ratio on adsorption in LTA 5A: a systematic molecular simulation study of alkane chain length

Recent adsorption isotherms of n-alkanes on Ca,Na-LTA-type zeolite afford development of a force field describing the interactions between calcium and n-alkanes in configurational-bias Monte Carlo simulations. The force field of Calero et al. (J. Am. Chem. Soc. 2004, 126, 11377-11386) is able to accurately describe the adsorption properties of linear alkanes in the sodium form of FAU-type zeolites. Here, we extend upon this type of force field by including calcium-type ions. The force field was fitted to reproduce the calcium and sodium cations positions on LTA 5A and the experimental adsorption properties of n-alkanes over all range of temperatures and pressures. This opens up a vast amount of experimental data on LTA 5A, both on adsorption and diffusion. Furthermore, evaluation of half a century of reported n-alkane adsorption data on LTA-type zeolites indicates that there are many inconsistencies between the various data sets, possibly as a result of (i) undisclosed calcium and sodium contents, (ii) less than perfect drying of the hygroscopic zeolite, and (iii) coadsorption of contaminants such as vacuum grease. Having obtained our force field, and confirmed its reliability on predictions outside the calibration set, we apply the force field on two "open" problems: (a) the heats of adsorption and Henry coefficient as a function of chain length and (b) the effect of cations in LTA-type zeolites. The molecular simulations shed new light on previous experimental findings, and we provide rationalizations on the molecular level that can be generalized to the class of cage/window-type nanoporous materials.     

适用于TATB,RDX,HMX含能材料的全原子力场的建立与验证

报道一个适用于三种常见的含能材料分子三硝基三氨基苯(TATB),环三亚甲基三硝胺(RDX),环四亚甲基四硝胺(HMX)的全原子力场.力场采用广泛使用的力场函数形式,其中键参数通过拟合量子化学密度泛函计算的数据获得,电荷参数和范德华参数通过拟合相应的分子晶体的物理性质(密度和升华焓)优化得到.通过计算分子和分子晶体的性质显示该力场可以用来准确地预测分子结构、分子振动频率和分子晶体的晶胞参数、密度和升华焓.进一步的验证显示该力场可用来较为准确地预测分子晶体的状态方程和机械模量.     

Ab initio parametrized polarizable force field for rutile-type SnO<Subscript>2</Subscript>

We report a new, polarizable classical force field for the rutile-type phase of SnO₂, casserite. This force field has been parametrized using results from ab initio (density functional theory) calculations as a basis for fitting. The force field was found to provide structural, dynamical and thermodynamic properties of tin oxide that compare well with both ab initio and experimental results at ambient and high pressures.     

Urea: an ab initio and force field study of the gas and solid phases

We have studied the gaseous and solid phases of urea using both quantum mechanics calculation and force field simulation methods. Our ab initio calculations confirmed experimental observations that urea structure is planar in the crystal, but nonplanar in the gas phase. Based on electron structure analysis, we suggest that the significant difference between these two structures in different environments can be qualitatively explained by two resonance structures. The planar structure is more polarized than the nonplanar one, and the former is stabilized in the solid phases due to strong electrostatic interactions. We found classical force field method is incapable to represent such strong polarization effect. Using molecular dynamics simulations with a force field optimized for condensed phases, we calculated the crystalline structures of urea in the temperature range of 12 to 293 K. The densities as well as cell parameters are within 2% deviation from the experimental data in the temperature range.     

Accurate prediction of absolute acidity constants in water with a polarizable force field: substituted phenols, methanol, and imidazole

Validity of a force field with explicit treatment of electrostatic polarization in a form of inducible point dipoles for computing acidity constants was tested by calculating absolute pK(a) values of substituted phenols, methanol, and imidazole in water with the molecular dynamics technique. The last two systems were selected as tyrosine and histidine side-chain analogues, respectively. The solvent was represented by an explicit polarizable water model. Similar calculations were also performed with a modified OPLS-AA nonpolarizable force field. The resulting pK(a) values were compared with available experimental data. While the nonpolarizable force field yields errors of about 5 units in the absolute pK(a) values for the phenols and methanol, the polarizable force field produces the acidity constant values within a ca. 0.8 units accuracy. For the case of imidazole, the fixed-charges force field was capable of reproducing the experimental value of pK(a) (6.4 versus the experimental 7.0 units), but only at a cost of dramatically underestimating dimerization energy for the imidazolium-water complex. At the same time, the polarizable force field yields an even more accurate result of pK(a) = 6.96 without any sacrifice of the accuracy in the dimerization energy. It has also been demonstrated that application of Ewald summation for the long-range electrostatics is important, and substitution of a simple cutoff procedure with Born correction for ions can lead to underestimation of absolute pK(a) values by more than 5 units. The accuracy of the absolute acidity constants computed with the polarizable force field is very encouraging and opens road for further tests on more diverse organic molecules sets, as well as on proteins.     

CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields

The widely used CHARMM additive all‐atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug‐like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug‐like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3‐phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform “all‐CHARMM” simulations on drug‐target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Evaluation of force fields for molecular simulation of polyhedral oligomeric silsesquioxanes

Polyhedral oligomeric silsesquioxanes (POSS) are nanometer-size molecules suitable for the production of organic-inorganic nanocomposite materials. These organic-inorganic nano-building blocks show promise for enabling the production of polymeric materials of exceptional mechanical properties as well as novel composite materials. While the experimental studies of these materials have rapidly evolved in the past decade, their theoretical investigation is still in its infancy. Toward the validation of force fields for the molecular simulation of POSS-containing systems, we present the charge-transfer reactive (CTR) force field for the molecular simulation of polyhedral oligomeric silsesquioxane (POSS) molecules and compare the ability of this, and several force fields taken from the literature, to predict the thermophysical properties of POSS-containing systems. The literature force fields compared include the universal force field (UFF) and the COMPASS and Hybrid-COMPASS force fields. Predictions from molecular dynamics simulations of the structural parameters (unit cell vectors), melting temperature, and FT-IR spectra of crystals of POSS monomers are presented. The POSS monomers investigated are octahydride, octamethyl, and octapropyl POSS. Predicted quantities are compared to experimental results where available and provide molecular-level physical insight into several aspects of the behavior of POSS molecules. While all the force fields tested perform reasonably well, our results indicate that the Hybrid-COMPASS and CTR force fields predict structural properties that are in good agreement with experimental data.     

New force field for molecular simulation of guanidinium-based ionic liquids

An all-atom force field was proposed for a new class of room temperature ionic liquids (RTILs), N,N,N',N'-tetramethylguanidinium (TMG) RTILs. The model is based on the AMBER force field with modifications on several parameters. The refinements include (1) fitting the vibration frequencies for obtaining force coefficients of bonds and angles against the data obtained by ab initio calculations and/or by experiments and (2) fitting the torsion energy profiles of dihedral angles for obtaining torsion parameters against the data obtained by ab initio calculations. To validate the force field, molecular dynamics (MD) simulations at different temperatures were performed for five kinds of RTILs, where TMG acts as a cation and formate, lactate, perchlorate, trifluoroacetate, and trifluoromethylsulfonate act as anions. The predicted densities were in good agreement with the experimental data. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) were investigated to depict the microscopic structures of the RTILs.     

A force field for molecular simulation of tetrabutylphosphonium amino acid ionic liquids

An all-atom force field was set up for a new class of ionic liquids (ILs), tetrabutylphosphonium amino acid, on the basis of the AMBER force field with determining parameters related to the phosphorus atom and modifying several parameters. Ab initio quantum chemical calculations were employed to obtain molecular geometries, infrared frequencies, and torsion energy profiles. Atom partial charges were obtained by using the one-conformation, two-step restraint electrostatic potential approach. Molecular dynamics simulation was carried out in the isothermal-isobaric ensemble for 14 tetrabutylphosphonium amino acid ILs at two temperatures to validate the force field against the experimental densities and heat capacities at constant pressure. Computed thermodynamic properties are in good agreement with available experimental values. Moreover, radial distribution functions were investigated to depict the microscopic structures of these ILs.     

Conformational analysis—129: Heats of formation and thermodynamic parameters for hydrocarbons,calculated by the molecular mechanics method including the effects of molecular vibrations,

A force field is presented in which the zero point energies and the effects of vibrational thermal excitation are included in the heat of formation calculations. The standard deviation of the difference between the calculated and experimental heats of formation of a diverse set of 44 hydrocarbons is 0.38 $\text{kcal}\text{mole}$, which is the same as the average experimental error for this set of compounds. This standard deviation is substantially lower than in previous similar force fields which did not include vibrational effects, and it is argued that much of this improvement is due to the inclusion of vibrational terms in the heat of formation. This force field also reproduces the experimental structures of compounds well.     

Classical Reactive Molecular Dynamics Implementations: State of the Art

Reactive molecular dynamics (RMD) implementations equipped with force field approaches to simulate both the time evolution as well as chemical reactions of a broad class of materials are reviewed herein. We subdivide the RMD approaches developed during the last decade as well as older ones already reviewed in 1995 by Srivastava and Garrison and in 2000 by Brenner into two classes. The methods in the first RMD class rely on the use of a reaction cutoff distance and employ a sudden transition from the educts to the products. Due to their simplicity these methods are well suited to generate equilibrated atomistic or material‐specific coarse‐grained polymer structures. In connection with generic models they offer useful qualitative insight into polymerization reactions. The methods in the second RMD class are based on empirical reactive force fields and implement a smooth and continuous transition from the educts to the products. In this RMD class, the reactive potentials are based on many‐body or bond‐order force fields as well as on empirical standard force fields, such as CHARMM, AMBER or MM3 that are modified to become reactive. The aim with the more sophisticated implementations of the second RMD class is the investigation of the reaction kinetics and mechanisms as well as the evaluation of transition state geometries. Pure or hybrid ab initio, density functional, semi‐empirical, molecular mechanics, and Monte Carlo methods for which no time evolution of the chemical systems is achieved are excluded from the present review. So are molecular dynamics techniques coupled with quantum chemical methods for the treatment of the reactive regions, such as Car–Parinello molecular dynamics.