1 / 24

Basics of molecular dynamics

Basics of molecular dynamics. Equations of motion for MD simulations. The classical MD simulations boil down to numerically integrating Newton’s equations of motion for the particles. Lennard-Jones potential.

betrys
Télécharger la présentation

Basics of molecular dynamics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Basics of molecular dynamics

  2. Equations of motion for MD simulations • The classical MD simulations boil down to numerically integrating Newton’s equations of motion for the particles

  3. Lennard-Jones potential • One of the most famous pair potentials for van der Waals systems is the Lennard-Jones potential

  4. Lennard-Jones Potential

  5. Dimensionless Units Advantages of using dimensionless units: • the possibility to work with numerical values of the order of unity, instead of the typically very small values associated with the atomic scale • the simplification of the equations of motion, due to the absorption of the parameters defining the model into the units • the possibility of scaling the results for a whole class of systems described by the same model.

  6. Dimensionless units • When using Lennard-Jones potentials in simulations, the most appropriate system of units adopts σ, m and ε as units of length, mass and energy,respectively, and implies making the replacements:

  7. Integration of the Newtonian Equation Classes of MD integrators: • low-order methods – leapfrog, Verlet, velocity Verlet – easy implementation, stability • predictor-corrector methods – high accuracy for large time-steps

  8. Boundary Condition

  9. Initial state • Atoms are placed in a BCC, FCC or Diamond lattice structure • Velocities are randomly assigned to the atoms. To achieve faster equilibration atoms can be assigned velocities with the expected equilibrium velocity, i.e. Maxwell distribution. • Temperature adjustment: Bringing the system to required average temperature requires velocity rescaling. Gradual energy drift depends on different factors- integration method, potential function, value of time step and ambient temperature.

  10. Conservation Laws • Momentum and energy are to be conserved throughout the simulation period • Momentum conservation is intrinsic to the algorithm and boundary condition • Energy conservation is sensitive to the choice of integration method and size of the time step • Angular momentum conservation is not taken into account

  11. Equilibration • For small systems whose property fluctuate considerably, characterizing equilibrium becomes difficult • Averaging over a series of timesteps reduces the fluctuation, but different quantities relax to their equilibrium averages at different rates • A simple measure of equilibration is the rate at which the velocity distribution converges to the expected Maxwell distribution

  12. Velocity distribution as a function of time

  13. Interaction computations • All pair method • Cell subdivision • Neighbor lists

  14. Thermodynamic properties at equilibrium

  15. Thermophysical and Dynamic Properties at equilibrium

  16. Example: Calculation of thermal conductivity at eqilibrium Following are the equations required to calculate thermal conductivity: Where Q is the heat flux Where k is the thermal conductivity

  17. Nonequilibrium dynamics • Homogeneous system: no presence of physical wall, all atoms perceive a similar environment • Nonhomogeneous system: presence of wall, perturbations to the structure and dynamics inevitable • Nonequlibrium more close to the real life experiments where to measure dynamic properties systems are in non equilibrium states like temperature, pressure or concentration gradient

  18. Example: Calculation of thermal equilibrium at nonequilibrium ( direct measurement) • To measure thermal conductivity of Silicon rod (1-D) heat energy is added at L/4 and heat energy is taken away at 3L/4 • After a steady state of heat current was reached, the heat current is given by • Using the Fourier’s law we can calculate the thermal conductivity as follows • Stillinger-Weber potential for Si has been used which takes care of two body and three body potential

  19. Example: continued..

  20. Molecular Dynamics Simulation of Thermal Transport at Nanometer Size Point Contact on a Planar Si Substrate Schematic diagram of the simulation box. Initial temperature is 300 K. Energy is added at the center of the top wall and removed from the bottom and side walls. Calculated temperature profile in the Si substrate for a 0.5 nm diameter contact radius.

  21. Thermal conductivity of nanofluids at Equilibrium Schematic diagram

  22. Results

  23. Results continued

  24. References • Rapaport D. C., “The Art of Molecular Dynamics Simulation”, 2nd Edition, Cambridge University Press, 2004 • W. J. Minkowycz and E. M. Sparrow (Eds), “Advances in Numerical Heat Transfer”,vol. 2, Chap. 6, pp. 189-226, Taylor & Francis, New York, 2000. • Koplik, J., Banavar, J. R. &Willemsen, J. F., “Molecular dynamics of Poiseuille flow and moving contact lines”, Phys. Rev. Lett. 60, 1282–1285 (1988); “Molecular dynamics of fluid flow at solid surfaces”, Phys.Fluids A 1, 781–794 (1989).

More Related