Hello LAMMPS users,
I have a question about some behavior I’m seeing with fix deform.
I’m trying to run a uniaxial tension simulation at constant strain rate for a bead-spring polymer model. The ensemble is periodic and is first equilibrated via fix npt at 0 pressure and 300K in a separate run. The tensile-test run starts by reading the resulting restart file.
I’m using fix deform to stretch the ensemble in the z direction. I’m also using fix npt for zero pressure conditions in the x and y directions. The commands I’m using are
fix STRETCH all deform 1 z vel ${vel} units box
fix NPT all npt temp 300 300 100 x 0 0 1000 y 0 0 1000
where ${vel} is a specified stretch velocity.
What I’m seeing is that if I run the simulation for two different lengths of time (e.g., 200,000 steps vs 400,000 steps), the simulation evolves slightly differently in each run. What might be causing this result? If I’m stretching the unit cell at the same velocity in each case, I would expect the simulations to exactly match at all the time points they have in common.
I’m running version April 12 2013. The issue occurs both on a single processor and in parallel.
Thanks,
Brian