Hi lammps users,
I am not very clear about the deformation method in Aidan's script for calculating elastic constants. Why change the box shape in xy and xz directions when dir=1 and in yz direction when dir=2? Why does the changing distance equal deltaxy, deltaxz, or deltayz, not delta?
Thank you for your help in advance.
Sincere,
Ann