计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 请教铁原子晶胞对甲烷分子的吸附，结构优化时分子间距离 ...

请教铁原子晶胞对甲烷分子的吸附，结构优化时分子间距离始终未改变并报错不收敛

查看数: 1132 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

limen

发布时间: 2023-4-25 23:11

正文摘要:

计算铁原子晶胞对甲烷分子的吸附时，Fe原子加SDD赝势基组，其它原子6-31g*基组，为什么算到129cycle时仍不收敛，而且14和19号原子距离从始至终都未改变。是否是因为电脑算不动？（没有工作站，电脑处理器为12t ...

回复 Reply

wzkchem5 发表于 Post on 2023-4-27 01:04:36

limen 发表于 2023-4-26 01:07
谢谢。
多原子体系多重度怎么确定？

查一下体相Fe每个Fe原子有几个单电子，由此可以算出你的团簇有几个单电子。因为Fe是铁磁性的，所以所有单电子都是平行的，进而得到你的团簇的自旋多重度。但是考虑到表面效应，实际自旋多重度可能不精确等于由此计算的值，所以在计算值附近试几个数，取能量最低的。
当然正如楼上所说，如果你算这个Fe团簇的目的是模拟体相铁，那么应该改为用周期性程序算，团簇模型对Fe来说太过粗糙，结果不能发表（起码当团簇这么小的情况下无法发表）。如果你算Fe团簇就是为了研究Fe团簇本身，那另说

zjxitcc 发表于 Post on 2023-4-26 10:08:16

本帖最后由 zjxitcc 于 2023-4-26 10:09 编辑

建模与实际体系天壤之别，无法体现Fe金属单胞，无需解决SCF不收敛，因为计算nonsense，要用周期性程序（例如QE, CP2K, VASP）算。

limen 发表于 Post on 2023-4-26 08:07:42

wzkchem5 发表于 2023-4-25 23:45
自旋多重度错误。铁是铁磁性的，所以这个团簇的自旋多重度显然不是1

谢谢。
多原子体系多重度怎么确定？

limen 发表于 Post on 2023-4-26 07:22:03

含光君发表于 2023-4-25 23:18
SCF不收敛还是结构优化不收敛？

SCF

sobereva 发表于 Post on 2023-4-26 06:59:09

如果你是想用簇模型代替无限延展的Fe晶体表面考察吸附，劝你断了这个念头，根本没戏。
应当用第一性原理程序当周期性体系算。不是什么晶体都适合用簇模型
B3LYP更是绝对不该用在这种问题上，下文说了
使用量子化学程序基于簇模型计算金属表面吸附问题
http://sobereva.com/540（http://bbs.keinsci.com/thread-16532-1-1.html）

wzkchem5 发表于 Post on 2023-4-25 23:45:50

自旋多重度错误。铁是铁磁性的，所以这个团簇的自旋多重度显然不是1

含光君 发表于 Post on 2023-4-25 23:18:04

SCF不收敛还是结构优化不收敛？

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