Enter a species (e.g. Al, Fe): Gd 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.812444350333819e-10, 6.225819716631428e-10 Calculating isothermal bulk modulus Warning: Maximum number of function evaluations has been exceeded. Finished calculation C11 = 1136.7108885868658 GPa C12 = 304.22089737557747 GPa B = 581.7175611126736 GPa