An embedded atom method potential of beryllium

We present an embedded atom method (EAM) potential for hexagonal beryllium, with a pair function in the form of a Morse potential and a Johnson embedding function with exponential electron density. The cohesive energy, elastic constants, lattice parameters and relaxed vacancy formation energy were used to fit the potential. The fitted-potential was validated by a comparison to first-principles and, wherever available, experimental results for the lattice energies of various crystal structures: vacancy cluster, interstitial formation and surface. Using a large cutoff distance of 5 angstrom, which includes interactions to approximately the eighth neighbor shell of beryllium, allows our potential to reproduce these quantities considerably better than previous EAM potentials. The accuracy obtained by our potential is similar to or in some cases even better than available modified EAM potentials, while being computationally less intensive.