计算化学公社

标题: 关于cp2k计算Fe表面磁矩设置 [打印本页]

---

作者
Author: 莫雨萱    时间: 2022-4-9 17:29
标题: 关于cp2k计算Fe表面磁矩设置
用cp2k计算Fe表面，初猜磁矩设为3，但为什么mulliken文件中显示磁矩几乎为零呀

---

作者
Author: 丁越    时间: 2022-4-9 18:03


UKS加了么？
参考http://bbs.keinsci.com/thread-22642-1-1.html

---

作者
Author: 莫雨萱    时间: 2022-4-9 19:54

丁越 发表于 2022-4-9 18:03
UKS加了么？
参考http://bbs.keinsci.com/thread-22642-1-1.html

UKS加了，然后测试体系的多重度，结果发现Fe磁矩几乎为零。老师我还有一个问题，关于多重度的测试，对于一个Fe来讲多重度是5，那如果我的体系存在23个Fe原子，那多重度测试上限是不是应该为115呀？

---

作者
Author: 丁越    时间: 2022-4-9 22:34
看12L的回复

---

作者
Author: 莫雨萱    时间: 2022-4-10 17:32

丁越 发表于 2022-4-9 22:34
看12L的回复

好的，谢谢

---

作者
Author: sobereva    时间: 2022-4-11 03:29
你的体系里有Pb，明显不光是Fe表面

除了KIND里设MAGNETIZATION，还得设整体的自旋多重度。
对于Fe表面，得开smearing，此时默认会根据自旋轨道能级自动判断轨道占据方式，而不严格满足整体的自旋多重度设置，此时也无所谓“多重度测试上限”。

---

作者
Author: 莫雨萱    时间: 2022-4-11 09:56

sobereva 发表于 2022-4-11 03:29
你的体系里有Pb，明显不光是Fe表面

除了KIND里设MAGNETIZATION，还得设整体的自旋多重度。

感谢老师的回答，我在开始计算时打开了smearing，没有设置自旋多重都，但是最后结果与文献结果相差很大，文献是用vasp计算的，几个能量差值（吸附能）为-0.75，我用cp2k计算出来，换算单位后，能量直接是-2.0--4.0左右，所以在想是不是因为没有设置多重都的原因，我把我input文件放在下面，麻烦老师帮我看看，是哪里的问题

---

作者
Author: 莫雨萱    时间: 2022-4-11 09:57
老师，这是我的input文件

---

作者
Author: jiejie    时间: 2022-7-20 12:26
请问这列Net charge是什么意思，是指电荷的得失，负的代表失去电荷吗

---

作者
Author: nusiew    时间: 2022-10-5 14:48

莫雨萱 发表于 2022-4-11 09:57
老师，这是我的input文件

请问楼主的问题解决了么

---

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

Powered by Discuz! X3.3