Category:Dielectric properties: Difference between revisions

From VASP Wiki
 
(29 intermediate revisions by 3 users not shown)
Line 1: Line 1:
==Dielectric function==
==Dielectric function==
[[File:exciton.png|thumbnail|400px|When light excites an electron from the valence to the conduction band, the created hole can bind with the electron to form an exciton.]]
When an external electric field <math>\mathbf E</math> acts on a medium, both the electronic and ionic charges will react to the perturbing field. For dielectric materials, in a very simplistic approach, one can think that the bound charges will create dipoles inside the medium leading to an induced polarization, <math>\mathbf P</math>. The combined effects of both fields are expressed in the electric displacement field <math>\mathbf D</math>, given by
When an external electric field <math>\mathbf E</math> acts on a medium, both the electronic and ionic charges will react to the perturbing field. For dielectric materials, in a very simplistic approach, one can think that the bound charges will create dipoles inside the medium leading to an induced polarization, <math>\mathbf P</math>. The combined effects of both fields are expressed in the electric displacement field <math>\mathbf D</math>, given by


Line 19: Line 21:
==== {{TAG|LEPSILON}}: density-functional-perturbation theory (DFPT) ====
==== {{TAG|LEPSILON}}: density-functional-perturbation theory (DFPT) ====


By setting {{TAG|LEPSILON}}=.True., VASP uses DFPT to compute the static ion-clamped dielectric matrix with or without local field effects ({{TAG|LRPA}}). 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 {{TAG|NBANDS}}.  
:By setting {{TAG|LEPSILON}}=.True., VASP uses DFPT to compute the static ion-clamped dielectric matrix with or without local field effects ({{TAG|LRPA}}). 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 {{TAG|NBANDS}}.  


At the end of the calculation, both the values of <math>\epsilon</math> including ({{TAG|LRPA}}=.True.) or excluding ({{TAG|LRPA}}=.False.) local-field effects are printed in the {{TAG|OUTCAR}} file. Perform a consistency check by comparing the values excluding local-field effects and static limit of <math>\epsilon</math> obtained with {{TAG|LOPTICS}}=.True., i.e., <math>\lim_{\omega\to0}\epsilon(\omega)</math>.  
:At the end of the calculation, both the values of <math>\epsilon</math> including ({{TAG|LRPA}}=.True.) or excluding ({{TAG|LRPA}}=.False.) local-field effects are printed in the {{TAG|OUTCAR}} file. Perform a consistency check by comparing the values excluding local-field effects and static limit of <math>\epsilon</math> obtained with {{TAG|LOPTICS}}=.True., i.e., <math>\lim_{\omega\to0}\epsilon(\omega)</math>.  


==== {{TAG|LCALCEPS}}: finite differences approach====
==== {{TAG|LCALCEPS}}: finite differences approach====
With {{TAG|LCALCEPS}}=.True., the dielectric tensor is computed from the derivative of the polarization, using
:With {{TAG|LCALCEPS}}=.True., the dielectric tensor is computed from the derivative of the polarization, using
:<math>
::<math>
\epsilon^\infty_{ij}=\delta_{ij}+
\epsilon^\infty_{ij}=\delta_{ij}+
\frac{4\pi}{\epsilon_0}\frac{\partial P_i}{\partial \mathcal{E}_j}
\frac{4\pi}{\epsilon_0}\frac{\partial P_i}{\partial \mathcal{E}_j}
Line 31: Line 33:
{i,j=x,y,z}.
{i,j=x,y,z}.
</math>
</math>
However, here the derivative is evaluated explicitly by employing finite differences. The direction and intensity of the perturbing electric field have to be specified in the {{FILE|INCAR}} file using the {{TAG|EFIELD_PEAD}} tag. As with DFPT, at the end of the calculation, VASP will write the dielectric tensor in the {{TAG|OUTCAR}} file. Control over the inclusion of local-field effects is done with the variable {{TAG|LRPA}}.
:However, here the derivative is evaluated explicitly by employing finite differences. The direction and intensity of the perturbing electric field have to be specified in the {{FILE|INCAR}} file using the {{TAG|EFIELD_PEAD}} tag. As with DFPT, at the end of the calculation, VASP will write the dielectric tensor in the {{TAG|OUTCAR}} file. Control over the inclusion of local-field effects is done with the variable {{TAG|LRPA}}.


