PSTRESS: Difference between revisions

From VASP Wiki
No edit summary
(Updated some mistakes in sign conventions. Improved wording. Explicitly discuss tensile and compressive stress and how this changes volume.)
Line 1: Line 1:
{{TAGDEF|PSTRESS|[real]}}
{{TAGDEF|PSTRESS|[real]}}


Description: Set the pressure in kB and add corrections to the stress tensor.
Description: Sets the external pressure in kB or add corrections to the stress tensor.


----
----


During [[Ionic minimization|ionic minimization]], an energy term <math> E= V * \mathrm{PSTRESS} </math> is added to the total energy and the pressure set by {{TAG|PSTRESS}} is added to the diagonals of the stress tensor. This allows performing structure optimization at a specific external pressure.  
The unit of {{TAG|PSTRESS}} is kB.


In [[Molecular dynamics|molecular-dynamics calculations]] within the {{TAG|NpT ensemble}}, {{TAG|PSTRESS}} controls the target pressure for the Parinello-Rahman{{cite|parrinello:prl:1980}}{{cite|parrinello:jap:1981}} barostat. The unit of {{TAG|PSTRESS}} is kB.
During [[Ionic minimization|ionic minimization]], an energy term <math> E= V \times \mathrm{PSTRESS} </math> is added to the total energy and the value of {{TAG|PSTRESS}} is subtracted from  the diagonals of the stress tensor. This allows to perform structure optimization at a specific external pressure.


{{TAG|PSTRESS}} can also be used to correct errors caused by Pulay stress, i.e., errors in the calculated stress tensor caused by the incomplete basis set. To correct for Pulay-stress errors, set {{TAG|PSTRESS}} to the negative value of Pulay stress. The Pulay stress is computed by taking the difference between the external pressure at a very large cutoff and the desired cutoff. Look for the line 'external pressure' in the {{TAG|OUTCAR}} file. Before using this tag in this way, please read the section {{TAG | Energy vs volume Volume relaxations and Pulay stress}}.
In [[Molecular dynamics|molecular-dynamics calculations]] within the {{TAG|NpT ensemble}}, {{TAG|PSTRESS}} controls the target pressure for the Parinello-Rahman{{cite|parrinello:prl:1980}}{{cite|parrinello:jap:1981}} barostat.
 
Generally, if a negative value is supplied, the system is under effective tensile strain and during relaxations (or molecular dynamics simulations) the volume
will increase. If a positive value is supplied, the system is under compressive strain; this will decrease the volume during relaxations and molecular dynamics simulations.
 
{{TAG|PSTRESS}} can also be used to correct errors caused by the Pulay stress, i.e., errors in the calculated stress tensor caused by the incomplete plane wave basis set. To correct for Pulay-stress errors, set {{TAG|PSTRESS}} to the negative value of the Pulay stress. The Pulay stress is computed by taking the difference between the external pressure at the desired cutoff and a very large energy cutoff (check the lines 'external pressure' in the {{TAG|OUTCAR}} file and calculate
pressure at cutoff you want to use <math> -  </math> pressure at huge cutoff; this must be a negative value). Before using this tag in this way, please read the following section carefully {{TAG | Energy vs volume Volume relaxations and Pulay stress}}.


{{sc|PSTRESS|Examples|Examples that use this tag}}
{{sc|PSTRESS|Examples|Examples that use this tag}}

Revision as of 18:39, 21 March 2023

PSTRESS = [real] 

Description: Sets the external pressure in kB or add corrections to the stress tensor.


The unit of PSTRESS is kB.

During ionic minimization, an energy term is added to the total energy and the value of PSTRESS is subtracted from the diagonals of the stress tensor. This allows to perform structure optimization at a specific external pressure.

In molecular-dynamics calculations within the NpT ensemble, PSTRESS controls the target pressure for the Parinello-Rahman[1][2] barostat.

Generally, if a negative value is supplied, the system is under effective tensile strain and during relaxations (or molecular dynamics simulations) the volume will increase. If a positive value is supplied, the system is under compressive strain; this will decrease the volume during relaxations and molecular dynamics simulations.

PSTRESS can also be used to correct errors caused by the Pulay stress, i.e., errors in the calculated stress tensor caused by the incomplete plane wave basis set. To correct for Pulay-stress errors, set PSTRESS to the negative value of the Pulay stress. The Pulay stress is computed by taking the difference between the external pressure at the desired cutoff and a very large energy cutoff (check the lines 'external pressure' in the OUTCAR file and calculate pressure at cutoff you want to use pressure at huge cutoff; this must be a negative value). Before using this tag in this way, please read the following section carefully Energy vs volume Volume relaxations and Pulay stress.

Examples that use this tag