Condensed-phase optimized molecular

Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies. 

COMPASS condensed phase optimized molecular potentials for atomistic simula-... 

A computational model based on molecular dynamics was developed to predict the miscibility of indomethacin in the carriers polyethylene oxide (PEO), glucose, and sucrose (Gupta et al. 2011). The cohesive energy density and the solubility parameters were determined by simulations using the condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field. The simulations predicted miscibility for indomethacin with PEO (A5 < 2), borderline miscibility with sucrose (A5 < 7), and immiscibility with glucose (A5 > 10 Table 2.2). 

A force field provides a mathematical description of the potential energy as a function of the atomistic configuration of the polymer. The amorphous packing of the PEEK-WC polymer have been constructed using the force field COMPASS (condensed-phase optimized molecular potentials for atomistic... 

Looking at the structure of HMR, it is clear that it occupies several sites on the surface. Some optimized adsorption structures on Pd(lll) surface modeled by molecular mechanics (MM) with the COMPASS (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies) force field are presented in Fig. 7.9. As the force fields cannot properly model chemisorption and, consequently, exact structures of the adsorbed molecules, the role of the metal surface in the MM calculations was mainly to act as a steric constraint. Depending on the adsorption structure, HMR covers 10-20 Pd surface atoms. 

The reason that atom-atom potentials are so popular, especially in the study of condensed phases [34] and more complex Van der Waals molecules [35], is that they contain few parameters and can be cheaply calculated, while they still describe (implicitly) the anisotropy of the intermolecular potential and they even model its dependence on the internal molecular coordinates. Moreover, they are often believed to be transferable, which implies that the same atom-atom interaction parameters in Eq. (6) can be used for the same types of atoms in different molecules. One should realize, however, that the accuracy of atom-atom potentials is limited by Eq. (5). Further inaccuracies are introduced when the atom-atom interaction parameters in Eq. (6) are transferred from one molecular environment to another. Furthermore, Eq. (6) does not include a term which represents the induction interactions and there is the intrinsic problem that these interactions are inherently not pairwise additive (see Sect. 1.4). Numerical experimentation on the C2H4-C2H4 and N2-N2 potentials, for example, has taught us [31, 33] that even when sufficient ab initio data are available, so that the terms in Eq. (6) can be fitted individually to the corresponding ab initio contributions and, moreover, the positions of the force centers for each term can be optimized, the average error in the best fit of each contribution still remains about 10%. Since the different contributions to the potential partly cancel each other... 

An important subset of many-body potentials shown to be important for simulating interfacial systems are those referred to as polarizable force fields.Various aspects of polarizable force fields, especially for use in biomolecular modeling, is explained by Ren et al. in Chapter 3 of this volume. If one treats the fixed charges in Eq. [3] as parameters to be fitted to obtain the best agreement of the condensed phase simulations with experiments, in many cases one finds that the optimal values are considerably different from those obtained from a fit to a molecular (gas phase) dipole moment or from quantum calculations on isolated molecules. This is because in a condensed medium, the local electric field E, (at the location of a particle i) is determined by all the fixed charges and by all the induced dipoles in the system ... 

---