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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » degauss的收敛
141/2返回列表12下一页
查看: 2636  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