Newtonian equations

In molecular dynamics applications there is a growing interest in mixed quantum-classical models various kinds of which have been proposed in the current literature. We will concentrate on two of these models the adiabatic or time-dependent Born-Oppenheimer (BO) model, [8, 13], and the so-called QCMD model. Both models describe most atoms of the molecular system by the means of classical mechanics but an important, small portion of the system by the means of a wavefunction. In the BO model this wavefunction is adiabatically coupled to the classical motion while the QCMD model consists of a singularly perturbed Schrddinger equation nonlinearly coupled to classical Newtonian equations, 2.2. 

Notice that the solution is not identical to J but an approximation of it. The evolution of a and S in time may conveniently be described via the following classical Newtonian equations of motion Given the initial values... 

Various kinds of mixed quantum-classical models have been introduced in the literature. We will concentrate on the so-called quantum-classical molecular dynamics (QCMD) model, which consists of a Schrodinger equation coupled to classical Newtonian equations (cf. Sec. 2). 

Application of the weighted residual method to the solution of incompressible non-Newtonian equations of continuity and motion can be based on a variety of different schemes. Tn what follows general outlines and the formulation of the working equations of these schemes are explained. In these formulations Cauchy s equation of motion, which includes the extra stress derivatives (Equation (1.4)), is used to preseiwe the generality of the derivations. However, velocity and pressure are the only field unknowns which are obtainable from the solution of the equations of continuity and motion. The extra stress in Cauchy s equation of motion is either substituted in terms of velocity gradients or calculated via a viscoelastic constitutive equation in a separate step. 

The temporal behavior of molecules, which are quantum mechanical entities, is best described by the quantum mechanical equation of motion, i.e., the time-dependent Schrdd-inger equation. However, because this equation is extremely difficult to solve for large systems, a simpler classical mechanical description is often used to approximate the motion executed by the molecule s heavy atoms. Thus, in most computational studies of biomolecules, it is the classical mechanics Newtonian equation of motion that is being solved rather than the quantum mechanical equation. 

For the model of free point particles the Newtonian equations present by far the simplest and most efficient analytical fonnalism. In contrast, for chains of rigid bodies, there are several different, but equally applicable, analytical methods in mechanics, with their spe-... 

Newtonian Approach Let Nb(t) and Nyj(t) represent the number of black and white balls at time t, respectively. Let Tib(t) and n (f) be the number of black and white balls having a marked site directly ahead of them at time t. The Newtonian equations of motion are then given by... 

When the fluid behaviour can be described by a power-law, the apparent viscosity for a shear-thinning fluid will be a minimum at the wall where the shear stress is a maximum, and will rise to a theoretical value of infinity at the pipe axis where the shear stress is zero. On the other hand, for a shear-thickening fluid the apparent viscosity will fall to zero at the pipe axis. It is apparent, therefore, that there will be some error in applying the power-law near the pipe axis since all real fluids have a limiting viscosity po at zero shear stress. The procedure is exactly analogous to that used for the Newtonian fluid, except that the power-law relation is used to relate shear stress to shear rate, as opposed to the simple Newtonian equation. 

The second step is the molecular dynamics (MD) calculation that is based on the solution of the Newtonian equations of motion. An arbitrary starting conformation is chosen and the atoms in the molecule can move under the restriction of a certain force field using the thermal energy, distributed via Boltzmann functions to the atoms in the molecule in a stochastic manner. The aim is to find the conformation with minimal energy when the experimental distances and sometimes simultaneously the bond angles as derived from vicinal or direct coupling constants are used as constraints. 

The Lagrangian equations contain nothing more than the original Newtonian equations, but have the advantage that the coordinates may be of any kind whatever. This is of particular importance when analyzing phenomena in which the motion of material particles is not observed directly, such... 

In the absence of friction, there are two forces acting on the mass m whose position vector at time t is denoted by the vector r[r] measured relative to the support point, which is the origin of a set of Cartesian axes with three-component k in the upward vertical direction. The first is the force of gravity on the mass, which acts downwards with a value —mgk. The second is the centripetal force, unknown for the moment, which is directed along the support towards the universal point. We denote this force by — Tr t, where Tis a scalar function of time to be found. The Newtonian equations of motion can then be written as... 

