24小时热门版块排行榜    

查看: 5945  |  回复: 8
【悬赏金币】回答本帖问题,作者大力2010将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

大力2010

金虫 (著名写手)

[求助] VASP结构优化以 writing wavefunctions 结束。 是不是表示不收敛? 已有2人参与

请教大神,我用vasp对掺杂了 v 和 ta 的异质结进行结构优化,incar文件时这样的:
global parameters
istart =  1            (read existing wavefunction; if there)
ispin  =  1            (non-spin polarised dft)
# icharg =  11         (non-self-consistent: gga/lda band structures)
lreal  = .false.       (projection operators: automatic)
# encut  =  400        (cut-off energy for plane wave basis set, in ev)
prec   =  normal       (precision level)
lwave  = .true.        (write wavecar or not)
lcharg = .true.        (write chgcar or not)
addgrid= .true.        (increase grid; helps gga convergence)
# lvtot  = .true.      (write total electrostatic potential into locpot or not)
# lvhar  = .true.      (write ionic + hartree electrostatic potential into locpot or not)
# nelect =             (no. of electrons: charged cells; be careful)
# lplane = .true.      (real space distribution; supercells)
# npar   = 4           (max is no. nodes; don't set for hybrids)
# nwrite = 2           (medium-level output)
# kpar   = 2           (divides k-grid into separate groups)
# ngx    = 500         (fft grid mesh density for nice charge/potential plots)
# ngy    = 500         (fft grid mesh density for nice charge/potential plots)
# ngz    = 500         (fft grid mesh density for nice charge/potential plots)
ivdw   = 11

electronic relaxation
ismear =  0            (gaussian smearing; metals:1)
sigma  =  0.05         (smearing value in ev; metals:0.2)
nelm   =  90           (max electronic scf steps)
nelmin =  6            (min electronic scf steps)
ediff  =  1e-05        (scf energy convergence; in ev)
# gga  =  ps           (pbesol exchange-correlation)

ionic relaxation
nsw    =  100          (max ionic steps)
ibrion =  2            (algorithm: 0-md; 1-quasi-new; 2-cg)
isif   =  2            (stress/relaxation: 2-ions, 3-shape/ions/v, 4-shape/ions)
ediffg = -2e-02        (ionic convergence; ev/aa)
# isym =  2            (symmetry: 0=none; 2=gga; 3=hybrids)




优化结束时,没出现“reached required accurancy”,以 writing wavefunctions 结束: 100 f= -.82977297e+03 e0= -.82954072e+03  d e =-.426037e-01
curvature:  -2.30 expect de=-0.187e+00 de for cont linesearch -0.376e-06
trial: gam= 1.72200 g(f)=  0.812e-01 g(s)=  0.000e+00 ort = 0.154e-03 (trialstep = 0.125e+00)
search vector abs. value=  0.514e+00
writing wavefunctions
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微笑述说忧伤

金虫 (小有名气)

【答案】应助回帖

主要离子步你总共就设了100步,到了100步还没收敛,就自动结束了,所以结果是并没有收敛。但是能不能收敛具体还要看你的输出文件。掺杂体系相比于本征体系而言,优化确实难一些,尤其是过渡金属掺杂,100个离子步有点少了。可以用低精度做个初步优化,再在这个基础上提高精度。另外V掺杂大概率会有磁性吧,最好加自旋。
8楼2021-12-07 20:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

物理类

新虫 (著名写手)

2楼2021-10-27 13:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大力2010

金虫 (著名写手)

引用回帖:
2楼: Originally posted by 物理类 at 2021-10-27 13:03:06
没收敛

怎样解决  求大师指点
3楼2021-10-27 17:19:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianyakl4

新虫 (小有名气)

4楼2021-10-28 00:16:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见