固体の比誘電率の計算.
液体の場合添付する式のように(周期境界条件下の)液体のシミュレーションで比誘電率を計算する論文を見かけます.が,結晶でこれに代入するとほぼepsilon_infになってしまいおかしい気がします.なぜなら〈M^2〉-〈M〉^2 は結晶だときれいにそろっているから0であるはず.〈M〉は系全体の電気双極子モーメントの和.でも実際に測定すると強誘電体の固体とかあるわけですよね.
固体の場合はどうやって分子論的情報をもとに比誘電率を計算できますか?
また,添付した式の適用範囲はどこまででしょうか?
また,理解の浅そうなところを教えて下さい