Lennard-Jones interactions molecular dynamics simulation

Ligand-Protein Interactions The energy of formation of ligand-protein contacts can be computed as a sum of non-bonded (Lennard-Jones and electrostatic) terms similar to those used in a molecular dynamics simulation. 

We now present results from molecular dynamics simulations in which all the chain monomers are coupled to a heat bath. The chains interact via the repiflsive portion of a shifted Lennard-Jones potential with a Lennard-Jones diameter a, which corresponds to a good solvent situation. For the bond potential between adjacent polymer segments we take a FENE (nonhnear bond) potential which gives an average nearest-neighbor monomer-monomer separation of typically a 0.97cr. In the simulation box with a volume LxL kLz there are 50 (if not stated otherwise) chains each of which consists of N -i-1... 

Extension to many dimensions provides insight into more sophisticated aspects of the method and into the nature of molecular interactions. In the second stage of this unit, the students perform molecular dynamics simulations of 3-D van der Waals clusters of 125 atoms (or molecules). The interactions between atoms are modeled using the Lennard-Jones potentials with tabulated parameters. Only pairwise interactions are included in the force field. This potential is physically realistic and permits straightforward programming in the Mathcad environment. The entire program is approximately 50 lines of code, with about half simply setting the initial parameters. Thus the method of calculation is transparent to the student. 

Fig. 2. Average mer distance from the host surface versus Npx for a brush in (a) a good solvent with x=l/3 (T=4.0e/kB) and (b) a 0 solvent with x=l/2 (T= 3.0 e/kB) from a molecular dynamics simulation in which the Lennard-Jones interaction between mers is truncated at rc=2.5cr. The results for the 0 solvent are shifted by 0.2 for clarity. The results are for chain length N= 25 (A), 50 ( ), 100 ( ), and 200 (O). From ref. [47].

Figure 1.38. Molecular dynamics simulation of the density profiles for spherical molecules in a cylinder, mimicking SFg in controlled pore glass (CPG-10). Fluid-fluid and fluid-wall interaction modelled by Lennard-Jones interactions. Reference A. de Keizer. T. Michalski and G.H. Findenegg, Pure Appl. Chem. 63(1991) 1495.

Figure 2.5 shows a Molecular Dynamics simulation, which is the counterpart of the previous example. Particle trajectories are shown for fluid particles, interacting via an law (as in the Lennard-Jones interaction), with a structureless "hard" wall. A striking feature is that the diffusion between the... 

In such molecular dynamic simulations, one starts with an array of atoms or molecules , initially on a lattice, interacting with one another via an interatomic potential. These interacting potentials were taken by Paskin et al (1980, 1981) to be the Lennard-Jones potential (l> rij) = e[(l/rij) — 2(l/rij) ], where e denotes the depth of the potential energy and rij denotes the interatomic separation of the atoms. This potential is assumed to have a sharp cut-off at an arbitrarily chosen value 1.6 (lattice constant) of the interatomic separation. The external stress or force is applied only at the boundary surface atoms. In order to investigate the Griffith fracture phenomena, one can consider for example a two-dimensional lattice of linear size L, remove a few I L) consecutive bonds along a horizontal row in the middle of the network, and apply tensile force on the upper and lower surface atoms in the vertical direction. 

Ionic monolayers can be, and have also been, analyzed theoretically either with advanced lattice theories or with Monte Carlo or molecular dynamics simulation. Basic principles and some illustrations of monolayer compositions have already been discussed in sec. 3.5. The step from Langmuir to Gibbs monolayers is theoretically realized through the choice of the adsoption energy. As before, the selection of the various parameters (x -interaction parameters in lattice theories, constants in the Lennard-Jones, interactions in MD, etc.) and approximations (choice of lattice, accounting for stereoisomery, or extent of truncation, respectively) remain a central issue. In view of the growing power of computers, increasingly better results may be expected in the near future. 

R0sjorde et al studied the phase transition in a pure fluid using non-equilibrium molecular dynamics simulations (NEMD). The NEMD method solves Newton s equations of motion for several thousand particles in an imaginary box see Hafskjold for a review. The particles interacted with a Lennard-Jones-type pair... 

Figure 1 An example of the pair distribution function g(r), the potential of mean force w(r), and the true interatomic potential u(r) based on our molecular dynamic simulations with a simple monatomic system. The only interatomic interactions in the system are the Lennard-Jones format potential.

In classical molecular dynamics simulations, atoms are generally considered to be points which interact with other atoms by some predehned potential form. The forms of the potential can be, for example, Lennard-Jones potentials or Coulomb potentials. The atoms are given velocities in random directions with magnitudes selected from a Maxwell-Boltzman distribution, and then they are allowed to propagate via Newton s equations of motion according to a finite-difference approximation. See the following references for much more detailed discussions Allen and Tildesley (1987) and Frenkel... 

These early MD studies were limited to hard sphere or Lennard-Jones interactions appropriate for simple fluids. In the early 1970s the first attempts to apply MD methods to molecular fluids were made. Harp and Berne [11] studied a model for carbon monoxide while Rahman and Stiflinger [12] carried out MD simulations of hquid water. The MD simulations of water were especially important because they addressed questions related to the dynamics and structure of an ubiquitous solvent that would figure prominently in later work on biological and condensed phase systems. Furthermore, because of the nature of the long-range Coulomb forces in this system and the dehcate nature... 

---