24小时热门版块排行榜    

查看: 4442  |  回复: 5
【悬赏金币】回答本帖问题,作者DDDDDTTT将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

DDDDDTTT

新虫 (初入文坛)

[求助] LAMMPS中模拟电场问题已有1人参与

lammps小白想给体系加电场,参考lammps官网上的 commands 后,决定用 fix efield command 命令

fix ID group-ID efield ex ey ez keyword value ...

页面最下面提示 :

This fix is part of the MISC package. It is only enabled if LAMMPS was built with that package. See the Making LAMMPS section for more info.

然后,我 make yes-MISC 命令行显示安装成功,然后用 make package-statue 查看,也是MISC 也是安装完成的。

但是,当我在In文件中用 fix 0 all efield 1.0 0.0 0.0 时显示错误 : ERROR : Unknown fix style (../modify.cpp:793)

我怀疑是不是包安装错了,我又 make yes-all 再试了一次,但还是 : ERROR : Unknown fix style (../modify.cpp:793)

这是为什么呀?求大神帮帮我.谢谢!


附代码:

# VARIABLES
variable fname index PE_cl100.txt
variable simname index PE_cl100

# Initialization
units                real
boundary        p p p
atom_style        molecular
log                 log.PE_cl100.txt
read_data        PE_cl100.txt

# Dreiding potential information
neighbor        0.4 bin
neigh_modify        every 10 one 10000
bond_style      harmonic
bond_coeff        1 350 1.53
angle_style     harmonic
angle_coeff        1 60 109.5
dihedral_style        multi/harmonic
dihedral_coeff        1 1.73 -4.49 0.776 6.99 0.0
pair_style        lj/cut 10.5
pair_coeff        1 1 0.112 4.01 10.5

compute csym all centro/atom fcc
compute peratom all pe/atom



#####################################################
# Equilibration (Langevin dynamics at 5000 K)

velocity         all create 5000.0 1231

fix           0 all efield 1.0 0.0 1.0
fix                1 all nve/limit 0.05
fix                2 all langevin 5000.0 5000.0 10.0 904297
thermo_style        custom step temp
thermo          10000
timestep        1
run                100000
unfix 1
unfix 2
write_restart         restart.PE_cl100.dreiding1


#####################################################
# Define Settings
compute eng all pe/atom
compute eatoms all reduce sum c_eng

#####################################################
# Minimization

dump                 1 all cfg 6 dump.comp_*.cfg mass type xs ys zs c_csym c_peratom fx fy fz

reset_timestep 0
fix 1 all nvt temp 500.0 500.0 100.0

thermo 20
thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms
min_style cg
minimize 1e-25 1e-25 500000 1000000


print "All done"
回复此楼

» 收录本帖的淘贴专辑推荐

分子动力学-初级问题资源

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 DDDDDTTT 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我标书代码变了 +60 学员NHuqdk 2024-06-16 81/4050 2024-06-18 22:39 by player2012
[考博] 这个博士要读吗 +10 Sea Breeze 2024-06-16 19/950 2024-06-18 22:23 by meier8
[基金申请] 2024国社科通讯评审 +17 qsd10086 2024-06-13 34/1700 2024-06-18 22:21 by LIU_V
[硕博家园] 硕博巨婴,也许才刚刚开始 +34 SNaiL1995 2024-06-12 87/4350 2024-06-18 20:49 by 大飞鱼鱼鱼
[考博] 34岁读博士晚吗 +33 emitdne 2024-06-13 33/1650 2024-06-18 19:03 by zdsxxy
[硕博家园] 关于硕博连读的一些疑问? +6 Lwenter 2024-06-14 7/350 2024-06-18 16:54 by cuiting
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+4 BruceChum 2024-06-15 21/1050 2024-06-18 13:03 by 冰蓝夜游神
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[基金申请] 有人中过人文社科类的博后特助吗? +3 outsider1986 2024-06-16 5/250 2024-06-18 11:10 by 袁天未然
[基金申请] 太卷了 +13 laoyuefubio 2024-06-17 26/1300 2024-06-18 10:56 by shuigubio
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[基金申请] 博后面上今天有bug可以看到是否资助? +22 lyfbangong 2024-06-12 33/1650 2024-06-17 16:17 by carolloo
[有机交流] 苯酚和炔醇反应,mitsunobu反应 5+3 TONGMEIMEI 2024-06-12 9/450 2024-06-17 15:00 by dschong
[找工作] 江西双非一本和四川双一流高校如何选择? 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
信息提示
请填处理意见