计算化学公社

标题: 求助：CP2K计算NMR屏蔽常数 [打印本页]

作者
Author: lwang2016 时间: 2023-9-28 09:51
标题: 求助：CP2K计算NMR屏蔽常数
本帖最后由 lwang2016 于 2023-9-28 09:51 编辑

各位老师好。

我使用CP2K进行了四甲基硅烷NMR屏蔽常数的计算。

1.几何优化使用了PBEsol+DZVP-MOLOPT-SR-GTH

2.NMR计算使用了revTPSS+pcsSeg-1

计算得到的C的屏蔽常数如下图，这些C原子的屏蔽常数应该是一样的，但结果却明显不是。
请问是哪里有问题呢，谢谢老师们。

作者
Author: sobereva 时间: 2023-9-28 09:53
晶胞明显太小
必须扩到足够大才能得到靠谱的结果

作者
Author: lwang2016 时间: 2023-9-28 09:54
本帖最后由 lwang2016 于 2023-9-28 11:28 编辑

sobereva 发表于 2023-9-28 09:53
晶胞明显太小
必须扩到足够大才能得到靠谱的结果

老师。

3*3*3 是否合适？

作者
Author: sobereva 时间: 2023-9-28 12:02

lwang2016 发表于 2023-9-28 09:54
老师。

3*3*3 是否合适？

2*3*4

作者
Author: lwang2016 时间: 2023-9-28 13:55

sobereva 发表于 2023-9-28 12:02
2*3*4

卢老师您好。

请问晶胞大小N1*N2*N3 的选择，一般是遵循怎样的原则呀。

谢谢。

作者
Author: sobereva 时间: 2023-9-28 16:53

lwang2016 发表于 2023-9-28 13:55
卢老师您好。

请问晶胞大小N1*N2*N3 的选择，一般是遵循怎样的原则呀。

几句话说不清楚，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里详细讲了
不放心就做结果随晶胞增大的收敛性测试。另外，对于有等价原子的情况，从原子的等价性满足的好坏上也可以一定程度判断当前超胞尺寸是否够大

作者
Author: lwang2016 时间: 2023-10-3 14:54

sobereva 发表于 2023-9-28 16:53
几句话说不清楚，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html ...

谢谢老师！

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