BCC鉄中に水素が入った場合の緩和を計算するのに、BFGS法てのが使えると先週知る。VASPで使われてる。Wikipedia や Numerical Recipie を見て速攻コーディングした。動くことは動くがうまくいかん。力がまだ働いてるのにそっちへ動かない。むむむ。わからん。もう Brain Fucking Goddamn Sonofabitch 法という感じだ。
追記：最適と判断された移動方向と、力の方向の内積を計算し、これが負の場合はHessianのアップデートをしないという改造で解決した。
そういえばBFG9000てのもあったな。映画DooMにも出てきた。
