24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 1440  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dfy50

新虫 (小有名气)

[求助] 单原子在金红石结构中的能量

我想算单原子例如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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rxd337

木虫 (正式写手)

养殖场经理

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
fzx2008: 金币+2, 鼓励交流 2012-04-05 11:45:16
dfy50: 金币+2, 有帮助, 谢谢您的回复 2012-04-05 14:51:41
你是要算金红石结构的SnO2的能量吗?你的计算参数有两点设的比较奇怪:1是你的ISMEAR=0,默认值是0.2,一般设成0.05就足够小了;2是你的K-POINTS设成了G点,你的晶胞很小,不需要K点取这么小,取个5x5x7就好了。
都是假设!
3楼2012-04-05 11:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

dfy50

新虫 (小有名气)

怎么木有人帮忙下呢,是我的算法有问题么
2楼2012-04-05 08:41:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfy50

新虫 (小有名气)

引用回帖:
3楼: Originally posted by rxd337 at 2012-04-05 11:07:00:
你是要算金红石结构的SnO2的能量吗?你的计算参数有两点设的比较奇怪:1是你的ISMEAR=0,默认值是0.2,一般设成0.05就足够小了;2是你的K-POINTS设成了G点,你的晶胞很小,不需要K点取这么小,取个5x5x7就好了。

您好,谢谢您的回复,我要计算的是Sn原子在SnO2结构中的单原子能量,因为我有看到计算单原子能量KPOINT就设成G点。就是建一个晶胞,把O原子从结构中移走只剩Sn原子,这样算可以么,谢谢您啦
4楼2012-04-05 14:51:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rxd337

木虫 (正式写手)

养殖场经理

【答案】应助回帖


liliangfang: 金币+1, 谢谢回帖交流 2012-04-06 18:06:02
引用回帖:
4楼: Originally posted by dfy50 at 2012-04-05 14:51:25:
您好,谢谢您的回复,我要计算的是Sn原子在SnO2结构中的单原子能量,因为我有看到计算单原子能量KPOINT就设成G点。就是建一个晶胞,把O原子从结构中移走只剩Sn原子,这样算可以么,谢谢您啦

我很是不理解,你要把O原子都移走那还是原来的体系吗?而且我也不懂你说的“Sn原子在SnO2结构中的单原子能量”到底指的是什么,听说去似乎没什么意义。请先说明白原委。
都是假设!
5楼2012-04-06 17:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-08 21:38 by 朱云虎202
[考研] 318求调剂 +13 ykyhsa 2026-04-05 15/750 2026-04-08 21:37 by wj165256
[考研] 286求调剂 +8 Faune 2026-04-08 8/400 2026-04-08 21:33 by sophia_93
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 材料调剂 +4 hzhahg 2026-04-06 4/200 2026-04-08 16:20 by luoyongfeng
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 求调剂 +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[考研] 农学,求调剂,314分 +4 访客记录可爱 2026-04-04 4/200 2026-04-07 21:07 by 等岸
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 13/650 2026-04-06 07:24 by houyaoxu
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 298求调剂 +5 zzz,,r 2026-04-02 8/400 2026-04-04 19:55 by 蓝云思雨
[考研] 085602 找调剂 +4 逆时针快乐 2026-04-02 4/200 2026-04-04 19:32 by 蓝云思雨
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
信息提示
请填处理意见