白皛白 发表于 2025-11-25 15:49 明白了,谢谢您, ![]() |
shinnashiori 发表于 2025-11-25 10:09 好的好的!谢谢谢谢! |
leixiaocui 发表于 2025-11-24 15:15 理论上这两个方法得到的能量是完全一致的,孤立的气体分子等于是放到了距离表面∞远的位置。不过气体分子进行自由能校正时需要考虑气体的气压,如果常温是液体的话需要用到饱和蒸汽压,这需要使用VASPKIT的502功能。不过体系2经过频率计算后是没法使用502功能进行自由能校正的,因此实际操作中使用方法1是更合理的。 |
leixiaocui 发表于 2025-11-24 15:15 你看你是以哪一个状态为最后一步了~如果你认为 CO 的脱附这个过程不需要另算那就算一个吸附 CO 的结构;如果你是以脱附后的 clean 表面,那就 clean 表面加上 CO 的能量。比如我在做 OER 的时候,我算的中间体就是 * → *OH → *O → *OOH → *OO → * |
白皛白 发表于 2024-6-27 10:09 您好,冒昧打扰,电催化在计算生成产物(以CO分子举例)的时候:1、是分开算把表面的能量加上co(孤立分子的能量)加和起来作为最后一步的自由能,2、还是优化一个把CO分子放在距离表面很远的POSCAR,整体优化直接得出能量呢? |
白皛白 发表于 2024-6-27 10:09 非常感谢,步骤很详细,帮助真心大 |
|
软件需要VASP,vaspkit,VESTA 1,固定一部分原子,以OOH*为例一般来说可以把除了OOH以外的原子全部固定,使用的POSCAR是之前优化好的结构,具体方法是使用vaspkit→402号功能→1(固定POSCAR)→1(通过原子列表选择)→输入要固定的原子的编号(编号可以在VESTA中查看)→all(三个方向都固定住)→生成POSCAR_FIX,确定POSCAR_FIX内容没问题后,mv POSCAR_FIX POSCAR。 2,修改INCAR,主要参数就是IBRION = 5以及NFREE = 2这两个,之后提交任务。 3,计算完成后,在提交任务的文件夹下,使用vaspkit→501号功能→输入温度(单位是k)。如果没问题的话会得到如图的输出,将其中的Thermal correction to G(T)这一栏的eV的能量加到结构的能量E上就可以了。 |
202406271002104822..png (13.12 KB, 下载次数 Times of downloads: 27)
501功能的输出
1.3 KB, 下载次数 Times of downloads: 11
INCAR样例
7.94 KB, 下载次数 Times of downloads: 4
固定了部分原子的POSCAR
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-24 08:13 , Processed in 0.175218 second(s), 25 queries , Gzip On.