关于在Gaussian16中结合BOMD和TD-DFT的问题

Kalinite

之前研究了一个金属有机配合物体系的非辐射过程的机理，其过程是S0激发至S1后，ISC至T2，IC至T1，在T1处与其他分子发生三重态能量转移。这个机理比较能够解释实验中的一些结果，暂且认为是合理的。但是不是很明白应该怎么研究IC过程，而且这一部分只是文章中的次要内容，于是用Gaussian16做了个BOMD和TD-DFT的结合，结果表明是比较容易到达T2 T1交叉点的，就这么写进文章了。
审稿人提了一个问题：
Won’t it be better to start with the S1 state populated post-excitation? Or possibly even start with the S0 state populated and study the excitation event under the pulse of light?
我的理解是：(1)从S1开始模拟，这涉及系间窜越必须考虑SOC；(2)从S0开始模拟，需要模拟光吸收过程，这似乎也没法实现，而且我的体系是通过多光子吸收到达激发态的。
故向大家提问，是否有程序能够实现这样的目的，比如Newton-X之类，但看了一些文献没见过这样的例子。