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

limen

计算铁原子晶胞对甲烷分子的吸附时，Fe原子加SDD赝势基组，其它原子6-31g*基组，为什么算到129cycle时仍不收敛，而且14和19号原子距离从始至终都未改变。
是否是因为电脑算不动？（没有工作站，电脑处理器为12th Gen Intel(R) Core(TM) i9-12900H   2.50 GHz）。
%chk=mx-Fe-CH4.chk
%nprocshared=8
#p opt b3lyp/genecp


Title Card Required


0 1
Fe         -1        0.00617446    0.01234901    0.00617446
Fe          -1       0.00617446    0.01234901    2.86917464
Fe          -1       2.86917464    0.01234901    0.00617446
Fe          -1       2.86917464    0.01234901    2.86917464
Fe         -1        1.43767455    0.01234901    1.43767455
Fe         -1        0.00617446    1.44384910    1.43767455
Fe          -1       1.43767455    1.44384910    2.86917464
Fe         -1       2.86917464    1.44384910    1.43767455
Fe         -1        1.43767455    1.44384910    0.00617446
Fe         -1        0.00617446    2.84077171    2.83459717
Fe         -1        2.83459717    2.84077171    2.83459717
Fe          -1       0.00617446    2.84077171    0.00617446
Fe          -1       2.83459717    2.84077171    0.00617446
Fe          0       1.42038581    2.84077171    1.42038581
C           0       1.91933806    5.17710389    1.63541561
H            0      2.92048068    5.07409142    1.27209499
H            0      1.94137634    5.42939687    2.67501292
H            0      1.42023512    5.95132383    1.09100791
H           0       1.39526049    4.25360336    1.50354646


Fe 0
SDD
****
C H 0
6-31g*
****


Fe 0
SDD

含光君

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

wzkchem5

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

sobereva

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

limen

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

SCF

limen

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

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

zjxitcc

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

wzkchem5

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

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