[CPMD-list] Energy conservation continued...

[Nick Cooper]

Firstly, thank you very much to those who replied to my initial question
regarding achieving energy conservation at elevated temperatures.

However, having increased the order of the Suzuki integrator and
decomposition ratio using

NOSE PARAMETERS
4 4 4 6 125 4

a constant downwards drift of the "conserved" quantity was observed.

Given this, I suspect that I have set the fictitious kinetic energy for
the electrons (the first number in the NOSE ELECTRONS keyword of CPMD)
too small. Having read "Adiabaticity in first-principles molecular dynamics"
by Blochl and Parrinello: PRB 45, 9413 (1992), I came upon their expression
for E_kin,ad (eqn 6 in the paper), which provides an estimate for the
transition point between two regimes. However, I have a couple of questions
regarding this expession:
1) I assume mu (EMASS in CPMD) is in electron mass units and M proton
mass units?
2) How do I determine the value of the of what's contained in the sigma_i?

Sorry for asking such a basic question but I haven't got anyone to ask
at present, so I was hoping someone here could enlighten me.

Thank you for your time.

Regards,

Nick Cooper