I am performing a uniaxial tensile test on a calcium silicate hydrate structure of 150 atoms. I am using the buckingham and stillenger and weber potentials. The structure is minimized before strainrates are appplied. I am applying a constant strainrate in the x direction via the erate command and applying zero applied pressure to the other dimensions via fix npt anisto command. The stress is computed using the following command. Are these correct commands to compute the stress in the my system?
compute peratom all stress/atom
compute p all reduce sum c_peratom c_peratom c_peratom
variable press equal -(c_p+c_p+c_p)/(3*vol)
thermo_style custom step temp etotal press v_press
I am currently using strainrates of 0.00002/ps and 0.002/ps. At these rates I get really low negative stresses (-3,000MPa) for tensile strains and high positive values for compressive strains. I even tried a strainrate of 2.96e-15/ps but the run was aborted.
Also, I wanted to know is there a minimum strainrate used for Md simulations / Lammps?