NELECT: Difference between revisions
| (2 intermediate revisions by the same user not shown) | |||
| Line 17: | Line 17: | ||
For atoms or molecules surrounded by vacuum, ε takes on the vacuum value ε=1. | For atoms or molecules surrounded by vacuum, ε takes on the vacuum value ε=1. | ||
VASP can automatically correct this leading error, see {{TAG|Electrostatic corrections}} for further information. | VASP can automatically correct this leading error, see {{TAG|Electrostatic corrections}} for further information. | ||
<!-- Consult the following pages for information about how to deal with charged materials | |||
* [[Dipole corrections for defects and charged defects| Dipole correction for defects and charged defects]] when dealing with bulk charged defects | |||
* [[Charged systems with density functional theory| Charged systems with density functional theory]] to understand why computing energies of charged materials requires special care. | |||
--> | |||
It is important to emphasize that the total energy cannot be corrected for charged slabs, since a charged slab results in an electrostatic potential that grows with the distance from the slab. | It is important to emphasize that the total energy cannot be corrected for charged slabs, since a charged slab results in an electrostatic potential that grows with the distance from the slab. | ||
This non-convergence of the potential is a result of the interaction between the charged slab and the compensating background. | This non-convergence of the potential is a result of the interaction between the charged slab and the compensating background. | ||
| Line 31: | Line 34: | ||
{{TAG|LDIPOL}}, | {{TAG|LDIPOL}}, | ||
{{TAG|EFIELD}} | {{TAG|EFIELD}} | ||
<!-- | |||
* [[Dipole corrections for defects and charged defects| Dipole correction for defects and charged defects]] when dealing with bulk charged defects | |||
* [[Charged systems with density functional theory| Charged systems with density functional theory]] to understand why computing energies of charged materials requires special care. | |||
--> | |||
{{sc|NELECT|Examples|Examples that use this tag}} | {{sc|NELECT|Examples|Examples that use this tag}} | ||
[[Category:INCAR tag]][[Category:Electronic occupancy]][[Category:Electrostatics]] | [[Category:INCAR tag]][[Category:Electronic occupancy]][[Category:Electrostatics]] | ||
Latest revision as of 11:57, 17 October 2024
NELECT = [real]
Default: NELECT = number of valence electrons
Description: NELECT sets the number of electrons.
| Warning: Unless you would like to perform a charged calculation, you should not set this line. |
The number of electrons is determined automatically from the POTCAR (ZVAL of the element) and POSCAR file (number of the atoms of the respective atom type) assuming the cell is charge-neutral. If the number of electrons is not compatible with the number derived from the valence and the number of atoms a homogeneous background charge is assumed. If the number of ions specified in the POSCAR file is 0 and NELECT=n, then the energy of a homogeneous electron gas is calculated.
Charged calculations
| Warning: Using the NELECT tag without the use of an appropriate correction leads to a very slow convergence of energies with respect to the size L of the supercell. |
The required first-order correction to the energy caused by an excess charge is given by
where q is the net charge of the system, α the Madelung constant of a point charge q placed in a homogeneous background charge -q, and ε the dielectric constant of the system. For atoms or molecules surrounded by vacuum, ε takes on the vacuum value ε=1. VASP can automatically correct this leading error, see Electrostatic corrections for further information. It is important to emphasize that the total energy cannot be corrected for charged slabs, since a charged slab results in an electrostatic potential that grows with the distance from the slab. This non-convergence of the potential is a result of the interaction between the charged slab and the compensating background. A practical consequence of this non-convergence is that the total energy depends linearly on the width of the vacuum.
| Warning: Presently, no simple a posteriori correction scheme is implemented in VASP for slab calculations. Total energies from charged slab calculations must be used with care. In certain cases, relative energies between two charged surface calculations may be useful |
Related tags and articles
Monopole Dipole and Quadrupole corrections, EPSILON, DIPOL, IDIPOL, LMONO, LDIPOL, EFIELD