| 查看: 836 | 回复: 1 | |||
wqg1806木虫 (著名写手)
科研农民工
|
[交流]
【求助】盼高手指点，NEB计算
|
|
我在计算过渡态的时候插入了4个点，计算收敛之后处理数据的时候发现找到的不是能量的鞍点。在这种情况下应该怎么办？是不是继续优化，找到鞍点？或者提高计算精度，去寻找鞍点？我把计算的精度设置以及得到的数据贴一下：SYSTEM =Nb O ENCUT =400 Ionic relaxation EDIFF = 1E-5 ! low accuracy EDIFFG =-0.01 !force #PREC = high #NELM = 100 ! maximum electronic steps #NELMIN = 4 ! minimum electronic steps #NELMDL = -8 ! update of charge for 3 steps (for continuation: 0) #IBRION = -1 ! use CG algorithm for ions #IBRION = 2 ! use CG algorithm for ions NSW =1550 ! No. of steps for ion #NSW =0 ! No. of steps for ion # new added lines: AMIX = 0.05 #AMIX = 0.1 BMIX = 0.0001 ! almost zero, but 0 will crash some versions AMIX_MAG = 0.3 #AMIX_MAG = 0.5 BMIX_MAG = 0.0001 ! almost zero, but 0 will crash some versions #LPARD=.TRUE. #NBMOD=-3 #EINT=-1 0 #ISIF=3 #ISYM=1 POTIM = 0.3 #add for NEB IMAGES=4 SPRING=-5.0 LCLIMB=.TRUE. LTANGENTOLD=.FALSE. LDNEB=.FALSE. IBRION=3 #IALGO = 48 ! RMM-DIIS algorithm for electrons ALGO = Fast LREAL =.TRUE. ! evaluate projection operators in real space #NPAR = 1 NSIM = 4 ! blocked algorithm update, four bands a time GGA = PE ! PBE DOS related values LVTOT = .TRUE. 计算收敛之后处理的数据：3 1.171351 -0.650848 1.218286 3.05 1.1930247 -0.677556873901334 1.24581005302295 3.1 1.2146984 -0.704827288335876 1.27010368993822 3.15 1.2363721 -0.732589228234031 1.29116691074581 3.2 1.2580458 -0.760772678526208 1.30899971544572 3.25 1.2797195 -0.789307624142813 1.32360210403796 3.3 1.3013932 -0.818124050014252 1.33497407652251 3.35 1.3230669 -0.847151941070934 1.34311563289939 3.4 1.3447406 -0.876321282243264 1.34802677316859 3.45 1.3664143 -0.90556205846165 1.34970749733011 3.5 1.388088 -0.9348042546565 1.34815780538395 3.55 1.4097617 -0.963977855758219 1.34337769733011 3.6 1.4314354 -0.993012846697216 1.33536717316859 3.65 1.4531091 -1.0218392124039 1.32412623289939 3.7 1.4747828 -1.05038693780867 1.30965487652251 3.75 1.4964565 -1.07858600784194 1.29195310403796 3.8 1.5181302 -1.10636640743411 1.27102091544572 3.85 1.5398039 -1.1336581215156 1.24685831074581 3.9 1.5614776 -1.1603911350168 1.21946528993822 3.95 1.5831513 -1.18649543286814 1.18884185302295 4 1.604825 -1.211901 1.154988 所以，3不是鞍点？CNEB不是可以自动找到能量鞍点的吗？ |
回复此楼» 猜你喜欢
315求调剂
已经有8人回复
-
312求调剂
已经有5人回复
-
一志愿郑州大学材料与化工085600，求调剂
已经有4人回复
-
315分 085602 求调剂
已经有7人回复
-
298分 070300求调剂
已经有5人回复
-
349求调剂
已经有9人回复
-
302求调剂
已经有4人回复
-
化工调剂303分，过四级
已经有10人回复
-
321求调剂
已经有14人回复
-
296材料专硕求调剂
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教关于本机的Dmol计算问题——望高手指点
已经有6人回复
- 加溶剂后计算，不收敛，请高手指点
已经有5人回复
- fluent计算发散，仅多了一个烟气回流的结构，求高手指点。
已经有5人回复
- MS discover 计算出错，请高手指点一下
已经有4人回复
- 新手请教用NEB方法计算energy barrier.
已经有15人回复
- 关于NEB计算的基本问题
已经有7人回复
- vasp用neb计算迁移势垒
已经有11人回复
- NEB计算的的收敛问题
已经有6人回复
- 用gaussian计算丙烷的标准生成热，结果与实验值相差很大，盼高手指导。
已经有7人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【求助】vasp用neb计算过渡态报错，错误提示是什么意思？
已经有9人回复
- 【其他】vasp neb 计算讨论
已经有20人回复
- 【求助成功】NEB计算的一个stdout输出怎么解读?
已经有6人回复
- 【求助成功】怎么样提高NEB计算速度?
已经有19人回复
- 【求助】vasp的NEB计算,输出需要哪些参数?
已经有7人回复
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
2楼2010-09-24 18:02:21
