有名なPseudopotentials that workの論文のパラメタ表から引っ張ってきて 擬ポテンシャル＋OPW法でSiとLiのエネルギーバンドを計算してみました。 ↑Liのバンド図 ↑Siのバンド図Siのバンドギャップがないのが気になりますが、大体形はあってそうです。 計算で…

一般化固有値問題をLAPACKを使って解いているのですが、行列が非対称になったときに たまに結果が収束しなかったりして困ってました。僕のLAPACKの使い方が駄目なだけかもしれませんが、 倍精度doubleでは駄目かもしれない。ちなみに、単精度(float)は全然話…

結晶のバンド計算をするとき、繰り返し計算の度に、全ｋ空間で状態密度を積分しないと いけなかったりするのですが、全k点を計算するのは大変なので、代表点の平均をとったり（特殊点法）、 金属なんかではテトラヘドロン法（四面体法）を使ったりして、近似…