I am trying to compute quadratic response properties after a GASCI calculation with .LUCITA. The CI converges and the CHECKPOINT.h5 is written correctly, but any attempt to activate *RESPONSE fails, either it is ignored or gives errors like "Illegal directory in PSIINP".
Is quadratic response from LUCITA wavefunctions currently implemented? If not, is it planned for future versions?
Thank you in advance!
Carles