Calculation update! New properties have been added to the website for dislocation monopole core structures, dynamic relaxes of both crystal and liquid phases, and melting temperatures! Currently, the results for these properties predominately focus on EAM-style potentials, but the results will be updated for other potentials as the associated calculations finish. Feel free to give us feedback on the new properties so we can improve their representations as needed.
Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: J. Wormald, J. Kempfert, J. March-Rico, and R. Smith (2025), "NNL.Fe.qSNAP-ZBL.2024.1: A Fe Spectral Neighbor Analysis Potential for Radiation Damage Simulations". DOI: 10.2172/3006967.
Abstract: The NNL.Fe.qSNAP-ZBL.2024.1 machine-learned potential (MLP) has been generated to support the development of an elemental body-centered cubic (BCC) Fe athermal recombination corrected neutron damage model and simulations of primary recoil atom (PRA) cascades in BCC Fe. This MLP is a quadratic spectral neighbor analysis potential (qSNAP) hybridized with the universal Ziegler-Beirsack-Littmark (ZBL) potential at short-range and is named according to Naval Nuclear Laboratory MLP naming conventions (NNL.material-system.MLP-type.year.version). Training set calculations for Fe are presented along with the subsequent MLP fitting procedure. A key criterion of the fitting procedure is that ZBL describes the short-range interaction with minimal impact on the MLP. The MLP is compared to density functional theory (DFT) predicted properties relevant to radiation damage simulation, including threshold displacement energies, for validation. The NNL.Fe.qSNAP-ZBL.2024.1 potential is considered suitable for molecular dynamics (MD) simulations of radiation defects up to 800 K and PRA cascades in BCC Fe up to around 10 keV. The potential can additionally be used on a limited basis for recoils of 10-20 keV, within which range the emergence of structures outside the training set in cascade simulations may cause system instabilities.
See Computed Properties Notes: These files were provided by Jonathan Wormald (Naval Nuclear Laboratory) on August 19, 2026. The .txt file includes example LAMMPS commands for defining the potential with the ZBL parameters. File(s):