| 查看: 2835 | 回复: 13 | |||
gemucai木虫 (正式写手)
|
[交流]
degauss的收敛 已有2人参与
|
|
大家好，我受到最近计算的声子谱虚频的启发，想看看一般体系取多大degauss会收敛。结果我就随便算了一下石墨烯，发现degauss不管取多大，计算出来的能量都是一直往下降的，我从0.01取到0.8（很大了吧？），能量就没有收敛的趋势。我最后两次测试所得的能量差居然还有0.02Ry，这也不能叫收敛吧？那么它如何才能算收敛呢？我在文献里看到一个方法，就是测gamma点的振动频率，并且是结合k点选取的。 我在这里想问大家，正常情况下degauss取多大范围？我觉得0.8都太大了。 |
回复此楼» 猜你喜欢
0703化学调剂 348分
已经有13人回复
-
327求调剂
已经有8人回复
-
材料307分求大佬组收留
已经有9人回复
-
求调剂
已经有33人回复
-
293调剂
已经有15人回复
-
求调剂求调剂
已经有11人回复
-
22408 266求调剂
已经有10人回复
-
一志愿211，化学310分，本科重点双非，求调剂
已经有8人回复
-
一志愿0703化学招61最终排名62化学求调剂
已经有10人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 最好的梦是一觉醒来，收敛了……
已经有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
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
5楼2012-05-31 21:43:06
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
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
7楼2012-06-01 12:17:59
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
8楼2012-06-01 12:24:04
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
