计算化学公社

标题: 请教铁原子晶胞对甲烷分子的吸附，结构优化时分子间距离始终未改变并报错不收敛 [打印本页]

作者
Author: limen 时间: 2023-4-25 23:11
标题: 请教铁原子晶胞对甲烷分子的吸附，结构优化时分子间距离始终未改变并报错不收敛
(, 下载次数 Times of downloads: 6)
计算铁原子晶胞对甲烷分子的吸附时，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

作者
Author: 含光君 时间: 2023-4-25 23:18
SCF不收敛还是结构优化不收敛？

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

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

作者
Author: limen 时间: 2023-4-26 07:22

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

SCF

作者
Author: limen 时间: 2023-4-26 08:07

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

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

作者
Author: zjxitcc 时间: 2023-4-26 10:08
本帖最后由 zjxitcc 于 2023-4-26 10:09 编辑

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

作者
Author: wzkchem5 时间: 2023-4-27 01:04

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

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

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