Intermolecular potentials, molecular modelling

Williams, D. E. J. Chem. Phys. 1967,47,4680-4684 Wdliams, D. E. Starr, T. L. Calculation of the crystal structures of hydrocarbons by molecular packing analysis, Comp. Chem. 1977,1,173-177 Hsu, L.-Y. Wdliams, D. E. Intermolecular potential-function models for crystalline perchlorohydrocarbons, Acta Cryst. 1980, A36, 277-281 Cox, S. R. Hsu, L.-Y. Williams, D. E. Nonbonded potential function models for crystalline oxohydrocarbons, Acta Cryst. 1981, A37, 293-301 Williams, D. E. Cox, S. R. Nonbonded potentials for aza-hydrocarbons the importance of coulombic interactions, Acta Cryst. 1984, B40, 404 17 Williams, D. E. Houpt, D. J. Fluorine nonbonded potentials derived from crystalline perfluorocarbons, Acta Cryst. 1986, B42,286-295. The Williams force field requires the calculation of atomic point and site charges for each molecule. The hydrocarbon force field was later improved (Williams, D. E. Improved intermolecular force field for crystalline hydrocarbons containing four- or three-coordinates carbon, J. Mol. Struct. 1999,485-486, 321-347)... 

Several colloidal systems, that are of practical importance, contain spherically symmetric particles the size of which changes continuously. Polydisperse fluid mixtures can be described by a continuous probability density of one or more particle attributes, such as particle size. Thus, they may be viewed as containing an infinite number of components. It has been several decades since the introduction of polydispersity as a model for molecular mixtures [73], but only recently has it received widespread attention [74-82]. Initially, work was concentrated on nearly monodisperse mixtures and the polydispersity was accounted for by the construction of perturbation expansions with a pure, monodispersive, component as the reference fluid [77,80]. Subsequently, Kofke and Glandt [79] have obtained the equation of state using a theory based on the distinction of particular species in a polydispersive mixture, not by their intermolecular potentials but by a specific form of the distribution of their chemical potentials. Quite recently, Lado [81,82] has generalized the usual OZ equation to the case of a polydispersive mixture. Recently, the latter theory has been also extended to the case of polydisperse quenched-annealed mixtures [83,84]. As this approach has not been reviewed previously, we shall consider it in some detail. 

This Chapter has outlined several different approaches to the computational determination of solution properties. Two of these address solute-solvent interactions directly, either treating the effects of individual solvent molecules upon the solute explicitly or by means of a reaction field due to a continuum model of the solvent. The other procedures establish correlations between properties of interest and certain features of the solute and/or solvent molecules. There are empirical elements in all of these methods, even the seemingly more rigorous ones, such as the parameters in the molecular dynamics/Monte Carlo intermolecular potentials, Eqs. (16) and (17), or in the continuum model s Gcavitation and Gvdw, Eqs. (40) and (41), etc. 

Spectroscopic measurement. Specifically, if the induced dipole moment and interaction potential are known as functions of the intermolecular separation, molecular orientations, vibrational excitations, etc., an absorption spectrum can in principle be computed potential and dipole surface determine the spectra. With some caution, one may also turn this argument around and argue that the knowledge of the spectra and the interaction potential defines an induced dipole function. While direct inversion procedures for the purpose may be possible, none are presently known and the empirical induced dipole models usually assume an analytical function like Eqs. 4.1 and 4.3, or combinations of Eqs. 4.1 through 4.3, with parameters po, J o, <32, etc., to be chosen such that certain measured spectral moments or profiles are reproduced computationally. 

In this paper, we make use of molecular modelling techniques, particularly the AMI semiempirical molecular orbital method, to study the intermolecular interactions that are important for determining the manner in which crystal formation takes place. We are particularly inter ested in compounds that can potentially exhibit nonlinear optical properties. The calculational techniques are directed towards providing insight into the manner in which the desired nonlinear optical properties can be op timized in the macromolecular crystal state.(1)... 

There are presently two main difficulties which handicap attempts at exact calculation. The first concerns the intermolecular potential, and the hazards of extrapolation from models derived from viscosity measurements have been discussed. Furthermore, such a method is of dubious validity for polyatomic molecules, because the intermolecular repulsive potential will generally appear to become progressively shallower with increasing molecular dimensions if the viscosity data are cast, for example, in the Lennard-Jones form. Energy transfer depends... 

In the main group of molecular models studied here the dielectric response of a linear molecule, characterized by a moment of inertia /, is examined. A molecule librates/rotates in a conservative intermolecular potential. Besides, in the harmonic oscillator model such a response is determined by the elastic constant k and by the masses of two vibrating particles. 

The examples shown above clearly demonstrate the potential of molecular modeling for structure analysis. Modeling helps to solve difficult structural problems and adds new information concerning long range intermolecular interactions. 

Other flexible molecular models of nitromethane were developed by Politzer et al. [131,132]. In these, parameters for classical force fields that describe intramolecular and intermolecular motion are adjusted at intervals during a condensed phase molecular dynamics simulation until experimental properties are reproduced. In their first study, these authors used quantum-mechanically calculated force constants for an isolated nitromethane molecule for the intramolecular interaction terms. Coulombic interactions were treated using partial charges centered on the nuclei of the atoms, and determined from fitting to the quantum mechanical electrostatic potential surrounding the molecule. After an equilibration trajectory in which the final temperature had been scaled to the desired value (300 K), a cluster of nine molecules was selected for a density function calculation from which... 

The above results apply to the ideal gas of molecules. The objects addressed in the context of molecular modeling of complex systems are known in the form of macroscopic samples, mostly in the condensed phase. Thus the intermolecular degrees of freedom significantly contribute to the thermodynamical and other properties due to intermolecular interactions. For taking these latter into account the Monte-Carlo (MC) or molecular dynamics (MD) techniques are applied to model systems containing from hundreds to thousands of molecules and correspondingly tens and hundreds of thousands of atoms. These two approaches represent two more modem contexts where a demand for efficient methods of calculation of molecular potential... 

As pointed out earlier, the present treatment attempts to clarify the connection between the sticking probability and the mutual forces of interaction between particles. The van der Waals attraction and Bom repulsion forces are included in the analysis of the relative motion between two electrically neutral aerosol particles. The overall interaction potential between two particles is calculated through the integration of the intermolecular potential, modelled as the Lennard-Jones 6-12 potential, under the assumption of pairwise additivity. The expression for the overall interaction potential in terms of the Hamaker constant and the molecular diameter can be found in Appendix I of (1). The Brownian motions of the two particles are no longer independent because of the interaction force between the two. It is, therefore, necessary to describe the relative motion between the two particles in order to predict the rate of collision and of subsequent coagulation. 

A key question about the use of any molecular theory or computer simulation is whether the intermolecular potential model is sufficiently accurate for the particular application of interest. For such simple fluids as argon or methane, we have accurate pair potentials with which we can calculate a wide variety of physical properties with good accuracy. For more complex polyatomic molecules, two approaches exist. The first is a full ab initio molecular orbital calculation based on a solution to the Schrddinger equation, and the second is the semiempirical method, in which a combination of approximate quantum mechanical results and experimental data (second virial coefficients, scattering, transport coefficients, solid properties, etc.) is used to arrive at an approximate and simple expression. 

Carbon tetrafluoride. Carbon tetra-fluoride, which undergoes a transition to a plastically crystalline (orientationally disordered) phase, has been investigated by the Parrinello-Rahman molecular dynamics method under constant-pressure conditions (6). A simple intermolecular potential model of the Lennard-Jones form was derived by taking into account the experimen-... 

---