ICORELEVEL: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 11: Line 11:
The core-level binding energies can be calculated either in the initial-state approximation or the final-state approximation. In the initial-state approximation a core electron is removed from the core states and added to the valence/conduction bands but no change of the potential (by e.g. relaxation of the valence electrons) is allowed. The core-level binding energy can then be directly calculated by the Kohn-Sham eigenvalues{{cite|lizzit:prb:2001}} of the core level <math>\epsilon_{c}</math> and the Fermi energy <math>\epsilon_{F}</math>
The core-level binding energies can be calculated either in the initial-state approximation or the final-state approximation. In the initial-state approximation a core electron is removed from the core states and added to the valence/conduction bands but no change of the potential (by e.g. relaxation of the valence electrons) is allowed. The core-level binding energy can then be directly calculated by the Kohn-Sham eigenvalues{{cite|lizzit:prb:2001}} of the core level <math>\epsilon_{c}</math> and the Fermi energy <math>\epsilon_{F}</math>
:<math>E_{CL}^{\mathrm{i}}=\epsilon_{c} - \epsilon_{F}</math>.
:<math>E_{CL}^{\mathrm{i}}=\epsilon_{c} - \epsilon_{F}</math>.
Th
In the final-state approximation the electrons (valence electrons in the frozen-core approximation) are allowed to relax, so that the thus created local hole is screened. In other words a fully self-consistent calculations is carried out with a core hole and an additional electron in the valence/conduction bands.
 
 
For {{TAG|ICORELEVEL}} = 0, which is the default, the core energies are not calculated.
 
VASP supports two options to calculate changes in the core level energies.
A detailed documentation of the implementation is presently missing, we
however refer to {{cite|koehler:prb:04}} and references therein for
an overview.
 
The simpler option ({{TAG|ICORELEVEL}} = 1) calculates the  core
levels in the initial state approximation,
which just involves recalculating the KS eigenvalues of the core states
after a self-consistent calculation of the valence charge density.
{{TAG|ICORELEVEL}} = 1 is a little bit more involved than the calculations using
{{TAG|LVTOT}} = ''.TRUE.'', since the Kohn-Sham energy of each core state
is recalculated.
This adds very little extra cost to the calculations; usually
the shifts correspond very closely to the change of
the electrostatic potential at the lattice sites (calculated
using {{TAG|LVTOT}} = ''.TRUE.'').
 
The second option  ({{TAG|ICORELEVEL}} = 2) is more involved.
In this case, electrons are removed from the core and placed
into the valence (effectively increasing {{TAG|NELECT}}).
The vasp implementation excites all selected core electrons for
all atoms of one species. The species as well as the
selected electrons are specified using


The following options are available in VASP:
*{{TAG|ICORELEVEL}}=0: The core energies are not calculated (default).
*{{TAG|ICORELEVEL}}=1: The initial-state approximation is used. This just involves recalculating the KS eigenvalues of the core states
after a self-consistent calculation of the valence charge density. {{TAG|ICORELEVEL}}=1 is a little bit more involved than the calculations using
{{TAG|LVTOT}}=''.TRUE.'', since the Kohn-Sham energy of each core state is recalculated. This adds very little extra cost to the calculations; usually
the shifts correspond very closely to the change of the electrostatic potential at the lattice sites (calculated using {{TAG|LVTOT}}=''.TRUE.'').
*{{TAG|ICORELEVEL}}=2: The final-state approximation is used. Electrons are removed from the core and placed into the valence (effectively increasing {{TAG|NELECT}}). The vasp implementation excites all selected core electrons for
all atoms of one species. The species as well as the selected electrons are specified using
  {{TAGBL|CLNT}} = species  
  {{TAGBL|CLNT}} = species  
  {{TAGBL|CLN}} =  main quantum number of excited core electron  
  {{TAGBL|CLN}} =  main quantum number of excited core electron  
  {{TAGBL|CLL}} =  l quantum number of excited core electron
  {{TAGBL|CLL}} =  l quantum number of excited core electron
  {{TAGBL|CLZ}} =  electron count
  {{TAGBL|CLZ}} =  electron count
 
