Category:Dielectric properties: Difference between revisions

Introduction

Optics - response to a varying E-field (absorption, reflectance, MOKE)

Methods based on screened interaction in solids (MBPT)

Most general definition of epsilon D=epsilon E

Methods for computing ${\displaystyle \epsilon}$

Static response: Density functional perturbation Theory (DFPT) and Finite differences based methods

LEPSILON

By setting LEPSILON=.True., VASP uses DFPT to compute the static ion-clamped dielectric matrix with or without local field effects. Derivatives are evaluated using Sternheimer equations, avoiding the explicit computation of derivatives of the periodic part of the wave function. This method does not require the inclusion of empty states via the NBANDS parameter.

At the end of the calculation the both the values of ${\displaystyle \epsilon}$ including (LRPA=.True.) or excluding (LRPA=.False.) local-field effects are printed in the OUTCAR file. Users can perform a consistency check by comparing the values with no local field to the zero frequency results for ${\displaystyle \epsilon}$ obtained from a calculation with LOPTICS=.True..

LCALCEPS

With LCALCEPS=.True., the dielectric tensor is computed from the derivative of the polarization, using

${\displaystyle \epsilon^\infty_{ij}=\delta_{ij}+ \frac{4\pi}{\epsilon_0}\frac{\partial P_i}{\partial \mathcal{E}_j} \qquad {i,j=x,y,z}. }$

However, here the derivative is evaluated explicitly by employing finite-differences. The direction and intensity of the perturbing electric field has to be specified in the INCAR using the EFIELD_PEAD variable. As with the previous method, at the end of the calculation VASP will write the dielectric tensor in the OUTCAR file. Control over the inclusion of local-field effects is done with the variable LRPA.

Dynamical response

LOPITCS

The variable LOPTICS allows for the calculation of the frequency dependent dielectric function once the ground state is computed. It uses the explicit expression to evaluate the imaginary part of ${\displaystyle \epsilon}$

${\displaystyle \epsilon^{(2)}_{\alpha \beta}\left(\omega\right) = \frac{4\pi^2 e^2}{\Omega} \mathrm{lim}_{q \rightarrow 0} \frac{1}{q^2} \sum_{c,v,\mathbf{k}} 2 w_\mathbf{k} \delta( \epsilon_{c\mathbf{k}} - \epsilon_{v\mathbf{k}} - \omega) \times \langle u_{c\mathbf{k}+\mathbf{e}_\alpha q} | u_{v\mathbf{k}} \rangle \langle u_{v\mathbf{k}} | u_{c\mathbf{k}+\mathbf{e}_\beta q} \rangle, }$

while the real part is evaluated using the Kramers-Kroning relation. At this level there are no effects coming from local fields.

This method requires a relatively large number of empty states, controlled by the variable NBANDS in the INCAR file and it should be checked for convergence.

Furthermore, the INCAR should also include values for CSHIFT (the broadening applied to the Lorentzian function which replaces the ${\displaystyle \delta}$-function), and NEDOS (the frequency grid for ${\displaystyle \omega}$).

ALGO = TDHF

This option performs a time-dependent Hartree-Fock or DFT calculation. It follows the Casida equation and uses a Fourier transform of the time-evolving dipoles to compute ${\displaystyle \epsilon}$.

The number of NBANDS controls how many bands are present in the time-evolution. This does not need to be as high as when LOPTICS is active, and only a few bands above the band gap need to be included.

The choice of time-dependent kernel is controlled by AEXX, HFSCREEN, and LFXC variables. For calculations using hybrid functionals, AEXX controls the fraction of exact exchange used in the exchange correlation potential, while HFSCREEN specifies the range-separation parameter. For a pure TDDFT calculation, LFXC uses the local exchange-correlation kernel in the time-evolution equations.

ALGO = TIMEEV

Uses a delta-pulse electric field to probe all transitions and calculate the dielectric function by following the evolution in time of the dipole momenta. This algorithm is able to fully reproduce the absorption spectra from standard Bethe-Salpeter calculations by setting the correct time-dependent kernel with LHARTREE=.True. and LFXC=.True.

