| 查看: 5348 | 回复: 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" |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学-初级问题资源 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
jluchemlib
铁杆木虫 (正式写手)
- 应助: 119 (高中生)
- 金币: 11613.4
- 散金: 200
- 红花: 19
- 帖子: 521
- 在线: 356.1小时
- 虫号: 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
28