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: E.L. Sikorski, M.A. Cusentino, M.J. McCarthy, J. Tranchida, M.A. Wood, and A.P. Thompson (2023), "Machine learned interatomic potential for dispersion strengthened plasma facing components", The Journal of Chemical Physics, 158(11), 114101. DOI: 10.1063/5.0135269.
Abstract: Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature, and at fusion reactor temperatures (≥1000 K), it may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperatures are still unknown. We present a machine learned Spectral Neighbor Analysis Potential for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that although the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500 K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500 K.
Notes: This potential was optimized for bulk and surface structures for both W and ZrC. No optimization was performed on pure C structures, and no physical performance should be expected for pure C simulations. We expect the potential to perform well in the temperature range of 300 - 2500K. Primary optimization was performed on bulk modulus, (100) and (110) surface energies, thermal expansion, and several stability checks detailed in the publication.
See Computed Properties Notes: These files were provided by Ember Sikorski on 28 March 2023. This potential can be used by adding "include in.pot_snapWZrC" to a LAMMPS input script. The zbl parameters in in.pot_snapWZrC must be included to achieve the accuracy and performance described in the publication. File(s):