I want to study the thermal properties of a mixed system which is constructed by virtual crystal approximation in VASP. When I try to run the ab initio Molecular Dynamics of this system in VASP, I found that there would be either nearly infinite Temperature (kinetic energy as well) or large ERROR output. So I am wondering if the system with VCA could be used to run the AIMD via vasp or any other packages?
Look forward your advice and thanks a lot!