Dear all,
I have to perform uni-axial deformation procedure to my polymer system which has rigid and non rigid atoms. There are 6660 number of atoms, out of which 800 atoms are rigid atoms. The rigid atoms are allowed to move only in a group. Below is the procedure
- Minimization
- NVT at 300 K
- NPT at 300 K
I am using the following script.
include “crystal_n.txt” # contains ids of rigid atoms
group all id 1:6660
group crystal union crys1
group amorphous subtract all crystal # contains ids of non-rigid atoms
neigh_modify exclude molecule crystal