HarmOsc_MD1_LJ

# HarmOsc_MD1_LJ - Introduction to Molecular Dynamics The...

This preview shows pages 1–7. Sign up to view the full content.

1 Introduction to Molecular Dynamics: The Harmonic Oscillator The LeapfrogAlgorithm The Lennard-Jones Potential BIOL 7110 / CHEM 8901 / BIOL 4105 / CHEM 4804 February 15, 2011

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
2 To be covered today Overview of molecular dynamics Newton’s Equations Differential equations converted to difference equations dt vs Δ t (difference equation) The MD algorithm: numerical integration of Newton’s equations The harmonic oscillator: theory The harmonic oscillator: treatment by MD (Pseudocode) The harmonic oscillator: results of the simulation The Lennard-Jones potential
3 Newton’s Equations: Energy, Force and Acceleration MD represents the numerical integration of Newton’s equations F i = [ ˆ k j ( E / r ij )] j = 1 3 We note that the force is just the negative gradient of the energy The trajectory of a cannonball. y x x ( t ) = v x t y ( t ) = v y t 1 2 gt 2 x (0) = 0 y (0) = 0 Analytic solution

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
4 Numerical Integration of Newton’s Equations (one-dimensional case) stepwise, replacing the differential equations By the difference equations Instead of solving Newton’s equations analytically, the computer uses numerical approximations x ( dt ) = x (0) + dx dx = vdt x ( Δ t ) = x (0) + Δ x Δ x = v Δ t We solve the integral equation x ( t ) = x (0) + v ( u ) du u = 0 u = t Numerical Integration of Newton’s Equations (two-dimensional case) stepwise, replacing the differential equations By the difference equations We solve the integral equations x ( t ) = x (0) + v x ( u ) du u = 0 u = t y ( t ) = y (0) + v y ( u ) du u = 0 u = t x ( dt ) = x (0) + dx y ( dt ) = y (0) + dy x ( Δ t ) = x (0) + Δ x y ( Δ t ) = y (0) + Δ y dx = v x dt dy = v y dt Δ x = v x Δ t Δ y = v y Δ t
5 The trajectory of a cannonball. y x Discretization of the trajectory The trajectory of a cannonball. y x The discretized trajectory closely matches the analytic solution

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
6 Molecular Dynamics Input to a molecular dynamics program: A set of molecular x,y,z coordinates E = f(x 1 ,y 1 ,z 1 ,…,x N ,y N ,z N ) An energy function x’ = x + v Δ t + a( Δ t) 2 /2 + … Algorithm: numerical integration of Newton’s equations of motion Just turn the crank Numerical Integration:
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}

### What students are saying

• As a current student on this bumpy collegiate pathway, I stumbled upon Course Hero, where I can find study resources for nearly all my courses, get online help from tutors 24/7, and even share my old projects, papers, and lecture notes with other students.

Kiran Temple University Fox School of Business ‘17, Course Hero Intern

• I cannot even describe how much Course Hero helped me this summer. It’s truly become something I can always rely on and help me. In the end, I was not only able to survive summer classes, but I was able to thrive thanks to Course Hero.

Dana University of Pennsylvania ‘17, Course Hero Intern

• The ability to access any university’s resources through Course Hero proved invaluable in my case. I was behind on Tulane coursework and actually used UCLA’s materials to help me move forward and get everything together on time.

Jill Tulane University ‘16, Course Hero Intern