The time step is controlled automatically by the CSHIFT and PREC variables. This means that the smaller the value of CSHIFT and the more accurate the level of precision chosen by the user, the higher the number of time steps that VASP will perform, and so the higher the cost of the calculation.

The number of valence and conduction bands involved in the time propagation are set by the NBANDSO and NBANDSV variables, respectively. Once again, users are advised to choose a small number of bands near the band gap if they wish to reproduce optical measurements.

Finally, the maximum energy used in both the Fourier transform and in calculating the frequency dependent dielectric function is set by OMEGAMAX, and the sampling of the frequency grid is controlled by NEDOS.

ALGO = CHI

Here the frequency dielectric function is computed at the independent particle level but starting from a GW calculation. VASP will compute the polarizability ${\displaystyle \chi}$ and then use Failed to parse (Conversion error. Server ("cli") reported: "[INVALID]"): {\displaystyle \epsilon_{\mathbf G\mathbf G'}(\mathbf q,\omega) = \delta_{\mathbf G\mathbf G'} - v(\mqthbf q+\mathbf G)\chi_{\mathbf G\mathbf G'}(\mathbf q,\omega)} .

within the Random-Phase Approximation (RPA) by setting ALGO=Chi.

ALGO = BSE

Setting ALGO=BSE computes the macroscopic dielectric function ${\displaystyle \epsilon_M}$ by solving the Bethe-Salpeter equations. Here the electron-hole pairs are treated a new quasi-particle, an exciton, and the dielectric function is built using the new eigenvectors (${\displaystyle X_{\lambda}^{cv\mathbf k}}$) and eigenvalues (${\displaystyle \omega_\lambda}$)

${\displaystyle \epsilon_M(\mathbf{q},\omega)= 1+v(\mathbf q)\sum_{\lambda\lambda'} \sum_{c,v,\mathbf k}\sum_{c',v',\mathbf k'}\langle c\mathbf{k}|e^{i\mathbf{qr}}|v\mathbf{k}\rangle X_\lambda^{cv\mathbf{k}} \langle c'\mathbf{k'}|e^{-i\mathbf{qr}}|v'\mathbf{k'}\rangle X_{\lambda'}^{c'v'\mathbf{k}',*}\times S^{-1}_{\lambda,\lambda'} \left(\frac{1}{\omega_\lambda - \omega - i\delta} + \frac{1}{\omega_\lambda+\omega + i\delta}\right)~. }$

where ${\displaystyle S_{\lambda\lambda'}}$ is the overlap between exciton states of indices ${\displaystyle \lambda}$ and ${\displaystyle \lambda'}$ (in general the BSE Hamiltonian is not hermitian, so eigenstates associated to different eigenvalues are not necessarily orthogonal).

The number of valence and conduction states which are included in the BSE Hamiltonian is controlled by the variables NBANDSO and NBANDSV, respectively. Note that normally only a few bands above and below the band gap are needed to converge the optical spectrum, so users should be careful in setting up these two variables. Otherwise the calculation might run out of memory.

For comparison with optical experiments (e.g. absorption, MOKE, reflectance), ${\displaystyle q}$ is the photon momentum and usually it is taken in the ${\displaystyle \mathbf q \to 0}$ limit. Furthermore, the coupling between the resonant and anti-resonant terms can be switched off, in what is called the Tamm-Dancoff approximation. This approximation can be activated with the variable ANTIRES set to 0. Setting this variable to 1 or 2 will include the coupling, but increase the computational cost.

Level of approximation

Micro-macro connection - Including local fields (GG' vs GG or 00, to confirm)

Inhomogeniety - Long wavelength vs short

Local fields in the Hamiltonian

Ion-clamped vs relaxed/dressed dielectric function

Static vs dynamic

Density-density versus current-current response functions

Relation to observables

Polarizability

Optical conductivity

Optical absorption

Reflectance

MOKE

Combination with other perturbations

Atomic displacements

Strain