应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 关于lammps中expected floating point parameter in input
141/2返回列表12下一页
查看: 8252  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
自己才是最靠得住的人！
1楼 2014-12-03 21:10:42
已阅   回复此楼   关注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 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