Enter a species  (e.g. Al, Fe): Tb
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: EAM_Dynamo_Mendelev_2018_Tb__MO_522239651961_000
Enter hexagonal lattice constants 'a' and 'c' (m): 3.62496330167414e-10, 5.657954530848218e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -8.381790
         Iterations: 71
         Function evaluations: 185

Finished calculation

C11 = 76.67653862273514 GPa
C12 = 18.857516939978215 GPa
  B = 38.13052416756385 GPa