Electric field response from density-functional-perturbation theory

Density-functional-perturbation theory provides a way to compute the second-order linear response to ionic displacement, strain, and electric fields. The equations are derived as follows.

In density-functional theory, we solve the Kohn-Sham (KS) equations

[math]\displaystyle{ H(\mathbf{k}) | \psi_{n\mathbf{k}} \rangle= e_{n\mathbf{k}}S(\mathbf{k}) | \psi_{n\mathbf{k}} \rangle, }[/math]

where [math]\displaystyle{ H(\mathbf{k}) }[/math] is the DFT Hamiltonian, [math]\displaystyle{ S(\mathbf{k}) }[/math] is the overlap operator and, [math]\displaystyle{ | \psi_{n\mathbf{k}} \rangle }[/math] and [math]\displaystyle{ e_{n\mathbf{k}} }[/math] are the KS eigenstates.