An auxiliary field approach for computing optical resonances in dispersive media