如何从gaussian log文件中读取跃迁偶极矩梯度（Electronic Transition Derivatives）

pendrangon

请教各位老师：我在尝试编写软件从gaussian的log文件中读取一些信息来计算共振拉曼光谱，在计算Herzberg-teller项时需要用到跃迁偶极矩的梯度：

以酞菁分子为例，我使用 #p freq wB97XD def2TZVP iop(7/33=1) td(root=1,nstat=40)对其在基态结构下进行了tddft计算，我目前知道的时跃迁偶极矩梯度信息在这一块：


但这一块共有16行3N+3列（N为原子数），我想知道这其中的每一行的和每一列的物理意义以及单位。请各位老师赐教。


附上结构：


此外如果方便的话请顺便告诉我跃迁偶极矩这一块的意义（这里40列应该就是nstat的值，16行应该和跃迁偶极矩块一样，但我不知道怎么对应的）：

感谢

sobereva

没注明单位几乎一定是原子单位

这种情况可以找个原子数很少的简单体系，自己通过有限差分求导数，和截图里的进行对照判断输出顺序。要么读源代码

pendrangon

谢老师指导，我用H2O（#p freq b3lyp/6-31g iop(7/33=1) td(root=1,nstat=20) nosymm）试了一下，取dx=0.01A的情况下做了跃迁偶极矩的有限差分，与” Electronic Transition Derivatives“部分的数据对比了一下， Electronic Transition Derivatives部分数据应该是按照如下方式摆放的：
该部分共有16行3N+3列，对于行，第一行占位，恒为零，第2，3，4行表示跃迁偶极矩的x,y,z分量，后面的行目前不清楚，应该是更高阶的梯度。对于列，前三列不知到是什么，从第4列开始依次为第一个原子的x,y,z，（即4,5,6列表示跃迁偶极矩对第一个原子的x,y,z坐标的梯度，以此类推后面共3N列）

附上差分计算的跃迁偶极矩梯度与 Electronic Transition Derivatives部分的数据对比：