| 查看: 2214 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者meizhixin将赠送您 20 个金币 | ||
[求助]
缺陷形成能计算结果跟文献出入很大，帮忙看一下问题，谢谢各位了
|
||
|
计算缺陷形成能，利用文献上的公式，但结果感觉很不靠谱。。。。 电中性的结果还好，带负电荷的缺陷其形成能计算结果显示在整个费米能级区间都是负的。。。负的十几电子伏特 修正项数值较小，对结果的影响小，我觉得是带电缺陷能量的计算结果影响比较大，但是我计算能量用的参数都一样，只是NELECT的设置不同 下面是一个空位缺陷的相关计算信息 电荷量 含缺陷体系能量 完美超胞能量 化学势 费米能级范围 -3 -457.2247659 -501.1467179 -9.285 0-2.25 -2 -467.2487878 -501.1467179 -9.285 0-2.25 -1 -477.1471905 -501.1467179 -9.285 0-2.25 0 -487.0161765 -501.1467179 -9.285 0-2.25 1 -496.8572304 -501.1467179 -9.285 0-2.25 2 -506.6024742 -501.1467179 -9.285 0-2.25 3 -514.48853 -501.1467179 -9.285 0-2.25 下面给出我计算能量的INCAR，以-1价缺陷为例 SYSTEM = VC ISTART = 0 # (Read existing wavefunction; if there) ISPIN = 2 # (Spin polarised DFT) LREAL = Auto # (Projection operators: automatic) ENCUT = 550 # (Cut-off energy for plane wave basis set, in eV) PREC = Accurate # (Precision level) LWAVE = .T. # (Write WAVECAR or not) LCHARG = .T. # (Write CHGCAR or not) LVHAR = .TRUE. # (Write ionic + Hartree electrostatic potential into LOCPOT or not) NELECT = 253 # (No. of electrons: charged cells; be careful) NPAR = 4 # (Max is no. nodes; don't set for hybrids) IVDW = 12 ISYM = 0 # Lattice Relaxation ISMEAR = 0 # (gaussian smearing method ) SIGMA = 0.05 # (please check the width of the smearing) EDIFFG = -0.001 # (Ionic convergence; eV/AA) EDIFF = 1E-06 GGA = PE 有做缺陷形成能计算的朋友帮我看一下吧，先谢谢各位了 |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有9人回复
-
今天务委会开完了，明天出结果吗
已经有18人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
2楼2021-02-01 23:50:02
3楼2021-02-01 23:52:35
4楼2021-04-12 13:34:25
28