格子からは、フォノンを計算して出す。
電子のヤツは、線形応答で出す。
CASTEPなら、こんな感じでは？