pressure calculations for liquid Pb
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 1
- Joined: Thu Jul 24, 2014 8:26 am
- License Nr.: I93blcf
pressure calculations for liquid Pb
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.
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.
-
- Global Moderator
- Posts: 1817
- Joined: Mon Nov 18, 2019 11:00 am
Re: pressure calculations for liquid Pb
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