24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2825  |  回复: 13

gemucai

专家顾问

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +8 勒布朗@ 2026-03-31 10/500 2026-03-31 22:53 by wwytracy
[考研] 0856求调剂 +9 楒桉 2026-03-28 9/450 2026-03-31 19:06 by 暮泽12
[考研] 材料科学与工程339求调剂 +3 hyz0119 2026-03-31 3/150 2026-03-31 18:31 by 1939136013狗壮
[考研] 086000调剂 +5 7901117076 2026-03-26 5/250 2026-03-31 17:45 by 544594351
[考研] 生物学296求调剂 +8 汤圆包 2026-03-29 12/600 2026-03-31 17:05 by 18828373951
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +10 XX小邓 2026-03-29 10/500 2026-03-31 16:46 by 不吃魚的貓
[考研] 求化学调剂 +12 wulanna 2026-03-28 12/600 2026-03-31 16:38 by 690616278
[考研] 085601 材料工程 313分 求调剂 +6 Ong3 2026-03-27 6/300 2026-03-31 15:18 by yanflower7133
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 295材料工程专硕求调剂 +10 1428151015 2026-03-27 10/500 2026-03-30 19:00 by 源_2020
[考研] 105500药学求调剂,一志愿山东大学药学,348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 298求调剂 +3 什么是胖头鱼 2026-03-30 5/250 2026-03-30 14:41 by 青海小西牛
[考研] 311求调剂 +10 lin0039 2026-03-26 10/500 2026-03-30 10:26 by herarysara
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
[考研] 0856求调剂 +13 zhn03 2026-03-25 14/700 2026-03-29 08:13 by fmesaito
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
信息提示
请填处理意见