Equation of motion technique

We now use the equation of motion technique (EOM) to obtain an expression for the GF Gj (t). We first remark that in the time domain two EOM can be written, depending on which time argument in the double-time GF the time derivative will act. One thus obtains in general ... 

Experimentally this has been found in dhcp Pr (Houmann et al. 1979, Jensen 1976) and PrAlj (Purwins et al. 1976). Various theoretical methods have been used to describe mixed-mode excitations pseudospon techniques, (Elliott et al. 1972), diagrammatic methods (Thalmeier and Fulde 1975) and equation of motion techniques (Morin et al. 1980b, Aksenov et al. 1981). In the RPA all these methods lead to mixed excitations of infinite lifetime. Lifetime effects can be included in a phenomenological manner, see, e.g., Becker et al. (1981). For a microscopic treatment of linewidth effects beyond RPA the projection operator technique is more suitable (Huggins et al. 1984). 

This edition, a completely rewritten and expanded version of the original, includes second quantization and diagrammatic perturbation theory, symmetric and unitary group methods, new forms of valence bond theory, dynamic properties and response, propagator and equation-of-motion techniques and the theory of intermolecular forces. 

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. 

The solution Xh(t) of the linearized equations of motion can be solved by standard NM techniques or, alternatively, by explicit integration. We have experimented with both and found the second approach to be far more efficient and to work equally well. Its handling of the random force discretization is also more straightforward (see below). For completeness, we describe both approaches here. 

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]. 

For a given potential energy function, one may take a variety of approaches to study the dynamics of macromolecules. The most exact and detailed information is provided by MD simulations in which one solves the equations of motion for the atoms constituting the macromolecule and any surrounding environment. With currently available techniques and methods it is possible... 

The simplest of the numerical techniques for the integration of equations of motion is leapfrog-Verlet algorithm (LFV), which is known to be symplectic and of second order. The name leapfrog steams from the fact that coordinates and velocities are calculated at different times. 

In this paper, we focus on numerical techniques for integrating the QCMD equations of motion. The aim of the paper is to systematize the discussion concerning numerical integrators for QCMD by ... 

A molecular dynamics simulation samples the phase space of a molecule (defined by the position of the atoms and their velocities) by integrating Newton s equations of motion. Because MD accounts for thermal motion, the molecules simulated may possess enough thermal energy to overcome potential barriers, which makes the technique suitable in principle for conformational analysis of especially large molecules. In the case of small molecules, other techniques such as systematic, random. Genetic Algorithm-based, or Monte Carlo searches may be better suited for effectively sampling conformational space. 

Elimination of the pressure term from the equation of motion does not automatically yield a robust scheme for incompressible flow and it is still necessary to satisfy the BB stability condition by a suitable technique in both forms of the penalty method. 

In the continuous penalty technique prior to the discretization of the governing equations, the pressure in the equation of motion is substituted from Fquation (3.6) to obtain... 

The use of selectively reduced integration to obtain accurate non-trivial solutions for incompressible flow problems by the continuous penalty method is not robust and failure may occur. An alternative method called the discrete penalty technique was therefore developed. In this technique separate discretizations for the equation of motion and the penalty relationship (3.6) are first obtained and then the pressure in the equation of motion is substituted using these discretized forms. Finite elements used in conjunction with the discrete penalty scheme must provide appropriate interpolation orders for velocity and pressure to satisfy the BB condition. This is in contrast to the continuous penalty method in which the satisfaction of the stability condition is achieved indirectly through... 

As described in the discrete penalty technique subsection in Chapter 3 in the discrete penalty method components of the equation of motion and the penalty relationship (i.e. the modified equation of continuity) are discretized separately and then used to eliminate the pressure term from the equation of motion. In order to illustrate this procedure we consider the following penalty relationship... 

Dissipative particle dynamics (DPD) is a technique for simulating the motion of mesoscale beads. The technique is superficially similar to a Brownian dynamics simulation in that it incorporates equations of motion, a dissipative (random) force, and a viscous drag between moving beads. However, the simulation uses a modified velocity Verlet algorithm to ensure that total momentum and force symmetries are conserved. This results in a simulation that obeys the Navier-Stokes equations and can thus predict flow. In order to set up these equations, there must be parameters to describe the interaction between beads, dissipative force, and drag. 

