Molecular Dynamics Algorithms: Difference between revisions
Jump to navigation
Jump to search
no edit summary
No edit summary |
No edit summary |
||
Line 36: | Line 36: | ||
ener%pot%tot = ener%pot%tot + pmfene | ener%pot%tot = ener%pot%tot + pmfene | ||
</pre> | </pre> | ||
<li>'''Update Velocities</li> | <li>'''Update Velocities''' + thermostat</li> | ||
<pre> | <pre> | ||
! backup old velocities | ! backup old velocities | ||
Line 43: | Line 43: | ||
v(t+dt/2) = v(t-dt/2) + a(t)*dt | v(t+dt/2) = v(t-dt/2) + a(t)*dt | ||
</pre> | </pre> | ||
<li>'''Update Positions</li> | <li>'''Update Positions''' + barostat</li> | ||
<pre> | <pre> | ||
! backup positions | ! backup positions |