LAECHG

LAECHG = .TRUE. | .FALSE.
Default: LAECHG = .FALSE. 

Description: when LAECHG=.TRUE. the all-electron charge density will be reconstructed explicitly and written out to file.

---

If LAECHG=.TRUE. is set VASP will reconstruct the all-electron charge density on the so-called "fine" FFT-grid. This "fine" FFT-grid consists of NGXF×NGYF×NGZF points in real space (i.e., the grid that is used to represent the augmented pseudo charge densities of the USPP and PAW methods).

In the language of the PAW method an "all-electron" density does not refer to the density of all electrons, instead it denotes a density that includes all the nodal features near the nucleus associated with the true (as opposed to the pseudized) one-electron orbitals. Within the PAW method, the all-electron density arising from a one-electron orbital ${\widetilde {\psi }}_{n\mathbf {k} }$ is given by:

$n(\mathbf {r} )=\langle {\widetilde {\psi }}_{n\mathbf {k} }|\mathbf {r} \rangle \langle \mathbf {r} |{\widetilde {\psi }}_{n\mathbf {k} }\rangle +\sum _{\alpha ,\beta }(\phi _{\alpha }^{\ast }(\mathbf {r} )\phi _{\beta }(\mathbf {r} )-{\widetilde {\phi }}_{\alpha }^{\ast }(\mathbf {r} ){\widetilde {\phi }}_{\beta }(\mathbf {r} ))\langle {\widetilde {\psi }}_{n\mathbf {k} }|{\widetilde {p}}_{\alpha }\rangle \langle {\widetilde {p}}_{\beta }|{\widetilde {\psi }}_{n\mathbf {k} }\rangle$

Normally one does not attempt to reconstruct this density explicitly since the second term on the right hand side varies rapidly near the nuclei and is too costly to expand in plane waves (see the bit about augmentation and compensation charges.

In fact, for LAECHG=.TRUE., VASP will reconstruct three distinct all-electron densities:

1. the core density.
2. the proto-atomic valence density (overlapping atomic densities).
3. the self-consistent valence density.

These are written to the files AECCAR0, AECCAR1, and AECCAR2, respectively.

Related Tags and Sections

AECCAR0, AECCAR1, AECCAR2

Examples that use this tag

---

Contents