| 查看: 2651 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
gemucai木虫 (正式写手)
|
[交流]
degauss的收敛 已有2人参与
|
||
|
大家好，我受到最近计算的声子谱虚频的启发，想看看一般体系取多大degauss会收敛。结果我就随便算了一下石墨烯，发现degauss不管取多大，计算出来的能量都是一直往下降的，我从0.01取到0.8（很大了吧？），能量就没有收敛的趋势。我最后两次测试所得的能量差居然还有0.02Ry，这也不能叫收敛吧？那么它如何才能算收敛呢？我在文献里看到一个方法，就是测gamma点的振动频率，并且是结合k点选取的。 我在这里想问大家，正常情况下degauss取多大范围？我觉得0.8都太大了。 |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 最好的梦是一觉醒来，收敛了……
已经有16人回复
- 自旋体系的声子计算应该用什么软件呢？
已经有5人回复
- PWSCF可以计算半导体材料费米能级变化吗？
已经有19人回复
- 【求助完结】Au原子的磁性计算为0，为什么不是1呢？求大侠解答
已经有9人回复
- 【求助】pwscf错误提示diagonalization(ZHEGV*) falied
已经有8人回复
- 【求助】求助pwscf一些参数的物理意义和设定经验
已经有13人回复
- 【求助】pwscf 能带计算
已经有16人回复
- 【求助】结构优化计算了几天，一直不停止，请大家帮忙！
已经有15人回复
- 【求助成功】结构优化过程中遇到了一些问题，请大家帮助！
已经有7人回复
- 【求助】PWSCF计算介电极化时为何出现如下提示？
已经有19人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
- 【求助】为什么计算出的PDOS自旋向上态存在负值？
已经有9人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
2楼2012-05-31 09:46:04
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
|
我把部分收敛结果贴出来： 固定k点变degauss: degauss efermi etot 0.01 -7.9121 -27.44587936 0.02 -7.7659 -27.44681480 0.03 -7.6153 -27.44800315 0.04 -7.4600 -27.44931083 0.05 -7.2994 -27.45069542 0.06 -7.1339 -27.45213457 可见，degauss越大，总能和费米能有发散趋势。 下面是固定degauss优化k点 kpoint efermi etot 200 -7.7659 -27.44681480 300 -7.7659 -27.44681451 100 -7.7658 -27.44681471 收敛的！ [ Last edited by gemucai on 2012-5-31 at 10:04 ] |
3楼2012-05-31 10:03:22
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
4楼2012-05-31 16:26:45
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
|
thank you very much. 但是我在定义文件中只看到了这个： var occupations -type CHARACTER { info { 'smearing': gaussian smearing for metals requires a value for degauss 'tetrahedra' : for metals and DOS calculation (see PRB49, 16223 (1994)) Requires uniform grid of k-points, automatically generated (see below) Not suitable (because not variational) for force/optimization/dynamics calculations 'fixed' : for insulators with a gap 'from_input' : The occupation are read from input file. Presently works only with one k-point (LSDA allowed). 没看见有冷冻法啊。还望前辈赐教，什么是冷冻法。谢谢！ |
6楼2012-06-01 08:06:37
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
9楼2012-06-01 14:16:55
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
10楼2012-08-16 17:00:15
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
12楼2012-08-17 08:08:35
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
14楼2012-08-17 14:56:57