24小时热门版块排行榜    

查看: 2347  |  回复: 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 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
Originally posted by kgdu at 2010-05-24 16:45:00:


请问优化ENCUT,K-point时,一般能量变化小于多少时就可以了,这个值是不是还体系原子的个数有关?还请问一下sigma怎么优化?谢谢了!!

候博士的指南中给出了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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

2楼2010-05-24 14:19:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)

cenwanglai(金币+2):谢谢参与讨论 2010-06-17 14:44:56
引用回帖:
Originally posted by cenwanglai at 2010-05-23 09:48:51:
在各位虫友的帮助启发下,现在VASP总能计算结果基本上都能与文献结果保持数量级上的一致.现在迫切需要提高计算精度,希望可以把偏差控制到0.02eV以下.我把我的计算流程列出来,希望大家提供改善或纠正的意见.

我做 ...

请问优化ENCUT,K-point时,一般能量变化小于多少时就可以了,这个值是不是还体系原子的个数有关?还请问一下sigma怎么优化?谢谢了!!
3楼2010-05-24 16:45:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rdw00

木虫 (著名写手)

cenwanglai(金币+2):谢谢回答. 2010-05-24 20:25:27
引用回帖:
Originally posted by kgdu at 2010-05-24 16:45:00:


请问优化ENCUT,K-point时,一般能量变化小于多少时就可以了,这个值是不是还体系原子的个数有关?还请问一下sigma怎么优化?谢谢了!!

不能看能量小于多少,是要求收敛。取最佳收敛点就可以
4楼2010-05-24 17:40:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天大材料与化工(085600)总分338 +3 蔡大美女 2026-03-13 3/150 2026-03-18 15:53 by jhhcooi
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +7 rare12345 2026-03-18 7/350 2026-03-18 14:31 by laoshidan
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +6 Liwangman 2026-03-15 6/300 2026-03-18 13:21 by 尽舜尧1
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见