非绝热振动耦合NAC与频率是什么关系，求助

adreamy · 发表于 Post on 2024-10-28 21:49:49

请问各位老师，像这种非绝热振动耦合NACME的值与频率关系的图可以怎么得到？一个分子的频率不是3N-5或者3N-6个吗？NACME是3N个吧，怎么能对应上的，求助。

wzkchem5 · 发表于 Post on 2024-10-29 10:08:01

这个是把NACME由笛卡尔基转到简正坐标基，然后把每个简正坐标对NACME的贡献对简正坐标的波数作图。
ORCA可以输出这个信息，如果默认的printlevel不够的话就把printlevel加大；MOMAP应该也可以做

adreamy · 发表于 Post on 2024-10-29 14:34:33

wzkchem5 发表于 2024-10-29 10:08
这个是把NACME由笛卡尔基转到简正坐标基，然后把每个简正坐标对NACME的贡献对简正坐标的波数作图。
ORCA可 ...

老师，要orca输出这个信息除了printlevel之外，还需要什么关键词吗

wzkchem5 · 发表于 Post on 2024-10-29 16:27:18

adreamy 发表于 2024-10-29 07:34
老师，要orca输出这个信息除了printlevel之外，还需要什么关键词吗

按计算IC速率的方法写输入文件，然后在ESD模块里调大printlevel

wal · 发表于 Post on 2024-10-29 16:54:11

楼主可以分享一下这篇文章嘛？没见过作这种图的，想看看

adreamy · 发表于 Post on 2024-10-29 17:35:32

wzkchem5 发表于 2024-10-29 16:27
按计算IC速率的方法写输入文件，然后在ESD模块里调大printlevel

谢谢老师

adreamy · 发表于 Post on 2024-10-29 17:37:14

wal 发表于 2024-10-29 16:54
楼主可以分享一下这篇文章嘛？没见过作这种图的，想看看

https://doi.org/10.1038/s41566-022-01079-8

adreamy · 发表于 Post on 2024-11-1 10:14:56

wzkchem5 发表于 2024-10-29 16:27
按计算IC速率的方法写输入文件，然后在ESD模块里调大printlevel

老师，再请教您一个问题，就是如果我想计算T2到T1的NACME（就用ESD（IC）模块计算），那么我需要怎么输入hess文件？我现在想的是，我以T1为参考态，以自旋多重度设为3，基态的hess文件就用T1计算得到的文件，那么T2的hess，我应该是以T1为参考态优化T2得到的hess，还是可以在S0坐标下用TDDFT优化计算得到的T2 hess文件？我见论坛有讨论自旋污染问题，我以T1为参考态自旋多重度设为3，关键词为td=(triplets,root=1）得到的第一个激发态能量是1.129eV，并且<S**2>=2.099和把坐标设为优化后的T1的坐标，自旋多重度设为1，关键词为td=(triplets,root=2,nstates=6)计算得到的第二激发态和第一激发态的能级差为1.0831eV，<S**2>=2.000，这样是不是我可以用TDDFT计算得到T2 hess去代替以T1为参考态优化计算得到的T2 hess？或者说我考虑的都不对，请老师指教。

wzkchem5 · 发表于 Post on 2024-11-1 11:18:34

adreamy 发表于 2024-11-1 03:14
老师，再请教您一个问题，就是如果我想计算T2到T1的NACME（就用ESD（IC）模块计算），那么我需要怎么输入 ...

这里分成两个问题，（1）以哪个结构为初始结构优化？（2）以哪个态（S0还是T1）为参考态？这两个问题是独立的。
先说（1），如果体系只存在一个T2构象，那么初始结构不重要；如果体系存在不止一个T2构象，需要估计不同T2构象之间的转变是否快于T2->T1 IC，如果是的话，需要搜索得到T2的最稳定构象；如果不是，需要知道T2是怎么来的，例如如果是从S1发生ISC而来，就以S1为初始结构优化。
再说（2），抛开程序限制的话，主要需要考虑的问题是哪个方法能把T2的电子结构描述得更准确。大体来说需要考虑两个问题：(a) S0和T1哪个参考态更容易用DFT描述，例如可以从多参考态性质的强弱来判断；(b) 以T1为参考态时，T2的自旋污染是否严重；(c) T2的双激发性质有多强。
如果限定用ORCA计算的话，因为ORCA不支持以S0为参考态算T1和T2的NACME，所以NACME只能以T1为参考态计算，这时T2的Hessian用哪个级别计算的问题还需要考虑一个因素：(d) 以S0为参考态的T2和以T1为参考态的T2的成分是否相符，如果不符的话NACME和Hessian的参考态应当一致。但如果用BDF计算的话，不受这个限制，可以以S0为参考态做所有T2计算

adreamy · 发表于 Post on 2024-11-1 12:43:12

wzkchem5 发表于 2024-11-1 11:18
这里分成两个问题，（1）以哪个结构为初始结构优化？（2）以哪个态（S0还是T1）为参考态？这两个问题是独 ...

谢谢老师，我在学习一下

phoenix · 发表于 Post on 2024-12-30 22:00:58

wzkchem5 发表于 2024-10-29 10:08
这个是把NACME由笛卡尔基转到简正坐标基，然后把每个简正坐标对NACME的贡献对简正坐标的波数作图。
ORCA可 ...

wzkchem5老师，您好，我想请教一下，我用VASP算频率得到Hessian矩阵之后，怎么将其转化为Normal mode呢？因为我想看一下不同频率之间的振动耦合。十分感谢老师的解答！

wzkchem5 · 发表于 Post on 2024-12-30 22:04:38

phoenix 发表于 2024-12-30 15:00
wzkchem5老师，您好，我想请教一下，我用VASP算频率得到Hessian矩阵之后，怎么将其转化为Normal mode呢？ ...

算phonon就行了，具体方法在本论坛第一性原理板块不难搜到。
但仅知道normal mode，还不能得出“不同频率之间的振动耦合”，因为在谐振近似下normal mode彼此之间没有耦合，需要做非谐振计算

phoenix · 发表于 Post on 2024-12-30 22:11:11

wzkchem5 发表于 2024-12-30 22:04
算phonon就行了，具体方法在本论坛第一性原理板块不难搜到。
但仅知道normal mode，还不能得出“不同频 ...

十分感谢w老师您的解答，那我先去论坛中找一下相关的计算方法。我还想请教一个问题：那非谐振计算之后的非对角元是不是就是不同振动模式之间的耦合呀？请w老师赐教

好兄弟 · 发表于 Post on 2025-1-13 10:48:24

wzkchem5 发表于 2024-10-29 16:27
按计算IC速率的方法写输入文件，然后在ESD模块里调大printlevel

老师好，ORCA6.0的 6.18.5 Internal Conversion Rates (unpublished)这部分没有ESD(IC)的写法，请问是怎么写kic的输入文件呢？

wzkchem5 · 发表于 Post on 2025-1-13 17:25:51

好兄弟发表于 2025-1-13 03:48
老师好，ORCA6.0的 6.18.5 Internal Conversion Rates (unpublished)这部分没有ESD(IC)的写法，请问是怎 ...

最新版手册写了
https://www.faccts.de/docs/orca/ ... ts/typical/esd.html

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

[综合交流] 非绝热振动耦合NAC与频率是什么关系，求助

浏览过的版块