| 查看: 1873 | 回复: 3 | |||
wujianxinfa铜虫 (初入文坛)
|
[求助]
vasp-b3lyp优化晶体表面缺陷结构电子步不收敛
|
|
TiO2 110表面4x2超胞点缺陷，想用b3lyp优化结构。 文献上是用crystal03做的，但我不会用。 尝试用vasp，但LDA可以收敛，b3lyp不能收敛（电子步），改了很多参数 IALGO = 38 -----> ALGO=Fast SIGMA=0.10-----> SIGMA=0.50 EDIFFG=1e-3 ----> EDIFFG=1e-2 仍然不能收敛。 貌似电子步的步长太大.输出文件会在后面附上。 那位高人能帮我解决这个问题？ 用什么程序包都可以，求一个完整的可靠的输入文件。 成功搞定后，可以拿走我全部的金币（虽然只有170，很少。谨表谢意） System = rutile110 ISTART=0 ICHARG=2 ISMEAR=0 NSW=100 IBRION=2 ISIF=2 ISYM=1 PREC=Medium ENCUT= 400 eV EDIFF=1E-3 EDIFFG=1E-2 SIGMA=0.10 ALGO = Fast LREAL=Auto LPLANE=.FALSE NSIM=1 LCHARG=F LWAVE=F POTIM=0.1 PRECFOCK=FAST LHFCALC=.TRUE. GGA=B3 AEXX=0.2 AGGAX=0.72 AGGAC=0.81 ALDAC=0.19 从第十步就这样了了，一直不能收敛。 RMM: 44 0.279821908089E+04 0.43378E-01 -0.12687E+02 2202 0.107E+01 0.263E-01 RMM: 45 0.279809779209E+04 -0.12129E+00 -0.12665E+02 2225 0.106E+01 0.356E-01 RMM: 46 0.279835902046E+04 0.26123E+00 -0.12476E+02 2208 0.107E+01 0.273E-01 RMM: 47 0.279833161472E+04 -0.27406E-01 -0.12697E+02 2214 0.106E+01 0.225E-01 RMM: 48 0.279810171254E+04 -0.22990E+00 -0.12634E+02 2222 0.107E+01 0.335E-01 RMM: 49 0.279813394587E+04 0.32233E-01 -0.12680E+02 2209 0.106E+01 0.594E-01 |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有8人回复
-
关于filecode
已经有4人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
2楼2017-01-07 23:45:50
走在路上i
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 95
- 散金: 30
- 帖子: 71
- 在线: 12小时
- 虫号: 10074762
- 注册: 2018-08-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
解决了吗？交流一下吧
3楼2019-06-08 19:33:42
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3563.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2019-06-20 04:56:51
100