Enter a species (e.g. Al, Fe): Pr 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.948607760006343e-10, 6.448178070383309e-10 Calculating isothermal bulk modulus Optimization terminated successfully. Current function value: -21.641405 Iterations: 73 Function evaluations: 188 Finished calculation C11 = 633.2324234696048 GPa C12 = 169.47364050368154 GPa B = 324.0599014923226 GPa