Force law

Molecular dynamics consists of the brute-force solution of Newton s equations of motion. It is necessary to encode in the program the potential energy and force law of interaction between molecules the equations of motion are solved numerically, by finite difference techniques. The system evolution corresponds closely to what happens in real life and allows us to calculate dynamical properties, as well as thennodynamic and structural fiinctions. For a range of molecular models, packaged routines are available, either connnercially or tlirough the academic conmuinity. 

Observe that, in principle, it is possible to introduce quaternions in the solution of the free rotational part of a Hamiltonian splitting, although there is no compelling reason to do so, since the rotation matrix is usually a more natural coordinatization in which to describe interbody force laws. 

We are going to be concerned with electrical and magnetic properties in this text, so I had better put on record the fundamental force laws for stationary charges and steady currents. These are as follows. 

This is a much more complicated force law, because the integrations have to be done around the complete electrical circuits A and B. The details of the integration do not matter, the point being this. Because cq and /aq are interrelated, we are... 

One way to do this is afforded by the predictor-corrector method. We ignore terms higher than those shown explicitly, and calculate the predicted terms starting with bP(t). However, this procedure will not give the correct trajectory because we have not included the force law. This is done at the corrector step. We calculate from the new position rP the force at time t + St and hence the correct acceleration a (t -f 5t). This can be compared with the predicted acceleration aP(f -I- St) to estimate the size of the error in the prediction step... 

The application of the L-J-D method to the present problem amounts to the assumption that the average contribution to the potential energy due to the interaction of a solute molecule with any of the elements constituting the wall of its cage can be described by the familiar force law... 

This determines the distance of closest approach in terms of the initial relative velocity, the impact parameter, and the dynamical quantities (masses and force law constants). The equation for the orbit of the relative motion is found from the first of Eqs. (1-8) and (1-9), using the identity (fj6) = (drjdd), as follows ... 

Angle of Deflection for Some Simple Cases.—If, as is often assumed for simplicity, the interparticle force law is given by... 

In the derivation of the Boltzmann equation, it was noted that the distribution function must not change significantly in times of the order of a collision time, nor in distances of the order of the maximum range of the interparticle force. For the usual interatomic force laws (but not the Coulomb force, which is of importance in ionized gases), this distance is less than about 10 T cm the corresponding collision times, which are of the order of the force range divided by a characteristic particle velocity (of the order of 10 cm/sec for hydrogen at 300° C), is about 10 12 seconds. 

While mathematically attractive, this force law is of limited interest physically it represents only the interaction between permanent quadrupoles, and even this with neglect of angles of orientation. However, although the details of the dependence of viscosity upon temperature are affected by the force law used, the general form of the hydrodynamic equation in the Navier-Stokes approximation is not affected. 

Further evaluation of the integrals for the general case of a v-power law will be discussed in Section 1.17 the calculation for v = 5, the Maxwell molecule force law, gives typical results, and will be completed here. 

Case of General Central Force Law.—The evaluation of the collision integrals for tire viscosity, Eq. (1-89), was simplified for... 

A similar approximation scheme, which gives successive approximations, in first order, to afl can be developed the thermal conductivity for the central force law is ... 

For hard sphere collisions, v(v) would be proportional to v9 and the mean free path independent of v A(v) is an equivalent mean free path for a- general force law. Cf. S. Chapman and T. G. Cowling, The Mathematical Theory of Non- Uniform Oases, pp. 91 and 348, Cambridge University Press, 1958. 

Campbell s Theorem, 174 Cartwright, M. L., 388 Caywood, T. E., 313 C-coefficients, 404 formulas for, 406 recursion relations, 406 relation to spherical harmonics, 408 tabulations of, 408 Wigner s formula, 408 Central field Dirac equation in, 629 Central force law... 

Evaluation for case of general central force law, 36 Lorentz approximation, 46 terms, 28,31,33 Collisions... 

Substituting Eq. (12) into Eq. (11) permits us to derive the Hookean spring force law, well-known in the classical theory of rubber elasticity ... 

It is worth recalling that any of the molecular force laws given by Eqs. (13-16) are derived within the framework of the freely-jointed model which considers the polymer chain as completely limp except for the spring force which resists stretching thus f(r) is purely entropic in nature and comes from the flexibility of the joints which permits the existence of a large number of conformations. With rodlike polymers, the statistical number of conformations is reduced to one and f(r) actually vanishes when the chain is in a fully extended state. 

This equation is useful in that it is applicable to any type of force law so long as the range of interaction and the separation are much less than the radius of the sphere. Thus the force to overcome the work of adhesion between a rigid sphere and a flat surface written in terms of the surface energy Ay is ... 

The basic principles are described in many textbooks [24, 26]. They are thus only sketchily presented here. In a conventional classical molecular dynamics calculation, a system of particles is placed within a cell of fixed volume, most frequently cubic in size. A set of velocities is also assigned, usually drawn from a Maxwell-Boltzmann distribution appropriate to the temperature of interest and selected in a way so as to make the net linear momentum zero. The subsequent trajectories of the particles are then calculated using the Newton equations of motion. Employing the finite difference method, this set of differential equations is transformed into a set of algebraic equations, which are solved by computer. The particles are assumed to interact through some prescribed force law. The dispersion, dipole-dipole, and polarization forces are typically included whenever possible, they are taken from the literature. 

Zimm [34] extended the bead-spring model by additionally taking hydrodynamic interactions into account. These interactions lead to changes in the medium velocity in the surroundings of each bead, by beads of the same chain. It is worth noting that neither the Rouse nor the Zimm model predicts a shear rate dependency of rj. Moreover, it is assumed that the beads are jointed by an ideally Hookean spring, i.e. they obey a strictly linear force law. 

First approaches to approximating the relaxation time on the basis of molecular parameters can be traced back to Rouse [33]. The model is based on a number of boundary assumptions (1) the solution is ideally dilute, i.e. intermolecular interactions are negligible (2) hydrodynamic interactions due to disturbance of the medium velocity by segments of the same chain are negligible and (3) the connector tension F(r) obeys an ideal Hookean force law. 

Attempts have been made to put the VSEPR model on a quantitative basis by describing the interaction between electron pairs in terms of force law of the type... 

We also determine both the gravitational force law perpendicular to the Galactic plane and the total surface mass density and thickness of the Galactic disk. The surface mass density of the Galactic disk within 800 pc derived from this analysis is E z <800 pc) = 76 Msun pc 2. 

In order to determine the force between plates as a function of their separation, one would have to perform a series of simulations with different wall separations and with the chemical potential of the fluid fixed at the bulk value. This is technically feasible, but very computationally intensive [42]. The qualitative behavior of the force law can, however, be estimated from the density profile of a fluid at a single wall using the wall sum rule and a superposition approximation [31,43]. The basic idea is that the density profile [denoted pH(z)] of a fluid between two walls at a separation H can be obtained from the density profile [denoted pj (z) of the same fluid at a single wall using... 

In this equation, cfo and q2 are the charges, r is the distance separating them, and e is the dielectric constant, which has a value of 1 for a vacuum or free space. This force law has no directional component 211... 

---