计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 如何正确使用cp2k计算MOF中碳原子的NMR化学位移？（目前 ...

如何正确使用cp2k计算MOF中碳原子的NMR化学位移？（目前遇到线性求解不收敛的问题）

查看数: 2779 | 评论数: 4 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Alset

发布时间: 2022-8-4 10:44

正文摘要:

老师们好，目前我打算计算一MOF（含C、H、Mg、O，晶胞大概含300个原子）中某linker上碳原子的NMR化学位移，由于先前所作的计算均在cp2k中进行（几何优化、REPEAT电荷计算等），所以我决定同样使用cp2k进行NMR计 ...

回复 Reply

Alset 发表于 Post on 2022-8-8 19:31:44

sobereva 发表于 2022-8-6 02:54
看不懂你的意思。
如果已经有一个结构，要计算某个部分原子的NMR，就扩展一部分，挖出来，边界恰当处理 ...

明白了，非常感谢您！

sobereva 发表于 Post on 2022-8-6 02:54:58

Alset 发表于 2022-8-5 14:37
非常感谢老师指教，
我在增大线性迭代次数后，便暂时没有出现先前的问题。但是还是希望能够算短一些耗时 ...

看不懂你的意思。
如果已经有一个结构，要计算某个部分原子的NMR，就扩展一部分，挖出来，边界恰当处理后算NMR就完了。不牵扯优化，也因此没有冻结的事

Alset 发表于 Post on 2022-8-5 14:37:09

本帖最后由 Alset 于 2022-8-5 15:01 编辑

sobereva 发表于 2022-8-5 10:06
revTPSS和PBE都是纯泛函，直接用revTPSS就完了，没必要先用PBE。pcSseg-1已经算很便宜的了。至于那个报错， ...

非常感谢老师指教，
我在增大线性迭代次数后，便暂时没有出现先前的问题。但是还是希望能够算短一些耗时，所以使用簇模型可能是比较好的方法。
按照我的理解，可否将我要计算的部分截取，之后在cp2k或gaussian冻结该部分所有原子（因为我使用了含pbc、整个优化后的结构进行了MD模拟，所以为了保证自洽，或许也应该保证计算NMR时的结构与先前相同），进行nmr的计算？

sobereva 发表于 Post on 2022-8-5 10:06:05

revTPSS和PBE都是纯泛函，直接用revTPSS就完了，没必要先用PBE。pcSseg-1已经算很便宜的了。至于那个报错，没具体输入输出文件难以判断。
Gaussian用的是簇模型算的，文献里那个图应该就是实际用的簇模型。

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