[lammps-users] trouble using rigid bodies in simulation (lost atoms)
surendra jain <jainsk.iitkgp@...24...>
Fri, 28 Jul 2017 17:32:25 +0530
Dear Lammps users,
I am studying the diffusion of CO2 molecules in a carbon model. In
my simulation I fix the carbon atoms (setforce 0.0 ..) and make the
CO2 molecules as a rigid body molecule (having 3 atoms). However I get
the following error : "lost atoms:". After debugging I found that one
of the atoms of CO2 molecule is very close to a carbon atom (rsq = 0.5
Angsgtrom). This causes huge force and hence huge velocity on the
atom. I have the following question :
(1) Does the rigid body moves as a whole or the coordinates of the
atoms are updated individually.
(2) If one atom faces huge force then the whole rigid body also
experiene the huge force.
(3) Since the molecule is rigid, is it possible for one of the atoms
to be detached from the molecule.
I ran a simulation of timestep 0.5fs and got the lost atoms error.
Should I use a smaller timestep. Can you please suggest me a suitable
timestep for rigid body simulations.