分享一个2维液态平衡模拟案例代码。
随机生成400*0.8个原子并进行温度初始化,原子在LJ力场作用下运动到一定位置后达到平衡,并输出rdf数据。
代码已经注释,以供参考。
当原子数量超过400*0.9时,体系接近固态,可自行更改参数并结合rdf图像进行验证。
in文件代码如下:
#设置原子个数
variable afrac equal 0.80
variable npart equal ${
afrac}*400
<
分享一个2维液态平衡模拟案例代码。
随机生成400*0.8个原子并进行温度初始化,原子在LJ力场作用下运动到一定位置后达到平衡,并输出rdf数据。
代码已经注释,以供参考。
当原子数量超过400*0.9时,体系接近固态,可自行更改参数并结合rdf图像进行验证。
in文件代码如下:
#设置原子个数
variable afrac equal 0.80
variable npart equal ${
afrac}*400
<