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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 采用murn.f计算体积弹性模量不正确？
51/1返回列表
查看: 3705  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

pgi1706

木虫 (小有名气)

[求助] 采用murn.f计算体积弹性模量不正确？

采用候老师的murn.f处理数据，计算体积弹性模量，发现与文献差别很大，不知是何原因？
输入文件如下：
2
0.03125
10.891 10.900 50
10
10.891 -.22071146E+03
10.892 -.22071168E+03
10.893 -.22071160E+03
10.894 -.22071159E+03
10.895 -.22071141E+03
10.896 -.22071119E+03
10.897 -.22071098E+03
10.898 -.22071078E+03
10.899 -.22071043E+03
10.900 -.22071035E+03
其中第一行 2 表示采用eV为能量单位，第二行是原胞与晶胞比例，结构为面心立方，我采用2＊2＊2超晶胞，是不是0.25／8＝0.03125？
后面行就不解释了，可是采用murn.f得到的结果部分如下：

FIT BY MURNAGHAN EQUATION EQUATION OF STATE
-------------------------------------------------------------------------------

ITER ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME      E0 (EV)    SUM OF SQUARES
                                                                           IERR

   20  20.5838    272.53892     72.82852   2131.16340   -220.71167      0.00000
                                                                              0
   20  20.5838    272.53892     72.82852   2131.16340   -220.71167      0.00000
                                                                              0
-------------------------------------------------------------------------------

      ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME       E0 (EV)      E0 (RY)

alat=  20.58383 b0=   72.8285 b0p=  2131.163 E0=  -8.110822
       20.5838    272.53892     72.82852   2131.16340   -220.71167    -16.22185
       10.8925
-------------------------------------------------------------------------------

        alat(A)          energy(Ryd)            alat(bohr)          ener gy(H)
-------------------------------------------------------------------------------

       10.89100         -16.2218344       20.58101          -8.1109172
       10.89118         -16.2218385       20.58136          -8.1109193
       10.89137         -16.2218418       20.58170          -8.1109209

发现B0＝72.8285Mbar，B0'=2131.163，而文献的结果却是B0＝89.67Gpa，B0'=4.44，相差几个数量级，以前在做弹性力学计算时，也总得不出正确的结果，不知是何原因？
能不能帮我验算一下？能量－晶格常数已经贴出来了，非常感谢！
回复此楼

» 猜你喜欢

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

bingmou

金虫 (著名写手)
★ ★ ★ ★
zzy870720z(金币+2): 谢谢指点 2011-09-20 12:27:22
zzy870720z(金币+2): 谢谢指点 2011-09-20 12:27:23
不要采用实验的方法验证，自己看看公式，从理论上进行解释。
你计算出来的所谓“总能量”，也只是含有若干个原子的总能量，不是实际材料的总原子个数
赞一下(2人) 回复此楼
9楼2011-09-19 23:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

pgi1706

木虫 (小有名气)
顶一下阿，在线等
回复此楼
2楼2011-09-15 14:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖

★ ★
pgi1706(金币+5): 谢谢版主，按你的意思，我再试试 2011-09-16 10:32:03
zzy870720z(金币+2): 谢谢指点 2011-09-16 18:02:17
一般体弹性模量的计算是很不准确的，但是数量级应该是正确的，而数量可以相差50%。不应该差距几个数量级的。
我记得好像应该采用原胞的长度而不是2×2×2的长度，其次是能量应该是平均到每个原子的能量
赞一下(1人) 回复此楼
4楼2011-09-16 05:03:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nright

铜虫 (小有名气)
★ ★
zzy870720z(金币+2): 谢谢指点，以后回答问题请选择“应助回帖”，呵呵 2011-09-16 18:02:52
晶格参数的变化的步长太小了，不同晶格参数对应的能量太接近，算得的弹性模量当然不准确。
这样算试一下：如果平衡体积为V0（应该接近实验V0)，算5个点，体积分别为1.10V0，1.05V0，V0，0.95V0，0.90V0，然后用这5个点的晶格参数和能量拟合，应可以得到正确的体模量。
赞一下(1人) 回复此楼
5楼2011-09-16 14:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