24小时热门版块排行榜    

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

前田庆次

木虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+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的回帖
查看全部 10 个回答

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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
请填处理意见