LVACPOTAV: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(21 intermediate revisions by 2 users not shown)
Line 1: Line 1:
'''VACPOTAV''' = .TRUE. | .FALSE.  
{{TAGDEF|LVACPOTAV|.True. {{!}} .False.|.False.}}


Default: '''VACPOTAV''' = .FALSE.  
Description: Switch on determining the vacuum potential by averaging the local potential over a field-free region.  
----


This tag switches on the computation of the vacuum potential; i.e. the average of the local potential in the vacuum region. It computes the average potential by iterating over regions of the 2D-averaged charge density (averaging done in the direction of [[IDIPOL]]) and tries to find a region in the cell where the charge density is nearly zero (corresponding to a field-free region in a cell without any compensating background charge). The potential is averaged over this region and is reported as the vacuum potential.
{{TAG|LVACPOTAV}} switches on the computation of the vacuum potential, i.e., the average of the local potential in the vacuum region.  
It computes the average potential by searching regions that are field-free ({{TAG|VACPOTFLAT}}), and the 2D-averaged charge density is nearly zero. The averaging is done in the direction of {{TAG|IDIPOL}} and is reported as the vacuum potential in the {{FILE|OUTCAR}}.
{{NB|tip| As {{TAG|LVACPOTAV}} performs a post-processing step, you may use it together with {{TAG|ALGO}} {{=}} None by restarting a converged calculation.}}
The vacuum potential is one of the quantities needed to [[Computing the work function|compute the work function]]. It can be extracted from the {{FILE|OUTCAR}} by the following bash command


{{NB|warning| This tag is available only for versions after 6.4.3. Prior to this version, the default algorithm used to determine the vacuum potential is used. The default algorithm simply reports the slice-averaged potential in all directions apart from {{TAG|IDIPOL}} a few grid points on either direction of the minimum 2D-averaged charge density.}}
  grep upper OUTCAR


The vacuum potential is one of the quantities needed to compute the workfunction (see [[Computing the work function|here]] for a step-by-step guide on how to compute the workfunction with VASP). The vacuum potentials are written to the OUTCAR file and can be accessed with the following command
Note that two vacuum potentials will be produced, one corresponding to the upper termination of the slab and one corresponding to the lower. Depending on the system, one might be more interesting than the other.
 
{{NB|tip|For determining the work function, we suggest using {{TAG|LVACPOTAV}} along with the {{TAG|LVHAR}} tag such that only the sum of the Hartree and ionic potentials are used in the calculation of the vacuum potential. This choice is because the exchange-correlation potential might be noisy in the vacuum region but should, in principle, be zero.}}
<code>grep upper OUTCAR</code>
{{NB|mind|{{TAG|LVACPOTAV}} is available only for versions after 6.4.3.}}
 
Before VASP 6.4.3, the default algorithm reports the 2D-averaged potential four grid points from the the minimum 2D-averaged charge density in the direction of {{TAG|IDIPOL}}, i.e., no averaging is performed along the surface normal of the 2D-averaged potential.
which will produce the vacuum potentials. Note that two vacuum potentials will be produced, one corresponding to the upper termination of the slab and one corresponding to the lower. Depending on the system you are studying, one might be more interesting than another.
=== Use in conjunction with the dipole correction ===
{{NB|tip| In case you would like to determine the workfunction, we suggest using '''VACPOTAV''' along with the [[LVHAR]] tag such that only the sum of the Hartree and Ionic potentials are used in the calculation of the vacuum potential}}
A typical use case for {{TAG|LVACPOTAV}} is together with the dipole correction (including tags {{TAG|LDIPOL}} and {{TAG|IDIPOL}}). Switching on the dipole correction is crucial for determining the vacuum potential; without it, there will be no field-free region for dipolar systems.
{{NB|mind| Note that {{TAG|LVACPOTAV}} is currently implemented for {{TAG|IDIPOL}} between 1 and 3.}}


=== Warnings ===
=== Warnings ===
In case '''VACPOTAV''' is not able to generate an accurate workfunction, the following warnings may be found in the [[OUTCAR]] file.  
In case {{TAG|LVACPOTAV}} is not able to generate an accurate work function, the following warnings may be found in the {{FILE|OUTCAR}} file.  


==== Vacuum region is likely too small ====
==== Vacuum region is likely too small ====
  -----------------------------------------------------------------------------
|                                                                            |
|          W    W    AA    RRRRR  N    N  II  N    N  GGGG  !!!          |
|          W    W  A  A  R    R  NN  N  II  NN  N  G    G  !!!          |
|          W    W  A    A  R    R  N N  N  II  N N  N  G      !!!          |
|          W WW W  AAAAAA  RRRRR  N  N N  II  N  N N  G  GGG  !            |
|          WW  WW  A    A  R  R  N  NN  II  N  NN  G    G                |
|          W    W  A    A  R    R  N    N  II  N    N  GGGG  !!!          |
|                                                                            |
  |    Did not find any points to average over, which means that no vacuum    |
  |    Did not find any points to average over, which means that no vacuum    |
  |    field-free region was found. Please increase the size of  your cell    |
  |    field-free region was found. Please increase the size of  your cell    |
  |    in the dimension of the dipole correction to obtain accurate            |
  |    in the dimension of the dipole correction to obtain accurate            |
  |    workfunction values.                                                    |
  |    workfunction values.                                                    |
