Dear LAMMPS user
I was trying to create a simulation of a solution with suspension in it. So my simulation consists of solvent molecules with nanoparticles suspended in it, with vapour on top of the suspension. I use LJ potential for potential between liquid-liquid and solvent-vapour, while using colloid for solid-liquid and solid-vapour.
Here is my code
units metal
dimension 3
boundary p p p
atom_style sphere
neighbor 2 bin
neigh_modify delay 0 every 1 check yes page 500000 one 50000