| 查看: 4423 | 回复: 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" |
回复此楼» 收录本帖的淘贴专辑推荐
 分子动力学-初级问题资源 |
» 猜你喜欢
蹲一个男朋友
已经有48人回复
-
执念
已经有10人回复
-
经常觉得挺累的
已经有9人回复
-
化学B口多少分能上会呀
已经有6人回复
-
答辩
已经有4人回复
-
TsCl保护羟基为什么不反应
已经有21人回复
-
期刊投稿进度
已经有8人回复
-
为什么有的影响因子高的期刊分区不高呢？
已经有5人回复
-
又忍不住想吐槽文科的学术圈
已经有9人回复
-
大扩招
已经有14人回复
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11200.4
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 354.5小时
- 虫号: 589656
- 注册: 2008-08-30
- 性别: GG
- 专业: 高分子物理与高分子物理化
2楼2016-09-21 16:13:07
3楼2016-09-25 04:23:57
4楼2017-03-10 15:32:20
5楼2018-06-10 12:51:42
6楼2018-06-10 12:53:26
15