Fig. presents the number of photons during the dynamics process. Note that â † â is not the "photon number" operator under the dipole gauge used in the PF Hamiltonian because the rigorous photon number operator should be obtained by applying the Power-Zienau-Woolley (PZW) Gauge transformation on the photon number operator â † â. The correct photon number operator under the dipole gauge is