计算化学公社

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

[CP2K] cp2k 使用CI-NEB计算过渡态报错Cholesky decomposition failed. Matrix ill condit...

[复制链接 Copy URL]

4

帖子

0

威望

27

eV
积分
31

Level 2 能力者

老师们好!我在用CP2K使用CI-NEB计算Xe原子在SiC相邻两个八面体间隙之间的过渡态,始末结构已经优化过了,用sobNEB生成了32个副本,每个副本分配6个线程,multiwfn生成的NEB.inp,但是刚开始运行的时候报了这个Cholesky decomposition failed. Matrix ill conditioned ?错误,用的192核的服务器跑的,试着用了不同版本的cp2k, popt ,psmp 都是这样,后来尝试改变计算级别由原来的PBE+DFT-D3(BJ)改成GFN1-xTB大概跑了17步还是出现了这个问题,不知道该如何解决?下边是我的输入文件,输出文件和始末结构,还有sobNEB生成的副本。

QQ图片20231108170701.png (57.9 KB, 下载次数 Times of downloads: 10)

QQ图片20231108170701.png

A_opted.xyz

11.14 KB, 下载次数 Times of downloads: 4

B_opted.xyz

11.14 KB, 下载次数 Times of downloads: 3

NEB.inp

18.53 KB, 下载次数 Times of downloads: 9

NEB.out

19.27 KB, 下载次数 Times of downloads: 3

副本.rar

122.44 KB, 下载次数 Times of downloads: 4

5万

帖子

99

威望

5万

eV
积分
112390

管理员

公社社长

2#
发表于 Post on 2023-11-9 19:08:45 | 只看该作者 Only view this author
改用对角化代替OT试试

插点没必要那么多,过于昂贵

随便看了下10.xyz,Xe与Si离得太近,本来Xe的半径就很大,这可能会造成数值问题。尝试重新产生插点,尽可能避免原子间过近距离出现。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

4

帖子

0

威望

27

eV
积分
31

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-11-16 14:15:29 | 只看该作者 Only view this author
sobereva 发表于 2023-11-9 19:08
改用对角化代替OT试试

插点没必要那么多,过于昂贵

社长,我用了对角化,最后跑出来的结果过渡态的能量比始末要低,是不是没有过渡态或者是路径上应该有中间态呀,而且我看了一下BAND.out,发现每个点都没有达到收敛的条件,是我反应路径选的不对么

8to8.png (58.68 KB, 下载次数 Times of downloads: 5)

8to8.png

QQ截图20231116140913.png (19.29 KB, 下载次数 Times of downloads: 6)

QQ截图20231116140913.png

NEB.out

126.5 KB, 下载次数 Times of downloads: 3

NEB_traj.xyz

152.02 KB, 下载次数 Times of downloads: 2

5万

帖子

99

威望

5万

eV
积分
112390

管理员

公社社长

4#
发表于 Post on 2023-11-17 01:07:08 | 只看该作者 Only view this author
shaozhuang 发表于 2023-11-16 14:15
社长,我用了对角化,最后跑出来的结果过渡态的能量比始末要低,是不是没有过渡态或者是路径上应该有中间 ...

不一定没过渡态,但至少当前NEB轨迹对应的路径没过渡态。可以用中间那个结构做个几何优化,看看能跑到什么结构,然后再考虑是否有其它路径、过渡态的可能
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

4

帖子

0

威望

27

eV
积分
31

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-11-17 11:12:11 | 只看该作者 Only view this author
sobereva 发表于 2023-11-17 01:07
不一定没过渡态,但至少当前NEB轨迹对应的路径没过渡态。可以用中间那个结构做个几何优化,看看能跑到什 ...

好的老师,是不是所有band的4个收敛条件都达到,结果才算准确呢

5万

帖子

99

威望

5万

eV
积分
112390

管理员

公社社长

6#
发表于 Post on 2023-11-17 12:43:19 | 只看该作者 Only view this author
shaozhuang 发表于 2023-11-17 11:12
好的老师,是不是所有band的4个收敛条件都达到,结果才算准确呢

也不是非得那么苛刻
从能量上看,变化很小时就基本可以当做收敛
而且从实际用途上来说,NEB的每个点也不需要像几何优化那样收敛得很精确,只要NEB轨迹能把反应路径正确描绘出来就够了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

4

帖子

0

威望

27

eV
积分
31

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-11-20 15:36:49 | 只看该作者 Only view this author
sobereva 发表于 2023-11-17 12:43
也不是非得那么苛刻
从能量上看,变化很小时就基本可以当做收敛
而且从实际用途上来说,NEB的每个点也 ...

谢谢老师,重新找了一个路径结果正常了,我想问老师,cp2k 用AIMD通过MSD算Xe间隙原子在SiC中的扩散系数是不是 等同于 相同情况下 用CI-NEB通过设置几个可能的扩散路径求得的几个相应的扩散系数的平均呢?AIMD是不是需要模拟够长的时间才能体现出来这个平均(对于那些扩散系数小的原子)

5万

帖子

99

威望

5万

eV
积分
112390

管理员

公社社长

8#
发表于 Post on 2023-11-21 01:48:57 | 只看该作者 Only view this author
shaozhuang 发表于 2023-11-20 15:36
谢谢老师,重新找了一个路径结果正常了,我想问老师,cp2k 用AIMD通过MSD算Xe间隙原子在SiC中的扩散系数 ...

CI-NEB并不给你扩散系数,这是静态的计算
顶多是用AIMD跑出来的扩散系数折算出来的活化能和CI-NEB计算的势垒相对比
想算准扩散系数,而且被计算的对象在体系中数目很少甚至就1个的情况、扩散势垒本身又相对较高的情况,需要相当相当长时间的MD才靠谱。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 17:57 , Processed in 0.426912 second(s), 29 queries , Gzip On.

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