Molecular Dynamics - Whats going wrong?
I have been writing a MATLAB code to reproduce the results of Verlets paper - Computer "experiments" on classical fluids, 1967.
Here is how I go about it:
- I initialize all the velocities (uniform distribution about [-1, 1])
- Initialize all positions (NOT random allocations).
- I use Verlet's algorithm to update the positions.