计算化学公社

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

[CP2K] 请教用CP2K计算吸附能的结果为正值

[复制链接 Copy URL]

32

帖子

0

威望

1618

eV
积分
1650

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本人使用CP2K计算在真空环境中,小分子在矿物表面的吸附能,结果却为正值。想请教各位问题出在哪里?
具体过程如下:1. 小分子和矿物的的复合结构在真空中优化, 计算级别为PBE/DZVP-MOLOPT-SR-GTH, 开启DFT-D3(BJ)色散校正; 2. 将优化后的结构在TZV2P-MOLOPT-GTH级别+色散校正下计算能量E(AB)=-4185.06875 Ha, 小分子和矿物也在相同级别下计算能量E(A)=-441.79871 Ha, E(B)=-3743.33368 Ha; 3. 将复合物的整体能量与小分子和矿物的单独计算的能量相减得到吸附能Eads=E(AB)-E(A)-E(B),结果为正0.06364Ha.不知道问题出在哪里?附上能量计算的inp文件,希望帮忙解答,谢谢!

Complex_SP.inp

17.54 KB, 下载次数 Times of downloads: 15

Mineral_SP.inp

15.71 KB, 下载次数 Times of downloads: 7

Molecule_SP.inp

4.88 KB, 下载次数 Times of downloads: 10

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2025-11-22 05:52:21 | 只看该作者 Only view this author
算单点仍使用DZVP-MOLOPT-SR-GTH看看什么情况
并且确保SCF充分收敛了,也可以尝试加大CUTOFF和Z方向盒子尺寸看看
北京科音自然科学研究中心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

32

帖子

0

威望

1618

eV
积分
1650

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2025-11-22 11:12:24 | 只看该作者 Only view this author
sobereva 发表于 2025-11-22 05:52
算单点仍使用DZVP-MOLOPT-SR-GTH看看什么情况
并且确保SCF充分收敛了,也可以尝试加大CUTOFF和Z方向盒子尺 ...

谢谢卢老师回复。我将计算级别改成DZVP-MOLOPT-SR-GTH后又重新计算了一遍,结果仍然是正值0.0777Ha. 之后我又将cutoff增大到500,结果是0.0814Ha. 把Z方向从20埃增加到25埃之后,结果为0.0136Ha.很奇怪,结果全是正值。此外,不同的盒子尺寸得到的吸附能差别很大,不知道怎样的设置才是最合理的。

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

4#
发表于 Post on 2025-11-23 01:03:35 | 只看该作者 Only view this author
418086011 发表于 2025-11-22 11:12
谢谢卢老师回复。我将计算级别改成DZVP-MOLOPT-SR-GTH后又重新计算了一遍,结果仍然是正值0.0777Ha. 之后 ...

做结果随Z尺寸的收敛性测试
北京科音自然科学研究中心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

32

帖子

0

威望

1618

eV
积分
1650

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2025-11-23 01:11:38 | 只看该作者 Only view this author
sobereva 发表于 2025-11-23 01:03
做结果随Z尺寸的收敛性测试

谢谢!但是目前还不知道计算的吸附能都是正值是什么原因导致的?会不会在输入文件中有某个校正没有考虑进来?

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

6#
发表于 Post on 2025-11-23 01:54:01 | 只看该作者 Only view this author
418086011 发表于 2025-11-23 01:11
谢谢!但是目前还不知道计算的吸附能都是正值是什么原因导致的?会不会在输入文件中有某个校正没有考虑进 ...

先把尺寸问题弄清楚了,之后再说其它的(几何优化收敛限、SCF收敛情况等,以及尝试PERIODIC的PSOLVER结合偶极校正)
北京科音自然科学研究中心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

52

帖子

0

威望

590

eV
积分
642

Level 4 (黑子)

7#
发表于 Post on 2025-11-24 09:26:20 | 只看该作者 Only view this author
请问楼主,这个是否和周期性有关系呢?看了3个输入文件,都使用的二维周期性XY。要不要考虑用三维周期性XYZ试试呢?

135

帖子

0

威望

1343

eV
积分
1478

Level 4 (黑子)

8#
发表于 Post on 2025-11-24 23:20:01 | 只看该作者 Only view this author
材料太薄了,应该至少3层;z轴真空有点小,至少27;让有机分子竖起来试试,即用羟基O与金属离子配位试试

32

帖子

0

威望

1618

eV
积分
1650

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2025-11-25 06:23:46 | 只看该作者 Only view this author
Assiduitas 发表于 2025-11-24 09:26
请问楼主,这个是否和周期性有关系呢?看了3个输入文件,都使用的二维周期性XY。要不要考虑用三维周期性XYZ ...

感觉跟Z轴的尺寸有关系,这个体系没有考虑Z方向的周期性。

32

帖子

0

威望

1618

eV
积分
1650

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2025-11-25 06:25:47 | 只看该作者 Only view this author
reid 发表于 2025-11-24 23:20
材料太薄了,应该至少3层;z轴真空有点小,至少27;让有机分子竖起来试试,即用羟基O与金属离子配位试试

待我做个收敛性测试。按理说不管竖着还是横着,能量相减应该都为负才对吧?

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 20:34 , Processed in 0.267216 second(s), 23 queries , Gzip On.

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