24小时热门版块排行榜    

查看: 965  |  回复: 0

fw7668172

捐助贵宾 (初入文坛)

[求助] 运行错误H原子种类不对 ,lj势语句错误

angle_coeff        1 2.60185 104.52
bond_coeff        1 10.84 0.9572

#angle_coeff        2  0 109.47
#bond_coeff        2  0  0.142
pair_style        lj/cut/tip4p/long 3 2 1 1 0.125 9 8
pair_coeff        1 1 0.0 0.0
pair_coeff        1 4 0.00324   3.19
pair_coeff        4 4 0.0070575 3.16435
pair_coeff        2 2 0.0 0.0
pair_coeff        1 2 0.0 0.0
pair_coeff        2 3 0.0 0.0
pair_coeff        2 4 0.00324   3.19
pair_coeff        3 4 0.0 0.0
pair_coeff        3 3 0.0 0.0
pair_coeff        1 3 0.0 0.0
kspace_style        pppm/tip4p 1.0e-5
group        g_A type 1
group        g_X type 2
group        g_C union g_A g_X
group        g_0 type 4
group        g_H type 3
group        g_W union g_0 g_H

#variable        field equal 0.04*cos(2*PI*20e12*1e-12*time)
#minimize        1e-5 1e-8 10000 100000
#reset_timestep        0

fix        1 g_C setforce 0.0 0.0 0.0
fix        2 g_W nvt temp 300 300 0.1
fix        3 g_W recenter INIT INIT NULL
fix        4 g_W shake 1e-5 500 0 b 1 a 1
#fix        5 g_W efield 0.0 0.0 v_field

velocity        g_W create 300.0 790316 dist gaussian
run        0
velocity        g_W scale 300.0
velocity        g_C set 0.0 0.0 0.0

timestep        0.001
compute 1 g_W temp
thermo        100
thermo_style        custom step c_1  epair  
#thermo_modify temp 1

#dump        1 g_W custom 200 c-h.lammpstrj type x y z
dump        1 all custom 200 move.lammpstrj mass id type x y z q element

run        200000
write_data data.*

#write_restart        restart.stable_singleC_water

#fix        5 g_W efield 0.0 0.0 0.08
#compute        1 g_W chunk/atom bin/3d x lower 1.5 y lower 1.5 z 0.0 1 nchunk every ids nfreq bound y 48.1425 49.6425 units box
#fix        5 g_W ave/chunk 100 5 10000 1 density/number ave running file angle.profile
#dump        1 g_W custom 500 angle.lammpstrj type x y z
#dump        2 all custom 200 single_C.lammpstrj id mol type x y z q



运行错误ERROR on proc 0:TIP4P hydrogen has incorrect atom type
模型是水滴和石墨炔
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fw7668172 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见