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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » degauss的收敛
141/2返回列表12下一页
查看: 2814  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考！ +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 277材料科学与工程080500求调剂 +6 自由煎饼果子 2026-03-16 6/300 2026-03-21 17:21 by 学员8dgXkO
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 261求B区调剂，科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 288求调剂，一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
信息提示
关闭
请填处理意见
关闭