XANES in Diamond: Difference between revisions

From VASP Wiki
Line 53: Line 53:
The super cell can be obtained either by taking the file POSCAR.3x3x3 provided with this tutorial. In the following we show how to get the super cell from the primitive cell using p4vasp:
The super cell can be obtained either by taking the file POSCAR.3x3x3 provided with this tutorial. In the following we show how to get the super cell from the primitive cell using p4vasp:
*Open p4vasp by typing ''p4v'' on the terminal.
*Open p4vasp by typing ''p4v'' on the terminal.
*Load the primitive cell by clicking on '''file'''<&rarr;>
*Load the primitive cell by clicking on '''file'''&rarr;'''Load system''':
[File:Fig XAS 1.png|200px]


== Download ==
== Download ==

Revision as of 15:01, 18 April 2018

(UNDER CONSTRUCTION)

Task

Calculation of the XANES K-edge in diamond using the supercell core-hole method.

Input

POSCAR

cubic diamond
 3.567
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
 2
direct
0.0 0.0 0.0
0.25 0.25 0.25
  • We will not use this structure as input but rather use it to construct a super cell from it that is actually used in the calculations.

INCAR

System = DIAMOND
ALGO = FAST
ISMEAR = 0; SIGMA = 0.1;
ICORELEVEL = 2
CLNT = 1
CLN = 1
CLL = 0
CLZ = 1.0
CH_LSPEC = .TRUE.
CH_SIGMA = 0.5
NBANDS = 300
LREAL = A
  • To promote a core electron into the conduction bands and hence create the core-hole ICORELEVEL=2 has to be set. This corresponds to the final state approximation
  • CLNT=1 selects the first atom species in the POSCAR file.
  • CLN=1 selects main quantum number 1 (hence K-edge).
  • CLL=0 selects angular quantum number 0 (s).
  • CLZ=1.0 selects the charge of the core hole. By setting this number to a fractional value we can mimick different screening of the electrons. Since this purely exploits error cancellation and the physical background of non-integer charges is not defined well, it should be only used with caution.
  • By setting CH_LSPEC=.TRUE. we enable the calculation of matrix elements between core and conduction states and the calculation of the core electron absorption spectrum.
  • The broadening of the core electron absorption spectrum is controlled by the tag CH_SIGMA. Usually it is good practice to set this value low and broaden the spectrum in post processing.
  • We have to set NBANDS to a larger value to consider enough conduction band states in the calculation.
  • Since super cells are used the calculation of the projection operators in real space is much faster, hence LREAL=A is set.

Calculation

Step 1 build a supercell

To perform calculations with a core-hole we need sufficiently large super cells to reduce the interaction of neighbouring core-holes. Usually the obtained spectrum with respect to the size of the super cell has to be converged. This means calculations are done with successively larger cells until no significant change is visible in the spectrum anymore. To save computational time we chose a cell for this tutorial, although for converged values one should use at least .

The super cell can be obtained either by taking the file POSCAR.3x3x3 provided with this tutorial. In the following we show how to get the super cell from the primitive cell using p4vasp:

  • Open p4vasp by typing p4v on the terminal.
  • Load the primitive cell by clicking on fileLoad system:

[File:Fig XAS 1.png|200px]

Download

Back to the main page.