| 查看: 2163 | 回复: 1 | ||
[求助]
在LAMMPS中使用compute voronoi/atom command命令输出错误
|
|
我在LAMMPS中使用了compute voronoi 1 voronoi/atom命令，并且用dump 2 all custom 100 voronoi.txt id c_voronoi 附上我的in文件： dimension 3 boundary p p p units real atom_style full read_data b8000.data neighbor 2 bin neigh_modify delay 0 every 1 pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 pair_coeff 1 1 0.1521 3.154000 8.5 pair_coeff * 2 0.000000 0.000000 8.5 bond_style harmonic bond_coeff 1 1000 0.9572 angle_style harmonic angle_coeff 1 3000 104.52 kspace_style pppm/tip4p 0.00001 thermo_style custom step temp density etotal press vol thermo 100 minimize 1.0e-12 1.0e-12 1000 1000 region 1 sphere 40 40 40 10 group 1 region 1 group 2 subtract all 1 fix 1 1 npt temp 298.15 298.15 $(100.0*dt) iso -2000.0 -2000.0 1000.0 fix 1 2 npt temp 298.15 298.15 $(100.0*dt) iso 1.0 1.0 1000.0 compute voronoi 1 voronoi/atom timestep 0.5 dump 1 all atom 100 dump.lammpstrj dump 2 all custom 100 voronoi.txt id c_voronoi run 500000 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
|
我找到原因了，我在dump时输出的只是group为1的原子，而compute voronoi/atom 命令只有该组的原子输出值其他组为0，所以我看到的大多数值都为0 发自小木虫Android客户端 |
2楼2019-04-10 09:54:36