| 查看: 4446 | 回复: 5 | |||
| 【悬赏金币】回答本帖问题，作者DDDDDTTT将赠送您 5 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
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" |
回复此楼» 收录本帖的淘贴专辑推荐
 分子动力学-初级问题资源 |
» 猜你喜欢
墨尔本大学博后offer要不要接
已经有3人回复
-
刚刚收到科研之友邮件
已经有5人回复
-
高校两个offer选择
已经有13人回复
-
面青地会评时间？？？
已经有5人回复
-
2024国社科通讯评审
已经有43人回复
-
2025年博士申请，可先做科研助理
已经有5人回复
-
Scientific Reports二审90天了，催稿也没反应怎么办？
已经有6人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有18人回复
-
青年基金会评专家到底是怎么会评的呀？主审专家是不是一般不会改动系统按函评给的顺序
已经有14人回复
-
有化学投稿的期刊推荐吗
已经有3人回复
4楼2017-03-10 15:32:20
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11208.9
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 354.5小时
- 虫号: 589656
- 注册: 2008-08-30
- 性别: GG
- 专业: 高分子物理与高分子物理化
2楼2016-09-21 16:13:07
3楼2016-09-25 04:23:57
5楼2018-06-10 12:51:42
5