Enter a species (e.g. Al, Fe): Xe Enter a hexagonal lattice type (e.g. hcp): hcp Enter a model name: LJ_Shifted_Bernardes_1958MedCutoff_Xe__MO_849320763277_003 Enter hexagonal lattice constants 'a' and 'c' (m): 4.344827670449861e-10, 7.095215419620568e-10 Calculating isothermal bulk modulus Optimization terminated successfully. Current function value: -0.333040 Iterations: 72 Function evaluations: 187 Finished calculation C11 = 7.315509203545162 GPa C12 = 1.9578687162120136 GPa B = 3.743748878656396 GPa