1 / 16

Molecular Dynamics

Molecular Dynamics. Classical trajectories and exact solutions Finite difference methods: the Verlet algorithm and variants Predictor-corrector methods Choosing an algorithm and time step. Classical Trajectory. Simulated atoms are modeled classically, by Newton’s second Law.

tao
Download Presentation

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. Molecular Dynamics • Classical trajectories and exact solutions • Finite difference methods: the Verlet algorithm and variants • Predictor-corrector methods • Choosing an algorithm and time step

  2. Classical Trajectory Simulated atoms are modeled classically, by Newton’s second Law For particle i with position xi, Forces are determined by the specific program used

  3. Simplest case: Hard Spheres In this, earliest model, interactions occur only in collision. (no Van der Walls forces) Acceleration only occurs at discrete times: when a collision occurs. The trajectories can be solved exactly More sophisticated models have continuous interactions, Each motion is coupled to all other motions; this results in a many-body problem with no analytic solution

  4. The time step dt The numeric solution to the dilemma of the coupling of motions is to use a finite time step, labeled dt. For this timeframe, the changes to the force on the particle due to its own motion are neglected, avoiding the many-body problem. Trajectories of the particles are approximated using a Taylor series:

  5. The Verlet algorithm The next step can be predicted from the Taylor expansion: Similarly, the previous step: Taking the difference cancels the velocity term: Velocity need not be explicitly calculated, but can be estimated, e.g.:

  6. Using Verlet The Verlet algorithm is straightforward and easy to implement However, it requires a previous step’s position, so the first step requires an additional position, which can be estimated simply: (from truncated Taylor expansion.) Numeric problems can arise from the small dt2a(t) being added to the much larger r(t) and r(t-dt) terms

  7. The Half-Step (‘Leapfrog’) Velocities are calculated, halfway between position steps: Explicitly includes the velocity, so energy scaling is easier Avoids a dt2 term, so numeric (rounding) errors are smaller Positions and velocities are not known simultaneously, complicating analysis of the results.

  8. Velocity Verlet This method also uses a half-step velocity calculation: After which a(t+dt) is computed, then the next v:

  9. (a): Verlet (b): Half-Step (Leapfrog) (c ): Velocity Verlet

  10. Beeman’s algorithm Somewhat similar to velocity Verlet: More complicated and requires more memory than velocity Verlet, but is more consistent in energetic consistency.

  11. Gear’s Predictor-Corrector methods Predict ac(t+dt) from the Taylor expansion at the starting point Begin with a simple prediction, as in any of the previous methods Initially step to r(t+dt), v (t+dt), a(t+dt),b(t+dt) at that point. The difference between the a(t+dt) and the predicted ac(t+dt): Estimates the error in the initial step, which is used to correct:

  12. More Predictor-Corrector The set of constants to use is parameterized. For the third-order example above Gear suggested: c0=1/6, c1=5/6, c2=1 (obviously), c3=1/3 These change with different order computations This method can be iterated over the same step to improve the step prediction. This is usually done two or three times. Gear methods are very accurate with short time steps, but are worse than simpler methods for long time steps.

  13. Error with respect to size of timestep Circles: Verlet Squares: Gear 4th order Triangles: Gear 5th order Diamonds: Gear 6th order (RMS Energy deviation) (log/log scale)

  14. Choosing a time step Too short - computation needlessly slow Too long - errors result from approximations Just right - errors acceptable, maximum speed

  15. Overlong Timesteps Particularly near collisions, The forces change quickly. Errors in these regions are compounded in subsequent steps. Simulation of the interatomic distance between two Argon atoms at two dts. The difference from the exact path is plotted.

  16. Choosing an Integration Algorithm Usual test of adequacy of the algorithm is the size of the RMS of the total energy variance over a number of steps. This is strongly dependent on the step size. While methods are better with smaller sizes, some scale up better than others. Often, the best algorithm to choose is the one which maximizes (maximum accurate step size / computation time per step) Other factors may be relevant (memory required, self-starting, compatibility with other algorithms)

More Related