I'm studying linear response theory from "Fundamentals of many body physics" by Nolting. We want to study electrical AC conductivity, so what happens when we apply a spatially homogeneous oscillating electric field to a system of charged particles.
Nolting starts with the assumption that the presence of the electric field $\vec E(t)$ causes an additional term in the hamiltonian given by $$V(t)=-\vec P\cdot\vec E(t)$$ where $\vec P$ is the electric dipole moment operator $$\vec P=\sum_iq_i\hat{\vec r_i}$$ where $q_i$ and $\hat{\vec r_i}$ are the charge and the position operator of the i-th particle.
I know in classical physics that this is the form of the interaction energy of an electric dipole with an external electric field, but I don't understand how to get this in quantum mechanics.
If the hamiltonian of the i-th particle is $$H_i=\dfrac{1}{2m}\left( \vec p_i-q_i\vec A(\vec r_i,t) \right)^2$$ to have an alternating omogeneous electric field we can choose the potential $$\phi=\text{electrostatic potential energy between the particles}$$ $$\vec A(\vec r_i,t)=\vec A_0\cos(\omega t)$$ if we expand we have, assuming all $q_i=e$ $$H_i=\dfrac{p_i^2}{2m}+\dfrac{e^2}{2m}A(\vec r_i)^2-\dfrac{e}{2m}\vec p_i\cdot\vec A(\vec r_i,t)-\dfrac{e}{m}\vec A(\vec r_i,t)\cdot \vec p_i$$
where I considered the fact that momentum and vector potential do not commute. The third term is zero if the vector potential is homogeneous and maybe we could neglect the second term if we are considering a small electric field, but the remaining term does not look like the one Nolting introduces. Here I find that the vector potential couples with the total momentum of the system, and not that the electric field couples with the positions.
How can we arrive to additional term in the hamiltonian given by $$V(t)=-\vec P\cdot\vec E(t)$$??