レナード-ジョーンズポテンシャル
Φ=4ε{(σ/r)^12-(σ/r)^6}
σ=0.20 nm
ε=100 k_B K

これに合わせるための、f.c.c.構造の
格子定数a_0は、
a_0=σ*(2)^(1/6)/(√2/2)=0.31748 nm (約)
になります。

温度はこのように 50 Kと低目がいいでしょう
（表面が多いのですぐ融解するので）。

(x長さ)=(y長さ)=(z長さ)の立方体形状、
単位胞が(2+2)×(2+2)×(2+2)=64個、
f.c.c.なので、原子数は(N=)64×4=256個になります。