Use: CAM-B3LYP/aug-cc-pVTZ TD(root=x) polar, where x is the index of the excited state of interest.

I'm trying to calculate the change in polarizability between the first excited state and ground state of an organic molecule.

My first thought was to calculate the polarizability of the molecule in the GS and S₁ using the polar keyword in Gaussian, and then compute the difference. However, Gaussian doesn't seem to support polar calculations for excited states with Density=(Check,Transition=1).

My question: Can I use Multiwfn to calculate the polarizability directly from the excited state density obtained from Gaussian TD-DFT? If so, what's the recommended workflow?

Thank you!
Dan