Connection between DFPT and cDFPT [1] ω ˜ qμν 2 = ω qμν 2 δ μν + 2 ω qμ ω qν Π qμν Phonon self-energy Π qμν = 2 N k g qkμ * f ( ε k+q ) - f ( ε k ) ε k+q - ε k g ˜ qkν
Add phonon self-energy
μ=  eV
kT=  eV
f(ε)=

[1] Nomura and Arita, PRB 92, 245108 (2015)