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: S. Malgope, M.K. Gupta, S. Bag, R. Mittal, S. Bhattacharya, A. Singh, and S.L. Chaplot (2024), "Untangling high-temperature thermal expansion and lattice thermal conductivity behavior of vanadium using machine-learned molecular dynamics", Physical Review B110(5), 054301. DOI: 10.1103/physrevb.110.054301.
Abstract: Vanadium metal is extensively used in modern technology, especially in the alloy and steel industry; it exhibits anomalous thermal expansion behavior across the entire temperature regime. Here, we extensively investigate the phonon anharmonicity contributed by volume change (implicit anharmonicity) and thermal amplitude (explicit anharmonicity) and their impact on thermal expansion and thermal transport in vanadium up to 2000 K, close to the melting temperature. We compared the different methods to evaluate the phonon anharmonicity, namely, the quasiharmonic approximation (QHA), temperature-dependent effective potential (TDEP) method, and machine-learned force-field-based molecular dynamics (MLMD) simulations. At 300 K, QHA overestimates the thermal expansion coefficient by ~20%, while TDEP provides an excellent description of the experimental data. This reveals a significant explicit anharmonicity at room temperature. At higher temperatures, the experimental thermal expansion coefficient continues to rise up to twice the QHA estimates, indicating significant anharmonicity. However, TDEP underestimates the experimental observations, as it only includes low orders of anharmonicity. MLMD, which includes all the anharmonic effects, successfully explains the anomalous expansion behavior over 500-2000 K. It is expected that the electronic entropy and the electron-phonon interaction would influence the thermal expansion, but their effect appears to be small. We used MLMD to calculate the spectral energy density of phonons up to 2000 K, which revealed small phonon shifts but large broadening. Above ~2000K, MLMD captures the melting and reproduces the experimental volume increase on melting. We have also calculated the lattice thermal conductivity using the TDEP-based third-order-perturbation method and the MLMD-based Green-Kubo method, over 300-1500 K, which includes a higher order of anharmonicity. This brings out the important contribution from the four and higher orders of phonon anharmonicity to the lattice thermal conductivity.
Notes: The potential was trained using first-principles ab initio molecular dynamics (AIMD) data generated at different temperatures. The training and testing details for this MLFF are discussed in the above citation.
See Computed Properties Notes: These files were provided by Samiran Malgope on Nov 3, 2025. The parameter file is in Deep Potential (DP) format and was generated using DeePMD-kit v2.0.1. It can be used with LAMMPS via the DeePMD plugin, which was tested by the authors using LAMMPS version 29 Oct 2020. File(s):