求助gaussian产生的*.fchk文件中的Polarizability Derivatives数据的解析

Hamilton

各位老师好，我想编写一个通过读取gaussian产生的fchk文件并用于计算Izz（拉曼张量Rzz的平方）的简单脚本，fchk文件中有Polarizability与Polarizability Derivativities等数据，我的体系共58个原子也就是58*3=174个位移自由度与174-6=168个振动自由度，fchk文件给出了1044=174*6个极化率导数（Polarizability Derivativities）数据，如最下图所示，我觉得他代表的是极化率张量对原子位移的导数。根据某个振动模式的拉曼张量等于极化率张量对该振动模式简正坐标的导数，同时其拉曼活性正比于其拉曼张量的平方和即 ，其中α为拉曼张量，Qk为振动模式k的简正坐标，我想请问各位老师如何将fchk文件中输出的Polarizability Derivativities转化为拉曼张量呢？