I have a Gromacs NVT simulation with a constant external electric field using the option "electric-field-z = 1.5 0 0 0".
When calling gmx energy -f (...).edr to calculate the potential energy of the system, is the energy due to interaction with the electric field included or not?
Particularly, I have a lot of molecular dipoles in the system which can align wrt to the electric field. I want to know if the dipole - field interaction energy is contained in the potential from gmx energy, and if not, how I can quantify it.