Question in Example/ELASTIC_T


Line 5 in states the value of elastic constant should not depend on the variable “up”.

I tried using different values for “up” and got the following values

  1. up = 2.0e-4 --> Elastic Constant C11all = 133.612 GPa
  2. up = 2.0e-4 --> Elastic Constant C11all = 359.269 GPa
  3. up = 2.0e-6 --> Elastic Constant C11all = 22932.335 GPa

It is stated in file in.elastic

It is also important to verify that the MD sampling of stress components is generating accurate statistical averages.

One indication of this is that the elastic constants are insensitive to the choice of the variable ${up} in init.mod.

Just wanted to verify is the script is valid even for the given example of Si.



This is a Q for Aidan (CCd).