鉄＋水素 - ita’s diary

鉄＋水素の分子動力学コードを２日で書き上げた。日本の職場にコードあったけど、持ってきたデータに入ってなかった。せっかくなんで最新のポテンシャルをリサーチして2003年のを見つけた。しかしまだバグあり。おかしいな。
全エネルギー
$\large \sum_{ij} V(R_{ij}) + \sum_i F(\rho_i)$
$\large \rho_i = \sum_j \phi(R_{ij})$
力
$\large \sum_{ij} ( V'(R_{ij})+(F'(\rho_i)+F'(\rho_j))\phi'(R_{ij})))\nabla_i R_{ij}$
でいいんだよなぁー。追記：r_{t+1} = r_t + dt f_t で dt が大き過ぎただけだった。正常に動作中。
詳細は Phil.Mag. 83 (35) 3977(2003)。関数形は下手な式でなく三次スプラインのテーブルで提供されてるのが面白い。Fe-H potential は PRB 54, 9765.