24小时热门版块排行榜    

查看: 1306  |  回复: 7

hwceng0816

木虫 (著名写手)


[交流] 求助 为什么计算了这么多次还是和文献相差挺大的,帮忙看下。

本人计算金属间化合物的弹性常数,结构换了很多次 参数也换了很多次,就是得不到文献的值,相差有点大,我发现同一个结构 用的参数不一样计算出来的结果相差很大。弹性常数计算:尤其我之前的一次K点取了8 8 6  和之后的6 6 4计算结果相差很大,为什么会这样,我已经快要绝望了,计算了一个多月了,没一次和文献接近的,我看到有的文献的精度比我还差都能得到比较准确的值,这是为什么?
我按照有些文献的参数计算 就是得不到一样的结果,哪怕是接近。

我补充说一句 我用的截断能和K点都比文献的高

[ Last edited by hwceng0816 on 2012-6-5 at 08:04 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

future_wl

木虫 (著名写手)


★ ★ ★
hwceng0816(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流。。。 2012-05-17 10:55:48
有时候K点和cutoff取的很大效果也不好,试试换赝势,我有一次换了一个原子的赝势,用默认的K点值就能优化很好了,另外你只算弹性常数的话用超软应该就行了
2楼2012-05-17 08:46:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)


引用回帖:
2楼: Originally posted by future_wl at 2012-05-17 08:46:49:
有时候K点和cutoff取的很大效果也不好,试试换赝势,我有一次换了一个原子的赝势,用默认的K点值就能优化很好了,另外你只算弹性常数的话用超软应该就行了

你说的换一个原子的赝势是怎么回事?这个我没听说过,一般都是整体选择一个赝势的。你说说的是怎么个弄法,我估计就是你说的还这个问题。 谢谢了。。
3楼2012-05-17 22:26:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)


★ ★
WDD880227: 金币+2, 多谢您的耐心指导 2012-05-18 08:55:19
引用回帖:
3楼: Originally posted by hwceng0816 at 2012-05-17 22:26:13:
你说的换一个原子的赝势是怎么回事?这个我没听说过,一般都是整体选择一个赝势的。你说说的是怎么个弄法,我估计就是你说的还这个问题。 谢谢了。。

在赝势的那个窗口中,可以给每一个原子选择不同的赝势,系统是给了一个默认的,但有时默认的会有问题
就像我这个图中那样,在norm conserving默认的赝势下,点击Ag元素后面的下拉框,就可以看到很多其它赝势,你随便换一个就行,但记住后缀名,如果你要用模守恒,必须选.recpot后缀的,超软对应.usp后缀
你肯定会问,选哪个最好呢,这个就没办法,只能试试了,一般来说选个价电子取的多,cutoff值高些的应该会好。
祝你好运!


4楼2012-05-18 08:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)


引用回帖:
4楼: Originally posted by future_wl at 2012-05-18 08:50:18
在赝势的那个窗口中,可以给每一个原子选择不同的赝势,系统是给了一个默认的,但有时默认的会有问题
就像我这个图中那样,在norm conserving默认的赝势下,点击Ag元素后面的下拉框,就可以看到很多其它赝势,你随 ...

你说的换一个赝势貌似没什么用吧  比如说超软赝势  能够选择的一般里面只有两个(后最为.usp) 我看了下两个的截断能和价电子都一样,这样有效果吗?指导一下,不胜感激!
5楼2012-06-16 19:48:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)


引用回帖:
2楼: Originally posted by future_wl at 2012-05-17 08:46:49
有时候K点和cutoff取的很大效果也不好,试试换赝势,我有一次换了一个原子的赝势,用默认的K点值就能优化很好了,另外你只算弹性常数的话用超软应该就行了

能不能加你QQ  我最近的计算出现问题  就是之前的几何优化和后面的弹性常数计算有关联吗  我想用PBE几何优化 用PW91计算弹性常数
6楼2012-06-18 16:22:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816(金币+1): 谢谢参与
你确定你的supercell和文献一样大?      

另外,886 kpoint对金属来说,不算高
7楼2012-06-18 18:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)


引用回帖:
7楼: Originally posted by franch at 2012-06-18 18:04:47
你确定你的supercell和文献一样大?      

另外,886 kpoint对金属来说,不算高

我用的是晶胞  不是超胞 有文献用VASP计算的 886
8楼2012-06-18 19:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hwceng0816 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程调剂 +4 大可digkids 2026-03-16 4/200 2026-03-16 14:00 by laoshidan
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 0703化学调剂,求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
请填处理意见