I've been reading the paper by J. R. Cheeseman, M. J. Frisch, “Basis Set Dependence of Vibrational Raman and Raman Optical Activity Intensities”, J. Chem. Theory and Comput., 7, (2011), 3323-3334. Moreover, as reported in the gaussian09 manual: "Raman and ROA intensities can be calculated separately from calculation of the force constants and normal modes, to facilitate using a larger basis for these properties as recommended."
I'm performing a frequency calculation with B3LYP/maug-cc-pVDZ and I want to perform a subsequent Raman computation with B3LYP/aug-cc-pVTZ using the previous checkpoint file. I've been tinkering with the keyword Polar=Raman but the job fails, complaining about missing input.
Is there anybody who can provide any suggestion and possibly a route card example to solve this problem?