| 查看: 969 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者蔡山万大鹏将赠送您 20 个金币 | ||
[求助]
CINEB重复文献的时候，比文献结果略高一点
|
||
|
体系：GaN 晶体P6_mc点群， 96原子方形超胞 问题：Ga 空位扩散势垒 方法：vasp， CINEB 计算过程： 1. 建立GaN超胞，设置好"初态"和"末态"Ga空位，PBE泛函进行原子位置优化得到CONGTCAR_i ; CONTCAR_f 。 2. 利用CINEB 插点 产生过程路径文件 “00 01 02 03 04 05 06” 3.利用CINEB优化计算到收敛，得到势垒 问题： 1. 利用了和文献一样的晶格常数和提到的参数，计算势垒为3.75eV，大部分文献为1.7~2eV ，不知道是INCAR不科学还是有什么地方没注意。 2. POSCAR中原子手动改成了一一对应，但是没有进行固定，不知是否需要固定部分原子。 计算文件 和 计算结如下： ##INCAR## ISTART = 0 ICHARG = 2 LREAL = F #K空间 ISPIN = 2 #自旋 ISYM = 0 #不对称 #PREC = A ENCUT = 450 #文献一致 #Ionic Relaxation ISIF = 2 NSW = 300 EDIFFG = -0.01 #文献一致 SYMPREC = 1E-5 ADDGRID = T #Electronic Relaxation ISMEAR = 0 SIGMA = 0.01 EDIFF = 1e-6 GGA = PS ALGO = V #优化算法 NELECT = 854 IDIPOL = 4 EPSILON = 9 LMONO = T #For CINEB IBRION = 3 POTIM = 0 IOPT = 1 ICHAIN = 0 LCLIMB = TRUE SPRING = -5 IMAGES = 5 ##KPOINTS## MK 222 ##计算结果## [grid269@lon12%tianhe2-B past]$ nebef.pl 0 0.000000 -603.044900 0.000000 1 0.002704 -600.672600 2.372300 2 0.002860 -599.723800 3.321100 3 0.005770 -599.293200 3.751700 （势垒） 4 0.002959 -599.728900 3.316000 5 0.004401 -600.697400 2.347500 6 0.000000 -603.044800 0.000100  1593828598.jpg |
回复此楼» 猜你喜欢
材料，纺织，生物（0856、0710），化学招生啦
已经有7人回复
-
301求调剂
已经有5人回复
-
281求调剂（0805）
已经有7人回复
-
材料与化工求调剂
已经有6人回复
-
299求调剂
已经有4人回复
-
277调剂
已经有6人回复
-
341求调剂
已经有7人回复
-
328求调剂，英语六级551，有科研经历
已经有8人回复
-
293求调剂
已经有12人回复
-
070300化学319求调剂
已经有4人回复
