24小时热门版块排行榜    

查看: 2073  |  回复: 8

gleerat

木虫 (正式写手)

[交流] 【求助成功】sigma优化,EENTRO的变化趋势对不对? 已有4人参与

这是我的sigma优化结果,像是sigma取值越小越好,这种趋势对不对?如果是这样的,那么当EENTRO的值取到什么量级时的SIGMA是比较好的?在sigma优化中好像没见过比0.01更小的sigma取值。
引用回帖:
0.001 -0.00000026
0.002 -0.00002343
0.003 -0.00007110
0.004 -0.00012185
0.005 -0.00017412

0.01 -0.00056390
0.02 -0.00222955
0.04 -0.00739715
0.06 -0.01604789
0.08 -0.03017025
0.10 -0.05057037
0.12 -0.07786173
0.14 -0.11236035
0.16 -0.15435267
0.18 -0.20452909
0.20 -0.26359441
0.22 -0.33240921
0.24 -0.41170476
0.26 -0.50093495
0.28 -0.60002319

引用回帖:
System = CaFe2As2 afm1-opt
PREC = Accurate
ISPIN = 2
MAGMOM = 4*0 3 3 -3 -3 3 -3 -3 3 8*0
ISMEAR = 0
SIGMA = 0.04
ISTART = 0
EDIFF = 10E-5
ICHARG = 2
ENCUT = 310

LORBIT = 11
LREAL = Auto

[ Last edited by gleerat on 2010-8-13 at 18:39 ]
回复此楼

» 猜你喜欢

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

明察、慎思、笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

★ ★
zzy870720z(金币+2):谢谢总结 2010-08-14 08:33:07
我问过候老师了,他是这样回复的:
引用回帖:
并不是要将考虑将SIGMA设置到使得entropy接近0。sigma的值,决定了电子的占有数,当它取不同值时,对金属体系的Fermi level处的电子占有数变化很大。
当对sigma在一个k-grid下进行测试时,然后再将k-grid,看entropy相对于前一个k-grid计算下的entropy的变化。
比如
sigma( s1,s2, .... s10) 在 k-grid: 8x8 x8 下面的entropy ( e1, e2, ...., e10), 然后对另一个大的k-grid: 10 x 10 x10下面的entropy( e'1, e'2, ...., e'10),其他的参数不变,看
e'1- e1,  e'2-e2,  ....., e'10- e10的变化情况。

e'_i - e_i的差值在 几个meV就够了。

再次谢谢候老师。
明察、慎思、笃行
2楼2010-08-13 18:39:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我想问一下楼主,怎么请教侯老师捏?
3楼2010-08-14 01:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by hn19870519 at 2010-08-14 01:23:32:
我想问一下楼主,怎么请教侯老师捏?

用邮箱问的。但是我不知道候老师是否愿意我散布他的邮箱,我是在论坛上看到地址的。你可以找找。
明察、慎思、笃行
4楼2010-08-14 09:23:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
好的,谢谢
5楼2010-08-14 10:04:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):感谢讨论 2010-08-14 21:38:32
引用回帖:
Originally posted by gleerat at 2010-08-10 18:26:25:
这是我的sigma优化结果,像是sigma取值越小越好,这种趋势对不对?如果是这样的,那么当EENTRO的值取到什么量级时的SIGMA是比较好的?在sigma优化中好像没见过比0.01更小的sigma取值。



[ Last edited by  ...

这个帖子是求助成功贴,但是问题还没有完。

SIGMA的选取与收敛速度有关,太小收敛慢,太大容易不收敛。但是最终要做到熵比较少才行,这样自由能和中能之间的差值就比较小。多小才算比较小呢?平均到每个原子上的能量差为1meV.也有的说是摊到每个原胞单位上自由能与总能的差在1meV。也就是小数点后第三位。因此对于原子比较多的超胞,比如100个原子,自由能与总能之间差0.1eV也是可以接受的。不过,大部分可以做到0.01eV.

个人理解,供参考。
6楼2010-08-14 14:21:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by cenwanglai at 2010-08-14 14:21:26:



这个帖子是求助成功贴,但是问题还没有完。

SIGMA的选取与收敛速度有关,太小收敛慢,太大容易不收敛。但是最终要做到熵比较少才行,这样自由能和中能之间的差值就比较小。多小才算比较小呢?平均到每个 ...

熵值最小,那EENTRO不就是熵么?
请教,到底那个是总能呢?我以前经常用TOTEN,而且见过最多的也是TOTEN(在候柱峰的《入门指南》中,提到总能也是用TOTEN)。有人说TOTEN就是总能,有人说总能在OSZICAR中的最后一行,有人说energy(sigma->0) 是总能,也有人说energy  without entropy是总能,这几个值经常不一样,到底要以那个为准呢?
自由能是什么呢?顺便给我解释一下吧:
引用回帖:
  free energy    TOTEN  =       -33.75078681 eV

  energy without entropy =      -33.72545217  energy(sigma->0) =      -33.74234193

这三个值是什么意思。
不胜感激。

[ Last edited by gleerat on 2010-8-14 at 18:26 ]
明察、慎思、笃行
7楼2010-08-14 18:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-08-14 21:39:00
计算总能是需要单独算一次
这时设定ISMEAR=-5
此时TOTEN(总能)是准确的
由天命,成人事!
8楼2010-08-14 18:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by lono75 at 2010-08-14 18:35:43:
计算总能是需要单独算一次
这时设定ISMEAR=-5
此时TOTEN(总能)是准确的

那么请问,自由能是什么呢?如何找到?
明察、慎思、笃行
9楼2010-08-15 10:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 338求调剂 +4 18162027187 2026-03-02 5/250 2026-03-02 22:18 by sunny81
[考研] 290分材料工程085601求调剂 数二英一 +8 llx0610 2026-03-02 9/450 2026-03-02 22:09 by 无际的草原
[考研] 085602化学工程350,调剂,有没有211的 +4 利好利好. 2026-03-02 7/350 2026-03-02 21:46 by sunny81
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +5 虫虫虫虫虫7 2026-03-01 9/450 2026-03-02 21:33 by sunny81
[考研] 085600求调剂 +4 LRZZZZZZ 2026-03-02 5/250 2026-03-02 20:33 by 云旗yunqi
[考研] 环境调剂 +4 chenhanheng 2026-03-02 4/200 2026-03-02 20:23 by hypershenger
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 282求调剂 +4 2103240126 2026-03-02 6/300 2026-03-02 18:07 by 2103240126
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 303求调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-02 15:01 by 向上的胖东
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 7/350 2026-03-02 10:38 by lature00
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
请填处理意见