24小时热门版块排行榜

返回列表

DDDDDTTT

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.1
帖子: 3
在线: 3.1小时
虫号: 3871280
注册: 2015-05-15
专业: 粘合剂

[求助] 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"

回复此楼