24小时热门版块排行榜    

查看: 4502  |  回复: 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娜 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见