Introduction

In this chapter we focus on a number of techniques used in MD simulation, primarily the methods for computing the interactions and integrating the equations of motion. The goal is to generate the atomic trajectories; subsequent chapters will deal with the all-important question of analyzing this raw ‘experimental’ data. We continue to work with the simplest atomic systems, in other words, with monatomic fluids based on the LJ potential, not only because we want to introduce the methodology gradually, but also because a lot of the actual qualitative (and even quantitative) behavior of many-body systems is already present in this simplest of models. Models of this kind are widely used in MD studies of basic many-body behavior, examples of which will be encountered in later chapters.

Equations of motion

While Newton's second law suffices for the dynamics of the simple atomic fluid discussed in this chapter, later chapters will require more complex forms of the equations of motion. The Lagrangian formulation of classical mechanics provides a general basis for dealing with these more advanced problems, and we begin with a brief summary of the relevant results. There are, of course, other ways of approaching the subject, and we will also make passing reference to Hamilton's equations. A full treatment of the subject can be found in textbooks on classical mechanics, for example [go180].