where the resolution of the identity is satisfied by the complete electronic Hilbert space. Following the discussion by Huo and co-workers, we can think of the truncated version of Eq. 1 as arising from transformation with a projected Power-Zienau-Woolley (PZW) operator from the minimal coupling Hamiltonian 57 . The projected PZW operator has the form