24小时热门版块排行榜    

查看: 1540  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见