A variety of techniques have been detailed for handling Newton s equations of motion, equation 5 (86—90). Integration techniques yield atomic... 

Langevin dynamics a technique to reduce the total number of equations of motion that are solved. Utilize the Coupled Heat Bath, wherein the method models the solvent effect by incorporating a friction constant into the overall expression for the force. 

In some cases, however, it is possible, by analysing the equations of motion, to determine the criteria by which one flow pattern becomes unstable in favor of another. The mathematical technique used most often is linearised stabiHty analysis, which starts from a known solution to the equations and then determines whether a small perturbation superimposed on this solution grows or decays as time passes. 

A variety of techniques have been introduced to increase the time step in molecular dynamics simulations in an attempt to surmount the strict time step limits in MD simulations so that long time scale simulations can be routinely undertaken. One such technique is to solve the equations of motion in the internal degree of freedom, so that bond stretching and angle bending can be treated as rigid. This technique is discussed in Chapter 6 of this book. Herein, a brief overview is presented of two approaches, constrained dynamics and multiple time step dynamics. 

There are basically two different computer simulation techniques known as molecular dynamics (MD) and Monte Carlo (MC) simulation. In MD molecular trajectories are computed by solving an equation of motion for equilibrium or nonequilibrium situations. Since the MD time scale is a physical one, this method permits investigations of time-dependent phenomena like, for example, transport processes [25,61-63]. In MC, on the other hand, trajectories are generated by a (biased) random walk in configuration space and, therefore, do not per se permit investigations of processes on a physical time scale (with the dynamics of spin lattices as an exception [64]). However, MC has the advantage that it can easily be applied to virtually all statistical-physical ensembles, which is of particular interest in the context of this chapter. On account of limitations of space and because excellent texts exist for the MD method [25,61-63,65], the present discussion will be restricted to the MC technique with particular emphasis on mixed stress-strain ensembles. 

By far the most common methods of studying aqueous interfaces by simulations are the Metropolis Monte Carlo (MC) technique and the classical molecular dynamics (MD) techniques. They will not be described here in detail, because several excellent textbooks and proceedings volumes (e.g., [2-8]) on the subject are available. In brief, the stochastic MC technique generates microscopic configurations of the system in the canonical (NYT) ensemble the deterministic MD method solves Newton s equations of motion and generates a time-correlated sequence of configurations in the microcanonical (NVE) ensemble. Structural and thermodynamic properties are accessible by both methods the MD method provides additional information about the microscopic dynamics of the system. 

In our calculations we make use of several standard techniques of molecular dynamics simulations. The integration of the equations of motions is done by the velocity form of the Verlet-algorithm with a time step of 1.5 The temperature is controlled... 

There are several attractive features of such a mesoscopic description. Because the dynamics is simple, it is both easy and efficient to simulate. The equations of motion are easily written and the techniques of nonequilibriun statistical mechanics can be used to derive macroscopic laws and correlation function expressions for the transport properties. Accurate analytical expressions for the transport coefficient can be derived. The mesoscopic description can be combined with full molecular dynamics in order to describe the properties of solute species, such as polymers or colloids, in solution. Because all of the conservation laws are satisfied, hydrodynamic interactions, which play an important role in the dynamical properties of such systems, are automatically taken into account. 

In the molecular dynamics (MD) [1, 2] technique, a system of particles evolves in time according to the equation of motion, E = nijXj, where L, is the net force acting on particle i, and m, and 3q are the mass and acceleration of particle i, respectively. In a molecular system, typical bond lengths are of the order of angstroms while bond vibrations take place at the time scale of 10-13 s. Therefore, the equations of motion for atoms have to be integrated with time steps on the order of 10 15 s. However, many important chemical and biological phenomena in macromolecules take place at much larger time scales, as shown in Table 8-1. 

So far, the solvent coordinate has not been defined. As noted at the beginning of this Section, the time dependent friction is to be found for the reacting solute fixed at the transition state value x of x. By (3.14), its dynamics were related to those of an (unspecified) solvent coordinate. v. One strategy to identify the solvent coordinate, its frequency, friction, etc., would be to derive an equation of motion for the relevant fluctuating force SF there. To this end, one can use a double-membered projection technique in terms of 8F and 8F. In particular, we define the projection operator... 

---