计算化学公社

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

[CASTEP/Dmol3/MS] 求助关于收敛性测试时k点、自旋极化参数、空带设置的一些疑问

[复制链接 Copy URL]

1

帖子

0

威望

21

eV
积分
22

Level 1 能力者

各位老师好,最近初学第一性原理计算,看了一些收敛性测试相关的帖子,目前有些疑问望解答:
1. 对于晶格参数a, b, c之间相差很大的晶体而言,k点要怎么设置才能做到精度和计算时间之间的平衡呢?我计算的晶胞大小为3.236×3.236×25.15埃,α=β=90°,γ=120°,如果要满足actual spacing数值几乎相等的话就需要非常密的k点设置,得从9×9×1开始。如若再往上增加k点精度(比如到36×36×4)以进行截断能的测试计算时间会非常长。
2. 对于六方晶系,有的文章中是以例如3×3×1、4×4×2、5×5×3这样每次对k值增加1的方式进行收敛性测试的,但是那些文章中的晶格参数a, b, c之间相差很小,不知道我这种情况是否适用?
3. 在做收敛性测试计算单点能时是否要设置自旋极化?还是说在确定了截断能以及k点后,做结构优化的时候再设置呢?我计算的材料叫铝酸镁钪(ScAlMgO4)里边有钪的话是要考虑自旋极化的吧?
4. 过去曾见到老师们有提到“初猜”,请问“初猜”是否就是初始自旋设置呢?如果是的话这个要怎么设置呢?是化学价的一半吗?
5. 之前对这个晶体做收敛性测试的时候发现自洽循环总是不收敛,就增加了10%的空带,展宽smearing设置为0.05。这样一来SCF就收敛了,请问这样做对DOS图什么的有影响吗?
6. 还有一点关于收敛标准的疑问,不论是k点还是截断能的收敛性测试,只要单点能计算前后相差在1mev每个原子就能认定为收敛吗?我一般是以50ev作为间隔增加截断能的。
望诸位老师不吝赐教,本人本科并非材料、数理相关方向,相关知识较为薄弱,如果犯了什么愚蠢的错误还请斧正。

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 16:28 , Processed in 3.910079 second(s), 29 queries , Gzip On.

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