Hi Naveen and Paul,
Thanks for your help. I tried to run the code using harmonic dihedral. Now, it starts the run but after around 50-60 time steps it says “ERROR on proc 0: Bad FENE bond”.
My bond_style is fene and my bond_coeffs 1 30 1.5 1.0 1.0.
I am not able to understand the problem. Are the bonds breaking? I bilieve I defined the bonds, angles and the dihedrals correctly (w/o leaving out any).
Thanks for your help. I really appreciate that as I am new to using LAMMPS.
Also, Can you explain how view the trajectories in a more detailed manner.
Thanks
Venkat