Potentials functions

The surface free energy can be regarded as the work of bringing a molecule from the interior of a liquid to the surface, and that this work arises from the fact that, although a molecule experiences no net forces while in the interior of the bulk phase, these forces become unbalanced as it moves toward the surface. As discussed in connection with Eq. Ill-IS and also in the next sections, a knowledge of the potential function for the interaction between molecules allows a calculation of the total surface energy if this can be written as a function of temperature, the surface free energy is also calculable. 

Another statistical mechanical approach makes use of the radial distribution function g(r), which gives the probability of finding a molecule at a distance r from a given one. This function may be obtained experimentally from x-ray or neutron scattering on a liquid or from computer simulation or statistical mechanical theories for model potential energies [56]. Kirkwood and Buff [38] showed that for a given potential function, U(r)... 

Two simulation methods—Monte Carlo and molecular dynamics—allow calculation of the density profile and pressure difference of Eq. III-44 across the vapor-liquid interface [64, 65]. In the former method, the initial system consists of N molecules in assumed positions. An intermolecule potential function is chosen, such as the Lennard-Jones potential, and the positions are randomly varied until the energy of the system is at a minimum. The resulting configuration is taken to be the equilibrium one. In the molecular dynamics approach, the N molecules are given initial positions and velocities and the equations of motion are solved to follow the ensuing collisions until the set shows constant time-average thermodynamic properties. Both methods are computer intensive yet widely used. 

There is, of course, a mass of rather direct evidence on orientation at the liquid-vapor interface, much of which is at least implicit in this chapter and in Chapter IV. The methods of statistical mechanics are applicable to the calculation of surface orientation of assymmetric molecules, usually by introducing an angular dependence to the inter-molecular potential function (see Refs. 67, 68, 77 as examples). Widom has applied a mean-held approximation to a lattice model to predict the tendency of AB molecules to adsorb and orient perpendicular to the interface between phases of AA and BB [78]. In the case of water, a molecular dynamics calculation concluded that the surface dipole density corresponded to a tendency for surface-OH groups to point toward the vapor phase [79]. 

The theoretical treatments of Section III-2B have been used to calculate interfacial tensions of solutions using suitable interaction potential functions. Thus Gubbins and co-workers [88] report a molecular dynamics calculation of the surface tension of a solution of A and B molecules obeying Eq. III-46 with o,bb/ o,aa = 0.4 and... 

One fascinating feature of the physical chemistry of surfaces is the direct influence of intermolecular forces on interfacial phenomena. The calculation of surface tension in section III-2B, for example, is based on the Lennard-Jones potential function illustrated in Fig. III-6. The wide use of this model potential is based in physical analysis of intermolecular forces that we summarize in this chapter. In this chapter, we briefly discuss the fundamental electromagnetic forces. The electrostatic forces between charged species are covered in Chapter V. 

Among several important developments, the potential function has been refined and one commonly used function takes the form [54]... 

The uncertainties in choice of potential function and in how to approximate the surface distortion contribution combine to make the calculated surface energies of ionic crystals rather uncertain. Some results are given in Table VII-2, but comparison between the various references cited will yield major discrepancies. Experimental verification is difficult (see Section VII-5). Qualitatively, one expects the surface energy of a solid to be distinctly higher than the surface tension of the liquid and, for example, the value of 212 ergs/cm for (100)... 

A fundamental approach by Steele [8] treats monolayer adsorption in terms of interatomic potential functions, and includes pair and higher order interactions. Young and Crowell [11] and Honig [20] give additional details on the general subject a recent treatment is by Rybolt [21]. 

A monolayer can be regarded as a special case in which the potential is a square well however, the potential well may take other forms. Of particular interest now is the case of multilayer adsorption, and a reasonable assumption is that the principal interaction between the solid and the adsorbate is of the dispersion type, so that for a plane solid surface the potential should decrease with the inverse cube of the distance (see Section VI-3A). To avoid having an infinite potential at the surface, the potential function may be written... 

Equation XVII-78 turns out to ht type II adsorption isotherms quite well—generally better than does the BET equation. Furthermore, the exact form of the potential function is not very critical if an inverse square dependence is used, the ht tends to be about as good as with the inverse-cube law, and the equation now resembles that for a condensed him in Table XVII-2. Here again, quite similar equations have resulted from deductions based on rather different models. 

Clearly, it is more desirable somehow to obtain detailed structural information on multilayer films so as perhaps to settle the problem of how properly to construct the potential function. Some attempts have been made to develop statistical mechanical other theoretical treatments of condensed layers in a potential field success has been reasonable (see Refs. 142, 143). 

In the presence of a potential function U(x,y), the density matrix in the high-temperature approximation has the fomi... 

