应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » degauss的收敛
141/2返回列表12下一页
查看: 2777  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gemucai

木虫 (正式写手)

[交流] degauss的收敛 已有2人参与

大家好，我受到最近计算的声子谱虚频的启发，想看看一般体系取多大degauss会收敛。结果我就随便算了一下石墨烯，发现degauss不管取多大，计算出来的能量都是一直往下降的，我从0.01取到0.8（很大了吧？），能量就没有收敛的趋势。我最后两次测试所得的能量差居然还有0.02Ry，这也不能叫收敛吧？那么它如何才能算收敛呢？我在文献里看到一个方法，就是测gamma点的振动频率，并且是结合k点选取的。

我在这里想问大家，正常情况下degauss取多大范围？我觉得0.8都太大了。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-31 08:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
并且还有一个奇怪的现象：
我固定K网格计算degauss的收敛，总是无法收敛。
但是我固定degauss来确定要取多少k点，却总是收敛的。
赞一下 回复此楼
2楼2012-05-31 09:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
我把部分收敛结果贴出来：
固定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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
那么，应该如何做收敛测试才是正确的？

是固定degauss，增加k点？还是固定k点，增加degauss？
赞一下 回复此楼
4楼2012-05-31 16:26:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
xueht987: 金币+2, 谢谢回帖，鼓励交流！ 2012-06-01 09:28:30
我想这取决于你采用的smearing的方法，你可以试试别的方法，低温下要采用冷冻法估计好一些。
原则上肯定是degauss小点好，尤其是当其没有真实的物理含义的时候。0.01-0.05是典型值了。一般是固定degauss测试k点的收敛性。你的100个K点，貌似不少了，所以收敛正常。
赞一下 回复此楼
5楼2012-05-31 21:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
5楼: Originally posted by goldenfisher at 2012-05-31 21:43:06
我想这取决于你采用的smearing的方法，你可以试试别的方法，低温下要采用冷冻法估计好一些。
原则上肯定是degauss小点好，尤其是当其没有真实的物理含义的时候。0.01-0.05是典型值了。一般是固定degauss测试k点的收 ...

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by gemucai at 2012-06-01 08:06:37
thank you very much.

但是我在定义文件中只看到了这个：
var occupations -type CHARACTER {
            info {
                'smearing':     gaussian smearing for metals
                                requires a value for degau ...

继续往下看。
应该是smearing=''
赞一下 回复此楼
7楼2012-06-01 12:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
sunyang1988: 金币+1, 专家考核, 谢谢交流 2012-06-01 18:42:03
引用回帖:
6楼: Originally posted by gemucai at 2012-06-01 08:06:37
thank you very much.

但是我在定义文件中只看到了这个：
var occupations -type CHARACTER {
            info {
                'smearing':     gaussian smearing for metals
                                requires a value for degau ...

smearing='mv'就是冷冻法。由Marzari和Vanderbilt搞出来的
赞一下 回复此楼
思想重于技巧，内涵重于表象
8楼2012-06-01 12:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
非常感谢！
回复此楼
9楼2012-06-01 14:16:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
5楼: Originally posted by goldenfisher at 2012-05-31 21:43:06
我想这取决于你采用的smearing的方法，你可以试试别的方法，低温下要采用冷冻法估计好一些。
原则上肯定是degauss小点好，尤其是当其没有真实的物理含义的时候。0.01-0.05是典型值了。一般是固定degauss测试k点的收 ...

前辈，倘若我的体系里有很明显的科恩反常，并且我要研究这种科恩反常的话，是不是degauss要尽量取小？
赞一下 回复此楼
10楼2012-08-16 17:00:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +9 大米饭！ 2026-03-15 12/600 2026-03-16 21:59 by 大米饭！
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[文学芳草园] 伙伴们，祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 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 余晖&
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
关闭
请填处理意见
关闭