[Pw_forum] question on the pressure in molecular dynamics
Davide Ceresoli
ceresoli at sissa.it
Wed Jan 17 18:36:33 CET 2007
Dear Eduardo,
you are totally right regarding the pressure. At finite temperature,
even at constant volume, you must add the kinetic contribution to the
pressure. This contribution always tends to "expand" the solid and for
high enough temperature, it may become important.
I've looked at the Andersen and Parrinello-Rahman papers, and the
kinetic term is always present in the formulas. The "external" pressure
you apply on the system must to balance both the kinetic term and the
virial term.
So, the kinetic term is correctly taken into account in the P-R
equations of motion, in the code. Could be that this contribution
(which averages to P_kin = (N/V) k_B T) is not printed in output.
Regards,
Davide
Eduardo Ariel Menendez P wrote:
> Hi,
> Let me insist in this question, as I have had no answer yet.
>
> >
> > Hi everybody,
> > I have a doubt concernig the pressure reported by pw or cp in a molecular
> > dynamics run. The stress tensor and its trace (the pressure) reported in
> > molecular
> > dynamics by pwscf and other
> > quantum codes is a static pressure. In a
> > molecular dinamics simulation,
> > if I want the "real" pressure, I must add
> > dp=(N/V) k_B T, N being the number of atoms and V the volume of the
> > cell, and T the instantaneous temperature or its meand value.
> >
> > In other words, the pressure of a non ideal gas (or liquid) is
> >
> > P=(N/V) k_B T + virial term (do not worry about the exact expression),
> >
> > then I understand that the pressure given by the codes is equivalent to
> > the virial term.
> > Is this correct? I found strange that nobody mentions it.
> >
> > Therefore, I understand that to have the "real" pressure I must add
> > (N/V)k_B T to the pressure reported by pw or cp. Please, correct me if I
> > am
> > wrong.
> >
> > Next level of question. If I do Parrinello-Rahman dynamics with a Nose
> > thermostast, what is the pressure specified by the keyword press ? Does
> > it include the thermal pressure?
> > I see that there is also a Nose thermostast for the cell. What is its
> > importance? Is it neccessary to use this cell thermostat for the NPT
> > ensemble?
> >
> Later I saw in the CP source that there is a subroutine
> 'add_thermal_stress', (which seems better than just 'add_thermal_press'),
>
> cpr_mod.f90: subroutine add_thermal_stress( stress, pmass, omega, h, vels, nsp, na )
>
> However, I found no CALL to this subroutine. Then my doubt increased.
>
>
> Eduardo
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
More information about the Pw_forum
mailing list