24小时热门版块排行榜    

查看: 1607  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 14:36 by winsaint
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
信息提示
请填处理意见