QMMM计算报错“CPASSERT failed”

liu256136 · 发表于 Post on 2024-12-25 14:09:44

本帖最后由 liu256136 于 2024-12-25 14:12 编辑

我的体系是Fe-Mn-S三元系，主要关注Mn和S，Fe更多的是起到溶剂的作用。我尝试过将Fe的基组从DZVP更换为SZV，但是还是偏慢，所以想尝试对Fe使用MM方法减小耗时。
但是一直报“CPASSERT failed”的错误，也尝试过CP2K官网上的faq:cholesky_decomp_failed [CP2K Open Source Molecular Dynamics ]方法，但是还是不行。
麻烦各位大佬答疑解惑

FeMnS-QMMM.inp (20.7 KB, 下载次数 Times of downloads: 3)

FeMnS-QMMM.out (7.77 KB, 下载次数 Times of downloads: 0)

FeC.eam (266.1 KB, 下载次数 Times of downloads: 1)

Uus/pMeC6H4-/キ · 发表于 Post on 2024-12-25 14:26:22

先不说计算如何，这个混沌无序的初始结构是怎么构建的，有没有化学意义……？

liu256136 · 发表于 Post on 2024-12-25 14:36:49

Uus/pMeC6H4-/キ发表于 2024-12-25 14:26
先不说计算如何，这个混沌无序的初始结构是怎么构建的，有没有化学意义……？

您好，这个从Fe里面随机替换的Mn和S，近似模拟一个钢液的状态，所以Mn和S就是偏向于无序的。另外其实这个结构我也用DFT优化过了，DFT是跑起来没问题的。

Uus/pMeC6H4-/キ · 发表于 Post on 2024-12-25 15:16:06

liu256136 发表于 2024-12-25 14:36
您好，这个从Fe里面随机替换的Mn和S，近似模拟一个钢液的状态，所以Mn和S就是偏向于无序的。另外其实这个 ...

这是要模拟高温钢液冷却形成无定形合金的过程？直接取个无序的结构做几何优化显然很难收敛，就算收敛了得到的结构多半也是能量过高的局部极小点，所以我会好奇有没有化学意义。另外如果原先是纯铁体系，直接替换原子种类肯定也得不到合理的电子状态（哪怕是作为几何优化初猜也很糟）。

或许更好的方法是构建掺杂的铁结构超胞，做升温熔化再缓慢降温的分子动力学模拟，再把退火完的结构拿去几何优化。DFT耗时过长先考虑解决SCF收敛问题（参考社长博文http://sobereva.com/665），不行的话还有GFN1-xTB可以尝试；做分子力学模拟得找适合Fe-Mn-S体系的EAM多体势凑齐所有参数，不比DFT计算简单多少，而且此时用LAMMPS还可能更合适。

liu256136 · 发表于 Post on 2024-12-25 15:44:22

Uus/pMeC6H4-/キ发表于 2024-12-25 15:16
这是要模拟高温钢液冷却形成无定形合金的过程？直接取个无序的结构做几何优化显然很难收敛，就算收敛了得 ...

首先感谢您的耐心回复，这个Fe-Mn-S其实仅仅只是一部分，实际上我的模型中还有别的物质作为固态基底，想去观察Mn和S与基底之间的相互作用。

关于Fe-Mn-S结构的优化，我最初确实是使用xTB方法跑过一段时间的高温分子动力学，然后切换为DFT方法开始结构优化。结构优化初期其实还是无法顺利收敛，这个可能确实是OT方法的限制，几乎是纯金属体系，所以按照社长的帖子尝试了DFT+U，顺利收敛后，使用波函数作为初猜，去掉DFT+U然后继续进行的结构优化。
当然其实这个结构是不是位于势能面最低其实对我的研究，影响不大。

有关这个经典分子动力学模拟，我最初其实是做经典分子动力学，但是确实是由于势函数限制的问题，不得已尝试从头算分子动力学。并且，我的研究其实有点类似于一个吸附的问题，在一些文献的对比中，使用DZVP的基组和TZV2P的基组，差异还是比较明显的。所以我只能尝试缩减Fe这部分的计算量。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] QMMM计算报错“CPASSERT failed”

浏览过的版块