计算化学公社

标题: 求助；关于铁原子计算级别的疑问 [打印本页]

作者
Author: Lisy 时间: 2026-1-29 18:19
标题: 求助；关于铁原子计算级别的疑问
我想计算一个锂氧化物对一个铁原子的吸附能，请问应该采用什么计算级别来结构优化比较好呢，用把b3lyp/def2svp可以吗

作者
Author: UW_0728. 时间: 2026-1-29 18:39
是孤立体系还是周期性材料？
前者推荐优化在Gaussian里用PBE0-D3(BJ)或者TPSSh-D3(BJ)，铁用SDD，其他用def2SVP或者6-31G(d)（想要更精确结果可以6-311G(d)或TZVP），单点在ORCA里用RI-PWPB95-D4/def2-TZVPP
后者推荐CP2K下优化用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH，单点用BEEF-vdW/TZV2P-MOLOPT-GTH（对Fe用带SR的，因为不带SR的对第四周期及以后的元素无定义）（RPBE或revPBE结合DFT-D3(BJ)也可以）

作者
Author: wzkchem5 时间: 2026-1-29 19:34
此外注意自旋多重度的恰当选取（简单地选最低的自旋多重度很可能是错的），以及检验波函数稳定性

作者
Author: Lisy 时间: 2026-1-29 19:46

UW_0728. 发表于 2026-1-29 18:39
是孤立体系还是周期性材料？
前者推荐优化在Gaussian里用PBE0-D3(BJ)或者TPSSh-D3(BJ)，铁用SDD，其他用de ...

好的，谢谢！

作者
Author: Lisy 时间: 2026-1-29 19:46

wzkchem5 发表于 2026-1-29 19:34
此外注意自旋多重度的恰当选取（简单地选最低的自旋多重度很可能是错的），以及检验波函数稳定性

好的，谢谢您，请问一个铁原子的话自旋多重度设置5应该没错吧

作者
Author: 122332 时间: 2026-1-29 19:55
应该没什么问题，你这里又不涉及到强场配体和弱场配体的事，实在不放心，对可能的自旋多重度都算一遍单点呗

作者
Author: sobereva 时间: 2026-2-2 06:50

Lisy 发表于 2026-1-29 19:46
好的，谢谢您，请问一个铁原子的话自旋多重度设置5应该没错吧

孤立的Fe原子就该这么设
要注意程序不一定能正确收敛到Fe的基态电子组态4s2 3d6，应当检查轨道占据情况。如果发现收敛到了错误的电子组态，通常需要自己调整初猜轨道占据方式

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