计算化学公社

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

[CP2K] cp2k对NO进行结构优化应该怎么做

[复制链接 Copy URL]

18

帖子

0

威望

77

eV
积分
95

Level 2 能力者

老师您好,我打算对一氧化氮NO做波函数的分析,对HOMO和LUMO进行分析,然后首先画了一个NO分子,用cp2k进行结构优化,但是不收敛,请问对于NO的结构优化,该怎么设置?
其中尝试过:自己画NO分子,保存为cif文件,a,b,c分别是10,α,β,γ分别是90°,对于电荷和自旋多重度我按照默认,分别是0和2,但是不收敛;其次,也尝试过用xyz文件进行结构优化,但也不收敛。
感觉结构不合理,但不知道应该要怎么做




25

帖子

0

威望

260

eV
积分
285

Level 3 能力者

2#
发表于 Post on 2025-12-15 23:09:40 | 只看该作者 Only view this author
给NO的N设置这么写
   &KIND N
      ELEMENT N
      BASIS_SET TZVP-MOLOPT-GTH-q5
      POTENTIAL GTH-PBE
      MAGNETIZATION  1.00
    &END KIND
然后设置MULTIPLICITY    2
打开UKS

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

3#
发表于 Post on 2025-12-16 00:11:42 | 只看该作者 Only view this author
至少给出你的完整输入文件,让别人弄清楚你的问题在哪、都有什么硬伤

另外强烈建议通过北京科音CP2K第一性原理计算培训班http://www.keinsci.com/KFP)好好系统学一遍怎么正确用CP2K,这么跌跌撞撞试错的方式试图掌握CP2K的正确用法极其困难,本来CP2K就是最难自学的程序之一。培训里的“能量的计算及相关问题”这一节超级详细讲了怎么正确算孤立体系,没这些背景知识不可能合理地做孤立体系计算:
北京科音自然科学研究中心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

18

帖子

0

威望

77

eV
积分
95

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2025-12-16 10:03:07 | 只看该作者 Only view this author
sobereva 发表于 2025-12-16 00:11
至少给出你的完整输入文件,让别人弄清楚你的问题在哪、都有什么硬伤

另外强烈建议通过北京科音CP2K第一 ...

好滴,谢谢老师解答,我昨天尝试调整N-O键距离,实现了收敛,优化出了结构,请老师帮忙看一下这个输入文件是否还有需要修改的地方,感谢老师
NO.inp (3.6 KB, 下载次数 Times of downloads: 1)


6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

5#
发表于 Post on 2025-12-17 02:25:21 | 只看该作者 Only view this author
CAOXI 发表于 2025-12-16 10:03
好滴,谢谢老师解答,我昨天尝试调整N-O键距离,实现了收敛,优化出了结构,请老师帮忙看一下这个输入文 ...

当前输入文件没显著问题,但建议用非周期边界条件结合PSOLVER wavelet,优雅得多,因为完全避免了分子与其镜像的静电相互作用,而且此时可以用小得多的盒子尺寸,计算耗时也低得多
北京科音自然科学研究中心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

18

帖子

0

威望

77

eV
积分
95

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2025-12-17 23:29:51 | 只看该作者 Only view this author
sobereva 发表于 2025-12-17 02:25
当前输入文件没显著问题,但建议用非周期边界条件结合PSOLVER wavelet,优雅得多,因为完全避免了分子与 ...

好滴好滴,谢谢老师,我再尝试一下

18

帖子

0

威望

77

eV
积分
95

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-12-18 09:29:38 | 只看该作者 Only view this author
sobereva 发表于 2025-12-17 02:25
当前输入文件没显著问题,但建议用非周期边界条件结合PSOLVER wavelet,优雅得多,因为完全避免了分子与 ...

感谢老师的解答,我按照非周期设置了inp文件,并且开启OT,但是不收敛;而不用非周期性的就能收敛,只是计算的结果和别人文献中的不一样,所以想再请问一下老师,这个inp文件还有什么可以改进的地方。周期性NONE: NO.inp (3.92 KB, 下载次数 Times of downloads: 1)

后面这两个文件是我不使用非周期设置进行结构优化和能量计算产生molden的两个inp文件,辛苦老师也帮我看看。周期性xyz: NO.inp (3.8 KB, 下载次数 Times of downloads: 0) NO.inp (3.01 KB, 下载次数 Times of downloads: 0)

我计算出来的LUMO能量结果为-3.5155 eV,大于文献里的-4.58 eV,而-4.58 eV这个结果与我的催化剂的HOMO的计算结果才是符合的,所以想计算出这个结果,请问老师有什么改进的方法,谢谢老师。



本版积分规则 Credits rule

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

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

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