Angular potential

３体項は次のStillinger-Weber (SW) 型の関数．

$${\varphi }^{3D}({\mathbf{r}}_i,{\mathbf{r}}_j,{\mathbf{r}}_k)={\alpha }_{ijk}\exp \left[\left(\frac{{\beta }_{ijk}}{r_{ij}-r_{\mathrm{c}}^{3D}}\right)+\left(\frac{{\beta }_{ijk}}{r_{jk}-r_{\mathrm{c}}^{3D}}\right)\right]\cdot \left[(\cos {\theta }_{ijk}-{\gamma }_{ijk}{)}^2+{\eta }_{ijk}\right]$$

次のような in.params.angular ファイルを読み込む必要がある．

# type,   cspi, cspj, cspk,  rc3,   alp,   bet,   gmm,   eta
  1   Zr    O     O      3.50    2.133   1.000  -0.677
  1   P     O     O      3.50    2.282   1.000  -0.553
  2   O     P     P      3.00    1.000   1.000  -0.333   1.00

• # や ! で始まる行はコメント行．
• type — angular potential のタイプ（1 or 2）．
  • type-1 の場合： ${\eta }_{ijk}=0$
  • type-2の場合： ${\eta }_{ijk}\ne 0$
• csp[i|j|k] — ３体を構成する元素．i が中心元素．
• rc3 — カットオフ関数．
• alp, bet, gmm — $\alpha ,\beta ,\gamma$