| 查看: 1378 | 回复: 2 | |||
[求助]
vasp几何优化离子收敛的问题 已有1人参与
|
|
我做的几何优化计算，我认为离子步没有达到收敛标准啊，输出结果怎么说精度已经被满足了呢？达到标准的含义是什么啊？ 我的INCAR是 SYSTEM = crystle in a box ISMEAR=-5 ISPIN=1 PREC=Accurate ENCUT=700 ISYM=0 NELMIN = 8 # do a minimum of ten electronic steps NELM=100 EDIFF = 4E-6 # high accuracy for electronic groundstate NSW =150 # 150 ionic steps should do IBRION = 2 # ISIF=3 #all are optimized NPAR=4 那么离子的精度默认应该为EDIFF = 4E-5 但是我最后的OUTCAR是这样： FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV) --------------------------------------------------- free energy TOTEN = -36.36260684 eV energy without entropy= -36.36260684 energy(sigma->0) = -36.36260684 d Force = 0.2275733E-03[-0.108E-04, 0.466E-03] d Energy = 0.3155840E-03-0.880E-04 d Force =-0.9954701E+00[-0.997E+00,-0.994E+00] d Ewald = 0.3477743E+01-0.447E+01 reached required accuracy - stopping structural energy minimisation OUTCAR前面是这样： Electronic Relaxation 1 ENCUT = 700.0 eV 51.45 Ry 7.17 a.u. 16.85 16.85 43.15*2*pi/ulx,y,z ENINI = 700.0 initial cutoff ENAUG = 531.4 eV augmentation charge cutoff NELM = 100; NELMIN= 8; NELMDL= -5 # of ELM steps EDIFF = 0.4E-05 stopping-criterion for ELM LREAL = F real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 0 Vosko Wilk Nusair interpolation ROPT = 0.00000 0.00000 Ionic relaxation EDIFFG = 0.4E-04 stopping-criterion for IOM NSW = 150 number of steps for IOM NBLOCK = 1; KBLOCK = 150 inner block; outer block IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG NFREE = 1 steps in history (QN), initial steepest desc. (CG) ISIF = 3 stress and relaxation IWAVPR = 11 prediction: 0-non 1-charg 2-wave 3-comb ISYM = 0 0-nonsym 1-usesym 2-fastsym LCORR = T Harris-Foulkes like correction to forces |
回复此楼» 猜你喜欢
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
帮忙看看fileCode
已经有5人回复
-
FileCode能看出啥？
已经有13人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP计算时，TOTAL-FORCE这么高是什么原因，求大神
已经有8人回复
- 懂ms和vasp 以及电荷布局分析的高手进来～
已经有6人回复
- 已经结构优化好，加了一个电子，算单点能，为啥就算不下来了？
已经有11人回复
- 用vasp优化Ni2团簇出现的问题！
已经有18人回复
- VASP优化几何时，ISIF=3时，如何只改变晶格，而不改变位矢方向.
已经有11人回复
- vasp结构优化出现问题：算完一个离子步后就停止
已经有12人回复
- castep 几何优化后的晶格常数在哪里(正交晶系)?
已经有6人回复
- 想请问一下，几何优化之后的晶格角度变化的问题
已经有12人回复
- 关于晶体优化问题
已经有6人回复
- vasp表面计算的结构优化问题
已经有7人回复
- vasp 几何优化在先？参数优化在先？
已经有17人回复
- vasp结构优化中EDIFF以及KPOINTS设置，以及氢吸附结合能问题！
已经有17人回复
- 如何做vasp掺杂
已经有15人回复
- 多目标优化经典测试问题
已经有12人回复
- VASP结构优化时出现的POTIM问题
已经有4人回复
- 求助如何固定原子坐标进行优化，快搞疯了
已经有16人回复
- VASP计算能带时是否读取WAVECAR
已经有16人回复
- vasp优化结构中出现的问题
已经有10人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助成功】晶格优化的结果很诡异
已经有13人回复
- 【求助】帮忙VASP写个小程序
已经有21人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
awdawdawdd
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 679
- 散金: 200
- 红花: 6
- 帖子: 111
- 在线: 67小时
- 虫号: 2237233
- 注册: 2013-01-12
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-05-16 06:59:21
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-05-16 06:59:21
|
因为你设置的ionic relaxation的criteria是energy啊，可以看到你最后给出的结果d Energy值是小于4E-5的。 还有一般ionic relaxation的criteria一般都设置为force condition吧。一般EDIFFG=-0.01 eV/A 就够了。之后如果你不放心的话在OUTCAR搜索优化最后一步的“TOTAL-FORCE” 然后可以看到每一个原子上的受力，每个力都小于你的EDIFFG就可以了。 例如：所有原子受力都小于0.01 eV/A POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- -0.87110 7.28732 6.55372 0.003283 -0.000509 0.002108 -1.18676 3.72336 6.55369 -0.001249 0.001833 0.002096 2.05755 5.23204 6.55361 -0.000977 -0.002545 0.002264 |
2楼2015-05-16 01:11:59
3楼2015-05-18 15:52:39
