24小时热门版块排行榜    

查看: 4595  |  回复: 9

xiewei娜

银虫 (小有名气)

[交流] 用vasp优化晶胞的几点疑问

在说我的问题前,我先讲讲我是怎么优化晶胞的:首先我的incar和kpoints中参数设置是来自于文献的,所以我从一开始就设置了文献给出的k点,而没有一步步提高k点。而ENCUT的值同一个体系,有的作者用的是600eV,有的用的是800eV,不知道用哪个好?我的优化步骤是按源资培训里面介绍的来做的,i:优化晶胞参数  ISIF=3,固定原子位置(即在POSCAR文件中原子坐标后加F)
ii:原子弛豫  ISIF=2 ,不固定原子位置(默认是不加F的,即T)
iii:用vasp算一般结构都是对称性降为P1,最后算完将CONTCAR用ms打开,通过find sysmery--impose sysmetry恢复对称性
以下是我的问题:
①我的这个优化过程可行吗
②第一步的时候我不确定到什么时候停止,就把contcar变为poscar重复第一步的计算,最近有一个我重复了20次,最后晶胞参数和总能量都不再变化了,然后我就进行原子弛豫,发现晶胞参数和总能量还是不变,我就不知道我这样做对不对,感觉重复太多次了。应该有一个什么标准来停止晶胞参数的优化这一步呢?我知道的就是看最护两步的总能量和晶胞参数基本一样,但有没有一个什么标准呢?跟incar中设置的EDIFF和EDIFFG有关吗?   
③excel中我会列入external pressure是看到文献的一句话 (见图)  ,不知是不是指的是最后优化总压力小于+-0.5KB,这个需不需要考虑到优化是否收敛中来?
④关于poscar也就是建模的问题,文献中有两种数据,一种是实验测得的晶胞参数和原子坐标,一种是别人理论计算优化好的晶胞参数和原子坐标,使用哪种数据建模更合适呢?

