24小时热门版块排行榜    

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

sy11052

新虫 (初入文坛)

[求助] 请问Al2O3的势函数单位的转换??

我想做Al2O3的模拟,调研了一些文献,发现现在大多的文献模拟Al2O3都是用的Matsui势函数。
但是有一个很大问题就是单位的转换问题;
我在lammps中用的单元是metal,能量单位是ev,现在Matsui势函数的单位是KJ/mol,我自己直观的转了一下,与所查的文献的结果查很大。
Matsui中能量值为719.48J/mol,而我查的的为12201ev
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sy11052

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by chaizhm at 2012-07-17 12:14:30
这两个有差别吗?你只要把单位统一不就行了吗...

差的很大,我将其在real下和将其在metal得到结果完全不一样。我感觉不是简单的除以阿伏伽德罗常数,A这样除完以后与我查的的Bucking势的参数相差太大,查好几个量级。不是简单的这样转换
10楼2012-07-24 00:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

chaizhm

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
zh1987hs: 金币+1, 鼓励交流 2012-07-15 11:49:40
你的能量是指什么能量呢?多少个原子?
壁立万仞,无欲则刚
2楼2012-07-13 18:01:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sy11052

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by chaizhm at 2012-07-13 18:01:51
你的能量是指什么能量呢?多少个原子?

Matsui势函数的表达式为:V=Qi*Qj/Rij+D(Bi+Bj)exp((Ai+Aj-rij)/(Bi+Bj))-Ci*Cj/Rij^6 ,其中D的值为4.184KJ/mol*A^3,Bi,Bj分别为0.034和0.138A。
而在lammps中的Bucking势函数的表达式为
V=Aexp(d-rij)/p)-C/Rij^6
我现在就想将上边的D*(Bi+Bj)如何等效到这边A,我的units为metal,故在lammps中的能量单位为ev
3楼2012-07-14 17:00:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaizhm

木虫 (著名写手)

【答案】应助回帖

★ ★
zh1987hs: 金币+2, 谢谢 2012-07-15 11:49:53
引用回帖:
3楼: Originally posted by sy11052 at 2012-07-14 17:00:46
Matsui势函数的表达式为:V=Qi*Qj/Rij+D(Bi+Bj)exp((Ai+Aj-rij)/(Bi+Bj))-Ci*Cj/Rij^6 ,其中D的值为4.184KJ/mol*A^3,Bi,Bj分别为0.034和0.138A。
而在lammps中的Bucking势函数的表达式为
V=Aexp(d-rij)/p)-C/Ri ...

两个单位都转换为标准能量单位,这个里面都转化为J

这个势能我用过,也用它自己编过程序,一点问题都没有
壁立万仞,无欲则刚
4楼2012-07-15 09:38:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见