Enter a species  (e.g. Al, Fe): Hf
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: LJ_ElliottAkerson_2015_Universal__MO_959249795837_003
Enter hexagonal lattice constants 'a' and 'c' (m): 3.403975975536272e-10, 5.55878036026327e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -55.972274
         Iterations: 71
         Function evaluations: 181

Finished calculation

C11 = 2556.366687135962 GPa
C12 = 684.1670558871812 GPa
  B = 1308.2335996367747 GPa