24小时热门版块排行榜    

查看: 7719  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

preach0519

木虫 (小有名气)

[求助] 关于lammps中expected floating point parameter in input 已有1人参与

如题 想要重复文献中的计算但是报错ERROR: Expected floating point parameter in input script or data file (../pair_buck_long_coul_long.cpp:203)
请大家帮忙看一下这是什么问题啊?
# alsicao
units metal
boundary p p p
atom_style atomic
dimension 3
region box block 0 38.7320 0 38.7320 0 38.7320
create_box 4 box
create_atoms 1 random 425 78913542 box
create_atoms 2 random 908 87577830 box
create_atoms 3 random 170 78686789 box
create_atoms 4 random 2496 96789413 box
mass 1 40.0798 #ca
mass 2 28.0860 #si
mass 3 26.9820 #al
mass 4 15.9994 #o
pair_style buck/long/coul/long cut long 10
pair_coeff 1 1 buck/long/coul/long 9684.976 1 1/3.448
pair_coeff 1 2 buck/long/coul/long 1362.401 1 1/4.492
pair_coeff 1 3 buck/long/coul/long 4879.785 1 1/3.448
pair_coeff 1 4 buck/long/coul/long 3718.745 1 1/3.448
pair_coeff 2 2 buck/long/coul/long 1.866e19 1 1/40.00
pair_coeff 2 3 buck/long/coul/long 2219.246 1 1/3.448
pair_coeff 2 4 buck/long/coul/long 223440.540 1 1/7.018
pair_coeff 3 3 buck/long/coul/long 2444.136 1 1/3.650
pair_coeff 3 4 buck/long/coul/long 1945.759 1 1/3.546
pair_coeff 4 4 buck/long/coul/long 15812.842 1 1/3.846
fix 1 all nvt temp 5000 5000 100
timestep 1
run 30000
fix 1 all nvt 5000 2000 100 drag 500
run 20000
fix 1 all nvt 2000 2000 100
run 30000
compute 3 all pe/atom
compute 4 all ke/atom
compute 5 coord/atom 3.0
dump 1 all custom 1 alsicao.lammpstrj id xs ys zs c_3 c_4 c_5
回复此楼
自己才是最靠得住的人!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

preach0519

木虫 (小有名气)

引用回帖:
4楼: Originally posted by sunep07 at 2015-05-01 19:41:51
请教楼主:使用random类型产生的原子很有可能是相互重叠的(译注:或者是距离较近,不太符合实际物理),从而导致计算出较大的力或能量。因此,在你开始进行正常的动力学计算时,最好先使用minimize对体系进行能量最 ...

请教谈不上 我也不是很懂 大家多多交流吧 嘿嘿 这个我之前没用过 刚看了一下手册里面说球形的原子默认的直径是1.0 但是并没有说能够避免重复的问题 所以感觉先对体系弛豫一下还是有必要的 我自己的理解是这样的 不对的地方多多指教啊
自己才是最靠得住的人!
6楼2015-05-10 21:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

preach0519

木虫 (小有名气)

log文件如下:
LAMMPS (29 Aug 2014-ICMS)
WARNING: OMP_NUM_THREADS environment is not set. (../comm.cpp:88)
  using 1 OpenMP thread(s) per MPI task
# alsicao
units metal
boundary p p p
atom_style atomic
dimension 3
region box block 0 38.7320 0 38.7320 0 38.7320
create_box 4 box
Created orthogonal box = (0 0 0) to (38.732 38.732 38.732)
  1 by 1 by 1 MPI processor grid
create_atoms 1 random 425 78913542 box
Created 425 atoms
create_atoms 2 random 908 87577830 box
Created 908 atoms
create_atoms 3 random 170 78686789 box
Created 170 atoms
create_atoms 4 random 2496 96789413 box
Created 2496 atoms
mass 1 40.0798 #ca
mass 2 28.0860 #si
mass 3 26.9820 #al
mass 4 15.9994 #o
pair_style buck/long/coul/long cut long 10
pair_coeff 1 1 buck/long/coul/long 9684.976 1 1/3.448
ERROR: Expected floating point parameter in input script or data file (../pair_buck_long_coul_long.cpp:203)
自己才是最靠得住的人!
2楼2014-12-03 21:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ysshj

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
preach0519(abinitio代发): 金币+2, 鼓励交流 2014-12-04 21:13:17
pair_coeff 里的参数部分不能有运算吧,比如 1/3.448 。
可以用variable equal 算 1/3.448 等,pair_coeff 里用v_或${},建议用${},v_可能也会有问题。
一个人要像一支队伍。
3楼2014-12-04 15:05:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunep07

禁虫 (小有名气)

本帖内容被屏蔽

4楼2015-05-01 19:41:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见