Hi!

I try to calculate soc-perturbed plasma frequency using VASP. In order to do it, I added the following keyword to my INCAR:

LSORBIT = .TRUE.

LEPSILON = .TRUE

LRPA = .TRUE

LOPTICS = .TRUE. ; CSHIFT = 0.1

However, it seems that this setup doesn't work, since no plasma_frequency_squared is produced in the output. Is there any way to calculate it? Or soc-perturbed calculations are currently not implemented to rpa optics?