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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】sigma优化，EENTRO的变化趋势对不对？
91/1返回列表
查看: 2012  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
明察、慎思、笃行
1楼 2010-08-10 18:26:25
已阅   回复此楼   关注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就够了。

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

hn19870519

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我想问一下楼主，怎么请教侯老师捏？
赞一下(1人) 回复此楼
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.

个人理解，供参考。
赞一下(2人) 回复此楼
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(总能)是准确的
赞一下(2人) 回复此楼
由天命,成人事！
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