24小时热门版块排行榜    

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

yangfangwi

新虫 (正式写手)

[交流] 【讨论】用vasp计算过弹性常数的请进!!更新中!!

大家好!!
我真是被计算弹性常数给搞怒了!!
为了方便大家共同探讨进步,今天特开此贴,大放金币!
希望能真的出现讨论的氛围!!
把这个小小的方向彻底攻破,理解!
不是简单的计算一堆貌似正确的数据~~~~~
(大家经常进来看看,我打算让它有生命力,互相讨论学习。你们也可以留下问题

[ Last edited by yangfangwi on 2009-12-10 at 10:55 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

我计算弹性常数大半年了,,,可是最近总是有种不祥的预感!
因为k点的选取对各个独立的弹性常数值的影响太大了,,而且如果没有实验值做依据。。那鬼知道那个是正确的??所以如何确定k是个很重要的问题,,不知道大家有什么妙招?
不是k点越大越精确,,,,我的测试结果
不知道你们遇到过这个问题没??12月6日问
嗯  一方面为了不让这么好的帖子沉下去,另一方面我又遇到了新的问题如下:
1、k点、encut、sigma是目前影响计算结果的几个明显因素。。不知道还有那些遗漏。。或者某些参数在计算弹性常数时有特殊要求的??
2、如果我的结果即独立弹性常数其他几个值(C1.C2.C3)都偏大了些。导致C4严重偏小(C4=E-C1-C2-C3)。。。。这种情况有没有存在的可能性。。。就是编辑会不会很在意??12月8日问

[ Last edited by yangfangwi on 2009-12-9 at 09:55 ]
2楼2009-12-06 15:37:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by zhxd1028 at 2009-12-6 16:00:

弹性常数计算是个蛮麻烦的事,我这学期就打算对vasp做弹性常数的影响因素做个详细的计算。有一大部分因素考虑了,k点是其中之一。
这个计算的起因是我在一片prb上看到有人讨论了计算得到的弹性常数与实验有差距 ...

呵呵,,把prb的文章分享一下啦,,,因为最近要开题自己的计算问题没找到原因很是着急!!谢谢了!
9楼2009-12-07 11:55:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by aijun0118 at 2009-12-6 18:41:
弹性常数对K点的取值确实很敏感,K点只要保持足够的大,那么计算弹性常数的时候,如果出现与实验值有偏差的话,一方面偏差不是太大的话是正常现象,可能大部分弹性常数是正常的,只是个别的不正常,尝试着与其他人 ...

恩 不知道你读过相关分析影响因素的文献没?如果有,可否把期刊号什么的发上来。。我想好好看看,找点突破口!!非常感谢!
10楼2009-12-07 15:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by 小虫009 at 2009-12-6 17:57:

我也遇到类似的问题 ,K点取大点是必须的,保证计算的精度 。如果K点取大后各个独立的弹性常数值的变化幅度还很大的话,可以尝试调大SIGMA的值,但也不要调太大,以稳定为主。

谢谢你的提示!只是你说的稳定为主??怎么满足!!
11楼2009-12-07 15:32:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by zhxd1028 at 2009-12-8 00:01:
k点肯定要测试啊,还是要计算很多k点的。因为考虑很多情况,所以计算量很大。不过做研究就要考虑完全一些。

Low-lying phonons in NaBH4 studied by inelastic scattering of synchrotron radiation  就是这篇吧??呵呵 我好好看看!
15楼2009-12-08 09:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by zhxd1028 at 2009-12-6 16:00:

弹性常数计算是个蛮麻烦的事,我这学期就打算对vasp做弹性常数的影响因素做个详细的计算。有一大部分因素考虑了,k点是其中之一。
这个计算的起因是我在一片prb上看到有人讨论了计算得到的弹性常数与实验有差距 ...

不好意思。。。我是真的很急。。。想解释目前数据出现的问题!!把你看到的那篇文献推荐一下吧???
16楼2009-12-08 09:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

牛人!赶紧来!帮帮我这只可怜的虫子
17楼2009-12-08 21:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by zhxd1028 at 2009-12-10 01:12:
Low-lying phonons in NaBH4 studied by inelastic scattering of synchrotron radiation  Phys. Rev B 78 172104 (2008)
这篇文章有部分论述,主要观点:DFT低估了晶格参数,导致计算得到的弹性常数偏离实验值; ...

嗯。。我读了。。好像只在conclusion部分分析了一下。。没有你理解的那么直观。。。好像是说
1、第一性原理高估了声子频率。(难道是这个影响了晶格参数?)
2、而高压衍射实验高估了体模量。。
另外提到DFT很难解释热力学的简谐振动和无序态。。。。因此否定了理论预测结构稳定性和热力学的准确性。。
可是我还是不能把他的结论和我们计算过程的参数联系起来?
20楼2009-12-10 10:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangfangwi

新虫 (正式写手)

引用回帖:
Originally posted by 小虫009 at 2009-12-6 17:57:

我也遇到类似的问题 ,K点取大点是必须的,保证计算的精度 。如果K点取大后各个独立的弹性常数值的变化幅度还很大的话,可以尝试调大SIGMA的值,但也不要调太大,以稳定为主。

我测试啦sigma值,,别的体系不确定,但是我计算的mg体系几乎没影响(我选择了0.2、0.3、0.4没变化)!!嗯,如果我计算的独立弹性常数普遍偏大,原因是什么呢?帮分析一下啊!谢谢
22楼2009-12-11 09:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangfangwi 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见