| 查看: 2848 | 回复: 9 | ||
Litter_horsr银虫 (小有名气)
|
[求助]
CLNEB计算过渡态的问题 已有1人参与
|
|
参照说明书，我做了在氧化镓表面吸附H2的过程过渡态 由于已经有文献给出H2吸附在表面时，它的稳定结构，所以就没有再去过多的考虑可能的吸附情况，直接采用最稳定结构作为终态进行优化。初始态是H2分子物理吸附在表面，距离表面高度为0.7A 优化和Clneb时候的IBCAR设置: 优化时： SYSTEM=GaO ENCUT=480 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.02 NSW=500 IBRION=2 ISIF=2 EDIFF=1E-5 EDIFFG=-0.01 NELM=180 LORBIT=0 POTIM=0.08 PREC=Accurate 初态和终态的结构优化之后，插入四个镜像 。00-05 SYSTEM=GaO ENCUT=480 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.02 NSW=500 ISIF=2 EDIFF=1E-5 EDIFFG=-0.1 NELM=180 LORBIT=0 PREC=Accurate #for Clneb IBRION=3 POTIM=0 IOPT=0 ICHAIN=0 LCLIMB=.TURE. SPRING=-5 IMAGES=4 结果能够运行，但是一直优化不完，最后一步的时候提示： 500 F= -.24183396E+03 E0= -.24183396E+03 d E =0.334694E-09 quench: g(F)= 0.000E+00 g(S)= 0.000E+00 dE (1.order)= 0.000E+00 writing wavefunctions real 28323.73 user 450211.62 sys 342.69 Sat Sep 5 17:00:42 CST 2015 并且查看原子间力的时候，有如下（nebef.pl脚本） 0 0.000000 -241.833900 0.000000 1 0.000000 -240.388700 1.445200 2 0.000000 -239.776100 2.057800 3 0.000000 -240.928700 0.905200 力全部为0…… 请问： 1、过渡态计算为什么不收敛？ 2、利用CLNEB主要采用力作为收敛条件，但是这个力全部为0.怎么解决？ |
回复此楼» 猜你喜欢
299求调剂
已经有7人回复
-
0703化学调剂，求导师收
已经有4人回复
-
300求调剂，材料科学英一数二
已经有3人回复
-
材料学硕，求调剂
已经有6人回复
-
307求调剂
已经有6人回复
-
材料专硕331求调剂
已经有3人回复
-
求调剂一志愿武汉理工大学材料工程（085601）
已经有7人回复
-
一志愿C9材料与化工专业总分300求调剂
已经有4人回复
-
321求调剂
已经有4人回复
-
307求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
阿克飞顿
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1426.6
- 散金: 37
- 红花: 2
- 帖子: 572
- 在线: 169.6小时
- 虫号: 2469736
- 注册: 2013-05-18
- 性别: GG
- 专业: 光学和光电子材料
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Litter_horsr: 金币+10, ★有帮助 2015-09-06 09:47:19
感谢参与，应助指数 +1
Litter_horsr: 金币+10, ★有帮助 2015-09-06 09:47:19
|
任务提交命令有问题，你是在超算中心算的还是自己的服务器上，VTSY脚本编译进去了吗 发自小木虫Android客户端 |
2楼2015-09-06 09:30:41
阿克飞顿
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1426.6
- 散金: 37
- 红花: 2
- 帖子: 572
- 在线: 169.6小时
- 虫号: 2469736
- 注册: 2013-05-18
- 性别: GG
- 专业: 光学和光电子材料
3楼2015-09-06 09:30:53
阿克飞顿
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1426.6
- 散金: 37
- 红花: 2
- 帖子: 572
- 在线: 169.6小时
- 虫号: 2469736
- 注册: 2013-05-18
- 性别: GG
- 专业: 光学和光电子材料
【答案】应助回帖
★ ★
cenwanglai: 金币+2, 谢谢~ 2015-09-06 10:03:02
cenwanglai: 金币+2, 谢谢~ 2015-09-06 10:03:02
|
这是我在超算中心的提交方法，你如果已经编译VTST脚本了，参照我的方法做一下改动试试 1）vasp和vasp-vtsts早加入到环境变量设置文件.bash_profile中了 2）提交vasp的计算脚本： bsub -J job_name -q q_x86_kanglh -n 48 -o log.txt?~/apps/vasp/5.3.2/x86_64/bin/vasp 3）提交vasp-vtst的计算脚本： bsub -J job_name -q q_x86_kanglh -n 48 -o log.txt ~/apps/vasp-vtst/5.3.2/x86_64/bin/vasp-vtst 发自小木虫Android客户端 |
4楼2015-09-06 09:33:46
Litter_horsr
银虫 (小有名气)
- 1ST强帖: 1
- 应助: 12 (小学生)
- 金币: 598.1
- 红花: 12
- 帖子: 162
- 在线: 100小时
- 虫号: 3160497
- 注册: 2014-04-24
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-09-06 09:44:41
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
|
1. VASP自身并不带climbing功能，需要把vtstcode加进vasp编译。 2. IOPT =0 是什么意思？如果加入了vtstcode代码，建议用更高效的优化器（调整IOPT=？）。请搜索neb tst 3. 反应路径计算，失败的几率很大，因此需要注意优化策略。最好不要一上来就高精度。 EDIFF=1E-5 ！作死 EDIFFG=-0.01 ！一上来就这样的精度，作死 NELM=180 ！正常情况下，离子步如果40步以内不收敛，说明结构及其不合理。默认即可。 LORBIT=0 POTIM=0.08 ！请用CNEB的优化器，POTIM = 0, IBRION = 3. PREC=Accurate !这个，看你的系统精度。如果对路径和计算技术本身不是很熟悉，可以考虑先Normal. 搜索出个大概来了，再Accurate. LCLIMB = .TRUE. !一开始的时候最好关掉。先搜索出路径的整体轮廓来（EDIFFG=-0.03~-0.05)，然后再打开。 4. 反应路径计算费神劳力，需要较多的低精度测试计算来加深对反应过程的理解和把握，最后在低精度计算筛选出来的少数可信的路径上提高精度才是有智慧的表现。 |
6楼2015-09-06 10:00:25
阿克飞顿
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1426.6
- 散金: 37
- 红花: 2
- 帖子: 572
- 在线: 169.6小时
- 虫号: 2469736
- 注册: 2013-05-18
- 性别: GG
- 专业: 光学和光电子材料
7楼2015-09-06 10:24:42
Litter_horsr
银虫 (小有名气)
- 1ST强帖: 1
- 应助: 12 (小学生)
- 金币: 598.1
- 红花: 12
- 帖子: 162
- 在线: 100小时
- 虫号: 3160497
- 注册: 2014-04-24
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2015-09-08 10:03:12
jackjiejl
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2490
- 散金: 127
- 帖子: 128
- 在线: 310小时
- 虫号: 2719699
- 注册: 2013-10-12
- 性别: GG
- 专业: 催化化学
9楼2016-05-06 14:43:18
10楼2016-05-10 09:46:17
