| 查看: 2261 | 回复: 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 有做缺陷形成能计算的朋友帮我看一下吧，先谢谢各位了 |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有11人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
2楼2021-02-01 23:50:02
3楼2021-02-01 23:52:35
4楼2021-04-12 13:34:25