|                                                                            |
 
  -----------------------------------------------------------------------------
A possible solution to this problem is to increase the size of the vacuum dimension in your cell.
A possible solution to this problem is to increase the size of the vacuum dimension in your cell.


==== The minimum charge density in your cell may be too large ====
==== The minimum charge density in your cell may be too large ====
  -----------------------------------------------------------------------------
|                                                                            |
|          W    W    AA    RRRRR  N    N  II  N    N  GGGG  !!!          |
|          W    W  A  A  R    R  NN  N  II  NN  N  G    G  !!!          |
|          W    W  A    A  R    R  N N  N  II  N N  N  G      !!!          |
|          W WW W  AAAAAA  RRRRR  N  N N  II  N  N N  G  GGG  !            |
|          WW  WW  A    A  R  R  N  NN  II  N  NN  G    G                |
|          W    W  A    A  R    R  N    N  II  N    N  GGGG  !!!          |
|                                                                            |
  |    The minimum charge density times volume of the cell along the axis      |
  |    The minimum charge density times volume of the cell along the axis      |
  |    of the dipole correction is larger 1E-1, which could mean that your    |
  |    of the dipole correction is larger 1E-1, which could mean that your    |
Line 51: Line 37:
  |    along the dipole correction (vacuum dimension) or perhaps              |
  |    along the dipole correction (vacuum dimension) or perhaps              |
  |    increasing the precision of your calculation.                          |
  |    increasing the precision of your calculation.                          |
|                                                                             |
 
  -----------------------------------------------------------------------------
Possible solutions include:
 
* Making sure you have a large enough vacuum dimension.
* Increasing the precision of your calculation by changing {{TAG|EDIFF}}.
 
== Related tags and articles ==
{{TAG|DIPOL}},
{{TAG|LDIPOL}},
{{TAG|IDIPOL}},
{{TAG|VACPOTFLAT}},
{{TAG|WRT_POTENTIAL}},
{{TAG|LVTOT}},
{{TAG|LVHAR}}
 
[[Computing the work function]]
 
[[Category:Electrostatics]][[Category:INCAR tag]]

Latest revision as of 09:24, 9 February 2024

LVACPOTAV = .True. | .False.
Default: LVACPOTAV = .False. 

Description: Switch on determining the vacuum potential by averaging the local potential over a field-free region.


LVACPOTAV switches on the computation of the vacuum potential, i.e., the average of the local potential in the vacuum region. It computes the average potential by searching regions that are field-free (VACPOTFLAT), and the 2D-averaged charge density is nearly zero. The averaging is done in the direction of IDIPOL and is reported as the vacuum potential in the OUTCAR.

Tip: As LVACPOTAV performs a post-processing step, you may use it together with ALGO = None by restarting a converged calculation.

The vacuum potential is one of the quantities needed to compute the work function. It can be extracted from the OUTCAR by the following bash command

 grep upper OUTCAR

Note that two vacuum potentials will be produced, one corresponding to the upper termination of the slab and one corresponding to the lower. Depending on the system, one might be more interesting than the other.

Tip: For determining the work function, we suggest using LVACPOTAV along with the LVHAR tag such that only the sum of the Hartree and ionic potentials are used in the calculation of the vacuum potential. This choice is because the exchange-correlation potential might be noisy in the vacuum region but should, in principle, be zero.
Mind: LVACPOTAV is available only for versions after 6.4.3.

Before VASP 6.4.3, the default algorithm reports the 2D-averaged potential four grid points from the the minimum 2D-averaged charge density in the direction of IDIPOL, i.e., no averaging is performed along the surface normal of the 2D-averaged potential.

Use in conjunction with the dipole correction

A typical use case for LVACPOTAV is together with the dipole correction (including tags LDIPOL and IDIPOL). Switching on the dipole correction is crucial for determining the vacuum potential; without it, there will be no field-free region for dipolar systems.

Mind: Note that LVACPOTAV is currently implemented for IDIPOL between 1 and 3.

Warnings

In case LVACPOTAV is not able to generate an accurate work function, the following warnings may be found in the OUTCAR file.

Vacuum region is likely too small

|     Did not find any points to average over, which means that no vacuum     |
|     field-free region was found. Please increase the size of  your cell     |
|     in the dimension of the dipole correction to obtain accurate            |
|     workfunction values.                                                    |

A possible solution to this problem is to increase the size of the vacuum dimension in your cell.

The minimum charge density in your cell may be too large

|     The minimum charge density times volume of the cell along the axis      |
|     of the dipole correction is larger 1E-1, which could mean that your     |
|     workfunction is not accurate as there is no field free region in        |
|     your cell. Please consider either increasing the size of your cell      |
|     along the dipole correction (vacuum dimension) or perhaps               |
|     increasing the precision of your calculation.                           |

Possible solutions include:

  • Making sure you have a large enough vacuum dimension.
  • Increasing the precision of your calculation by changing EDIFF.

Related tags and articles

DIPOL, LDIPOL, IDIPOL, VACPOTFLAT, WRT_POTENTIAL, LVTOT, LVHAR

Computing the work function