24小时热门版块排行榜    

查看: 4596  |  回复: 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娜 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 8/400 2026-03-16 13:42 by Dream007008
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 333求调剂 +3 球球古力 2026-03-09 3/150 2026-03-14 01:57 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见