Fe crystal 鉄の結晶(FSポテンシャル)
FSポテンシャルを用いて、bccのFe結晶を構成しました。a=0.287nmです。
FeのFinnis-Sinclair potentialのパラメータ:(他のFSポテンシャルとφ(r)の形が異なります)
Φ=V(r)-2AF'φ(r)、F'=∂F/∂ρ、F=Sqrt(ρ)、ρi=sum(φ(rij),j=1,..n, j!=i)
φ(r)=(r-d)^2+b*(r-d)^3/d, (r<d else 0),V(r)=(r-c)^2*(c0+c1*r+c2*r^2) (r<c else 0)
A=1.828905eV,d=3.569745A,b=1.8,c=3.4A,c0=1.2371147,c1=-0.3592185,c2=-0.0385607
(出典:北川他著、「初心者のための分子動力学法」養賢堂1997、表3.3より)
inserted by FC2 system