24小时热门版块排行榜    

Znn3bq.jpeg
查看: 8094  |  回复: 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

木虫 (小有名气)

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的回帖

sunep07

禁虫 (小有名气)

本帖内容被屏蔽

5楼2015-05-09 11:52:45
已阅   回复此楼   关注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的回帖

preach0519

木虫 (小有名气)

引用回帖:
5楼: Originally posted by sunep07 at 2015-05-09 11:52:45
你好,请教一下:
我模拟磷酸盐玻璃,势函数编写如下,结果运行几步之后E_pair 大到10^86,密度,buck参数是完全按照文献输入的,不知道错出自哪里:
set       group   o     charge  -1.2
set       group   p ...

交流啊 请教谈不上 哈哈 我没太明白 你前面的pair_coeff      * *    0      1   0是什么的势函数参数啊?还有就是 后面的参数里面只有1与其它的?那2 2 2 3等其它原子类型之间就没有势函数参数的关系了吗?
自己才是最靠得住的人!
7楼2015-05-10 21:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunep07

禁虫 (小有名气)

本帖内容被屏蔽

8楼2015-05-11 14:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunep07

禁虫 (小有名气)

本帖内容被屏蔽

9楼2015-05-11 14:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

preach0519

木虫 (小有名气)

引用回帖:
9楼: Originally posted by sunep07 at 2015-05-11 14:20:55
开始设定所有原子对之间的buck 势能为0,但是这样的话库伦力也为零了不是?我的本意是让buck为零,22,23之间只考虑库伦力。...

能这样用吗?我没试过 觉得这块是有问题的 你可以把你想要为零的地方设为零 我个人观点 你再考虑一下
自己才是最靠得住的人!
10楼2015-05-11 22:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 preach0519 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 6/300 2026-05-21 01:37 by 17520222914
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +11 1234567wang 2026-05-17 13/650 2026-05-20 18:55 by tanfancy
[基金申请] 提交了我也来说说感想 +6 fummck 2026-05-20 6/300 2026-05-20 17:29 by 姜太公烤鱼
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 8/400 2026-05-20 09:32 by xtlilibin
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
信息提示
请填处理意见