| 查看: 2343 | 回复: 12 | |||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[交流]
【求助成功】从原胞到表面吸附结构,一套完整的计算,怎么规划优化? 已有5人参与
|
||
|
在各位虫友的帮助启发下,现在VASP总能计算结果基本上都能与文献结果保持数量级上的一致.现在迫切需要提高计算精度,希望可以把偏差控制到0.02eV以下.我把我的计算流程列出来,希望大家提供改善或纠正的意见. 我做计算的流程如下: (1)原胞优化,建模. (2)ISTART=0,ICHARG=2,EDIFF=1E-3,EDIFFG=-0.1,IBRION=2,LREAL=Auto粗结构优化; (3)rm WAVECAR.EDIFF=1E-4,EDIFFG=-0.01,IBRION=2,LREAL=Auto结构优化; (4)静态计算,优化ENCUT,K-points,SIGMA (5)rm WAVECAR,用优化好的参数和结构,做静态自恰计算,得到总能和电荷密度数据. (6)能带等其它性质的计算. (7)数据处理. 问题: (1)FFT的网格优化需要做吗?在哪个阶段做比较合适? (2)对于含有f电子的模型,带U优化结构有必要吗?先不带U优化结构,计算总能和电子结构的时候再把U加上满足精度要求吗? (3)同种材料不同结构的模型,SIGMA要分别优化吗?比如原胞总能计算,需要优化SIGMA;做出表面后SIGMA又要优化一次;做出表面吸附结构后,SIGMA需要再优化一次?提出这样的问题是因为我检查同样的SIGMA设置,不同模型中EENTRO差别较大,只有少数情况是较小的(<0.01) (4)过了一段时间,计算任务积累的越来越多了,各位怎么管理和保存这些结算结果呢? 这些还是入门性的问题,很多可以自己慢慢摸索.但是希望可以和大家交流一下,以节约时间.也希望这个帖子对一些人有帮助. [ Last edited by cenwanglai on 2010-6-17 at 14:46 ] |
回复此楼» 猜你喜欢
085600材料与化工调剂 324分
已经有5人回复
-
344求调剂
已经有6人回复
-
266求调剂
已经有9人回复
-
化学工程321分求调剂
已经有18人回复
-
314求调剂
已经有8人回复
-
08工科 320总分 求调剂
已经有5人回复
-
085601专硕，总分342求调剂，地区不限
已经有5人回复
-
311求调剂
已经有6人回复
-
化工学硕306求调剂
已经有10人回复
-
286求调剂
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 表面吸附体系的结构优化（弛豫）无法收敛 求分析求指导！
已经有5人回复
- 离子吸附结构优化后太奇怪了！？？？
已经有9人回复
- vasp表面和吸附分子的结构分别优化后，如何将它们放到一起
已经有11人回复
- fluent 计算收敛后，看不到结果
已经有4人回复
- 关于在构建表面吸附结构时吸附物的位置问题
已经有5人回复
- EDIFF EDIFFG之间设置关系
已经有12人回复
- vasp结构优化中EDIFF以及KPOINTS设置，以及氢吸附结合能问题！
已经有17人回复
- 文献中通常所说的结构优化是单胞的优化还是对超胞进行的优化呢？
已经有5人回复
- 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
已经有4人回复
- Encut和K点个数之间有没有联系？
已经有6人回复
- 【求助】在用MS搭建表面模型时，怎么看出来撘了几层模型呢？？
已经有15人回复
- 不算完美的2010（我的年终总结和展望）
已经有10人回复
- 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
已经有5人回复
- 【求助】结构优化时如何将原子位置固定在原胞内
已经有9人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助完毕】硫酸根离子的结构优化，偶极修正及势能修正。
已经有20人回复
- 【求助】ENCUT的意思
已经有12人回复
- 【求助完毕】计算吸附能的时候要除以原胞的个数,有什么物理意义呢?
已经有10人回复
7楼2010-05-24 18:28:07
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
|
候博士的指南中给出了SIGMA优化的一个例子.去找来看看.意思是在0.1-0.3之间设置几个SIGMA,分别做静态计算,然后看EENTRO和SIGMA的关系,一般要控制EENTRO较小,比如0.01(这个我不肯定,请高手给予判断)就可以了. 不过很多文献直接给出所有计算SIGMA=0.2不变化,或者算某个特别的值是指定SIGMA=0.05或者0.07这样很小的数值. 但是可以看到,SIGMA对OUTCAR中的三个能量数值是有影响的: freen energy TOTEN energy without entropy energy (sigma->0) 到底选什么能量作为计算能量差计算基准?不知道.因为有文献交代所有能量都插值sigma->0.候博士的指南中好象用的是energy without entropy.但是经常看到用grep "TOTEN" OUTCAR |tail -1 来看能量的. 希望高手可以给个结论. |
9楼2010-05-24 20:43:23
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
2楼2010-05-24 14:19:12
3楼2010-05-24 16:45:00
4楼2010-05-24 17:40:38
5楼2010-05-24 17:56:42
lono75
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6963.1
- 散金: 100
- 红花: 2
- 帖子: 2433
- 在线: 536.4小时
- 虫号: 932458
- 注册: 2009-12-23
- 专业: 凝聚态物性 II ：电子结构
cenwanglai(金币+2):谢谢回复. 2010-05-24 20:26:20
|
1meV/atom [ Last edited by lono75 on 2010-5-24 at 20:13 ] |
6楼2010-05-24 18:19:45
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
8楼2010-05-24 20:23:49
lono75
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6963.1
- 散金: 100
- 红花: 2
- 帖子: 2433
- 在线: 536.4小时
- 虫号: 932458
- 注册: 2009-12-23
- 专业: 凝聚态物性 II ：电子结构
10楼2010-05-24 21:08:58
50