24小时热门版块排行榜    

查看: 2801  |  回复: 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的回帖

lono75

铁杆木虫 (著名写手)

cenwanglai(金币+2):谢谢回复. 2010-05-24 20:26:20
引用回帖:
Originally posted by kgdu at 2010-05-24 17:56:42:
对于每个给定的ENCUT,计算收敛后,比较不同ENCUT下的系统能量,看其差值,小于某一值时,就可以取那个较小ENCUT。这个差值一般是多少呢?

1meV/atom

[ Last edited by lono75 on 2010-5-24 at 20:13 ]
由天命,成人事!
6楼2010-05-24 18:19:45
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 cHIDTWCTtLlG 2026-08-19 5/250 2026-08-19 17:30 by oCMhZAsgNFoO
[基金申请] filecode,4个jtjc了 +10 ziyangfang 2026-08-19 11/550 2026-08-19 17:19 by wdc28
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 117/5850 2026-08-19 17:12 by 医学老男孩
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见