24小时热门版块排行榜    

查看: 8429  |  回复: 51
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 yangfangwi 的 22 个金币

yangfangwi

新虫 (正式写手)

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

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

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

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

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

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

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

zhxd1028

禁虫 (小有名气)

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

3楼2009-12-06 16:00:41
已阅   回复此楼   关注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的回帖

aijun0118

木虫 (正式写手)

博士

★ ★ ★ ★
yangfangwi(金币+2,VIP+0):恩,差别很大的。。所以才让我很想不通。。甚至怀疑vasp的指导意义啦! 12-7 11:40
qasd(金币+2,VIP+0):xiexie~ 12-7 20:05
弹性常数对K点的取值确实很敏感,K点只要保持足够的大,那么计算弹性常数的时候,如果出现与实验值有偏差的话,一方面偏差不是太大的话是正常现象,可能大部分弹性常数是正常的,只是个别的不正常,尝试着与其他人的计算结果做个比较,或许你的结果是个新发现。另一方面,找其他原因,不要拘泥于K点的问题,这一点你要确信。
资源共享
6楼2009-12-06 18:41:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanglianli136

铁杆木虫 (著名写手)


yangfangwi(金币+1,VIP+0):轻轻的你走了,留下了一片云彩~~~~~ 12-10 10:59
引用回帖:
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低估了晶格参数,导致计算得到的弹性常数偏离实验值; ...

那篇文章被引用两次,都是你引的!哈哈
19楼2009-12-10 09:03:20
已阅   回复此楼   关注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的回帖

wangzg9292

木虫 (正式写手)


yangfangwi(金币+1,VIP+0):大体系时就必须考虑了!! 12-7 11:45
大家讲到K点,我这部分做错了,我没有进行K点测试,只是用到经验数值计算,不过结果跟别人的计算结果还算吻合,我就没再测试,可是计算新物质时,我也就直接用了。
8楼2009-12-07 09:22:08
已阅   回复此楼   关注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的回帖

zhxd1028

禁虫 (小有名气)

★ ★
qasd(金币+2,VIP+0):xiexie~ 12-7 20:05
本帖内容被屏蔽

12楼2009-12-07 16:50:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhxd1028

禁虫 (小有名气)

★ ★
zdhlover(金币+2,VIP+0):谢谢,欢迎常来交流 12-11 00:47
本帖内容被屏蔽

18楼2009-12-10 01:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhxd1028

禁虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★
yangfangwi(金币+5,VIP+0):谢谢了!!嗯 我对vasp的认识不够深入。。之把它当机械的工具了。。所以遇到问题就蒙了!可否留个其他的联系方式。。希望不嫌麻烦!! 12-10 15:19
wuchenwf(金币+3,VIP+0):谢谢 12-10 20:42
本帖内容被屏蔽

21楼2009-12-10 13:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖仅楼主可见
34楼2015-05-14 16:22:47
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

sunshinerose

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我现在也在计算弹性常数,ISYM=0取消对称性后,计算总是异常结束,没有错误提示,OUTCAR里面也没有弹性矩阵。不知道大家有遇到过这样的问题吗?该怎么解决呢?

发自小木虫IOS客户端
做一只科研girl~
46楼2016-05-27 10:28:34
已阅   回复此楼   关注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的回帖

小虫009

金虫 (初入文坛)

★ ★ ★ ★ ★
yangfangwi(金币+3,VIP+0):非常感谢!我一直忽略了sigma的问题!恩0.2算正常。。你可以象征性的给个建议吗?? 12-7 11:38
qasd(金币+2,VIP+0):xiexie~ 12-7 20:04
引用回帖:
Originally posted by yangfangwi at 2009-12-6 15:37:
我计算弹性常数大半年了,,,可是最近总是有种不祥的预感!
因为k点的选取对各个独立的弹性常数值的影响太大了,,而且如果没有实验值做依据。。那鬼知道那个是正确的??所以如何确定k是个很重要的问题,,不知 ...

我也遇到类似的问题 ,K点取大点是必须的,保证计算的精度 。如果K点取大后各个独立的弹性常数值的变化幅度还很大的话,可以尝试调大SIGMA的值,但也不要调太大,以稳定为主。
5楼2009-12-06 17:57:13
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 应届大龄土博教职选择,石河子大学和广西大学 10+4 茅山老道士 2024-06-10 34/1700 2024-06-16 22:05 by yugege2009
[找工作] 江西双非一本和四川双一流高校如何选择? 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[考博] 2025考博 +6 自强不息a?a 2024-06-15 8/400 2024-06-16 19:39 by 鱼翔浅底1
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 17/850 2024-06-16 19:14 by onroad24
[考博] 34岁读博士晚吗 +27 emitdne 2024-06-13 27/1350 2024-06-16 18:57 by yugege2009
[基金申请] 面青地会评时间 +6 tanjydd 2024-06-15 6/300 2024-06-16 15:42 by 我是王小帅
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[论文投稿] 二审返修送审10天了,原来一审的3个审稿人只有2个接受了审稿,会邀请新审稿人么? 50+3 huanpo116 2024-06-15 5/250 2024-06-16 10:27 by bobvan
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[基金申请] 面上基金有一个(两个)C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[硕博家园] 申博 +3 悦悦小小鱼 2024-06-12 3/150 2024-06-12 15:11 by chen5805
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
信息提示
请填处理意见