===Dynamical response: Green-Kubo and many-body perturbation theory===
===Dynamical response: Green-Kubo and many-body perturbation theory===
===={{TAG|LOPTICS}}: Green-Kubo formula====
===={{TAG|LOPTICS}}: Green-Kubo formula====
{{TAG|LOPTICS}} allows for the evaluation of the frequency-dependent dielectric function once the ground state is computed. It uses the explicit expression to evaluate the imaginary part of <math>\epsilon</math>:  
:{{TAG|LOPTICS}} allows for the evaluation of the frequency-dependent dielectric function once the ground state is computed. It uses the explicit expression to evaluate the imaginary part of <math>\epsilon</math>:  
:<math>
::<math>
\epsilon^{(2)}_{\alpha \beta}\left(\omega\right) = \frac{4\pi^2 e^2}{\Omega}  
\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)
\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)
Line 42: Line 44:
             \langle u_{v\mathbf{k}} | u_{c\mathbf{k}+\mathbf{e}_\beta q} \rangle,
             \langle u_{v\mathbf{k}} | u_{c\mathbf{k}+\mathbf{e}_\beta q} \rangle,
</math>
</math>
while the real part is evaluated using the Kramers-Kroing relation. At this level, there are no effects coming from local fields.  
:while the real part is evaluated using the Kramers-Kroing relation. At this level, there are no effects coming from local fields.  


This method requires two steps: First, obtain the electronic ground state. Secondly, increase the value of {{TAG|NBANDS}} in the {{TAG|INCAR}} file to include unoccupied states. Always check for convergence w.r.t. the number of unoccupied states.  
:This method requires two steps: First, obtain the electronic ground state. Secondly, increase the value of {{TAG|NBANDS}} in the {{TAG|INCAR}} file to include unoccupied states. Always check for convergence w.r.t. the number of unoccupied states.  


Furthermore, the {{TAG|INCAR}} should also include values for {{TAG|CSHIFT}} (the broadening applied to the Lorentzian function which replaces the <math>\delta</math>-function), and {{TAG|NEDOS}} (the frequency grid for <math>\omega</math>).
:Furthermore, the {{TAG|INCAR}} should also include values for {{TAG|CSHIFT}} (the broadening applied to the Lorentzian function which replaces the <math>\delta</math>-function), and {{TAG|NEDOS}} (the frequency grid for <math>\omega</math>).


