24小时热门版块排行榜    

查看: 2365  |  回复: 1

放逐孤岛

铜虫 (小有名气)

[求助] 在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
  • 回复此楼
    已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

    放逐孤岛

    铜虫 (小有名气)

    我找到原因了,我在dump时输出的只是group为1的原子,而compute voronoi/atom 命令只有该组的原子输出值其他组为0,所以我看到的大多数值都为0

    发自小木虫Android客户端
    2楼2019-04-10 09:54:36
    已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
    相关版块跳转 我要订阅楼主 放逐孤岛 的主题更新
    最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
    [论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:10 by xSrQNSSjcgoX
    [硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
    [论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 uKOSCocCb6Fm 2026-11-12 10/500 2026-11-15 02:15 by xSrQNSSjcgoX
    [硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 10/500 2026-11-15 00:39 by Equinoxhua
    [考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 7/350 2026-11-14 22:44 by Equinoxhua
    [考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-14 22:39 by Equinoxhua
    [考博] 考博求助 10+3 社会实践对别 2026-11-08 3/150 2026-11-13 14:04 by 沐小虫252
    [公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-13 11:37 by Equinoxhua
    [论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
    [考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
    [考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
    [考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 11:22 by HajghKOMynmz
    [教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:14 by HajghKOMynmz
    [考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
    [考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:49 by HajghKOMynmz
    [硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
    [硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
    [公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 10:03 by HajghKOMynmz
    [公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
    [考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
    信息提示
    请填处理意见