24小时热门版块排行榜    

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

yangfangwi

新虫 (正式写手)

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

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

[ Last edited by yangfangwi on 2009-12-10 at 10:55 ]
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 vasp 第一性原理

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
yangfangwi(金币+5,VIP+0):受益匪浅啊!谢谢牛人!可是如果是个很大的体系如96个原子的mg六角体系,,它的精度要如何取值才相对合适呢?? 12-7 11:49
qasd(金币+3,VIP+0): 12-7 20:05
弹性性能的计算确实要求k点精度比较高,因为计算中涉及到晶格的改变。
其他的计算,如电子结构计算以及原子位置的弛豫,都是在相同晶格下计算的,布里渊区形状没有发生改变,因此即使k点收敛测试结果不好,对结果的影响也不是很大,定性的结论还是正确的。
而计算弹性性能,需要晶格发生变形,计算变形前后的能量差,而这会导致布里渊区形状发生改变,相当于k点取样发生了改变,如果能量的计算对于k点取样敏感,将对计算得到的能量差产生很大的影响。因此必须保证k点取样的改变不会对能量计算造成太大影响,也就是通常说的收敛性测试一定要满足。
同样截断能也要达到收敛

» 本帖已获得的红花(最新10朵)

7楼2009-12-06 19:26:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 52 个回答

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的回帖

zhxd1028

禁虫 (小有名气)

★ ★ ★ ★ ★
yangfangwi(金币+5,VIP+0):谢谢啊,一定拿出来晒晒。。可否把你那篇prb文章上传or给个期刊号??万分感谢!! 12-7 11:29
本帖内容被屏蔽

3楼2009-12-06 16:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhxd1028

禁虫 (小有名气)

★ ★ ★ ★ ★
yangfangwi(金币+3,VIP+0):嗯,我也这么考虑过,如果一味的靠近实验值就失去计算的指导意义啦!其他原因?不知道你有什么独到见解? 12-7 11:33
qasd(金币+2,VIP+0):xiexie~ 12-7 20:04
本帖内容被屏蔽

4楼2009-12-06 16:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见