pressure calculations for liquid Pb

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Locked
Message
Author
dkosov
Newbie
Newbie
Posts: 1
Joined: Thu Jul 24, 2014 8:26 am
License Nr.: I93blcf

pressure calculations for liquid Pb

#1 Post by dkosov » Thu Jul 24, 2014 8:46 am

Hello

I am trying to compute pressure for liquid Pb running MD. I have 216 Pb atoms in the cubic box os side 19.25 Ang - which gives the correct liquid lead density at T=800K. The system is well equilibrated (20 ps VASP equilibration starting from the melted structure obtained via classical MD). I run Nose dynamics at T=800 and compute the average pressure along MD run.

When I use LDA the pressure is too low ~ -14 kB (should be physically zero for this T and density), if I use the same structure but employ PBE functional the pressure becomes ~+15kB (way too large). It looks wrong or I missing some physics here.

I would appreciate any help and advices,
Many thanks, Daniel.
Last edited by dkosov on Thu Jul 24, 2014 8:46 am, edited 1 time in total.

support_vasp
Global Moderator
Global Moderator
Posts: 1817
Joined: Mon Nov 18, 2019 11:00 am

Re: pressure calculations for liquid Pb

#2 Post by support_vasp » Thu Sep 12, 2024 8:40 am

Hi,

We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.

Best wishes,

VASP


Locked