计算化学公社

标题: CP2K跑周期性体系TDDFT/激发态效果如何 [打印本页]

作者
Author: kkyle 时间: 2021-7-2 21:00
标题: CP2K跑周期性体系TDDFT/激发态效果如何
各位大神：不知道大家有没有用CP2K跑过激发态吗，效果如何，我想看光激发下电子转移的方向，不知道可否实现。
谢谢！

作者
Author: sobereva 时间: 2022-3-1 15:21
按此文的做法，非常容易
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
http://sobereva.com/634（http://bbs.keinsci.com/thread-28006-1-1.html）

作者
Author: flowerstudy 时间: 2023-2-26 20:22

sobereva 发表于 2022-3-1 15:21
按此文的做法，非常容易
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
http://sober ...

卢老师，请问cp2k能做rt-tddft吗？cp2k能做tddft，又能做AIMD，两个结合起来不就是rt-tddft吗？

作者
Author: reid 时间: 2023-2-26 23:24

flowerstudy 发表于 2023-2-26 20:22
卢老师，请问cp2k能做rt-tddft吗？cp2k能做tddft，又能做AIMD，两个结合起来不就是rt-tddft吗？

我用orca跑过非周期性体系的激发态下的MD，cp2k没试过

作者
Author: sobereva 时间: 2023-2-27 03:25

flowerstudy 发表于 2023-2-26 20:22
卢老师，请问cp2k能做rt-tddft吗？cp2k能做tddft，又能做AIMD，两个结合起来不就是rt-tddft吗？

不是一码事
RT-TDDFT强调的是电子波函数随时间的演化，是对含时薛定谔方程积分，而且CP2K的RT-TDDFT是原子核坐标不变的情况下做的
一般意义的TDDFT结合AIMD是在激发态势能面下做的绝热动力学。CP2K结合Newton-X也可以做非绝热动力学。

作者
Author: flowerstudy 时间: 2023-2-28 08:48

sobereva 发表于 2023-2-27 03:25
不是一码事
RT-TDDFT强调的是电子波函数随时间的演化，是对含时薛定谔方程积分，而且CP2K的RT-TDDFT是原 ...

好的，谢谢sob老师的回复

作者
Author: flowerstudy 时间: 2023-2-28 08:49

reid 发表于 2023-2-26 23:24
我用orca跑过非周期性体系的激发态下的MD，cp2k没试过

谢谢您的回复。

作者
Author: cdw66666666 时间: 2023-3-3 20:04

sobereva 发表于 2022-3-1 15:21
按此文的做法，非常容易
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
http://sober ...

社长您好，我用cp2k跑周期性体系的TDDFT老是因为内存不够而结束进程，请问我该从哪个方面改进呢？之前计算一个周期性体系成功了，这次相比只是加了个34个原子的有机分子进到之前的体系了，应该不会计算量突然增大那么多吧

作者
Author: sobereva 时间: 2023-3-8 22:31

cdw66666666 发表于 2023-3-3 20:04
社长您好，我用cp2k跑周期性体系的TDDFT老是因为内存不够而结束进程，请问我该从哪个方面改进呢？之前计 ...

减少MPI进程数，或者保持并行核数相同情况下用psmp版结合一些OpenMP线程以省内存再试

作者
Author: cdw66666666 时间: 2023-3-11 09:42

sobereva 发表于 2023-3-8 22:31
减少MPI进程数，或者保持并行核数相同情况下用psmp版结合一些OpenMP线程以省内存再试

好的，谢谢sob老师

作者
Author: Cara 时间: 2023-7-14 17:32

sobereva 发表于 2023-2-27 03:25
不是一码事
RT-TDDFT强调的是电子波函数随时间的演化，是对含时薛定谔方程积分，而且CP2K的RT-TDDFT是原 ...

sob老师，请问我想用cp2k做TDDFT结合AIMD，也就是激发态的动力学，可以帮我介绍一些您已知的资料么？

作者
Author: sobereva 时间: 2023-7-15 08:50

Cara 发表于 2023-7-14 17:32
sob老师，请问我想用cp2k做TDDFT结合AIMD，也就是激发态的动力学，可以帮我介绍一些您已知的资料么？

直接做就完了

作者
Author: Cara 时间: 2023-7-28 21:09

sobereva 发表于 2023-7-15 08:50
直接做就完了

谢谢sob老师，我之前的描述有问题。确实使用multiwfn建立cp2k输入文件之后，MD和TDDFT一起是对每一步瞬态构型都分析TDDFT。我重新组织了一下问题，是想请教对一有机发光分子的S1单重激发态进行动力学模拟，来获取S1单重激发态在一段时长内随机畸变下的构型合集。整个过程都保持S1激发态的电子状态，不考虑surface hopping的问题。请问sob老师这个用cp2k是否可行？

作者
Author: sobereva 时间: 2023-7-28 22:41

Cara 发表于 2023-7-28 21:09
谢谢sob老师，我之前的描述有问题。确实使用multiwfn建立cp2k输入文件之后，MD和TDDFT一起是对每一步瞬态 ...

原则上可以

但对有机小分子这种问题更建议用Gaussian或ORCA做BOMD，效率比CP2K更高，数值精度更好

作者
Author: Cara 时间: 2023-7-28 23:09

sobereva 发表于 2023-7-28 22:41
原则上可以

但对有机小分子这种问题更建议用Gaussian或ORCA做BOMD，效率比CP2K更高，数值精度更好

谢谢sob老师，我去学习一下BOMD的相关。
我的分子是并四苯的三聚体分子，加上真实溶剂分子后原子数能达500，倾向使用cp2k进行S1激发态的动力学模拟，理论上是可以的，但是作为初学者查看cp2k官网的手册里MD部分，没有查到能对激发态设限的关键词，multiwfn构建输入文件也没有相关项。请教sob老师有没有办法限定对S1激发态进行动力学模拟？
此外，在吧内帖子中认识到Newton-X，擅长surface-hopping问题，想请教sob老师，如果cp2k不能实现，那这个软件能否满足我只对S1激发态进行动力学模拟获取畸变构型集（轨迹）的需求？

作者
Author: sobereva 时间: 2023-7-29 12:00

Cara 发表于 2023-7-28 23:09
谢谢sob老师，我去学习一下BOMD的相关。
我的分子是并四苯的三聚体分子，加上真实溶剂分子后原子数能达5 ...

打开TDDFT计算设置和动力学设置就完了

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