计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1954|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] Fe2P的slab优化后结构变化很大,阴离子向外跑,阳离子向内跑

[复制链接 Copy URL]

38

帖子

0

威望

554

eV
积分
592

Level 4 (黑子)

构建了Fe2P的slab,如果不扩胞直接进行结构优化可以正常收敛并且结构也没问题,但是扩胞成2x2的超胞后再进行结构优化,slab结构变化都非常大,P向层外跑,Fe往层内压缩的趋势(尝试固定和不固定部分层数、更改暴露表面均是该结果)。
INCAR在附件中,扩胞后K点只算111,赝势是正常的赝势,用的vasp_gam版去运行的。
想请教一下各位问题出在哪里呢?

202211102234051981..png (83.26 KB, 下载次数 Times of downloads: 9)

优化但未收敛时的结构

优化但未收敛时的结构

202211102233234488..png (89.45 KB, 下载次数 Times of downloads: 10)

初始结构

初始结构

CONTCAR

11.9 KB, 下载次数 Times of downloads: 2

POSCAR

6.37 KB, 下载次数 Times of downloads: 4

INCAR

2.38 KB, 下载次数 Times of downloads: 4

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

2#
发表于 Post on 2022-11-11 09:33:11 | 只看该作者 Only view this author
以当前POSCAR里的盒子尺寸来说,只考虑gamma点不合适,至少2*2*2

slab太薄
并且考虑实际可能出现的表面重构的可能性
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

38

帖子

0

威望

554

eV
积分
592

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-11-11 10:02:55 | 只看该作者 Only view this author
sobereva 发表于 2022-11-11 09:33
以当前POSCAR里的盒子尺寸来说,只考虑gamma点不合适,至少2*2*2

slab太薄

谢谢sob大佬的回复!如此简化确实是出于无奈,之前slab是4层,依旧有P逸出slab的情况。
我有一点比较好奇,扩胞前后真空层厚度都是一样的,为什么扩胞后才出现这么明显的slab变化呢?是因为自由度改变了吗?

255

帖子

1

威望

2682

eV
积分
2957

Level 5 (御坂)

4#
发表于 Post on 2022-11-11 10:33:03 | 只看该作者 Only view this author
优化slab模型的时候至少也要固定最下面三层的原子,哪能不固定就优化呢。主要是为了尽可能的模拟真实的晶体内部结构,而对表面的原子进行弛豫。还要在当前的结构下对k点、截断能进行测试,最好也设置一下磁性的初猜。表面重构也很正常,不过要你根据经验或者查找文献看看重构后的表面长啥样。
我唯一知道的就是我一无所知,但我是化学小迷弟

38

帖子

0

威望

554

eV
积分
592

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-11-11 11:01:47 | 只看该作者 Only view this author
ChemG 发表于 2022-11-11 10:33
优化slab模型的时候至少也要固定最下面三层的原子,哪能不固定就优化呢。主要是为了尽可能的模拟真实的晶体 ...

谢谢大佬的回复!固定3层放开4层的话,以我那个小破机器估计得跑一年 :-(,K点、截断能以及色散校正等对于表面重构影响大不大呢?我这个重构太夸张了,继续结构优化,一些P跑到距离slab7-8A的位置上了感觉实在不合理。

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

6#
发表于 Post on 2022-11-11 11:14:58 | 只看该作者 Only view this author
wang5360848 发表于 2022-11-11 10:02
谢谢sob大佬的回复!如此简化确实是出于无奈,之前slab是4层,依旧有P逸出slab的情况。
我有一点比较好 ...

没准之前的结果本来就是不合理的,过小的晶胞造成的周期性约束抑制了理应有的表面重构
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

7#
发表于 Post on 2022-11-11 11:16:12 | 只看该作者 Only view this author
wang5360848 发表于 2022-11-11 11:01
谢谢大佬的回复!固定3层放开4层的话,以我那个小破机器估计得跑一年 :-(,K点、截断能以及色散校正等对 ...

扩胞后,k点2*2*2够了
色散校正跟这个没直接联系
截断能做个收敛性测试便知
嫌VASP慢不如用又免费又极快的CP2K
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

38

帖子

0

威望

554

eV
积分
592

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-11-11 11:27:32 | 只看该作者 Only view this author
sobereva 发表于 2022-11-11 11:16
扩胞后,k点2*2*2够了
色散校正跟这个没直接联系
截断能做个收敛性测试便知

好的,谢谢sob大佬!最近是想转CP2K来着,只是刚学会vasp[捂脸]

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 22:36 , Processed in 0.168238 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list