计算化学公社

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

作者
Author: lwang2016 时间: 2023-9-13 09:37
标题: 求助：CP2K计算NMR屏蔽常数
各位老师，

使用CP2K计算有机晶体NMR屏蔽常数是否为1.几何优化；2.单点能计算；3.NMR屏蔽常数计算这3个步骤。

请问有没有进行过NMR屏蔽常数计算的老师，分享一下教程。

谢谢

作者
Author: sobereva 时间: 2023-9-13 23:16
不需要单点计算，优化过就可以算NMR

基本设置方式：

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

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里有NMR计算涉及的各个选项的详细解读和完整的计算例子

作者
Author: lwang2016 时间: 2023-9-15 09:05

sobereva 发表于 2023-9-13 23:16
不需要单点计算，优化过就可以算NMR

基本设置方式：

您好卢老师

进行NMR计算提示“Method not implemented for k-points ”
是由于K点设置有问题么，请问该如何解决。
谢谢。

作者
Author: sigma-jlu 时间: 2023-9-15 15:40

lwang2016 发表于 2023-9-15 09:05
您好卢老师

进行NMR计算提示“Method not implemented for k-points ”

程序已经提示的很清楚，k点功能在cp2k此版本中尚未实现，最新2023.2也没有实现（别说你用的2023.1了）；不是你的问题。

作者
Author: sobereva 时间: 2023-9-16 01:49

lwang2016 发表于 2023-9-15 09:05
您好卢老师

进行NMR计算提示“Method not implemented for k-points ”

NMR计算不支持k点

作者
Author: lwang2016 时间: 2023-9-16 15:39

sobereva 发表于 2023-9-16 01:49
NMR计算不支持k点

卢老师，您好！

去掉了K点重新进行了计算。
不知道任务是否正常结束，另外NMR屏蔽常数的数据是否也是在out文件中查询到。

谢谢

作者
Author: sobereva 时间: 2023-9-18 09:07

lwang2016 发表于 2023-9-16 15:39
卢老师，您好！

去掉了K点重新进行了计算。

正常结束了
.data文件里看

作者
Author: lwang2016 时间: 2023-9-18 14:50

sobereva 发表于 2023-9-18 09:07
正常结束了
.data文件里看

您好卢老师

屏蔽常数不知道如何去看。请指教

作者
Author: sobereva 时间: 2023-9-18 18:44

lwang2016 发表于 2023-9-18 14:50
您好卢老师

屏蔽常数不知道如何去看。请指教

各向同性磁屏蔽值明显是isotropy后面的

作者
Author: lwang2016 时间: 2023-9-18 20:05

sobereva 发表于 2023-9-18 18:44
各向同性磁屏蔽值明显是isotropy后面的

多谢老师。

现在计算NMR用哪种理论方法+基组赝势较多？

作者
Author: sobereva 时间: 2023-9-18 20:27

lwang2016 发表于 2023-9-18 20:05
多谢老师。

现在计算NMR用哪种理论方法+基组赝势较多？

看元素
C、H用下文相同的级别
revTPSS泛函结合pcSseg-1基组是计算NMR很好的选择
http://sobereva.com/623（http://bbs.keinsci.com/thread-25918-1-1.html）

作者
Author: lwang2016 时间: 2023-9-19 11:10

sobereva 发表于 2023-9-18 20:27
看元素
C、H用下文相同的级别
revTPSS泛函结合pcSseg-1基组是计算NMR很好的选择

多谢卢老师。

我是使用multiwfn工具进行cp2k输入文件制作的，进行基组赝势选择的时候没有您提到的选项，不知道如何设置。

请问您有现成的inp文件，我用来参照么。谢谢。

作者
Author: lwang2016 时间: 2023-9-19 17:38

sobereva 发表于 2023-9-18 20:27
看元素
C、H用下文相同的级别
revTPSS泛函结合pcSseg-1基组是计算NMR很好的选择

您好卢老师

在multiwfn软件里进行基组的设置，没有 pcSseg-1这个基组，
请问我该如何设置呀

谢谢老师。

作者
Author: sobereva 时间: 2023-9-20 00:17

lwang2016 发表于 2023-9-19 17:38
您好卢老师

在multiwfn软件里进行基组的设置，没有 pcSseg-1这个基组，

自己从BSE基组库上导出CP2K格式的此基组的定义的文件，自己手动设

作者
Author: lwang2016 时间: 2023-9-20 11:01
本帖最后由 lwang2016 于 2023-9-20 15:49 编辑

sobereva 发表于 2023-9-20 00:17
自己从BSE基组库上导出CP2K格式的此基组的定义的文件，自己手动设

您好卢老师，

下载了pcsseg-1.1.cp2k文件，请问具体如何使用。

是否将文件内容全部复制到Linux中cp2k/data/其它基组文件中（如，EMSL_BASIS_SETS）

谢谢

作者
Author: sobereva 时间: 2023-9-21 01:47

lwang2016 发表于 2023-9-20 11:01
您好卢老师，

下载了pcsseg-1.1.cp2k文件，请问具体如何使用。

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt

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

作者
Author: lwang2016 时间: 2023-9-21 08:10
本帖最后由 lwang2016 于 2023-9-21 10:02 编辑

sobereva 发表于 2023-9-21 01:47
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）的ppt

...

卢老师您好

赝势有下载的么？还是使用程序自带potential文件里的 ALL？

除去周期性是否要用超胞 2*2*2？

谢谢

作者
Author: sobereva 时间: 2023-9-21 10:07

lwang2016 发表于 2023-9-21 08:10
卢老师您好

赝势有下载的么？还是使用程序自带potential文件里的 ALL？

GAPW做全电子基组计算一律用ALL

“除去周期性是否要用超胞 2*2*2” 完全不知道你在说什么

作者
Author: lwang2016 时间: 2023-9-21 10:18

sobereva 发表于 2023-9-21 10:07
GAPW做全电子基组计算一律用ALL

“除去周期性是否要用超胞 2*2*2” 完全不知道你在说什么

多谢卢老师，

另外，请问几何优化的K点如何进行设置。
如，a=3.8669, b=16.6497,c=21.7468,alpha=90.0000,beta=99.5215,gamma =90.0000 这样的晶胞，进行几何优化，取6 1 1么？

作者
Author: sobereva 时间: 2023-9-22 02:45

lwang2016 发表于 2023-9-21 10:18
多谢卢老师，

另外，请问几何优化的K点如何进行设置。

k点怎么设，没有体系的具体信息没法说

作者
Author: lwang2016 时间: 2023-9-22 08:48

sobereva 发表于 2023-9-22 02:45
k点怎么设，没有体系的具体信息没法说

卢老师您好

对附件中的两个结构进行几何优化，K点如何设置

谢谢。

作者
Author: sobereva 时间: 2023-9-22 09:03

lwang2016 发表于 2023-9-22 08:48
卢老师您好

对附件中的两个结构进行几何优化，K点如何设置

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）这页ppt判断

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

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