BCC 金属中のH原子の拡散をやるMCコードについて作成前に文献調査中。簡単だと思ってたけどいろいろ複雑だ。まずH原子は小さいので格子の隙間に入る。こればBCCの場合 Octahedral と Tetrahedral の二種類がある。鉄などでは T に入りやすいが引張り応力かけて拡がった場合もしかしたら O にも入るかもしれん。T から T へ移動する場合の障壁や、T サイトでのエネルギーの応力依存性なんかも妥当な形(数値は適当、符号は慎重に)を決めないと。そんで拡散係数は液体中とほとんど同じなんだそうな。速すぎー。あとはHが隣接した場合のエネルギーとかも要る。たぶん反発だな。
目的
