gaussian和qchem如何计算非绝热电子耦合

杲杲出日

我想计算基态和激发态之间的非绝热电子耦合然后进一步计算非辐射复合速率，遇到了以下问题：

1. 在使用Gaussian计算非绝热耦合时，关键词设置为：# wb97xd/6-31G* NoSymm tda(NAC,nstates=3)iop(3/107=0002500000,3/108=0002500000,9/40=4) scrf(pcm,read)。log文件在最后出现了图1的结果，请问这就是S0和S1态之间的NAC吗？假如我想计算S0和S2之间的，关键词改为了 tda(NAC,root=2)，不知道这样是否正确。

2. 使用qchem-5.0.1计算时，设置如图2所示，想要在TDDFT下计算基态和激发态之间的NAC，但为何报错：Q-Chem fatal error occurred in module /scratch/svnadmin/18824_zgan/qchem/qparser/ParseQInput.C, line 1985:Please set CALC_NAC = 2 when SET_QUADRATIC is set to TRUE。明明我已经改了相关设置

3. 非绝热耦合元是波函数相对核坐标的一阶导数，请问应该如何计算才是耦合元的值呢？ 是Σ(xi2+yi2+zi2)^2 对所有原子坐标平方求和开根号，还是(Σxi)^2+(Σyi)^2+(Σzi)^2先对所有x、y、z的值进行加和再平方后开根号呢？

wal

理论上td(nac,root=2)算的就是S2/S0的nac，但是跟S1/S0的一样有点奇怪。你要不把输出传上来看看？
对xyz加和再平方听上去有点离谱，这是个向量，把xyz加一起感觉没啥物理意义，应该是平方和

杲杲出日

谢谢老师，我检查了下闹了乌龙，root=1和2的结果不一样。gaussian有没有办法像qchem一样一次性把S0和计算的所有激发态之间的NAC结果都输出呢？