The above PF Hamiltonian has the advantage as a pure real Hamiltonian and the photonic DOF can be viewed and computationally treated as an additional "nuclear coordinate". We emphasize that both the operators as well as the wavefunctions should be Gauge transformed (though Û ), in order to have a Gauge invariant expectation values. This means that