Electron-phonon coupling

At the heart of the electron-phonon coupling is the calculation of the gradient of the effective potential, which is done using finite displacements just like the phonons. Those two calculations can run simultaneous, if the required set of parameters coincide. The electron-phonon matrix is quite sensitive to self-interaction of a displaced atom with its periodic images, so a sufficiently large supercell needs to be used. The (3x3x3) supercell used in this example is a bit too small. When using a supercell for the calculation you have to consider, that the atoms object needs to contain the primitive cell, not the supercell, while the parameters for the calculator object need to be good for the supercell, not the primitive cell. (effective_potential.py)

from ase.build import bulk
from gpaw import GPAW, FermiDirac
from gpaw.elph import DisplacementRunner

atoms = bulk('Si', 'diamond', a=5.431)

calc = GPAW(mode='lcao', h=0.18, basis='dzp',
            kpts=(4, 4, 4),
            xc='PBE',
            occupations=FermiDirac(0.01),
            symmetry={'point_group': False},
            convergence={'energy': 2e-5, 'density': 1e-5},
            txt='displacement.txt')

elph = DisplacementRunner(atoms=atoms, calc=calc,
                          supercell=(3, 3, 3), name='elph',
                          calculate_forces=True)
elph.run()

This script executes \(2*3*N\) displacements and saves the change in total energy and effective potential into a file cache in the directory \(elph\). The phonon/effective potential calculation can take quite some time, but can be distributed over several images. The DisplacementRunner() class is based on ASEs ase.phonon.Displacement class, which allows to select atoms to be displaced using the set_atoms function.

In the second step we map the gradient of the effective potential to the LCAO orbitals of the supercell (supercell.py). For this we first need to create the supercell

atoms = bulk('Si', 'diamond', a=5.431)
atoms_N = atoms * (3, 3, 3)

and run GPAW to obtain the wave function. This step currently currently only works for k-point parallelization, so you need to include parallel={'domain': 1, 'band': 1} in your script.

sc = Supercell(atoms, supercell=(3, 3, 3))
sc.calculate_supercell_matrix(calc, fd_name='elph')

invokes the calculation of the supercell matrix, which is stored in the \(supercell\) file cache.

Note

The real space grids used in the finite displacement calculations needs to be the same as the one used in the supercell matrix calculation. If you use planewave mode for the finite displacement calculation you should set the required grid manually, for example by adding gpts=(nx, ny, nz) where nx, ny, nz need be substituted with the required number of grid points in each direction. You can use python3 -m gpaw.elph.gpts to get help with this.

After both calculations are finished the final electron-phonon matrix can be constructed. (gmatrix.py)

from gpaw import GPAW
from gpaw.elph import ElectronPhononMatrix

calc = GPAW("scf.gpw")
atoms = calc.atoms

elph = ElectronPhononMatrix(atoms, 'supercell', 'elph')
q = [[0., 0., 0.], [1. / 11., 1. / 11., 1. / 11.]]
g_sqklnn = elph.bloch_matrix(calc, k_qc=q,
                             savetofile=True, prefactor=False)

For this we need to perform another ground state calculation (scf.py) to obtain the wave function used to project the electron-phonon coupling matrix into. The electron-phonon matrix is computed for a list of q values, which need to be commensurate with the k-point mesh chosen. The ElectronPhononMatrix class doesn’t like parallelization so much and should be done separately from the rest.

Exercise

The electron-phonon matrix can be used to calculate acoustic and optical deformation potentials without much effort. The optical deformation potential is identical to the bare electron-phonon coupling matrix Ref. [1]:

\[D_{ODP} = \langle m \vert \nabla_u V_{eff} \cdot \mathbf e_l \vert n \rangle .\]

For the silicon VBM at Gamma, for the LO phonons at Gamma, we expect a deformation potential of about \(\vert M \vert \approx 3.6\) eV/Angstrom.

Converge the deformation potential of Si with respect to the supercell size, k-points, grid spacing and SCF convergence parameters.

References