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

adreamy

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

wzkchem5

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

adreamy

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

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

wzkchem5

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

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

wal

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

adreamy

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

谢谢老师

adreamy

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

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

adreamy

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

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

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

谢谢老师，我在学习一下