POTIM=3 时间步长,单位 fs, 通常 1 到 3. ISIF=2 计算外界的压力. NBLOCK= 1 多少个时间步长,写一次 CONTCAR,CHG 和CHGCAR,PCDAT. KBLOCK=50 NBLOCK*KBLOCK 个步长写一次 XDATCAR.(个离子步写一次 PCDAT.) ISMEAR=-1 费米迪拉克分布. SIGMA =0.05 单位:...
6.依次点击Edit-> Edit Data-> Volumetric data -> Import,分别导入a.vasp和b.vasp,即相应原子的...
ISMEAR = -5 ; SIGMA = 0.20 PREC=Accurate ISIF=2 IBRION=-1 LWAVE=.TRUE. LCHARG=.TRUE. NSW=0 LREAL=.FALSE. EDIFFG=1E-3 LORBIT=10 不知道是哪里设置有问题还是什么原因,请求大家帮助! He-dos.png He原子位于超包内DOS图.png 返回小木虫查看更多分享...
此时energy without entropy是Free energy TOTEN在SIGMA趋于0的极限。 注意:(1)有人在算单个原子的能量时会发现单个原子的能量虽然很小但并不是0,但是按我上面的推导,固体中的结合能是相对孤立体系的能量而来的,所以单个原子得到的TOTEN肯定是0啊,原因在于我们的POTCAR不可能绝对合理,而且我们也知道计算单个原子的...
INCAR :PREC=Accurate EDIFF=1e-6 EDIFFG = 1e-4 ENCUT = 500 ISMEAR = 0 SIGMA = 0.05 ISYM...
这里的Tao_sigma是个动能密度,就是所有波函数梯度模的平方求和,Rho_sigma就是电荷密度啦,C_deltaV是个常数,他们写了篇论文扔去了JCP——又见JCP。这么个好东东发表了之后貌似没引起太大的动静,或许是大家一直在蠢蠢欲动,反正看起来故事似乎要结束了。
ISMEAR = 0; SIGMA = 0.1 IALGO = 48 GGA=91 ISPIN=2 VOSKOWN = 1 AMIX=0.2 BMIX=0.0001 AMIX_MAG=0.8 BMIX_MAG=0.0001 MAGMOM= 8*0 LREAL = .FALSE. #option for system contain less than 20 atoms NPAR = 1 LWAVE = .F. LORBIT = 1 #非自洽计算选项,计算spd分波态密度 RWIGS = 1.312 非...
ISMEAR = 0; SIGMA = 0.1 IALGO = 48 GGA=91 ISPIN=2 VOSKOWN = 1 AMIX=0.2 BMIX=0.0001 AMIX_MAG=0.8 BMIX_MAG=0.0001 MAGMOM= 8*0 LREAL = .FALSE. #option for system contain less than 20 atoms NPAR = 1 LWAVE = .F. LORBIT = 1 #非自洽计算选项,计算spd分波态密度 RWIGS = 1.312 非...
这里的Tao_sigma是个动能密度,就是所有波函数梯度模的平方求和,Rho_sigma就是电荷密度啦,C_deltaV是个常数,他们写了篇论文扔去了JCP——又见JCP。这么个好东东发表了之后貌似没引起太大的动静,或许是大家一直在蠢蠢欲动,反正看起来故事似乎要结束了。