In quantum mechanics, the spatial variables are constituted by generalized coordinates (, ), which replace the individual Cartesian coordinates of all single particles in the set. The Lagrangian equations of motion are the Newtonian equations transposed to the generalized coordinate system. 

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 N 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 Newtonian 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. 

However, before embarking on this analysis, a brief excursion is of interest [285]. Let us specifically ignore the spatial dependence of the distribution of the N particles. It could have been calculated from the deterministic Newtonian equations of motion. Now considering, in particular, the motion of particles 1 and 2 as above, average the distribution over all position of the N particles to give the velocity distribution function... 

Equation (7) is the famous Hellmann-Feynman theorem which allows the full set of quantum-many-body forces to be calculated which can then be used to optimize the atomic geometry or to study the dynamics of the atoms by integrating the Newtonian equations of motion,... 

In 1985, Car and Parrinello introduced a very useful concept that the optimization with respect to the ionic and electronic degrees of freedom need not be done separately. In this concept or theory, one can use a coupled set of pseudo-Newtonian equations of... 

For almost steady flows one can expand yl1 or y about t t and obtain second-order fluid constitution equations in the co-deforming frame. When steady shear flows are considered, the CEF equation is obtained, which, in turn, reduces to the GNF equation for T i = 2 = 0 and to a Newtonian equation if, additionally, the viscosity is constant. 

In principle, one can write down all of these forces and formulate the Newtonian equations of motion for the fluid this yields a complicated differential equation known as the Navier-Stokes equation [1-3]. A complete solution of the Navier-Stokes equation gives the exact trajectory and velocity of each fluid element. In practice, the calculations are often difficult because one must simultaneously account for all fluid elements and the interactions between these elements caused by the viscous drag forces. (The simultaneous motion of many interacting fluid elements is analogous to the simultaneous motion of many interacting mechanical objects, the latter being so complicated that it is described as the many body problem. ) However, in certain cases, the Navier-Stokes equation is reduced to a tractable form by the existence of steady low-velocity flow and high symmetry in the flow conduit (e.g., capillary tubes of circular cross section). We will examine such simple cases shortly. 

A simple equation that shows much of viscoelastic behaviour is due to Maxwell (Figure C4-5). You see that it is the same as the Newtonian equation, but with one additional stress term containing the relaxation time. We can see two extreme cases ... 

Differential Equations for Fluidized Bed Gasifier Model. In a hydrodynamical sense, the processes in fluidized bed gasifiers involve the interaction of a system of particles with flowing gas. The motion of these particles and gas is, at least in principle, completely described by the Navier-Stokes equations for the gas and by the Newtonian equations of motion for the particles. Solution of these equations together with... 

Equation (47) has been derived previously by the rheologists from the so-ealled power law G = hr" which assumes a constant n at any chosen applied pressure (k being a constant). This derivation is not exact since n varies with the position of the fluid inside the capillary. It is unity at the center of the capillary where the flow is Newtonian and becomes greater than one toward the wall when the flow is non-Newtonian. Equation (47) can now be derived without this unrealistic assumption. 

When we evaluate the Green-Kubo relations for the transport coefficients we solve the equations of motion for the molecules. They are often modelled as rigid bodies. Therefore we review some of basic definition of rigid body dynamics [10]. The centres of mass of the molecules evolve according to the ordinary Newtonian equations of motion. The motion in angular space is more complicated. Three independent coordinates a, =(a,a,-2, ,3), i = 1, 2,. . N where N is the number of molecules, are needed to describe the orientation of a rigid body. (Note that a, is not a vector because it does not transform like a vector when the coordinate system is rotated.) The rate of change of a, is... 

For later understanding and interpretation of the dynamics of the ion (see in particular Sect. 3.3) it is also instructive to examine the Newtonian equations of motion. They take on the following appearance ... 

---