计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2000|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 用Multiwfn辅助生成CP2K输入脚本做S1和T1激发态动力遇到无法产生SOMO轨道的问题

[复制链接 Copy URL]

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Cara 于 2024-3-7 00:25 编辑

众所周知,S0态和S1态由于自选多重度均为1,在进行AIMD计算时仅设置自旋多重度为1会回落S0,所以无法进行S1激发态动力学获取一条瞬态结构均为S1态可能构型合集的xyz轨迹和相关MO的cube文件们。
之前请教过sob老师能够算有机分子光激发后S1态的激发态动力学,老师指导我说:“动力学中用的受力是S1的自然就是S1态的AIMD,显然让AIMD的程序计算激发态、指定算哪个态的受力就完了。但凡程序支持BOMD、支持TDDFT计算受力,都能干这个。Gaussian、ORCA、CP2K全都能干”,
于是我用Multiwfn辅助得到了CP2K的输入脚本,对于T1激发态,在自旋多重度为3的情况下,用EXCITED_STATES关键词设定用第1个激发态的受力(如图所示),并且加入了TDDFPT部分,其中RKS_TRIPLETS 为T,按备注含义是计算三重态。需要输出结构轨迹xyz和3个HOMO3个LUMO的cube,理论上应该包含SOMO。
但这个脚本计算完成后,MO的cube并不包含SOMO以及变为开壳层,而是类似闭壳层基态的HOMO LUMO状态。
同理对自旋多重度为1也设计了S1态激发态动力学的输入文件,结果MO也是没有SOMO,为闭壳层的HOMO LUMO。
这种情况是否说明我的受力设置是失效的?
我的课题急需得到有机分子光激发后大量处在S1态的瞬态结构们,以及他们的SOMO轨道。在这里我诚心请教sob老师和各位老师们,我的错误在哪以及正确的操作方法。输入脚本已附上。
自知学识太少,果断报名了本届CP2K的培训班,奈何个人原因只能选择不到场的方式。这个课题问题也比较紧迫,如果需要别的软件才能有效解决这个需求,也感谢各位老师不吝赐教,我愿意去学习。



Excited.png (52.89 KB, 下载次数 Times of downloads: 9)

Excited.png

TDDFPT.png (74.71 KB, 下载次数 Times of downloads: 8)

TDDFPT.png

S1.inp

4.46 KB, 下载次数 Times of downloads: 8

T1.inp

4.46 KB, 下载次数 Times of downloads: 5

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2024-3-8 09:02:48 | 只看该作者 Only view this author
sobereva 发表于 2024-3-8 05:31
我用2024.1跑了几步,没任何问题

谢谢sob老师百忙中测试,我去换用2024.1版本

6万

帖子

99

威望

6万

eV
积分
125130

管理员

公社社长

9#
发表于 Post on 2024-3-8 05:31:25 | 只看该作者 Only view this author
Cara 发表于 2024-3-8 00:34
想再请教一下sob老师,顺利用Multiwfn编写CP2K输入文件后,在计算过程中,对角化和OT收敛方式出现了如图 ...

我用2024.1跑了几步,没任何问题

DIAGONALIZATION.out (39.23 KB, 下载次数 Times of downloads: 7)

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2024-3-8 00:34:54 | 只看该作者 Only view this author
sobereva 发表于 2024-3-7 14:16
Multiwfn创建CP2K输入文件的界面里,把任务设为MD,打开TDDFT计算设置,20 Choose the state to evaluate ...

想再请教一下sob老师,顺利用Multiwfn编写CP2K输入文件后,在计算过程中,对角化和OT收敛方式出现了如图两种报错,查阅了部分报错分享贴仍没有真正理解,请问老师如何去解决这两种情况?

DIAGONALIZATION.png (36.77 KB, 下载次数 Times of downloads: 5)

DIAGONALIZATION.png

OT.png (28.92 KB, 下载次数 Times of downloads: 5)

OT.png

OT.inp

12.38 KB, 下载次数 Times of downloads: 12

OT.out

25.08 KB, 下载次数 Times of downloads: 2

DIAGONALIZATION.out

22.8 KB, 下载次数 Times of downloads: 1

DIAGONALIZATION.inp

12.24 KB, 下载次数 Times of downloads: 8

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2024-3-7 15:11:05 | 只看该作者 Only view this author
sobereva 发表于 2024-3-7 14:16
Multiwfn创建CP2K输入文件的界面里,把任务设为MD,打开TDDFT计算设置,20 Choose the state to evaluate ...

好的,谢谢sob老师,我明白了。

6万

帖子

99

威望

6万

eV
积分
125130

管理员

公社社长

6#
发表于 Post on 2024-3-7 14:16:47 | 只看该作者 Only view this author
Cara 发表于 2024-3-7 14:09
谢谢sob老师指教,我会注意措辞严谨。
sob老师能提示我Multiwfn编写对S1态跑AIMD的CP2K输入脚本的关键词 ...

