24小时热门版块排行榜

返回列表

放逐孤岛

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 1.2
帖子: 95
在线: 24.4小时
虫号: 5225683
注册: 2016-11-11
性别: GG
专业: 能源化工

[求助] 在LAMMPS中使用compute voronoi/atom command命令输出错误

我在LAMMPS中使用了compute voronoi 1 voronoi/atom命令，并且用dump 2 all custom 100 voronoi.txt id c_voronoi

命令进行输出，但是得到的结果却显示体积与面数全是0，请问这是什么情况？
附上我的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

回复此楼