Enter a species  (e.g. Al, Fe): Kr
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: LJ_Shifted_Bernardes_1958LowCutoff_Kr__MO_995724792024_003
Enter hexagonal lattice constants 'a' and 'c' (m): 3.997291501543885e-10, 6.531054338208935e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -0.215780
         Iterations: 72
         Function evaluations: 186

Finished calculation

C11 = 6.317103006725966 GPa
C12 = 1.6873112903551923 GPa
  B = 3.2305751958121167 GPa