24小时热门版块排行榜    

查看: 4601  |  回复: 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娜 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见