A3.11.2 QUANTUM SCATTERING THEORY FOR A ONEDIMENSIONAL POTENTIAL FUNCTION... 

In this chapter, we look at the techniques known as direct, or on-the-fly, molecular dynamics and their application to non-adiabatic processes in photochemistry. In contrast to standard techniques that require a predefined potential energy surface (PES) over which the nuclei move, the PES is provided here by explicit evaluation of the electronic wave function for the states of interest. This makes the method very general and powerful, particularly for the study of polyatomic systems where the calculation of a multidimensional potential function is an impossible task. For a recent review of standard non-adiabatic dynamics methods using analytical PES functions see [1]. 

Here, t is the nuclear kinetic energy operator, and so all terms describing the electronic kinetic energy, electron-electron and electron-nuclear interactions, as well as the nuclear-nuclear interaction potential function, are collected together. This sum of terms is often called the clamped nuclei Hamiltonian as it describes the electrons moving around the nuclei at a particular configrrration R. 

In this picture, the nuclei are moving over a PES provided by the function V(R), driven by the nuclear kinetic energy operator, 7. More details on the derivation of this equation and its validity are given in Appendix A. The potential function is provided by the solutions to the electronic Schrddinger equation. 

This situation, despite the fact that reliability is increasing, is very undesirable. A considerable effort will be needed to revise the shape of the potential functions such that transferability is greatly enhanced and the number of atom types can be reduced. After all, there is only one type of carbon it has mass 12 and charge 6 and that is all that matters. What is obviously most needed is to incorporate essential many-body interactions in a proper way. In all present non-polarisable force fields many-body interactions are incorporated in an average way into pair-additive terms. In general, errors in one term are compensated by parameter adjustments in other terms, and the resulting force field is only valid for a limited range of environments. 

Grubmiiller described a method to induce conformational transitions in proteins and derived rate constants for these ([Grubmiiller 1994]). The method employs subsequent modifications of the original potential function based on a principal component analysis of a short MD simulation. It is discussed in more detail in the chapter of Eichinger et al. in this volume. 

By applying a pulling force at a portion of the solute molecule in a specific direction (see chapters of Eichinger et al. and Schulten in this volume), conformational transitions can be induced in specific directions. In order to reconstruct information about the underlying potential function governing protein motion, the irreversible work performed on the system by these forces must be discounted ([Balsera et al. 1997]). 

An ideal empirical potential function on the residue level is a function V that assigns to each sequence-coordinates pair s,x) an energy V s,x) such that... 

While this is disappointing, the nonuniqueness theorem also shows that if some empirical potential is able to predict correct protein folds then many other empirical potentials will do so, too. Thus, the construction of empirical potentials for fold prediction is much less constrained than one might think initially, and one is justified in using additional qualitative theoretical assumptions in the derivation of an appropriate empirical potential function. 

M. Oobatake and G.M. Crippen, Residue-residue potential function for conformational analysis of proteins, J.Phys. Chem. 85 (1981), 1187-1197. 

W, g potential functions, k 1, has been discussed in various papers (see, for example, [6, 11, 9, 16, 3]). It has been pointed out that, for step-sizes /j > e = 1/ /k, the midpoint method can become unstable due to resonances [9, 16], i.e., for specific values of k. However, generic instabilities arise if the step-size k is chosen such that is not small [3, 6, 18], For systems with a rotational symmetry this has been shown rigorously in [6j. This effect is generic for highly oscillatory Hamiltonian systems, as argued for in [3] in terms of decoupling transformations and proved for a linear time varying system without symmetry. 

Among the main theoretical methods of investigation of the dynamic properties of macromolecules are molecular dynamics (MD) simulations and harmonic analysis. MD simulation is a technique in which the classical equation of motion for all atoms of a molecule is integrated over a finite period of time. Harmonic analysis is a direct way of analyzing vibrational motions. Harmonicity of the potential function is a basic assumption in the normal mode approximation used in harmonic analysis. This is known to be inadequate in the case of biological macromolecules, such as proteins, because anharmonic effects, which MD has shown to be important in protein motion, are neglected [1, 2, 3]. 

To determine the vibrational motions of the system, the eigenvalues and eigenvectors of a mass-weighted matrix of the second derivatives of potential function has to be calculated. Using the standard normal mode procedure, the secular equation... 

The accuracy of the CSP approximation is, as test calculations for model. systems show, typically very similar to that of the TDSCF. The reason for this is that for atomic scale masses, the classical mean potentials are very similar to the quantum mechanical ones. CSP may deviate significantly from TDSCF in cases where, e.g., the dynamics is strongly influenced by classically forbidden regions of phase space. However, for simple tunneling cases it seems not hard to fix CSP, by running the classical trajectories slightly above the barrier. In any case, for typical systems the classical estimate for the mean potential functions works extremely well. 

---