The electron count {{TAGBL|CLZ}} specifies how many electrons are excited
The electron count specifies how many electrons are excited
from the core. Usually 1 or 0.5 (Slaters transition state) are sensible choices.
from the core. Usually 1 or 0.5 (Slaters transition state) are sensible choices.
{{TAG|CLNT}} selects for which species in the {{TAG|POTCAR}} file the electrons
{{TAG|CLNT}} selects for which species in the {{TAG|POTCAR}} file the electrons

Revision as of 10:27, 8 April 2022

ICORELEVEL = 0 | 1 | 2
Default: ICORELEVEL = 0 

Description: ICORELEVEL controls whether the core energies are explicitely calculated or not and how they are calculated.


The binding energy of core electrons is given as

.

Here, is the energy from a standard density-functional calculation in which the number of core electrons corresponds to the unexcited ground state. ist the energy of a calculation where one electron is removed from the core of one particular atom and added to the valence or conduction band.

The core-level binding energies can be calculated either in the initial-state approximation or the final-state approximation. In the initial-state approximation a core electron is removed from the core states and added to the valence/conduction bands but no change of the potential (by e.g. relaxation of the valence electrons) is allowed. The core-level binding energy can then be directly calculated by the Kohn-Sham eigenvalues[1] of the core level and the Fermi energy

.

In the final-state approximation the electrons (valence electrons in the frozen-core approximation) are allowed to relax, so that the thus created local hole is screened. In other words a fully self-consistent calculations is carried out with a core hole and an additional electron in the valence/conduction bands.

The following options are available in VASP:

  • ICORELEVEL=0: The core energies are not calculated (default).
  • ICORELEVEL=1: The initial-state approximation is used. This just involves recalculating the KS eigenvalues of the core states

after a self-consistent calculation of the valence charge density. ICORELEVEL=1 is a little bit more involved than the calculations using LVTOT=.TRUE., since the Kohn-Sham energy of each core state is recalculated. This adds very little extra cost to the calculations; usually the shifts correspond very closely to the change of the electrostatic potential at the lattice sites (calculated using LVTOT=.TRUE.).

  • ICORELEVEL=2: The final-state approximation is used. Electrons are removed from the core and placed into the valence (effectively increasing NELECT). The vasp implementation excites all selected core electrons for

all atoms of one species. The species as well as the selected electrons are specified using

CLNT = species 
CLN =  main quantum number of excited core electron 
CLL =  l quantum number of excited core electron
CLZ =  electron count

The electron count CLZ specifies how many electrons are excited from the core. Usually 1 or 0.5 (Slaters transition state) are sensible choices. CLNT selects for which species in the POTCAR file the electrons are excited. Usually one would like to excite the electrons for only one atom, this requires to change the POSCAR and POTCAR file, such that the selected atom corresponds to one species in the POTCAR file. i.e. if the calculation invokes a supercell with 64 atoms of one type, the selected atom needs to be singled out, and the POSCAR file will than contain 63 "standard" atoms as well as one special species, at which the excited core hole will be placed (the POTCAR file will hold two identical PAW datasets in this case).

Several caveats apply to this mode. First the excited electron is always spherical, multipole splittings are not available. Second, the other core electrons are not allowed to relax, which might cause a slight error in the calculated energies. Third, absolute energies are not meaningful, since VASP usually reports valence energies only. Only relative shifts of the core electron binding energies are relevant (in some cases, the VASP total energies might become even positive).

Supercell core-hole method

ICORELEVEL=2 and it's related tags are necessary for the calculation of X-ray absorption spectra (XAS) using the super-cell core-hole method.

A description how to set up super-cell core-hole calculations is given here.

A tutorial for the calculation of XAS is here.

Related tags and articles

CLNT, CLN, CLL, CLZ, CH_LSPEC, CH_SIGMA, CH_NEDOS

Examples that use this tag

References