用vasp优化晶胞的几点疑问
DN}RHRJ@JIY6)14~{[AFJUF.png


用vasp优化晶胞的几点疑问-1
RKRX]Q1Q2F`%IC}0{H{@`NB.png
回复此楼

» 猜你喜欢

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

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

xrdcvw73

新虫 (正式写手)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiewei娜: 金币+5, 谢谢! 2015-08-24 21:09:58
楼主问题太详细了
从优化的角度来说,楼主的思路是对的
一下是我的看法,个人观点:
首先 优化晶胞参数,原子位置(这个一起弄 ISIF=3 对单胞优化)——建超胞(ISIF=2)只优化位置——设计性质
不过ENCUT 建议优化前还是测试哈
KPOINT可以参考文献或者经验值
2楼2015-08-24 20:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
2楼: Originally posted by xrdcvw73 at 2015-08-24 20:10:11
楼主问题太详细了
从优化的角度来说,楼主的思路是对的
一下是我的看法,个人观点:
首先 优化晶胞参数,原子位置(这个一起弄 ISIF=3 对单胞优化)——建超胞(ISIF=2)只优化位置——设计性质
不过ENCUT 建议 ...

意思是优化晶胞参数和原子弛豫合在一起,在单胞优化好的基础上建超胞的话就只进行原子弛豫?
还有看有没有优化好事看什么标准呢?
3楼2015-08-24 21:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mlanqiang

木虫之王 (文学泰斗)

蓝博士

blessing
蓝精灵
4楼2015-08-24 21:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrdcvw73

新虫 (正式写手)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiewei娜: 金币+5, 好的,谢谢哈 2015-08-25 11:06:53
引用回帖:
3楼: Originally posted by xiewei娜 at 2015-08-24 21:08:56
意思是优化晶胞参数和原子弛豫合在一起,在单胞优化好的基础上建超胞的话就只进行原子弛豫?
还有看有没有优化好事看什么标准呢?...

你把CONTCAR复制为POSCAR 然后重新拿之前的INCAR优化一哈,基本上一步就好了,就说明打到所需的精度了
经验:如果优化在几步或者十几步结束,基本上不用上面的过程,可以认为CONTCAR就是达到精度的结构...如果不放心,也可以试试.....
5楼2015-08-24 21:40:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiewei娜: 金币+5 2015-08-25 11:04:53
在说我的问题前,我先讲讲我是怎么优化晶胞的:首先我的incar和kpoints中参数设置是来自于文献的,所以我从一开始就设置了文献给出的k点,而没有一步步提高k点。而ENCUT的值同一个体系,有的作者用的是600eV,有的用的是800eV,不知道用哪个好?
"那肯定是800eV好!但是这东西要测试,如果600达到你的精度要求了,就没必要800了……"
我的优化步骤是按源资培训里面介绍的来做的,i:优化晶胞参数  ISIF=3,固定原子位置(即在POSCAR文件中原子坐标后加F)
ii:原子弛豫  ISIF=2 ,不固定原子位置(默认是不加F的,即T)
iii:用vasp算一般结构都是对称性降为P1,最后算完将CONTCAR用ms打开,通过find sysmery--impose sysmetry恢复对称性
以下是我的问题:
①我的这个优化过程可行吗②第一步的时候我不确定到什么时候停止,就把contcar变为poscar重复第一步的计算,最近有一个我重复了20次,最后晶胞参数和总能量都不再变化了,然后我就进行原子弛豫,发现晶胞参数和总能量还是不变,我就不知道我这样做对不对,感觉重复太多次了。应该有一个什么标准来停止晶胞参数的优化这一步呢?我知道的就是看最护两步的总能量和晶胞参数基本一样,但有没有一个什么标准呢?跟incar中设置的EDIFF和EDIFFG有关吗?   
“姐,你这步骤真详细,但是我想,没有哪个方法,是通用的方法,对于不同的体系,是有不同的方法的,当然,你这个步骤,应该没有问题,但是真麻烦”
③excel中我会列入external pressure是看到文献的一句话 (见图)  ,不知是不是指的是最后优化总压力小于+-0.5KB,这个需不需要考虑到优化是否收敛中来?
“这个不懂……我都是利用EDIFFG控制”
④关于poscar也就是建模的问题,文献中有两种数据,一种是实验测得的晶胞参数和原子坐标,一种是别人理论计算优化好的晶胞参数和原子坐标,使用哪种数据建模更合适呢?
“如果你没有文献参考,当然是用实验的了。但是有文献参考,如果方法和你一样,那当然是用文献的了,但是前提是,文献一定要可靠,比如PRB一类的,或者是知名学者的文献”
以上个人愚见,仅供参考
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
6楼2015-08-25 08:45:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福,,,,
7楼2015-08-25 08:59:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
7楼: Originally posted by huzhang01 at 2015-08-25 08:59:50
祝福,,,,

呵呵
8楼2015-08-25 11:02:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
6楼: Originally posted by 前田庆次 at 2015-08-25 08:45:19
在说我的问题前,我先讲讲我是怎么优化晶胞的:首先我的incar和kpoints中参数设置是来自于文献的,所以我从一开始就设置了文献给出的k点,而没有一步步提高k点。而ENCUT的值同一个体系,有的作者用的是600eV,有的用 ...

如果我打算自己测试encut和kpoint,那还可以参考别人优化好的结构建模吗,会不会矛盾
9楼2015-08-25 11:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g_wcai

木虫 (正式写手)

开心

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiewei娜: 金币+5, 嗯嗯,还是测试一下吧 2015-08-25 19:30:30
ENCUT参数要看POTCAR中各个原子ENMAX  =  208.733和ENMIN  =  156.550 eV,一般取最大值的1.5倍!当然也可以测试一下!!!
哈哈,笑口常开!
10楼2015-08-25 11:18:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiewei娜 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们,祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见