模型中含有O，必须进行自选极化计算吗？

kaikai · 发表于 Post on 2024-5-19 09:31:41

各位老师好，使用Dmol3计算Al2O3模型，表面存在一定数量的羟基，在计算时，模型中含有O，必须进行自选极化计算吗？Dmol3好像不能进行自旋极化吧，勾选Spin unrestricted是否能够起到类似自旋极化的作用？
各位老师，目前不是很清楚什么情况下要用自旋极化，希望各位老师给予解答，谢谢

sobereva · 发表于 Post on 2024-5-21 23:33:48

我不用Dmol3，如果本身是闭壳层体系，当限制性闭壳层计算和非限制性开壳层结果有差别，属于程序自己的bug。至少CP2K、Gaussian等各种正经的主流程序都没这个问题，原理上也绝对不会有这种问题。

Observer-Aya · 发表于 Post on 2024-5-21 15:46:35

虽然sob老师说的没错，你的属于开壳层体系理应不加，但按照B站上孙成华老师的建议，如果是计算闭壳层体系的单原子能，开与不开自旋极化依然会个位数ev的差距，所以建议加上，算是软件小经验一类的。

kaikai · 发表于 Post on 2024-5-19 14:38:38

sobereva 发表于 2024-5-19 11:56
你这是典型的闭壳层体系，没自旋极化的事

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.c ...

谢谢老师的解答

sobereva · 发表于 Post on 2024-5-19 11:56:08

你这是典型的闭壳层体系，没自旋极化的事

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt：

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CASTEP/Dmol3/MS] 模型中含有O，必须进行自选极化计算吗？