求助：如何利用CP2K计算拉曼光谱（例子或者参考资料）

HuangL · 发表于 Post on 2021-12-7 11:26:09

最近了解了一些CP2K相关的内容，想请问各位老师如果要用CP2K计算拉曼光谱的话有相关的参考资料或者例子推荐吗？另外，了解到CP2K一般用来算周期性体系，如果是算聚合物的拉曼光谱的话，在拉曼位移和强度方面跟实验做比较相差大吗？

sobereva · 发表于 Post on 2021-12-7 13:55:46

Google搜Computing Bulk Phase Vibrational Spectra with CP2k and TRAVIS文档

sigma-jlu · 发表于 Post on 2023-9-7 17:57:32

sobereva 发表于 2021-12-7 13:55
Google搜Computing Bulk Phase Vibrational Spectra with CP2k and TRAVIS文档

google给出的解答，Martin Brehm教授的讲稿与tutorial，给出的cp2k的方法均是aimd和md runtype。
但在runtype vibration_analysis下面，不是也有一个Raman选项吗，不知为何无法正常计算结束，写信骚扰过Brehm教授，但尚未获得回复。

不知卢老师有没有用vibration_analysis跑成过Raman的input例子？

youyno · 发表于 Post on 2023-9-7 19:46:07

我们课题组一个已毕业博士写的一些脚本比TRAVIS好用，
https://github.com/AlexanderHoff ... ipts/PostProcessing

sobereva · 发表于 Post on 2023-9-7 21:20:32

sigma-jlu 发表于 2023-9-7 17:57
google给出的解答，Martin Brehm教授的讲稿与tutorial，给出的cp2k的方法均是aimd和md runtype。
但在ru ...

Multiwfn直接就能基于CP2K的振动分析的输出绘制拉曼光谱。绘制操作同下文
使用Multiwfn计算特定方向的UV-Vis吸收光谱
http://sobereva.com/648（http://bbs.keinsci.com/thread-30873-1-1.html）

需要用的CP2K输入文件，在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）培训里都讲了

仅当必须用MD的时候才需要用travis，此时非常昂贵

sigma-jlu · 发表于 Post on 2023-9-8 16:33:02

sobereva 发表于 2023-9-7 21:20
Multiwfn直接就能基于CP2K的振动分析的输出绘制拉曼光谱。绘制操作同下文
使用Multiwfn计算特定方向的UV ...

多谢卢老师解答。
原本想问的是，您有没有vibration_analysis成功运行过的inp文件可以分享？
我这里IR分析OK，但Raman分析在中间中断退出，out中看像是系统掉电一样。均附如下。

note：在Raman分析中，均根据您cp2k培训班教材指导填写--
INTENSITIES T
……

DO_RAMAN T
PERIODIC_DIPOLE_OPERATOR T

……
其中基组用到pcSseg-2，在rar文件中都附了。

如您时间有限而无法对Raman-break中的inp文件进行挑错，能否分享一个运行正常结束的Raman输入文件？

sobereva · 发表于 Post on 2023-9-8 21:46:22

sigma-jlu 发表于 2023-9-8 16:33
多谢卢老师解答。
原本想问的是，您有没有vibration_analysis成功运行过的inp文件可以分享？
我这里IR ...

freq.inp (4.33 KB, 下载次数 Times of downloads: 83)

freq.out (159.4 KB, 下载次数 Times of downloads: 40)

sigma-jlu · 发表于 Post on 2023-9-13 16:38:24

sobereva 发表于 2023-9-8 21:46

收到，测试OK。
发现Raman计算break的原因主要是其对泛函的依赖，在使用Molopt基组情况下，只是把TPSS换成PBE，就能顺利计算出Raman结果。原想有过渡元素首选TPSS，看来用在Raman分析上，必须例外对待。

谢谢卢老师。

sobereva · 发表于 Post on 2023-9-13 21:58:15

sigma-jlu 发表于 2023-9-13 16:38
收到，测试OK。
发现Raman计算break的原因主要是其对泛函的依赖，在使用Molopt基组情况下，只是把TPSS换 ...

TPSS是meta-GGA泛函，PBE是GGA，有可能在对meta-GGA的支持上有限制。另外，TPSS在CP2K里有自身的计算代码，也可以基于libXC接口利用libXC的代码计算，一种不行的话可以试另一种。
在计算Raman上，没有证据TPSS比PBE有任何优势。不需要特意用TPSS。

不会扣篮的后卫 · 发表于 Post on 2024-2-11 10:58:53

youyno 发表于 2023-9-7 19:46
我们课题组一个已毕业博士写的一些脚本比TRAVIS好用，
https://github.com/AlexanderHoffman/supporting-i ...

老师，请问可以简单介绍一下画CP2K输出的拉曼光谱代码吗（我在PostProcessing文件夹里，只看到了针对于QE，VASP的输出拉曼光谱绘制代码)?

sobereva · 发表于 Post on 2024-2-19 18:02:07

不会扣篮的后卫发表于 2024-2-11 10:58
老师，请问可以简单介绍一下画CP2K输出的拉曼光谱代码吗（我在PostProcessing文件夹里，只看到了针对于QE ...

如果是对单个结构的话，Multiwfn直接就能基于CP2K振动分析输出文件（计算了拉曼活性的情况）绘制拉曼光谱，按下文操作就完了
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD、VCD和ROA光谱图
http://sobereva.com/224

不会扣篮的后卫 · 发表于 Post on 2024-2-19 21:53:07

sobereva 发表于 2024-2-19 18:02
如果是对单个结构的话，Multiwfn直接就能基于CP2K振动分析输出文件（计算了拉曼活性的情况）绘制拉曼光谱 ...

谢谢卢老师！

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

[CP2K] 求助：如何利用CP2K计算拉曼光谱（例子或者参考资料）