I'm trying to calculate of the Duschinsky rotation matrix between two oxidation states of a cobalt (II and III) complex. I'm using G09. I used freq=READFCHT in comment route and PrtMat=1 and JDusch in the end of input file but I think G09 didn't read this comments and I just have a regular freq output file without any Duschinsky rotation matrix.