24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1405  |  回复: 3

Alchemy4eM

新虫 (初入文坛)

[求助] 求助lammps定义group出现错误Illegal group command 已有1人参与

定义group的程序码:
variable        temp1 equal 298.15               
variable        temp2 equal 358.15                       
variable        r index 3                       
units                real
atom_style        full
boundary        p p p
pair_style        lj/cut/coul/long 10 12
bond_style        harmonic
angle_style        harmonic
dihedral_style        harmonic
improper_style        cvff
read_data        water-.data
kspace_style        pppm 1.0e-6
neighbor        2.0 bin
neigh_modify        delay 5
velocity        all create ${temp1} 123456
timestep        1
write_dump        all atom dump.setup_Water_$rnm_${temp2}K.lammpstrj

group                O type 1
group                H type 2
group                water union O H

fix                1 water nvt temp ${temp1} ${temp1} 100 drag 0.2

thermo                1000
restart                5000000 restartWater_$rnm_${temp2}K.*
dump                1 all atom 10000 dump_Water_$rnm_equil_${temp2}K.lammpstrj                               

run                100000
unfix                1
undump                1

fix                1 water npt temp ${temp2} ${temp2} 100 x 1.0 1.0 0.2 y 1.0 1.0 0.2 z 1.0 1.0 0.2 drag 0.2

dump                1 all atom 100000 dump_Water_$rnm_${temp2}K.lammpstrj                                                #1000000


variable        step1 equal 100000                               
variable        step2 equal 1000000                               
variable        current equal "step"                               
variable        step3 equal 10000                               

compute                neighbor O coord/atom 5.1 1
compute                neighborH H coord/atom 5.1 2

variable        vapor atom " c_neighbor < 4 && c_neighborH == 0"
group                vapor dynamic all var vapor every ${step1}
variable        liquid atom " c_neighbor >= 4 && c_neighborH == 0"
group                liquid dynamic all var liquid every ${step1}



恳请各位老师帮助!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

★ ★
Alchemy4eM: 金币+2 2016-05-04 11:13:10
本帖仅楼主可见
2楼2016-05-02 20:12:19
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)

本帖仅楼主可见
3楼2016-05-02 20:18:53
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

Alchemy4eM

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by philipchan at 2016-05-02 20:18:53
抱歉,没注意,搞错了。。。没发现哪边有错。。。也许”group                liquid dynamic all var liquid every ${step1}“ 中" all" 该去掉?或者group中有些ID和特指名词冲突了?...

谢谢了~是因为我的LAMMPS版本过低,改过版本和优化步骤之后就好了
4楼2016-05-04 11:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Alchemy4eM 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26年申博自荐-计算机视觉 20+3 光亮我自己 2026-05-01 4/200 2026-05-03 09:01 by lfy8008
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-03 07:22 by jmr1kibhvw
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 07:19 by jmr1kibhvw
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 06:52 by jmr1kibhvw
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-03 06:34 by jmr1kibhvw
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 d8tbypao1q 2026-05-01 4/200 2026-05-03 05:49 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-03 05:07 by jmr1kibhvw
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 04:34 by w5idbsc3mz
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 03:49 by w5idbsc3mz
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 02:36 by w5idbsc3mz
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k17zww5nna 2026-05-01 4/200 2026-05-03 02:06 by w5idbsc3mz
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k17zww5nna 2026-05-01 5/250 2026-05-03 02:04 by w5idbsc3mz
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 01:34 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 00:00 by 火炎焱燚炙烤
[考博] 考博 +4 Eternitf 2026-04-29 4/200 2026-05-02 16:10 by yejenny
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[文学芳草园] 太白金星有点烦 +3 yexuqing 2026-04-27 3/150 2026-04-30 21:56 by myrtle
[考博] 中国地质大学(北京)博士招生补录,数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
[论文投稿] 有没有快的中文核心比较快录用的,纳米材料光催化 5+3 2914439618 2026-04-27 4/200 2026-04-29 08:59 by 北京莱茵润色
信息提示
请填处理意见