| 查看: 1437 | 回复: 4 | ||
[求助]
单原子在金红石结构中的能量
|
|
我想算单原子例如Sn或者O原子在金红石型SnO2中的能量，各位大牛看看小弟这样的输入文件对不对，例如算Sn原子的时候就把金红石里面的O原子删了再计算，这样算成立么，输入文件如下 INCAR system=TiO2-Anatase-gga ISIF=3 NSW=15 ISTART=0 IBRION=2 ISMEAR=0 !0 ! for less K-point, set 0. try -5 (tetrahedron method valid for NKPT>=4) with multi-KPOINTS NFREE =15 NELM =120 NELMIN= 4 NELMDL = -12 ENCUT=250 NBANDS=60 Prec=Normal ROPT=-5E-4 -5E-4 -5E-4 ! * Making the calculations faster IALGO = 48 ! RMM-DIIS algorithm for electrons LREAL =.FALSE. ! evaluate projection operators in real space NSIM = 4 ! blocked algorithm update, four bands at a time LELF=TURE ICHARG=1 KPOINTS Gamma-point only 1 ! one k-point Rec ! in units of the reciprocal lattice vector 0 0 0 1 ! 3 coordinates and weight POSCAR RT-TiO2 ! FROM NRL 1.0 4.59373000 0.00000000 0.00000000 0.00000000 4.59373000 0.00000000 0.00000000 0.00000000 2.95812000 2 DIRECT 0.00000000 0.00000000 0.00000000 0.50000000 0.50000000 0.50000000 结果输出如下 vasp.4.6.28 25Jul05 complex POSCAR found : 1 types and 2 ions LDA part: xc-table for Ceperly-Alder, standard interpolation POSCAR, INCAR and KPOINTS ok, starting setup WARNING: wrap around errors must be expected FFT: planning ... 4 reading WAVECAR WARNING: chargedensity file is incomplete entering main loop N E dE d eps ncg rms rms(c) RMM: 1 -0.683565056167E+01 -0.68357E+01 -0.92993E+02 60 0.563E+01 RMM: 2 -0.200878692289E+02 -0.13252E+02 -0.76576E+01 60 0.118E+01 RMM: 3 -0.207429613605E+02 -0.65509E+00 -0.36739E+00 60 0.314E+00 RMM: 4 -0.207883692799E+02 -0.45408E-01 -0.34644E-01 60 0.105E+00 RMM: 5 -0.207935292095E+02 -0.51599E-02 -0.48420E-02 60 0.406E-01 RMM: 6 -0.207943351652E+02 -0.80596E-03 -0.77193E-03 60 0.165E-01 RMM: 7 -0.207944659283E+02 -0.13076E-03 -0.12731E-03 60 0.660E-02 RMM: 8 -0.207944879625E+02 -0.22034E-04 -0.21513E-04 60 0.273E-02 RMM: 9 -0.207944919941E+02 -0.40316E-05 -0.39754E-05 110 0.109E-02 RMM: 10 -0.207944923987E+02 -0.40467E-06 -0.39817E-06 96 0.361E-03 RMM: 11 -0.207944924384E+02 -0.39669E-07 -0.39119E-07 89 0.643E-04 RMM: 12 -0.207944924376E+02 0.77387E-09 0.00000E+00 73 0.000E+00 0.295E+00 RMM: 13 -0.206605450558E+02 0.13395E+00 -0.20905E-02 120 0.214E-01 0.166E+00 RMM: 14 -0.206030344905E+02 0.57511E-01 -0.42321E-02 120 0.307E-01 0.206E-01 RMM: 15 -0.206026093126E+02 0.42518E-03 -0.17013E-03 120 0.681E-02 0.691E-02 RMM: 16 -0.206025996737E+02 0.96389E-05 -0.11676E-04 109 0.175E-02 1 F= -.20602600E+02 E0= -.20432338E+02 d E =-.206026E+02 curvature: 0.00 expect dE= 0.000E+00 dE for cont linesearch 0.000E+00 trial: gam= 0.00000 g(F)= 0.352E-63 g(S)= 0.657E+02 ort = 0.000E+00 (trialstep = 0.100E+01) search vector abs. value= 0.657E+02 bond charge predicted LAPACK: Routine ZPOTRF failed! 4 是什么问题，谢谢各位大牛帮忙瞧瞧，我这样算单原子可以么 [ Last edited by dfy50 on 2012-4-4 at 16:46 ] |
回复此楼» 猜你喜欢
070300化学279求调剂
已经有8人回复
-
土木水利专硕276分求调剂
已经有9人回复
-
328求调剂
已经有16人回复
-
化学工程调剂289
已经有26人回复
-
331求调剂
已经有3人回复
-
材料与化工专硕306分找合适调剂
已经有24人回复
-
材料工程调剂
已经有8人回复
-
材料专硕调剂
已经有14人回复
-
334求调剂
已经有7人回复
-
材料工程085601，270求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算单个氧原子和分子的能量
已经有20人回复
- MS单原子能量
已经有19人回复
- [求助完毕]单原子掺杂金属氧化物，不同位置氧空位形成能不一样，这是为什么呢？
已经有10人回复
- 关于高对称性和P1对称性
已经有5人回复
- 大家用castep计算单原子O或N有没有出错？
已经有11人回复
- 【求助】单原子的总能量怎么算？
已经有37人回复
- 【求助】计算杂质形成能时 如何算单个原子的能量
已经有11人回复
2楼2012-04-05 08:41:55
rxd337
木虫 (正式写手)
养殖场经理
- 应助: 13 (小学生)
- 金币: 2475.2
- 红花: 1
- 帖子: 372
- 在线: 117.9小时
- 虫号: 256580
- 注册: 2006-06-03
- 性别: GG
- 专业: 材料物理化学
3楼2012-04-05 11:07:00
4楼2012-04-05 14:51:25
rxd337
木虫 (正式写手)
养殖场经理
- 应助: 13 (小学生)
- 金币: 2475.2
- 红花: 1
- 帖子: 372
- 在线: 117.9小时
- 虫号: 256580
- 注册: 2006-06-03
- 性别: GG
- 专业: 材料物理化学
5楼2012-04-06 17:01:06
