Author Topic: Presure calculations  (Read 2968 times)

0 Members and 1 Guest are viewing this topic.

Offline GJK

  • QuantumATK Guru
  • ****
  • Posts: 198
  • Country: in
  • Reputation: 0
    • View Profile
Presure calculations
« on: April 12, 2015, 12:33 »
Hi,

It is well known that the enthalpy is given by H=U-pV , we can get the enthalpy from total energy and volume from lattice parameters  but how to get the pressure so that we can easily find out internal internal energy and hence calculate the behavior of two interacting systems.

Kindly  provide the codes for calculating pressure if possible?

Thanks & Regards
GJK

Offline Julian Schneider

  • QuantumATK Staff
  • QuantumATK Guru
  • *****
  • Posts: 164
  • Country: dk
  • Reputation: 25
    • View Profile
Re: Presure calculations
« Reply #1 on: April 14, 2015, 09:36 »
For a static system (i.e. obtained from a geometry optimization), the pressure can be calculated from the stress tensor via P = -1/3 (s_xx + s_yy + s_zz) (see attached script).
At finite temperature (i.e. in a molecular dynamics run) there would be an additional contribution of the kinetic energy: P = 2/3 E_kin/V - 1/3 (s_xx + s_yy + s_zz).

If you want to compare the enthalpy of two systems, you have to simulate them at the desired target pressure. In a molecular dynamics simulation, this can easily be achieved by employing one the NPT-methods. If you intend to obtain the structures via a geometry optimization, you have to optimize the cell until the desired target stress is reached. In ATK2015, this can easily be achieved via the target_stress parameter.
« Last Edit: April 14, 2015, 09:40 by Julian Schneider »