I was wondering if there is a ready-to-use option (i.e. by setting some rem keywords) to let Q-Chem print out the E-field at all nuclear positions. The underlying functionalities must be available since those derivatives are needed for geometry optimization. I found ESP_EFIELD from the manual but that only seems to work when a plotting grid is given (not for nuclear positions).
If there is no such a ready-to-use option, I think I can enable that with a tiny modification in the code.