计算化学公社

标题: Gaussian算吸附能结构优化不出自己想要的构型怎么解决 [打印本页]

作者
Author: LIberate 时间: 2024-12-4 13:08
标题: Gaussian算吸附能结构优化不出自己想要的构型怎么解决
本人最近刚刚接触计算，想请教一下各位老师同学，我想参考的文献（图1）算一下BN不同氧气吸附类型的吸附能，但是pauiling吸附构型（图2）结构优化后就变成了另外两种图（3），pauiling吸附构型的结构可能的确不稳定，那我如何结构优化能得到这种构型来计算他的吸附能呢？

目前已经尝试了多种冻结优化的方法，只有pauling吸附的构型，我把氧原子固定，能得到大约符合的构型。但我看到文献，标准做法是将固体的原子冻结，于是我把BN冻结，只做氧气吸附，得到的构型依旧变成如图3的样式。

我所使用的优化方法如下：

# opt=(cartesian, MaxStep=5, NoTrustUpdate, GDIIS) b3lyp/6-31g(d) nosymm empiricaldispersion=gd3bj

(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1)

作者
Author: Uus/pMeC6H4-/キ 时间: 2024-12-4 13:16
这截的BN基底又小又薄还不符合晶体表面结构。虽然对簇模型做量子化学计算也是研究吸附的一种方法，但刚开始接触的话还是改用CP2K之类的做周期性体系计算吧。

作者
Author: cokie 时间: 2024-12-4 13:51
本帖最后由 cokie 于 2024-12-4 13:56 编辑

你取片段之后，边缘的BN有用H封边饱和吗。

作者
Author: sobereva 时间: 2024-12-4 17:27
片段取得太小，而且不加氢饱和边缘根本不行，当前模型虚假的边界效应太强
什么是合理的模型，参考下文
18碳环（cyclo[18]carbon）与石墨烯的相互作用：基于簇模型的研究一例
http://sobereva.com/615（http://bbs.keinsci.com/thread-25306-1-1.html）

体系这么小还用6-31G*差点意思，起码O2用6-311G*

整个slab绝对不能都冻结，只能冻结边缘的

作者
Author: LIberate 时间: 2024-12-4 20:52

sobereva 发表于 2024-12-4 17:27
片段取得太小，而且不加氢饱和边缘根本不行，当前模型虚假的边界效应太强
什么是合理的模型，参考下文
18 ...

非常感谢sob老师，学习修改后我再进行计算

作者
Author: LIberate 时间: 2024-12-4 20:54

cokie 发表于 2024-12-4 13:51
你取片段之后，边缘的BN有用H封边饱和吗。

没有，非常感谢

作者
Author: LIberate 时间: 2024-12-4 20:57

Uus/pMeC6H4-/キ发表于 2024-12-4 13:16
这截的BN基底又小又薄还不符合晶体表面结构。虽然对簇模型做量子化学计算也是研究吸附的一种方法，但刚开始 ...

非常感谢，我学习一下

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