Molecular dynamics classical

Ciccotti G and Ferrario M 1998 Constrained and nonequilibrium molecular dynamics Classical and Quantum Dynamics In Condensed Phase Simulations ed B J Berne, G Ciccotti and D F Coker (Singapore World Scientific) pp 157-77... 

Tab. 1.1 Comparison of the properties of quantum chemical electronic structure calculations (QC methods), classical molecular dynamics (Classical MD) based on empirical force fields and first-principles molecular dynamics (ab initio MD) simulations.

Molecular dynamics itself can be further divided into two classes classical molecular dynamics and ab initio molecular dynamics. Classical molecular dynamics treats molecules as point masses and the interactions between molecules are represented by simple potential functions, which are based on empirical data or fi om independent quantum mechanical calculations. The so-called ab initio molecular dynamics unifies classical molecular dynamics and density-function theory and takes into account the electronic structure when calculating the forces on atomic nuclei. In this entry, we only present a brief summary of the classical molecular dynamics simulation method. Readers interested in Monte Carlo or ab initio molecular dynamics simulation methods is referred to other entry such as Monte Carlo Method. ... 

Bartell and co-workers have made significant progress by combining electron diffraction studies from beams of molecular clusters with molecular dynamics simulations [14, 51, 52]. Due to their small volumes, deep supercoolings can be attained in cluster beams however, the temperature is not easily controlled. The rapid nucleation that ensues can produce new phases not observed in the bulk [14]. Despite the concern about the appropriateness of the classic model for small clusters, its application appears to be valid in several cases [51]. 

The summation of pair-wise potentials is a good approximation for molecular dynamics calculations for simple classical many-body problems [27], It has been widely used to simulate hyperthennal energy (>1 eV) atom-surface scattering ... 

As in any field, it is usefiil to clarify tenninology. Tliroughout this section an atom more specifically refers to its nuclear centre. Also, for most of the section the /)= 1 convention is used. Finally, it should be noted that in the literature the label quantum molecular dynamics is also sometimes used for a purely classical description of atomic motion under the potential created by tlie electronic distribution. 

The method of molecular dynamics (MD), described earlier in this book, is a powerful approach for simulating the dynamics and predicting the rates of chemical reactions. In the MD approach most commonly used, the potential of interaction is specified between atoms participating in the reaction, and the time evolution of their positions is obtained by solving Hamilton s equations for the classical motions of the nuclei. Because MD simulations of etching reactions must include a significant number of atoms from the substrate as well as the gaseous etchant species, the calculations become computationally intensive, and the time scale of the simulation is limited to the... 

Equation (C3.5.3) shows tire VER lifetime can be detennined if tire quantum mechanical force-correlation Emotion is computed. However, it is at present impossible to compute tliis Emotion accurately for complex systems. It is straightforward to compute tire classical force-correlation Emotion using classical molecular dynamics (MD) simulations. Witli tire classical force-correlation function, a quantum correction factor Q is needed 5,... 

G. D. Billing, Quanium-Classical Methods in Reaction and Molecular Dynamics, Lecture Notes in Chemistry, A. Lagana and A. Riganelli, eds,. Springer-Verlag, Berlin, 2000,... 

In this section, the basic theory of molecular dynamics is presented. Starting from the BO approximation to the nuclear Schrddinger equation, the picture of nuclear dynamics is that of an evolving wavepacket. As this picture may be unusual to readers used to thinking about nuclei as classical particles, a few prototypical examples are shown. 

Verlet, L. Computer experiments on classical fluids. I. Thermodynamical properties of Lennard-Jones molecules. Phys. Rev. 165 (1967) 98-103. Ryckaert, J.-P., Ciccotti,G., Berendsen, H.J.C. Numerical integration of the cartesian equations of motion of a system with constraints Molecular dynamics of n-alkanes. Comput. Phys. 23 (1977) 327-341. 

Bala, R, Lesyng, B., McCammon, J.A. Extended Hellmann-Feynman theorem for non-stationary states and its application in quantum-classical molecular dynamics simulations. Chem. Phys. Lett. 219 (1994) 259-266. 

Bornemann, F.A., Nettesheim, R, Schiitte, C. Quantum-classical molecular dynamics as an approximation to full quantum dynamics. J. Chem. Rhys. 105 (1996) 1074-1083. 

Since 5 is a function of all the intermediate coordinates, a large scale optimization problem is to be expected. For illustration purposes consider a molecular system of 100 degrees of freedom. To account for 1000 time points we need to optimize 5 as a function of 100,000 independent variables ( ). As a result, the use of a large time step is not only a computational benefit but is also a necessity for the proposed approach. The use of a small time step to obtain a trajectory with accuracy comparable to that of Molecular Dynamics is not practical for systems with more than a few degrees of freedom. Fbr small time steps, ordinary solution of classical trajectories is the method of choice. 

In this paper, we discuss semi-implicit/implicit integration methods for highly oscillatory Hamiltonian systems. Such systems arise, for example, in molecular dynamics [1] and in the finite dimensional truncation of Hamiltonian partial differential equations. Classical discretization methods, such as the Verlet method [19], require step-sizes k smaller than the period e of the fast oscillations. Then these methods find pointwise accurate approximate solutions. But the time-step restriction implies an enormous computational burden. Furthermore, in many cases the high-frequency responses are of little or no interest. Consequently, various researchers have considered the use of scini-implicit/implicit methods, e.g. [6, 11, 9, 16, 18, 12, 13, 8, 17, 3]. 

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

Verlet, L. Computer Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules. Physical Review 159 (1967) 98-103 Janezic, D., Merzel, F. Split Integration Symplectic Method for Molecular Dynamics Integration. J. Chem. Inf. Comput. Sci. 37 (1997) 1048-1054 McLachlan, R. I. On the Numerical Integration of Ordinary Differential Equations by Symplectic Composition Methods. SIAM J. Sci. Comput. 16 (1995) 151-168... 

The principal idea behind the CSP approach is to use input from Classical Molecular Dynamics simulations, carried out for the process of interest as a first preliminary step, in order to simplify a quantum mechanical calculation, implemented in a subsequent, second step. This takes advantage of the fact that classical dynamics offers a reasonable description of many properties of molecular systems, in particular of average quantities. More specifically, the method uses classical MD simulations in order to determine effective... 

Approximation Properties and Limits of the Quantum-Classical Molecular Dynamics Model... 

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. 

This paper is meant as a contribution to systematize the quantum-classical modeling of molecular dynamics. Hence, we are interested in an extended theoretical understanding of the models rather than to further contribute to the bunch of numerical experiments which have been performed on certain models by applying them to particular molecular systems. Thus, we will carefully review the assumptions under which our models are known to approximate the full quantum dynamical (QD) evolution of the system. This knowledge... 

P. Bala, P. Grochowski, B. Lesyng, and J. A. McCammon Quantum-classical molecular dynamics. Models and applications. In Quantum Mechanical Simulation Methods for Studying Biological Systems (M. Fields, ed.). Les Houches, France (1995)... 

Bornemann, F. A., Schiitte, Ch. On the Singular Limit of the Quantum-Classical Molecular Dynamics Model. Preprint SC 97-07 (1997) Konrad-Zuse-Zentrum Berlin. SIAM J. Appl. Math, (submitted)... 

Numerical Integrators for Quantum-Classical Molecular Dynamics... 

---