|
发布时间: 2024-10-28 21:49
正文摘要:请问各位老师，像这种非绝热振动耦合NACME的值与频率关系的图可以怎么得到？一个分子的频率不是3N-5或者3N-6个吗？NACME是3N个吧，怎么能对应上的，求助。 |
wzkchem5 发表于 2024-11-1 11:18 谢谢老师，我在学习一下 |
adreamy 发表于 2024-11-1 03:14 这里分成两个问题，（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计算 |
wzkchem5 发表于 2024-10-29 16:27 老师，再请教您一个问题，就是如果我想计算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？或者说我考虑的都不对，请老师指教。 |
wal 发表于 2024-10-29 16:54 https://doi.org/10.1038/s41566-022-01079-8 |
wzkchem5 发表于 2024-10-29 16:27 谢谢老师 |
| 楼主可以分享一下这篇文章嘛？没见过作这种图的，想看看 |
adreamy 发表于 2024-10-29 07:34 按计算IC速率的方法写输入文件，然后在ESD模块里调大printlevel |
wzkchem5 发表于 2024-10-29 10:08 老师，要orca输出这个信息除了printlevel之外，还需要什么关键词吗 |
|
这个是把NACME由笛卡尔基转到简正坐标基，然后把每个简正坐标对NACME的贡献对简正坐标的波数作图。 ORCA可以输出这个信息，如果默认的printlevel不够的话就把printlevel加大；MOMAP应该也可以做 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 05:31 , Processed in 0.438346 second(s), 27 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites