Electron-phonon interactions theory: Difference between revisions

From VASP Wiki
(Created page with "== Electron-phonon interactions from statistical sampling == The probability distribution of finding an atom within the coordinates <math>\kappa+d\kappa</math> (where <math>\...")
 
No edit summary
 
(9 intermediate revisions by 3 users not shown)
Line 1: Line 1:
== Electron-phonon interactions from statistical sampling ==
== Electron-phonon interactions from statistical sampling ==


The probability distribution of finding an atom within the coordinates <math>\kappa+d\kappa</math> (where <math>\kappa</math> denotes the Cartesian coordinates as well as the atom number) at temperature <math>T</math> in the harmonic approximation is given by the following expression
The probability distribution of finding an atom within the coordinates <math>\kappa+d\kappa</math> (where <math>\kappa</math> denotes the Cartesian coordinates as well as the atom number) at temperature <math>T</math> in the harmonic approximation is given by the following expression{{cite|bloch:zfp:1932}}{{cite|landau:ctp:1959}}


<math>
<math>
Line 15: Line 15:
Here <math>M_{\kappa}</math>, <math>\nu</math> and <math>\omega_{\nu}</math> denote the mass, phonon eigenmode and phonon eigenfrequency, respectively.
Here <math>M_{\kappa}</math>, <math>\nu</math> and <math>\omega_{\nu}</math> denote the mass, phonon eigenmode and phonon eigenfrequency, respectively.
The equation for <math>dW</math> is valid at any temperature and the high (Maxwell--Boltzmann distribution) and low temperature limits are easily regained.
The equation for <math>dW</math> is valid at any temperature and the high (Maxwell--Boltzmann distribution) and low temperature limits are easily regained.
In order to obtain an observable <math>O(T)</math> at a given temperature <math>T</math, the average of the observable sampled at different coordinate sets <math>x_{T}^{\textrm{MC},i}</math> with sample size <math>n</math> is taken
In order to obtain an observable <math>O(T)</math> at a given temperature <math>T</math>, the average of the observable sampled at different coordinate sets <math>x_{T}^{\textrm{MC},i}</math> with sample size <math>n</math> is taken


<math>
<math>
\langle O(T)\rangle = \frac{1}{n} \sum\limits_{i=1}^{n} O(x_{T}^{\textrm{MC,i}}).
\langle O(T)\rangle = \frac{1}{n} \sum\limits_{i=1}^{n} O(x_{T}^{\textrm{MC,i}}).
</math>
</math>
=== Full Monte-Carlo sampling ===


Each set <math>i</math> is obtained from the equilibrium atomic positions <math>x_{\textrm{eq}}</math> as
Each set <math>i</math> is obtained from the equilibrium atomic positions <math>x_{\textrm{eq}}</math> as
Line 34: Line 37:


Here <math>\varepsilon_{\kappa,\nu}</math> denotes the unit vector of eigenmode <math>\nu</math> on atom <math>\kappa</math>. The magnitude of the displacement in each Cartesian direction is obtained from the normal-distributed random variable <math>\mathcal{N}</math> with a probability distribution according to <math>dW</math>.
Here <math>\varepsilon_{\kappa,\nu}</math> denotes the unit vector of eigenmode <math>\nu</math> on atom <math>\kappa</math>. The magnitude of the displacement in each Cartesian direction is obtained from the normal-distributed random variable <math>\mathcal{N}</math> with a probability distribution according to <math>dW</math>.
=== ZG configuration (one-shot method) ===
Motivated by the empirical observation that for increasing super-cell sizes the number of required structures in the MC method can be decreased, M. Zacharias and F. Giustino{{cite|zacharias:prb:2016}} proposed a one-shot method where only a single set of displacements is used
<math>
\Delta \tau^{\textrm{OS}} = \sqrt{\frac{1}{M_{\kappa}}} \sum\limits_{\nu}^{3(N-1)} (-1)^{\nu-1} \varepsilon_{\kappa,\nu} \sigma_{\nu,T},
</math>
where the summation over the eigenmodes runs in an ascending order with respect to the values of the eigenfrequencies, and the magnitude of each displacement is given by
<math>
\sigma_{\nu,T} = \sqrt{(2 n_{\nu,T}+1) \frac{\hbar}{2 \omega_{\nu}}}.
</math>
Here <math>n_{\nu,T}=[\mathrm{exp}(\hbar \omega_{\nu} /k_{B}T)-1]^{-1}</math> denotes the Bose-Einstein occupation number.
In this way the sum for the observable <math>\langle O(T)\rangle</math>  is reduced to a single calculation. In Ref. {{cite|zacharias:prb:2016}} it was shown that for super-cell sizes <math>N\rightarrow \infty</math> the structural configuration obtained using the ZG configuration should lead to equivalent results as fully converged MC calculations. In practice, it was shown that already relatively small super-cell sizes are sufficient to achieve good accuracy, but the convergence with respect to the cell size can vary between different materials. In Ref. {{cite|karsai:njp:2018}} we have also used a slightly modified approach, in which the signs of the displacements are chosen randomly instead of <math>\pm 1</math>. This was only necessary, when calculating volume dependent ZPR, since the modes sometimes change the order as the volume changes. Using alternating signs for the displacement then causes small discontinuities in the ZPR volume curve of the order of 5 meV for carbon diamond. By averaging over many random phases this problem can be eliminated. Nevertheless 5 meV difference between the ZG configuration and full MC is considered as accurate enough in most calculations.
== References ==
<references/>


----
----
[[Category:Lattice Vibrations]][[Category:Electron-phonon interactions]][[Category:Theory]][[Category:VASP6]]
[[Category:Phonons]][[Category:Electron-phonon interactions]][[Category:Theory]]

Latest revision as of 14:44, 20 February 2024

Electron-phonon interactions from statistical sampling

The probability distribution of finding an atom within the coordinates (where denotes the Cartesian coordinates as well as the atom number) at temperature in the harmonic approximation is given by the following expression[1][2]

where the mean-square displacement of the harmonic oscillator is given as

Here , and denote the mass, phonon eigenmode and phonon eigenfrequency, respectively. The equation for is valid at any temperature and the high (Maxwell--Boltzmann distribution) and low temperature limits are easily regained. In order to obtain an observable at a given temperature , the average of the observable sampled at different coordinate sets with sample size is taken


Full Monte-Carlo sampling

Each set is obtained from the equilibrium atomic positions as

with the displacement

Here denotes the unit vector of eigenmode on atom . The magnitude of the displacement in each Cartesian direction is obtained from the normal-distributed random variable with a probability distribution according to .


ZG configuration (one-shot method)

Motivated by the empirical observation that for increasing super-cell sizes the number of required structures in the MC method can be decreased, M. Zacharias and F. Giustino[3] proposed a one-shot method where only a single set of displacements is used

where the summation over the eigenmodes runs in an ascending order with respect to the values of the eigenfrequencies, and the magnitude of each displacement is given by

Here denotes the Bose-Einstein occupation number. In this way the sum for the observable is reduced to a single calculation. In Ref. [3] it was shown that for super-cell sizes the structural configuration obtained using the ZG configuration should lead to equivalent results as fully converged MC calculations. In practice, it was shown that already relatively small super-cell sizes are sufficient to achieve good accuracy, but the convergence with respect to the cell size can vary between different materials. In Ref. [4] we have also used a slightly modified approach, in which the signs of the displacements are chosen randomly instead of . This was only necessary, when calculating volume dependent ZPR, since the modes sometimes change the order as the volume changes. Using alternating signs for the displacement then causes small discontinuities in the ZPR volume curve of the order of 5 meV for carbon diamond. By averaging over many random phases this problem can be eliminated. Nevertheless 5 meV difference between the ZG configuration and full MC is considered as accurate enough in most calculations.

References