应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 单原子在金红石结构中的能量
51/1返回列表
查看: 1344  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-04 16:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfy50

新虫 (小有名气)
怎么木有人帮忙下呢，是我的算法有问题么
赞一下 回复此楼
2楼2012-04-05 08:41:55
已阅   回复此楼   关注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的回帖

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的回帖
相关版块跳转 我要订阅楼主 dfy50 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