24小时热门版块排行榜    

查看: 2198  |  回复: 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的回帖
    相关版块跳转 我要订阅楼主 放逐孤岛 的主题更新
    最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
    [考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
    [考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
    [考研] 0703化学调剂 ,六级已过,有科研经历 +14 曦熙兮 2026-03-15 14/700 2026-03-21 19:12 by ColorlessPI
    [考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
    [考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
    [考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
    [基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
    [考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
    [考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
    [考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
    [考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
    [考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
    [考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
    [考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
    [考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
    [考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
    [考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
    [考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
    [考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
    [考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
    信息提示
    请填处理意见