LAMMPS is a really popular code for doing classical molecular dynamics simulations. In this post, we will dive deep into the heart of molecular dynamics and also how to construct the input file for LAMMPS.
It is really important to know that when you try to learn any softwares, it is really important to understand the underlying physics. That means you should be able to match the input file to a particular physical system that you want to simulate, and you need to be able to tell all the assumptions in that physical model in order to determine its accuracy compared to the real system. You need to know pros and cons of all the assumptions and how they influence the accuracy of your simulation, that's the key in learning any software. Without understanding, all you get is just numbers, without meaning. It is this understanding that is crucial in computational chemistry, otherwise you can be replaced by the AI and automated systems.
Fundamentals of Molecular Dynamics
Molecular dynamics is really easy to explain: just treating all the nuclei as point particles, and using the Newtonian equation of motion to solve the time evolution of the system, In mathematical term, that means:
$$m_i\frac{d^2\vec{r}}{dt^2}=-\vec{f}_i$$
where: \(\vec{f}_i=-\frac{\partial U}{\partial \vec{r}_i}\) and \(U\) is the many-body potential of the whole system. This potential energy will be approximated by using classical force field (in classical MD) or ab-initio simulation (in Ab-Initio MD).
So in order to create the input file, we need to know:
- Structure of the system
- Types of atoms in the system
- The box of the system (periodic boundary condition)
- The coordinates of all the atoms (whether it is cartesian or ractional)
- The initial velocity of all the atoms
- Choice of the force field
- In Classical MD, there are many different possibilies
- But what I'm interested in is how to use DeePMD as the potential to conduct MD simulations
- Settings for the time-evolution
- time step
- how many iterations we want to conduct (run)
- Environment
- thermostate
Based on these knowledge, we can begin constructing our system.
Building the input files for LAMMPS
In general, there are two input files which we need for LAMMPS: (1) the input file for specifying all the parameters that we need (usually called in.lammps
(2) the structural file, which contains all the structural information about the system (usually called some_name.data
.
Here is an example file of in.lammps
:
boundary p p p
units metal
atom_style atomic
neighbor 2.0 bin
neigh_modify every 10 delay 0 check no
read_data lmp_pbWater.data
mass 1 1
mass 2 16
mass 3 207
pair_style deepmd graph.pb
pair_coeff
velocity all create 300.0 23456778
fix 1 all nvt temp 300.0 0.5
timestep 1
thermo 100
#thermo_modify lost warn
dump myDump2 all custom 1000 geometry.xyz type x y z
run 100000
write_data lmp_pbWater.dat
write_restart lmp_pbWater.rest
And in below is the example file of pbWater.data
:
lmp_pbWater.data (written by ASE)
162 atoms
3 atom types
0.0 10.649027999999999 xlo xhi
0.0 12.296438999999999 ylo yhi
0.0 20 zlo zhi
Atoms
1 3 10.611976730389317 12.281261505342298 0.14421639999999999
2 3 3.40217883344196 0.18510936306209999 0.25271199999999999
3 3 7.0938796064978398 0.53783837213903996 19.963010999999998
4 3 1.6040585085579595 3.1612538754066595 0.025070400000000003
5 3 5.2367321453709597 3.3484584287099701 0.3163494
6 3 8.7861238873168794 3.3167923930683902 0.23802740000000003
7 3 0.035780521099439994 6.4401701622453 19.889976599999997
8 3 3.23187414666168 6.3574113158792098 19.974749800000001
9 3 6.9956113345263597 6.3705788346179695 19.955363999999999
10 3 1.8640909468537195 9.6805229719140904 19.898860800000001
In the above I only show the coordinates of first 10 atoms, the lammps structure file can be generated by using ASE (Atomic Simulation Environment).