Multi-principal-component alloys have attracted great interest as a novel paradigm in alloy design, with often unique properties and a vast compositional space auspicious for materials discovery. High entropy alloys (HEAs) belong to this class and are being investigated for prospective nuclear applications with reported superior mechanical properties including high-temperature strength and stability compared to conventional alloys. Computational materials design has the potential to play a key role in screening such alloys, yet for high-temperature properties, challenges remain in finding an appropriate balance between accuracy and computational cost. Here we develop an approach based on density-functional theory (DFT) and thermodynamic integration aided by machine learning based interatomic potential models to address this challenge. We systematically evaluate and compare the efficiency of computing the full free energy surface and thermodynamic properties up to the melting point at different stages of the thermodynamic integration scheme. Our new approach provides a ×4 speed-up with respect to comparable free energy approaches at the level of DFT, with errors on high-temperature free energy predictions less than 1 meV/atom. Calculations are performed on an equiatomic HEA, TaVCrW - a low-activation composition and therefore of potential interest for next generation fission and fusion reactors.

The unique and unanticipated properties of multiple principal component alloys have reinvigorated the field of alloy design and drawn strong interest across scientific disciplines. The vast compositional parameter space makes these alloys a unique area of exploration by means of computational design. However, as of now a method to compute efficiently, yet with high accuracy the thermodynamic properties of such alloys has been missing. One of the underlying reasons is the lack of accurate and efficient approaches to compute vibrational free energies—including anharmonicity—for these chemically complex multicomponent alloys. In this work, a density-functional-theory based approach to overcome this issue is developed based on a combination of thermodynamic integration and a machine-learning potential. We demonstrate the performance of the approach by computing the anharmonic free energy of the prototypical five-component VNbMoTaW refractory high entropy alloy.

One of the effective potentials that has proven to be very versatile and useful for describing metals is the modified embedded atom method (MEAM) potential. The reference-free version of the MEAM (RF-MEAM) potential provides more flexibility for fitting than the 2NN-MEAM because it also describes the pair potential as an explicit function. In this work, we present a methodology to fit RF-MEAM potentials to DFT data. We then evaluate the performance of the fitted potential by comparing MD simulations with experimental and DFT data. As an example, the methodology is applied to a binary and a quaternary alloy, namely NiTi and NbMoTaW. In the case of the equi-atomic NiTi shape memory alloy, our attention focuses on designing a potential that properly captures its mechanical behavior, given that the existing potentials fail to predict elastic constants in agreement with experiments. To reach our aim, we included the stress tensors of different high temperature NiTi configurations in the fitting database. The obtained RF-MEAM potential outperforms existing EAM and MEAM potentials in predicting the lattice and elastic constants of austenitic and martensitic phases as well as the corresponding transformation temperatures. To demonstrate the suitability of this methodology also for more complex systems, a RF-MEAM potential is fitted to model the multi-component NbMoTaW high-entropy alloy. Validation is achieved through comparison between observables obtained through the MD output and ab initio data. The article also reports key improvements to the optimization code MEAMfit v2 and the freely-available LAMMPS implementation of the RF-MEAM formalism. Most notably, resorting to analytic derivatives of the objective function with respect to the potential parameters rather than derivatives through finite differences, the time necessary for fitting has decreased by an order of magnitude.