In Q-chem version 5.4, when attempting to perform a calculation with $external_charges in NEO-DFT, the calculation was reflected in the general SCF, but ignored in the NEO-SCF. Since jobtype=polarizability allowed the calculation of polarisation rates reading the NEO code, I thought it would be possible to include external charges, but is this not possible? If this is possible through a higher version, I would like to know about it too.