| 查看: 2307 | 回复: 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 |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有4人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有3人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
|
我找到原因了，我在dump时输出的只是group为1的原子，而compute voronoi/atom 命令只有该组的原子输出值其他组为0，所以我看到的大多数值都为0 发自小木虫Android客户端 |
2楼2019-04-10 09:54:36
20