计算化学公社

标题: 新手求教一个关于晶体表面吸附能计算的问题 [打印本页]

作者
Author: 冰皮麻薯麻 时间: 2023-3-28 11:10
标题: 新手求教一个关于晶体表面吸附能计算的问题
对于晶体表面吸附吸附能计算的公式：Ead=Esm-Es-Em。其中Esm表示几何优化后的吸附模型能量，Es表示底面能量，Em表示吸附分子能量。方法一：之前看B站上的一个教学视频里获取Es和Em的方式是：把优化好的吸附模型复制两份，一份删除吸附分子，一份删除底面；然后分别计算能量。
方法二：Es和Em分别用开始吸附之前，他们刚分别独立完成几何优化后的模型进行能量计算。
两种方法算出来的能量有明显的差异。
想问下大佬们是怎么做的。

作者
Author: wzkchem5 时间: 2023-3-28 21:38
取决于计算目的。如果是为了和实验的吸附能对比，那吸附分子和基底必须分别优化，但至于是用吸附之前优化的结构，还是用吸附之后重新优化的结构，影响不大。对于极个别其他目的，例如能量分解分析，可能需要用第一种方法。

作者
Author: sobereva 时间: 2023-3-29 01:20
第二种是常规做法

第一种只适合计算吸附结构下的相互作用能，或者计算作为第二种方法的近似（仅当吸附导致的变形效应可忽略不及时才能接受）

作者
Author: 冰皮麻薯麻 时间: 2023-3-29 18:00
谢谢大佬们

作者
Author: Hubery_fei 时间: 2023-4-12 17:53
请教一下，我在计算吸附能时发现，若使用未优化的Slab模型构建吸附模型进行优化和使用结构优化完成的Slab模型构建的吸附模型进行计算后，两者的数值差距也有点大。请问大佬，该趋向于哪种呢？

作者
Author: sobereva 时间: 2023-4-12 18:09

Hubery_fei 发表于 2023-4-12 17:53
请教一下，我在计算吸附能时发现，若使用未优化的Slab模型构建吸附模型进行优化和使用结构优化完成的Slab模 ...

根据结构差异，弄清楚吸附能差异的原因

作者
Author: Hubery_fei 时间: 2023-4-12 18:22

sobereva 发表于 2023-4-12 18:09
根据结构差异，弄清楚吸附能差异的原因

感谢大佬的回复，吸附能的计算是金属表面吸附H原子。吸附模型的并没有大的差异。另外自选极化这个参数是否有必要添加呢？添加前后对于吸附能的数值也有一定的影响。麻烦大佬解惑。感谢

作者
Author: sobereva 时间: 2023-4-12 18:27

Hubery_fei 发表于 2023-4-12 18:22
感谢大佬的回复，吸附能的计算是金属表面吸附H原子。吸附模型的并没有大的差异。另外自选极化这个参数是 ...

假设你用的优化收敛限足够严，两个结构确实是不同极小点，可以取能量低的模型

如果体系是开壳层的（包括但不限于体系电子数是奇数的情况）就开开自旋极化

作者
Author: Hubery_fei 时间: 2023-4-12 19:50

sobereva 发表于 2023-4-12 18:27
假设你用的优化收敛限足够严，两个结构确实是不同极小点，可以取能量低的模型

如果体系是开壳层的（包 ...

感谢大佬的回复。结构优化已然明了。有关开闭壳层的判断想再请教下，我目前在做的Pt(111)面上吸附H原子，体系的电子数是Pt的核外电子数外加H的电子数吗？本身并未施加外场。麻烦大佬解惑，十分感谢

作者
Author: sobereva 时间: 2023-4-12 22:26

Hubery_fei 发表于 2023-4-12 19:50
感谢大佬的回复。结构优化已然明了。有关开闭壳层的判断想再请教下，我目前在做的Pt(111)面上吸附H原子， ...

是

作者
Author: Hubery_fei 时间: 2023-4-13 09:16

sobereva 发表于 2023-4-12 22:26
是

感谢您的回复 sob老师

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