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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 势函数的构件
51/1返回列表
查看: 1776  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yongning88

铜虫 (小有名气)

[求助] 势函数的构件 已有1人参与

我按照文献中的公式和参数计算了一下pair-potentiial, 算出的值和文献中给出的图一样（如图3所示），可为什么当 r 小于1.5 A 时，对势就很大（如图5所示），文献中没给出这部分的值，而LAMMPS中给出的势文件却是有个峰值的（如图4所示）。跪求大牛们的帮助。

势函数的构件
function.PNG


势函数的构件-1
parmeters.PNG


势函数的构件-2
pair_potential_in-paper.PNG


势函数的构件-3
pair_potential_for_LAMMPS.PNG


势函数的构件-4
pair_potenial_in_paper.PNG
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-29 09:51:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
月只蓝: 金币+3, 感谢指导！ 2014-05-16 20:11:11
引用回帖:
8楼: Originally posted by yongning88 at 2014-05-05 09:42:33
谢谢您的解释，既然r很小是可以忽略不用拟合，那为什么在有些已有的势函数文件中，对势值是从零开始增长，直到出现峰值，他这有事怎么处理的呢？？...

他既然提出了一个函数，这个函数的参数拟合出来之后，自然会有一个函数值随r变化的过程，r从零开始增长，不久出现峰值，这是这个函数的性质所决定的，是一个数学问题。而这个函数在r很小的时候的性质，是不影响它的使用的。一个势函数的可用性和势函数的好坏，不是看这个函数在r很小的时候的情况，而是看这个函数在原子最近邻距离的附近开始之后的情况，比如r＝1.0ang 或者1.5ang以后的情况。
赞一下 回复此楼
9楼2014-05-06 09:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

lsloneil

专家顾问 (正式写手)
你是如何用lammps得到图4和图5的？
赞一下 回复此楼
2楼2014-04-29 10:31:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongning88

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by lsloneil at 2014-04-29 10:31:08
你是如何用lammps得到图4和图5的？

这个不是用LAMMPS得到的，是用MATLAB 算的，我是想自己构架那个可以在LAMMPS中用的势函数
赞一下 回复此楼
3楼2014-04-29 10:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
★
月只蓝: 金币+1, 鼓励交流！ 2014-04-30 13:22:50
mark一下。
楼主的那张图跟文献一样？我觉得从这个尺度看哪有一样的，或者楼主标清楚一下图几。。。
最好把文献贴出来，帖子里面只放自己的图。
赞一下 回复此楼
4楼2014-04-30 09:37:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 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文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