Enter a species  (e.g. Al, Fe): Tb
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.773548335011991e-10, 6.162301447658867e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -23.315636
         Iterations: 74
         Function evaluations: 193

Finished calculation

C11 = 781.6411422801003 GPa
C12 = 209.19266490324034 GPa
  B = 400.0088240288603 GPa