Multiwfn创建CP2K输入文件的界面里,把任务设为MD,打开TDDFT计算设置,20 Choose the state to evaluate force设1就完了。受力是S1态的,自然MD就是在S1态势能面上采样。即便你用Gaussian、ORCA,也都是同样的理,对哪个电子态算受力,(绝热)动力学就是对哪个态跑的。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-3-7 14:09:06 | 只看该作者 Only view this author
sobereva 发表于 2024-3-7 12:34
没有所谓的“S0的结构合集”的概念
S1和S0的势能面不同,你对S1态跑AIMD,自然是对S1态的势能面进行采样 ...

谢谢sob老师指教,我会注意措辞严谨。
sob老师能提示我Multiwfn编写对S1态跑AIMD的CP2K输入脚本的关键词么?
之前阅读Multiwfn写CP2K脚本的博文,加上自己的探索,仍然不能确定EXITED STATE和TDDFPT部分的加入是否实现了对S1态势能面采样。得到的xyz轨迹受导师质疑,而我又无法从能量或者其他角度验证这是S1势能面的采样,不是S0势能面的采样。
加上CP2K可能不是进行S1的AIMD的主流软件,也找不到相关的文献证明我的计算方法和结果是合理。
这个实操技术的问题困扰我很久。sob老师可能很少进行这么低级的细节操作,不过还是很期待老师的知识储备能点醒我问题的关键

6万

帖子

99

威望

6万

eV
积分
125130

管理员

公社社长

4#
发表于 Post on 2024-3-7 12:34:23 | 只看该作者 Only view this author
Cara 发表于 2024-3-7 09:25
谢谢Sob老师,我明白了空穴、电子分布对我更有实际意义。已经开始尝试把博文中CP2K结合Multiwfn用TDD ...

没有所谓的“S0的结构合集”的概念
S1和S0的势能面不同,你对S1态跑AIMD,自然是对S1态的势能面进行采样,轨迹会较为倾向于在S1态的势阱里晃悠

“单纯用于统计式的分析S1态和S0态在结构上弛豫前后的差异” 病句
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

24

帖子

0

威望

287

eV
积分
311

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-3-7 09:25:13 | 只看该作者 Only view this author
sobereva 发表于 2024-3-7 03:49
做T1态AIMD一般应当用UKS方式计算,不牵扯到TDDFT,直接把自旋多重度设为3就完了。UKS计算没有严格的SOMO的 ...

    谢谢Sob老师,我明白了空穴、电子分布对我更有实际意义。已经开始尝试把博文中CP2K结合Multiwfn用TDDFT计算考察电子激发态的电子-空穴时的单点计算,改为MD计算了,期望对瞬态构型们均做TDDFT计算,获取多个molden文件。后期继续补习初中级量子培训班的基础。
    目前仍然有两个疑问想要请教:(1)CP2K中用TDDFT算S1的AIMD不额外产生一套轨道,仍是基态S0的轨道这一点是很好接受的,那么xyz轨迹文件其实也是S0的结构合集对么?(2)我如何才能只获取S1态势能面下的结构轨迹,单纯用于统计式的分析S1态和S0态在结构上弛豫前后的差异?(不想优化破坏瞬态构型,那样失去了动力学模拟的真实性)

6万

帖子

99

威望

6万

eV
积分
125130

管理员

公社社长

2#
发表于 Post on 2024-3-7 03:49:47 | 只看该作者 Only view this author
做T1态AIMD一般应当用UKS方式计算,不牵扯到TDDFT,直接把自旋多重度设为3就完了。UKS计算没有严格的SOMO的概念(除非之后用Multiwfn做双正交化,http://sobereva.com/448),ROKS才有,但ROKS在CP2K里很不好用。看下面的北京科音初级量子化学培训班(http://www.keinsci.com/workshop/KEQC_content.html)的ppt了解常识(PS:这些最基本程度的量化知识都是初级班讲的,不在CP2K培训班里重复)



CP2K里S1态AIMD一般需要用TDDFT算,而TDDFT算激发态过程并不额外产生一套轨道。激发态轨道这事老有人问,所以在北京科音基础(中级)量子化学培训班(http://www.keinsci.com/workshop/KBQC_content.html)里我专门做了张ppt:



另外,TDDFT算的激发态也绝对不是仅仅能简单描述为一个占据轨道的电子跃迁到一个空轨道上的,看下文了解常识
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230

总之,你对SOMO的认识存在明显误区。很可能按下文得到各个结构下的空穴、电子分布对你更有实际意义
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434http://bbs.keinsci.com/thread-10775-1-1.html
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
http://sobereva.com/634http://bbs.keinsci.com/thread-28006-1-1.html




北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 05:46 , Processed in 0.178021 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list