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

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

Originally posted by cenwanglai at 2010-08-14 14:21:26:



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

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

  free energy    TOTEN  =       -33.75078681 eV

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

并不是要将考虑将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就够了。

Originally posted by hn19870519 at 2010-08-14 01:23:32:
我想问一下楼主，怎么请教侯老师捏？