My Community

Hi,
VASP does geometry optimization at zero temperature.
However, experimental measurements may be performed at more higher temperatures.
Is it possible to combine geometry optimization and molecular dynamics
to obtain lattice constants and atomic positions at, for example, T=300K ?
And if yes, how?

1) Run the MD for a suitable period of time.
2) Time-average the coordinates!

(For lattice constants, applying the quasiharmonic approximation---calculating phonon density of states and minimising the free energy---can also be applied. See PRB 73 (2006) 174302. But that approach is much harder.)