计算化学公社

标题: cp2k对NO进行结构优化应该怎么做 [打印本页]

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

(, 下载次数 Times of downloads: 0)

作者
Author: 喝酸奶的巧乐兹 时间: 2025-12-15 23:09
给NO的N设置这么写
&KIND N
   ELEMENT N
   BASIS_SET TZVP-MOLOPT-GTH-q5
   POTENTIAL GTH-PBE
   MAGNETIZATION  1.00
&END KIND
然后设置MULTIPLICITY 2
打开UKS

作者
Author: sobereva 时间: 2025-12-16 00:11
至少给出你的完整输入文件，让别人弄清楚你的问题在哪、都有什么硬伤

另外强烈建议通过北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）好好系统学一遍怎么正确用CP2K，这么跌跌撞撞试错的方式试图掌握CP2K的正确用法极其困难，本来CP2K就是最难自学的程序之一。培训里的“能量的计算及相关问题”这一节超级详细讲了怎么正确算孤立体系，没这些背景知识不可能合理地做孤立体系计算：
(, 下载次数 Times of downloads: 1)

作者
Author: CAOXI 时间: 2025-12-16 10:03

sobereva 发表于 2025-12-16 00:11
至少给出你的完整输入文件，让别人弄清楚你的问题在哪、都有什么硬伤

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

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

(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 1)

作者
Author: sobereva 时间: 2025-12-17 02:25

CAOXI 发表于 2025-12-16 10:03
好滴，谢谢老师解答，我昨天尝试调整N-O键距离，实现了收敛，优化出了结构，请老师帮忙看一下这个输入文 ...

当前输入文件没显著问题，但建议用非周期边界条件结合PSOLVER wavelet，优雅得多，因为完全避免了分子与其镜像的静电相互作用，而且此时可以用小得多的盒子尺寸，计算耗时也低得多

作者
Author: CAOXI 时间: 2025-12-17 23:29

sobereva 发表于 2025-12-17 02:25
当前输入文件没显著问题，但建议用非周期边界条件结合PSOLVER wavelet，优雅得多，因为完全避免了分子与 ...

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

作者
Author: CAOXI 时间: 2025-12-18 09:29

sobereva 发表于 2025-12-17 02:25
当前输入文件没显著问题，但建议用非周期边界条件结合PSOLVER wavelet，优雅得多，因为完全避免了分子与 ...

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

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

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

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