计算化学公社

标题: cp2k计算偶极矩，并绘制光谱求助 [打印本页]

作者
Author: @Bruce 时间: 2022-2-16 17:15
标题: cp2k计算偶极矩，并绘制光谱求助
请问cp2k输出偶极矩文件后（RTTDDFT-output-moments.dat），怎么提取数据并画光谱图？谢谢

作者
Author: sobereva 时间: 2022-2-16 18:41
为什么不用线性响应TDDFT？用最新版Multiwfn基于CP2K的TDDFT输出文件直接就能绘制UV-Vis谱

作者
Author: @Bruce 时间: 2022-2-16 19:15

sobereva 发表于 2022-2-16 18:41
为什么不用线性响应TDDFT？用最新版Multiwfn基于CP2K的TDDFT输出文件直接就能绘制UV-Vis谱

卢老师，您好。因为我看一些相关的文献，是用cp2k的RT-TDDFT计算的。所以我尝试用了REAL_TIME_PROPAGATION。但是后处理画图一直没搞对

作者
Author: sobereva 时间: 2022-3-25 19:42

@Bruce 发表于 2022-2-16 19:15
卢老师，您好。因为我看一些相关的文献，是用cp2k的RT-TDDFT计算的。所以我尝试用了REAL_TIME_PROPAGATIO ...

看https://groups.google.com/g/cp2k/c/HWVSO5quEPQ/m/gE_M98ULAQAJ回复里有人提供的spectrum.py。自己用origin之类程序做傅里叶变换也可以实现

我感觉通过CP2K的real time propagation算电子激发问题相当难用，跑轨迹花的时间很长，结果还可能一团糟，手册写得也过于简单。只适合纯理论层面上的鼓捣，不适合作为应用层面的研究。

作者
Author: @Bruce 时间: 2022-3-28 11:31
本帖最后由 @Bruce 于 2022-3-28 11:33 编辑

sobereva 发表于 2022-3-25 19:42
看https://groups.google.com/g/cp2k/c/HWVSO5quEPQ/m/gE_M98ULAQAJ回复里有人提供的spectrum.py。自己用 ...

谢谢卢老师。我很赞同您说的手册该部分介绍简单

。我发帖之前也看到了这个python脚本，画出来的图很奇怪，后面我改了nwchem的RT-TDDFT例子的脚本，画出来的图稍微好了点，但还是不很好。另外，感觉算的时候 $MD 部分没生效似的，我就觉得输入文件有问题，可对比了google组的一些例子还是没找到问题。

作者
Author: sobereva 时间: 2022-3-28 11:33

@Bruce 发表于 2022-3-28 11:31
谢谢卢老师。我很赞同手册该部分介绍简单。我发帖之前也看到了这个，画出来的图很奇怪，后面我改了nw ...

官网的H2CO例子居然连个结果都不给，就给个输入文件，实在不友好。而且那个例子的盒子偏小，跑到后来都出现密度截断了。用NWChem、Octopus跑这个大抵更好，虽然我没试过，但应该至少不至于比这更差。

作者
Author: @Bruce 时间: 2022-3-28 11:38

sobereva 发表于 2022-3-28 11:33
官网的H2CO例子居然连个结果都不给，就给个输入文件，实在不友好。而且那个例子的盒子偏小，跑到后来都出 ...

谢谢卢老师。试过Octopus，结果和别人已发表的类似分子的结果偏差甚大，其通过CP2K得到。NWchem没接触过

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