[Pw_forum] CP simulation with ensemble-dft algorithm

Giovanni La Penna lapenna.giovanni at gmail.com
Sun Jun 19 11:51:53 CEST 2011


Dear Pascal and QE-users,

 In the *.evp (and in stdout) of CP, ekinc (Kel) and actual T are printed,
as, respectively, columns 2 and 4 (in *.evp). Kel is the fictitious
kinetic energy
of electronic variables.
Empirically, I found that:

Kel/Kion < 1/10

for a reliable CP-MD run, with Kion= 3/2 Ndeg R T . Since Kel is
(please confirm) in Ha,
R must be converted in Ha/K.
Even though the word reliable may be an opinion, when the ratio
becomes larger, then Kel starts diverging very quickly (while in the
other case only a slow drift occurs).
I would be glad if other CP users could confirm this empirical observation.

                      Giovanni La Penna (Cnr-Iccom, Italy)



More information about the users mailing list