This is a simple program calculating the Verlet algorithm for molecular dynamics (MD) taking into account only kinetic energy and the Lennard-Jones (LJ) potential. Therefore it is most suitable for noble gas simulations, like Argon which is in the example. Inert gasses (e.g., N2 ) can also be simulated using this approximation but covalent bonds are not yet implemented.
The LJ potential is defined for same-species atoms/molecules therefore mixed systems are not supported.
Stadard libraries that this repository depend on:
- numpy
- scipy
- Verlet algorithm now produces output file according to tutorial file
- Added parameters file, currently only for masses
- changed angstrom to nm
- fixed get_acceleration
- finished verlet algorithm
- formatting now produces .xyz files readable by VMD
- minor function changes (like get velocity)
- mass is now read from dictionary and not from input file
- Add parameters dictionaries: masses using #n and name, LJ
- Change verlet function to output steps into memeory and then save as file instead of std
- Add Columb interactions
- Check TMPChem youtube tutorials
- Random seed for veloctiy distribution, temperature based
- Random positioning
- Parallelization
- PBC
- Covalent bond option