===={{TAG|ALGO}} = TDHF: Casida equation====
===={{TAG|ALGO}} = TDHF: Casida equation====
This option performs a [[Bethe-Salpeter-equations calculations#Time-dependent Hartree-Fock calculation|time-dependent Hartree-Fock]] or [[Bethe-Salpeter-equations calculations#Time-dependent DFT calculation|time-dependent density-functional-theory (TDDFT)]] calculation. It follows the Casida equation and uses a Fourier transform of the time-evolving dipoles to compute <math>\epsilon</math>.
:This option performs a [[Bethe-Salpeter-equations calculations#Time-dependent Hartree-Fock calculation|time-dependent Hartree-Fock]] or [[Bethe-Salpeter-equations calculations#Time-dependent DFT calculation|time-dependent density-functional-theory (TDDFT)]] calculation. It follows the Casida equation and uses a Fourier transform of the time-evolving dipoles to compute <math>\epsilon</math>.


The number of {{TAG|NBANDS}} controls how many bands are present in the time evolution. This can be fewer empty states compared to {{TAG|LOPTICS}}.  
:The number of {{TAG|NBANDS}} controls how many bands are present in the time evolution. This can be fewer empty states compared to {{TAG|LOPTICS}}.  


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


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


The time step is controlled automatically by the {{TAG|CSHIFT}} and {{TAG|PREC}}. This means that the smaller the value of {{TAG|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 the higher the cost of the calculation.
:The time step is controlled automatically by the {{TAG|CSHIFT}} and {{TAG|PREC}}. This means that the smaller the value of {{TAG|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 the higher the cost of the calculation.


The number of valence and conduction bands involved in the time propagation is set by the {{TAG|NBANDSO}} and {{TAG|NBANDSV}} tags, respectively. Choose a small number of bands near the band gap to reproduce optical measurements.  
:The number of valence and conduction bands involved in the time propagation is set by the {{TAG|NBANDSO}} and {{TAG|NBANDSV}} tags, respectively. Choose a small number of bands near the band gap 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 {{TAG|OMEGAMAX}}, and the sampling of the frequency grid is controlled by {{TAG|NEDOS}}.
:Finally, the maximum energy used in both the Fourier transform and in calculating the frequency-dependent dielectric function is set by {{TAG|OMEGAMAX}}, and the sampling of the frequency grid is controlled by {{TAG|NEDOS}}.


===={{TAG|ALGO}} = CHI: polarizability within RPA approximation====
===={{TAG|ALGO}} = CHI: polarizability within RPA approximation====
Here, the frequency dielectric function is computed within the Random-Phase approximation. VASP will compute the polarizability <math>\chi</math> by setting {{TAG|ALGO}}=Chi in the {{TAG|INCAR}} file and then use  
:Here, the frequency dielectric function is computed within the Random-Phase approximation. VASP will compute the polarizability <math>\chi</math> by setting {{TAG|ALGO}}=Chi in the {{TAG|INCAR}} file and then use  


::<math>
::<math>
Line 71: Line 73:
</math>
</math>


to compute the dielectric function. Here, <math>v</math> is the bare Coulomb potential describing electron-electron interaction. This requires increasing {{TAG|NBANDS}} to include unoccupied states as generally the case for [[GW_and_dielectric_matrix|GW calculations]]. Effects coming from local fields in the random-phase approximation (RPA) can be activated by setting {{TAG|LRPA}}=.True..  
:to compute the dielectric function. Here, <math>v</math> is the bare Coulomb potential describing electron-electron interaction. This requires increasing {{TAG|NBANDS}} to include unoccupied states as generally the case for [[GW_and_dielectric_matrix|GW calculations]].  


Two methods for computing the polarizability are available: For {{TAG|LSPECTRAL}}=.True., VASP will avoid direct computation of <math>\chi</math> and use a fast matrix-vector product. However, this can introduce spurious peaks at low frequencies for some values of {{TAG|CSHIFT}} and {{TAG|NOMEGA}}. The second method computes <math>\chi</math> directly and is activated by setting {{TAG|LSPECTRAL}}=.False.. However, it is much slower than the former method.
:Two methods for computing the polarizability are available: For {{TAG|LSPECTRAL}}=.True., VASP will avoid direct computation of <math>\chi</math> and use a fast matrix-vector product. However, this can introduce spurious peaks at low frequencies for some values of {{TAG|CSHIFT}} and {{TAG|NOMEGA}}. The second method computes <math>\chi</math> directly and is activated by setting {{TAG|LSPECTRAL}}=.False.. However, it is much slower than the former method.


===={{TAG|ALGO}} = BSE: macroscopic dielectric function including excitons====
===={{TAG|ALGO}} = BSE: macroscopic dielectric function including excitons====
Setting {{TAG|ALGO}}=BSE computes the macroscopic dielectric function <math>\epsilon_M</math> by solving the [[Bethe-Salpeter equations|Bethe-Salpeter equations]]. The electron-hole pairs are treated as a new quasi-particle called exciton, and the dielectric function is built using the eigenvectors (<math>X_{\lambda}^{cv\mathbf k}</math>) and eigenvalues (<math>\omega_\lambda</math>):  
:Setting {{TAG|ALGO}}=BSE computes the macroscopic dielectric function <math>\epsilon_M</math> by solving the [[Bethe-Salpeter equations|Bethe-Salpeter equations]]. The electron-hole pairs are treated as a new quasi-particle called exciton, and the dielectric function is built using the eigenvectors (<math>X_{\lambda}^{cv\mathbf k}</math>) and eigenvalues (<math>\omega_\lambda</math>):  


::<math>
::<math>
Line 86: Line 88:
</math>
</math>


Here, <math>S_{\lambda\lambda'}</math> is the overlap between exciton states of indices <math>\lambda</math> and <math>\lambda'</math> (in general, the BSE Hamiltonian is not hermitian, so eigenstates associated to different eigenvalues are not necessarily orthogonal).
:Here, <math>S_{\lambda\lambda'}</math> is the overlap between exciton states of indices <math>\lambda</math> and <math>\lambda'</math> (in general, the BSE Hamiltonian is not hermitian, so eigenstates associated to different eigenvalues are not necessarily orthogonal).


The number of occupied and unoccupied states that are included in the BSE Hamiltonian is controlled by the {{TAG|NBANDSO}} and {{TAG|NBANDSV}}, respectively. Normally only a few bands above and below the band gap are required to converge the optical spectrum, and the memory requirements increase quickly with the number of bands. Thus, be careful in setting up these two tags.  
:The number of occupied and unoccupied states that are included in the BSE Hamiltonian is controlled by the {{TAG|NBANDSO}} and {{TAG|NBANDSV}}, respectively. Normally only a few bands above and below the band gap are required to converge the optical spectrum, and the memory requirements increase quickly with the number of bands. Thus, be careful in setting up these two tags.  


Regarding the comparison with optical experiments, (e.g., absorption, MOKE, reflectance), <math>q</math> is the photon momentum. Often, the <math>\mathbf q \to 0</math> limit is considered. Furthermore, the coupling between the resonant and anti-resonant terms can be switched off, in what is called the [[Bethe-Salpeter equations#Theory#Tamm-Dancoff approximation|Tamm-Dancoff approximation]]. This approximation can be activated with the variable {{TAG|ANTIRES}} set to 0. Setting this variable to 1 or 2 will include the coupling, but increase the computational cost.
:Regarding the comparison with optical experiments, (e.g., absorption, MOKE, reflectance), <math>q</math> is the photon momentum. Often, the <math>\mathbf q \to 0</math> limit is considered. Furthermore, the coupling between the resonant and anti-resonant terms can be switched off, in what is called the [[Bethe-Salpeter equations#Theory#Tamm-Dancoff approximation|Tamm-Dancoff approximation]]. This approximation can be activated with the variable {{TAG|ANTIRES}} set to 0. Setting this variable to 1 or 2 will include the coupling, but increase the computational cost.


==Level of approximation==
==Level of approximation==
Line 111: Line 113:


::<math>
::<math>
\lim_{\mathbf q\to0}\frac{\langle c\mathbf k + \mathbf q|e^{\mathrm i\mathbf q\cdot\mathbf r}|v\mathbf k\rangle}{q} \approx \lim_{\mathbf q\to0}\frac{\langle c\mathbf k+\mathbf q|1 + \mathrm i\mathbf q\cdot\mathbf r|v\mathbf k\rangle}{q} = \hat\mathbf q\cdot \langle c\mathbf k+\mathbf q|\mathbf r|v\mathbf k\rangle.
\lim_{\mathbf q\to0}\frac{\langle c\mathbf k + \mathbf q|e^{\mathrm i\mathbf q\cdot\mathbf r}|v\mathbf k\rangle}{q} \approx \lim_{\mathbf q\to0}\frac{\langle c\mathbf k+\mathbf q|1 + \mathrm i\mathbf q\cdot\mathbf r|v\mathbf k\rangle}{q} = \hat{\mathbf{q}}\cdot \langle c\mathbf k+\mathbf q|\mathbf r|v\mathbf k\rangle.
</math>
</math>


VASP can also analyze the effects of finite momentum excitons. This is important, e.g., in the case of the optical absorption of bulk hexagonal BN. To calculate the absorption spectrum at finite momentum, set {{TAG|KPOINT_BSE}} to the index of the desired '''q''' point.  
VASP can also analyze the effects of finite momentum excitons. This is important, e.g., in the case of the optical absorption of bulk hexagonal BN. To calculate the absorption spectrum at finite momentum, set {{TAG|KPOINT_BSE}} to the index of the desired '''q''' point.


===Local fields in the Hamiltonian===
===Local fields in the Hamiltonian===
Line 160: Line 162:


==Other dielectric properties==
==Other dielectric properties==
===Polarizability===
===Electron energy loss spectroscopy (EELS)===
In EELS experiments a narrow beam of electrons with a well defined energy is shot at the sample. These electrons then lose energy to the sample by exciting plasmons, electron-hole pairs, or other higher-order quasiparticles. The loss function can then be expressed as
 
::<math>
\mathrm{EELS} = -\mathrm{Im}\left[\epsilon^{-1}(\omega)\right].
</math>
 
===Optical conductivity===
===Optical conductivity===
From Maxwell's equations and the microscopic form of Ohm's law it is possible to arrive at the following relation between the tensorial dielectric function and the optical conductivity <math>\sigma(\omega)</math>
From Maxwell's equations and the microscopic form of Ohm's law it is possible to arrive at the following relation between the tensorial dielectric function and the optical conductivity <math>\sigma(\omega)</math>


::<math>
::<math>
\epsilon_{\alpha\beta}(\omega) = \delta_{\alpha\beta} + \mathrm i \frac{4\pi}{\omega}\sigma_{\alpha\beta}(\omega).
\sigma_{\alpha\beta}(\omega) = \mathrm i\frac{\omega}{4\pi}\left[\delta_{\alpha\beta} - \epsilon_{\alpha\beta}(\omega)\right].
</math>
</math>


Line 199: Line 207:


==Electric response combined with perturbations of the ionic degrees of freedom==
==Electric response combined with perturbations of the ionic degrees of freedom==
<!---
{{TAG|IBRION}} = 5,6 (DFPT) and 7,8 (finite differences)


===Atomic displacements===
===Low-frequency corrections from atomic displacements===
The corrections from ionic motion to the low frequency regime can be added to <math>\epsilon_{\alpha\beta}^\infty</math> following{{cite|gonze:prb:1997}}
::<math>
\epsilon_{\alpha\beta}(\omega) = \epsilon_{\alpha\beta}^\infty + \frac{4\pi e^2}{\Omega_0}\sum_\nu\frac{S_{\alpha\beta,\nu}}{\omega_\nu^2 - (\omega+\mathrm i\eta)^2},
</math>
where <math>\omega_\nu</math> is the phonon frequency of mode <math>\nu</math>, and <math>S_{\alpha\beta,\nu}</math> is the mode-oscillator strength, defined by
::<math>
S_{\alpha\beta,\nu} = \left(\sum_{I,\delta}Z^*_{I\alpha\delta}\varepsilon^*_{I\delta,\nu}(\mathbf{q = 0})\right)\left(\sum_{J,\delta'}Z^*_{J\beta\delta'}\varepsilon_{J\delta',\nu}(\mathbf{q = 0})\right).
</math>
Here <math>Z^*_{J\beta\delta'}</math> are the Born effective charges and <math>\varepsilon_{J\delta',\nu}(\mathbf{q = 0})</math> are the eigendisplacements associated with the vibration mode <math>\nu</math> for atom <math>J</math> along direction <math>\delta'</math>. More information on the theory and methods behind the computation of phonon frequencies and eigendisplacements can be found in the [[Phonons:_Theory|phonons]] dedicated page.
 
Inclusion of the low-frequency corrections can be activated in the {{TAG|INCAR}} file with either {{TAG|IBRION}} = 5,6 (DFPT) or 7,8 (finite differences), and by setting to .True. the variables {{TAG|LEPSILON}} or {{TAG|LCALCEPS}}.


Keywords:
====Polar materials====
Born dynamical effective charge tensor aka Born effective charges.  
For polar materials it is important to recall that there is a discontinuity near <math>\Gamma</math>, i.e. <math>\omega^2_\nu(\mathbf{q\to 0}) \neq \omega^2_\nu(\mathbf{q= 0})</math>, and in fact, for a given unitary directional vector <math>\mathbf q</math>, it can be shown that the Lyddane-Sachs-Teller relationship holds{{cite|gonze:prb:1997}}
Polar materials
::<math>
which experiment measures this?
\prod_\nu\frac{\omega^2_\nu(\mathbf{q \to 0})-\omega^2}{\omega^2_\nu(\mathbf{q= 0}) - \omega^2} = \frac{\sum_{\alpha\beta}q_\alpha\epsilon_{\alpha\beta}(\omega)q_\beta}{\sum_{\alpha\beta}q_\alpha\epsilon_{\alpha\beta}^\infty q_\beta},
--->
</math>
meaning that the splitting in frequencies between the LO and TO modes at zero momentum carries over the evaluation of the dielectric function.


===Strain===
In order to obtain smooth phonon dispersions and to properly account for the LO-TO splitting in the evaluation of the optical limit of the dielectric function, users should read the dedicated page on [[Phonons:_Theory#Long-range_interatomic_force_constants_(LO-TO_splitting)|LO-TO splitting]], where it is explained how to set the variables {{TAG|LPHON_POLAR}}, {{TAG|PHON_DIELECTRIC}}, and {{TAG| PHON_BORN_CHARGES}} in the {{TAG|INCAR}} file.


piezoelectric tensor
===Corrections from strain===


The dielectric tensor can also be included in the evaluation of the elastic tensor, <math>C_{jk}</math>, (see [[Static_linear_response:_theory|theory of static linear response]] for more information on derived quantities from the static linear response). While this quantity is normally evaluated at fixed <math>\mathbf E</math>-field, in cases where a thin film is placed between layers of insulating materials, it is more convenient to evaluate the elastic tensor for fixed displacement field <math>\mathbf D</math>, since the boundary conditions fix the components of this vector in the direction normal to the surface.  
The dielectric tensor can also be included in the evaluation of the elastic tensor, <math>C_{jk}</math>, (see [[Static_linear_response:_theory|theory of static linear response]] for more information on derived quantities from the static linear response). While this quantity is normally evaluated at fixed <math>\mathbf E</math>-field, in cases where a thin film is placed between layers of insulating materials, it is more convenient to evaluate the elastic tensor for fixed displacement field <math>\mathbf D</math>, since the boundary conditions fix the components of this vector in the direction normal to the surface.  
Line 225: Line 243:
== References ==
== References ==


[[Category:VASP|Dielectric Functions]] [[Category:Linear response|Linear response]]
[[Category:VASP|Dielectric Functions]] [[Category:Linear response]]

Latest revision as of 13:23, 19 December 2024

Dielectric function

When light excites an electron from the valence to the conduction band, the created hole can bind with the electron to form an exciton.

When an external electric field acts on a medium, both the electronic and ionic charges will react to the perturbing field. For dielectric materials, in a very simplistic approach, one can think that the bound charges will create dipoles inside the medium leading to an induced polarization, . The combined effects of both fields are expressed in the electric displacement field , given by

.

If the external field is not strong enough to significantly change the properties of the dielectric medium, one can treat the induced polarization within the so-called linear response regime. Here, the information on how the dielectric reacts on the external field is given by the dielectric function:

,

which (assuming that the system has time-reversal symmetry) leads to

.

Depending on the nature of the external field, there are different approaches for the calculation of . If is static, then one can rely on perturbative methods based on finite differences or density functional perturbation theory (DFPT). However, if the perturbation is a time-dependent -field, (e.g., in measurements of the optical absorption, reflectance, magneto-optical Kerr effect (MOKE), etc.), the response will depend on the frequency of the external field. For these cases, one must employ methods based on time-dependent linear response, (e.g., Green-Kubo) or many-body perturbation theory.

Below we present an overview of all possible cases where VASP employs either one of such methods for the calculation of .

Static response

LEPSILON: density-functional-perturbation theory (DFPT)

By setting LEPSILON=.True., VASP uses DFPT to compute the static ion-clamped dielectric matrix with or without local field effects (LRPA). 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 NBANDS.
At the end of the calculation, both the values of including (LRPA=.True.) or excluding (LRPA=.False.) local-field effects are printed in the OUTCAR file. Perform a consistency check by comparing the values excluding local-field effects and static limit of obtained with LOPTICS=.True., i.e., .

LCALCEPS: finite differences approach

With LCALCEPS=.True., the dielectric tensor is computed from the derivative of the polarization, using
However, here the derivative is evaluated explicitly by employing finite differences. The direction and intensity of the perturbing electric field have to be specified in the INCAR file using the EFIELD_PEAD tag. As with DFPT, 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: Green-Kubo and many-body perturbation theory

LOPTICS: Green-Kubo formula

LOPTICS allows for the evaluation of the frequency-dependent dielectric function once the ground state is computed. It uses the explicit expression to evaluate the imaginary part of :
while the real part is evaluated using the Kramers-Kroing relation. At this level, there are no effects coming from local fields.
This method requires two steps: First, obtain the electronic ground state. Secondly, increase the value of NBANDS in the INCAR file to include unoccupied states. Always check for convergence w.r.t. the number of unoccupied states.
Furthermore, the INCAR should also include values for CSHIFT (the broadening applied to the Lorentzian function which replaces the -function), and NEDOS (the frequency grid for ).

ALGO = TDHF: Casida equation

This option performs a time-dependent Hartree-Fock or time-dependent density-functional-theory (TDDFT) calculation. It follows the Casida equation and uses a Fourier transform of the time-evolving dipoles to compute .
The number of NBANDS controls how many bands are present in the time evolution. This can be fewer empty states compared to LOPTICS.
The choice of time-dependent kernel is controlled by AEXX, HFSCREEN, and LFXC tags. 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: delta-pulse electric field

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. 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 the higher the cost of the calculation.
The number of valence and conduction bands involved in the time propagation is set by the NBANDSO and NBANDSV tags, respectively. Choose a small number of bands near the band gap 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: polarizability within RPA approximation

Here, the frequency dielectric function is computed within the Random-Phase approximation. VASP will compute the polarizability by setting ALGO=Chi in the INCAR file and then use
to compute the dielectric function. Here, is the bare Coulomb potential describing electron-electron interaction. This requires increasing NBANDS to include unoccupied states as generally the case for GW calculations.
Two methods for computing the polarizability are available: For LSPECTRAL=.True., VASP will avoid direct computation of and use a fast matrix-vector product. However, this can introduce spurious peaks at low frequencies for some values of CSHIFT and NOMEGA. The second method computes directly and is activated by setting LSPECTRAL=.False.. However, it is much slower than the former method.

ALGO = BSE: macroscopic dielectric function including excitons

Setting ALGO=BSE computes the macroscopic dielectric function by solving the Bethe-Salpeter equations. The electron-hole pairs are treated as a new quasi-particle called exciton, and the dielectric function is built using the eigenvectors () and eigenvalues ():
Here, is the overlap between exciton states of indices and (in general, the BSE Hamiltonian is not hermitian, so eigenstates associated to different eigenvalues are not necessarily orthogonal).
The number of occupied and unoccupied states that are included in the BSE Hamiltonian is controlled by the NBANDSO and NBANDSV, respectively. Normally only a few bands above and below the band gap are required to converge the optical spectrum, and the memory requirements increase quickly with the number of bands. Thus, be careful in setting up these two tags.
Regarding the comparison with optical experiments, (e.g., absorption, MOKE, reflectance), is the photon momentum. Often, the limit is considered. 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

Microscopic and macroscopic quantities

It is important to distinguish between macroscopic quantities, measured over several repetitions of the unit cell, and microscopic quantities, which include fields that change rapidly in all regions of the unit cell.

When measuring a property experimentally, it is the macroscopic version that will be represented in the experimental data. On the other hand, computationally, the microscopic quantities are more accessible. In other words, in order to compare experimental and computational results, the microscopic quantities, e.g., the dielectric function, must be averaged over several repetitions of the unit cell. It is possible to show that the macroscopic dielectric function, is related to the microscopic one via

where is the inverse dielectric function at .

Note that this does not mean that ! The full matrix has to be inverted and it is the component at that is used to calculate .

Finite momentum dielectric function

In the optical limit, the momentum of the incoming photon, , is almost zero, since the wavelength of the electric field is several times larger than the dimensions of the unit cell. Since the Coulomb potential diverges at very small momenta, the optical limit of the dielectric function must be obtained by taking with the limit of instead of setting . For instance, in the case of the independent particle approximation of full BSE, one yields

VASP can also analyze the effects of finite momentum excitons. This is important, e.g., in the case of the optical absorption of bulk hexagonal BN. To calculate the absorption spectrum at finite momentum, set KPOINT_BSE to the index of the desired q point.

Local fields in the Hamiltonian

Local fields, i.e., terms with finite , can be turned on or off in the Coulomb potential when evaluating the polarizability. VASP will distinguish the results in the OUTCAR file with:

MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects) 
BORN EFFECTIVE CHARGES (including local field effects) 
PIEZOELECTRIC TENSOR (including local field effects) 

and

MACROSCOPIC STATIC DIELECTRIC TENSOR (excluding local field effects) 
BORN EFFECTIVE CHARGES (excluding local field effects) 
PIEZOELECTRIC TENSOR (excluding local field effects)

for both cases.

Another approximation can also be taken, where the contributions from the exchange-correlation kernel are neglected when evaluating the polarizability. This is equivalent to the so-called random-phase approximation (RPA) and can be activated by setting LRPA=.True. in the INCAR.

Ion-clamped vs relaxed-ion/dressed dielectric function

The dielectric function computed either in the static or dynamical response regimes does not consider the effects coming from changes of the ionic positions due to the incoming electric field. This can be corrected by computing the relaxed-ion (or dressed) dielectric function [1]

where is the volume of the unit cell, is the Born effective charge, and is the force constants matrix.

Density-density versus current-current response functions

The inclusion of an electromagnetic field in the Hamiltonian is subject to a gauge choice. For instance, a classical electric field can be described by either a scalar potential or a longitudinal vector potential in the incomplete Weyl gauge ( = 0). The former means that the perturbing potential couples to the electronic density, while the second, implies a vector potential couples to a current. The fundamental consequence is that one can define two different response functions: a density-density response function for the first, ; and a current-current response function, . This circumstance is a common source of error when comparing experimental and computational optical properties of periodic systems, as discussed by Sangalli et al.[2].

Infact, perturbations associated with longitudinal fields will be described by the density-density polarisability function , (e.g., laser fields taken in the classical limit), while transverse fields will be described by the current-current polarizability , which is in fact a 3x3 tensor, (e.g., required to obtain the MOKE). Fundamentally, the time-dependent density is associated only with the longitudinal part of the current via the continuity equation. This also links both response functions via

It guarantees that the dielectric functions obtained from either approach match at finite momentum and frequency: .

The current-current dielectric function is exact at both and at . Especially for metals, it will reproduce the proper behavior of the Drude tail at . However, is more prone to numerical instabilities.

Other dielectric properties

Electron energy loss spectroscopy (EELS)

In EELS experiments a narrow beam of electrons with a well defined energy is shot at the sample. These electrons then lose energy to the sample by exciting plasmons, electron-hole pairs, or other higher-order quasiparticles. The loss function can then be expressed as

Optical conductivity

From Maxwell's equations and the microscopic form of Ohm's law it is possible to arrive at the following relation between the tensorial dielectric function and the optical conductivity

Optical absorption

For an electromagnetic wave traveling through a medium, one can express the electric field as , and the effects of the medium in the wave propagation are contained inside the dispersion relation . Using Maxwell's equations, one can arrive at

If the magnetic permeability of the material is assumed to be equal to that of vacuum, the equation above implies that the refractive index can be written as . Since is complex, the exponential factor in will have a dampening factor, , which accounts for the absorption of electromagnetic energy by the medium. With this relation one can define the absorption coefficient, as

Reflectance

From the previous subsection, one can also define the reflectivity coefficient at normal incidence as

This equation can be generalized for any angle of incidence , resulting in the general form of Fresnel equations.

Magneto-optical Kerr effect (MOKE)

The incoming electromagnetic wave interacts with the finite magnetic moment of the material. Usually, the interaction is with the magnetization of the medium, but there are also antiferromagnetic systems that can observe a finite MOKE. Generally, the reflected wave will gain an extra complex phase with respect to the incident -field. For a surface or two-dimensional material, (e.g., hexagonal BN, MoS), this phase can be computed using the off-diagonal components of the current-current dielectric tensor:

Electric response combined with perturbations of the ionic degrees of freedom

Low-frequency corrections from atomic displacements

The corrections from ionic motion to the low frequency regime can be added to following[3]

where is the phonon frequency of mode , and is the mode-oscillator strength, defined by

Here are the Born effective charges and are the eigendisplacements associated with the vibration mode for atom along direction . More information on the theory and methods behind the computation of phonon frequencies and eigendisplacements can be found in the phonons dedicated page.

Inclusion of the low-frequency corrections can be activated in the INCAR file with either IBRION = 5,6 (DFPT) or 7,8 (finite differences), and by setting to .True. the variables LEPSILON or LCALCEPS.

Polar materials

For polar materials it is important to recall that there is a discontinuity near , i.e. , and in fact, for a given unitary directional vector , it can be shown that the Lyddane-Sachs-Teller relationship holds[3]

meaning that the splitting in frequencies between the LO and TO modes at zero momentum carries over the evaluation of the dielectric function.

In order to obtain smooth phonon dispersions and to properly account for the LO-TO splitting in the evaluation of the optical limit of the dielectric function, users should read the dedicated page on LO-TO splitting, where it is explained how to set the variables LPHON_POLAR, PHON_DIELECTRIC, and PHON_BORN_CHARGES in the INCAR file.

Corrections from strain

The dielectric tensor can also be included in the evaluation of the elastic tensor, , (see theory of static linear response for more information on derived quantities from the static linear response). While this quantity is normally evaluated at fixed -field, in cases where a thin film is placed between layers of insulating materials, it is more convenient to evaluate the elastic tensor for fixed displacement field , since the boundary conditions fix the components of this vector in the direction normal to the surface.

If is the elastic tensor defined at fixed -field and the elastic tensor defined at fixed -field, then they are related by

where is the ion-relaxed piezoelectric tensor.

References

Subcategories

This category has the following 2 subcategories, out of 2 total.