24小时热门版块排行榜    

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见