计算化学公社

标题: vasp 程序如何处理非周期性带电荷体系 [打印本页]

作者
Author: zhangshu 时间: 2023-9-8 13:19
标题: vasp 程序如何处理非周期性带电荷体系
我研究的体系为[CoW12O40]6-，体系带6个负电荷（Co为二价），根据POSCAR和POTCAR文件计算可得，[CoW12O40]中应该含有417个外层电子，
1. 我想研究[CoW12O40]6-体系，是不是只要设置NELECT=423就行了？
2. 如果我想研究[CoW12O40]6-体系得一个电子，是不是NELECT要设置为424？
3. 如果设置NELECT=424，ISPIN是不是要设置为1，还是说处理含过渡金属的体系，必须要考虑磁性，ISPIN必须设置为2？

作者
Author: 卡开发发 时间: 2023-9-9 05:04
对于NELECT的设置来说，你的1、2没问题。3而言，如果体系有磁性那么就应当设置，仅仅是ISPIN=2是不够的，要设置的初始磁矩要比实际磁矩稍大一点但比其更高的磁矩要小（因为存在高低自旋态都存在的情况）。

另外带电荷的体系需要考虑LMONO，如果分子有偶极还应该使用LDIPOL配合IDIPOL=4，并且要把盒子尺寸测试收敛。带这么多的电荷体系不一定能稳定，也许实际情况实在溶剂当中或其他环境，你可能需要根据你自己的实际情况进行调整。

作者
Author: zhangshu 时间: 2023-9-9 11:21

卡开发发发表于 2023-9-9 05:04
对于NELECT的设置来说，你的1、2没问题。3而言，如果体系有磁性那么就应当设置，仅仅是ISPIN=2是不够的，要 ...

感谢老师您百忙之中的回复，我还想请教一下：
对于多重度不唯一的体系，ispin应该设置为2；但对于闭壳层体系，ispin是否可以直接设置成1？有没有什么方法是直接可以判断材料磁性的？（ps：之前用的都是高斯程序，不太懂磁性与多重度之间的联系）

作者
Author: 卡开发发 时间: 2023-9-9 19:34

zhangshu 发表于 2023-9-9 11:21
感谢老师您百忙之中的回复，我还想请教一下：
对于多重度不唯一的体系，ispin应该设置为2；但对于闭壳层 ...

如果确实是闭壳层，按1设置没问题。材料磁性稍微简单一些的可能通过结构化学或者固体物理知识判断出来，但太复杂的也没啥办法，可以去一些晶体数据库或者文献去查，如果都没办法确定，但尝试在有限次的可以试试，这些要求都不满足就很棘手了。

作者
Author: zhangshu 时间: 2023-9-9 20:42
谢谢老师，我准备把ispin=1和=2的都算一下，看看哪个能量低一点。感谢老师您百忙之中的解答

作者
Author: xexlalalan 时间: 2023-9-10 06:59
本帖最后由 xexlalalan 于 2023-9-10 07:01 编辑

zhangshu 发表于 2023-9-9 11:21
感谢老师您百忙之中的回复，我还想请教一下：
对于多重度不唯一的体系，ispin应该设置为2；但对于闭壳层 ...

磁性计算确实棘手。即使是单重态，材料如果是反铁磁性也得用上ispin=2。反铁磁性可以理解成某种自旋极化单重态，一个单胞里6个钴，3个自旋朝上，3个自旋朝下

材料的磁性很难根据经验判断他的磁性一定是什么样，而且磁性复杂的材料也不少，都是超出人们经验预测范围的，要不然人们也不会累死累活的去打中子衍射。如果这个材料有实验报导过，就根据他们中子衍射的结果去算。如果没有的话，而且你们组不长单晶不打衍射，就把这六个钴按照可能的反铁磁性的排列组合都算一遍，看看哪个能量最低。

作者
Author: xexlalalan 时间: 2023-9-10 07:00

xexlalalan 发表于 2023-9-10 06:59
磁性计算确实棘手。即使是闭壳层，材料如果是反铁磁性也得用上ispin=2。反铁磁性可以理解成某种自旋极化 ...

才发现你这里面还有钨，如果他们也有磁性的话，排列组合会更加复杂。

祝好运

作者
Author: zhangshu 时间: 2023-9-11 14:37

xexlalalan 发表于 2023-9-10 07:00
才发现你这里面还有钨，如果他们也有磁性的话，排列组合会更加复杂。

祝好运

好的好的，感谢您的回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)