计算化学公社

标题: 关于使用CP2K进行硅表面悬键键能计算的求助 [打印本页]

作者
Author: xinxin19 时间: 2025-10-31 21:48
标题: 关于使用CP2K进行硅表面悬键键能计算的求助
各位老师好，我最近需要计算Si表面悬键Si-H键、Si-OH、Si-F、Si-Cl、Si-Br的键能，但得到的结论与文献或平常认知的不同，比如计算得到的Si-H键键能大于Si-OH的键能，Si-Cl、Si-Br的键能小于计算得到的Si-OH的键能，这令我十分困惑。所以想求助各位老师，辛苦帮我看一下是哪里出现了问题。我以计算的Si-H键为例，将计算的过程详细列在了下面。

（1）结构优化。我根据文献所述，以每个硅原子带一个氢的表面进行结构优化（即表面重构），具体的输入文件以附件形式上传，名字为Si_H_opt_chonggou.inp，计算得到的.restart与.out文件也一并上传；并用同样基组、k点和周期性下计算了这个体系的能量，具体的输入文件和输出文件也一并以附件形式上传，名字为Si_H_opt_chonggou_energy.inp/out。
（2）删掉一个氢原子后的能量计算。接下来我将硅表面的一个氢原子删掉，计算剩下的体系的能量，k点与周期性的设置与结构优化的条件是一致的，其中体系电荷设置如下：CHARGE 0 #Net charge MULTIPLICITY 2 #Spin multiplicity UKS。具体的输入文件和输出文件也一并以附件形式上传，名字为Si_H_opt_chonggou_deleted_H.inp/out
（3）计算氢原子的能量。最后我将第（2）步中被删掉的氢原子留下，计算H的能量，不设置k点与周期性，体系的电荷设置如下：CHARGE 0 #Net charge MULTIPLICITY 2 #Spin multiplicity UKS。具体的输入文件和输出文件也一并以附件形式上传，名字为H_ENERGY.inp/out
（4）键能计算公式。最后键能照样按E(AB)-E(A)-E(B)得到。
上述就是我计算硅表面键能的详细方法，还请大家帮忙看一下哪里有误，这个过程应该怎样是合理准确的。万分感谢！！！

作者
Author: xinxin19 时间: 2025-11-16 17:37
还请各位可能知道什么原因的老师停步来帮助指导一下哦，至今还是没有解决这个问题

作者
Author: 15939474177 时间: 2026-1-5 10:11
楼主请问这个问题解决了吗

作者
Author: xinxin19 时间: 2026-2-11 14:20

15939474177 发表于 2026-1-5 10:11
楼主请问这个问题解决了吗

没有呢，如果你有想法怎么操作，我们可以继续交流哦

作者
Author: northern_mar 时间: 2026-2-11 21:24
我也是新手，楼主方便把计算得到的键能数据发一下嘛？如果数值相差很小，比如只有10kJ/mol左右的差别的话，提高计算精度可能有效。
我刚刚看了下你输入文件，计算精度不够。卢老师上课讲的能量计算至少要用TZVP等级，并且你这个体系较小，使用TZV2P也没问题。k点收敛性测试不知道你有没有做过，你的晶胞大小建议使用441。
同时周期性条件可以试试调整为一致的。

作者
Author: xinxin19 时间: 2026-2-12 15:07

northern_mar 发表于 2026-2-11 21:24
我也是新手，楼主方便把计算得到的键能数据发一下嘛？如果数值相差很小，比如只有10kJ/mol左右的差别的话， ...

方便，这是我计算得到的结果，文献实验值Si-H键能为3.2ev，差异还挺大的

断键前能量（hartree） Si- H F OH Cl Br等原子 键能（Hartree） 键能（ev）
Si-H -171.0754527 -170.451457 -0.464070509 0.159925137 4.351786885
Si-F -360.7544296 -336.4246936 -24.12720648 0.202529502 5.5111113
Si-OH -299.636761 -282.9417077 -16.52859259 0.166460677 4.529628061
Si-Cl -286.5809919 -271.5203861 -14.91379348 0.146812331 3.99496905
Si-Br -274.0053303 -260.518306 -13.36676549 0.120258813 3.272410674

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