Dear Jixuan Wu,
For applying E-field in DFT optimization (or in ab initio MD) one would also use the BornEffectiveCharges (always computed from DFT) to compute the impact of E-field on DFT energy and forces. This is something we are working on, but we do not yet have a working prototype.
Regards,
Jess