24小时热门版块排行榜    

查看: 2346  |  回复: 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文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:42 by XiYP04It3Evb
    [考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
    [考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:20 by XiYP04It3Evb
    [博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:07 by XiYP04It3Evb
    [考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:28 by XiYP04It3Evb
    [硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
    [找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
    [找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
    [博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:49 by XiYP04It3Evb
    [找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:36 by XiYP04It3Evb
    [考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:32 by XiYP04It3Evb
    [教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
    [公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GVGRDgNfv7P4 2026-10-23 9/450 2026-10-25 12:10 by XiYP04It3Evb
    [博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
    [考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
    [公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-24 23:35 by XcOYG53d3x4p
    [硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-24 22:56 by XcOYG53d3x4p
    [公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
    [教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
    [文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
    信息提示
    请填处理意见