Topological invariants based on generalized position operators and application to the interacting Rice-Mele model

We discuss different properties and the ability of several topological invariants based on position operators to identify phase transitions, and compare with more accurate methods, like crossing of excited energy levels and jumps in Berry phases. The invariants have the form $\text{Im} \text{ln} \left\langle \exp \left[ i(2\pi /L)\Sigma _{j}x_{j} \left( m_{_{\uparrow }}\hat{n}_{j\uparrow } +m_{\downarrow }\hat{n}_{j\downarrow }\right) \right] \right\rangle$, where $L$ is the length of the system, $x_{j}$ the position of the site $j$, $\hat{n}_{j\sigma }$ the operator of the number of particles at site $j$ with spin $\sigma$. We show that $m_{\sigma }$ should be integers, and in some cases of magnitude larger than 1 to lead to well defined expectation values. For the interacting Rice-Mele model (which contains the interacting Su-Schrieffer-Heeger and the Ionic Hubbard model as specific cases), we show that three different invariants give complementary information and are necessary and sufficient to construct the phase diagrams in the regions where the invariants are protected by inversion symmetry. We also discuss the consequences for pumping of charge and spin, and the effect of an Ising spin-spin interaction or a staggered magnetic field.