计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1281|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 一个有机物吸附在同一催化剂的不同晶面上,计算出来的吸附能为正值合理吗

[复制链接 Copy URL]

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

各位老师好,我将同一个有机物吸附在同一催化剂的不同晶面上时计算得出的吸附能为正值,请问这正常吗?(自行谷歌了一下,众说纷纭)如果正常的话,正的吸附能是越大越稳定还是越小越稳定呢?由于非化学类专业出身,不太能分析这种情况,望各位懂的老师指点一下,非常感谢!

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

2#
发表于 Post on 2024-3-21 10:55:58 | 只看该作者 Only view this author
如果是电子能量求差,明显算得不对,除非你算的表面和被吸附物带的净电荷相同。如果是自由能求差,由于熵罚效应,温度不很低时吸附自由能为正很正常,说明热力学不自发

吸附能越负越稳定
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-3-21 20:51:33 | 只看该作者 Only view this author
sobereva 发表于 2024-3-21 10:55
如果是电子能量求差,明显算得不对,除非你算的表面和被吸附物带的净电荷相同。如果是自由能求差,由于熵罚 ...

sob老师,我算的吸附能就是用CP2K结构优化完成后得到的能量(Total FORCE_EVAL ( QS ) energy)算的,根据我对电子能量和自由能的理解,我这个应该是自由能求差,所以正值是正常的吧?还有吸附能越负越稳定对于吸附能为正值时是不是就是越小越稳定?

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

4#
发表于 Post on 2024-3-21 21:47:24 | 只看该作者 Only view this author
科研好难 发表于 2024-3-21 20:51
sob老师,我算的吸附能就是用CP2K结构优化完成后得到的能量(Total FORCE_EVAL ( QS ) energy)算的,根 ...

你用电子能量求差得到正值,这已经是错的了,必须先解决这个问题,用自由能求差是后话
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-3-22 10:13:52 | 只看该作者 Only view this author
sobereva 发表于 2024-3-21 21:47
你用电子能量求差得到正值,这已经是错的了,必须先解决这个问题,用自由能求差是后话

sob老师,CP2K计算什么是电子能量求差?计算什么又是自由能求差呢?我不是很懂,我以为Total FORCE_EVAL ( QS ) energy是自由能所以我这是自由能求差呢。还有我的错误是什么导致的呢,结构错误?计算方法错误?设置错误?麻烦老师帮我解释一下并告诉我该怎么修改一下重新计算吸附能。

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

6#
发表于 Post on 2024-3-22 17:02:15 | 只看该作者 Only view this author
科研好难 发表于 2024-3-22 10:13
sob老师,CP2K计算什么是电子能量求差?计算什么又是自由能求差呢?我不是很懂,我以为Total FORCE_EVAL  ...

你不做振动分析怎么可能给你自由能。先把此文看了(PS:亟需恶补基本常识)
谈谈该从Gaussian输出文件中的什么地方读电子能量
http://sobereva.com/488http://bbs.keinsci.com/thread-13450-1-1.html
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552http://bbs.keinsci.com/thread-17494-1-1.html

参考北京科音基础(中级)量子化学培训班(http://www.keinsci.com/workshop/KBQC_content.html)的ppt:
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-3-25 21:58:47 | 只看该作者 Only view this author
sobereva 发表于 2024-3-21 21:47
你用电子能量求差得到正值,这已经是错的了,必须先解决这个问题,用自由能求差是后话

sob老师,我分子是B3LYP/6-31G*计算的,催化剂用的PBEsol/pob-TZVP-rev2,吸附用的PBE/DZVP-MOLOPT-SR-GTH,电子能量求差为正值是不是因为泛函和机组不统一的原因?

6万

帖子

99

威望

5万

eV
积分
120141

管理员

公社社长

8#
发表于 Post on 2024-3-26 00:47:59 | 只看该作者 Only view this author
科研好难 发表于 2024-3-25 21:58
sob老师,我分子是B3LYP/6-31G*计算的,催化剂用的PBEsol/pob-TZVP-rev2,吸附用的PBE/DZVP-MOLOPT-SR-GT ...

显然必须相同计算级别得到的能量才能求差,否则系统性误差没法充分抵消,结果没任何意义
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

148

帖子

0

威望

737

eV
积分
885

Level 4 (黑子)

9#
发表于 Post on 2024-3-26 09:09:18 | 只看该作者 Only view this author
可能还要注意一下吸附能的定义式,有的文献中吸附能的定义是和一般的做法反着来得

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-3-26 10:29:26 | 只看该作者 Only view this author
sobereva 发表于 2024-3-26 00:47
显然必须相同计算级别得到的能量才能求差,否则系统性误差没法充分抵消,结果没任何意义

好的,谢谢sob老师的解答

130

帖子

0

威望

703

eV
积分
833

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-3-26 10:30:34 | 只看该作者 Only view this author
123wo 发表于 2024-3-26 09:09
可能还要注意一下吸附能的定义式,有的文献中吸附能的定义是和一般的做法反着来得

好的,我再好好看一下,谢谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 22:10 , Processed in 0.155294 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list