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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » degauss的收敛
141/2返回列表12下一页
查看: 2801  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂 +7 隔壁陈先生 2026-03-12 7/350 2026-03-18 16:14 by 枫桥ZL
[考研] 344求调剂 +5 knight344 2026-03-16 6/300 2026-03-18 15:17 by knight344
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 0703化学调剂 ，六级已过，有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 331求调剂（0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
关闭
请填处理意见
关闭