一个适用于高温高压二氧化硅MD模拟的table

1. 势函数表达

短距离内用一个协调项代替BKS
At small distances, we replace the BKS form with a harmonic term to prevent unphysical particle fusion that we find can infrequently occur at extreme high temperatures

Wij是一个常数，来保证在rm处势能和力的连续，但是按照文献中计算出来的并不连续，下边代码中对Wij进行了修改。

Shen, Yuan, et al. “Nanosecond homogeneous nucleation and crystal growth in shock-compressed SiO 2.” Nature materials 15.1 (2016): 60-65.