Quantum Integrability of 1D Ionic Hubbard Model

The quantum integrability of the 1D ionic Hubbard model (IHM) is established using two independent numerical methods, namely i) energy level spacing statistics and ii) occupation profile of one-particle density matrix (OPDM) eigen-values. Both methods suggest that the 1D IHM is integrable. The calculations of energy level statistics reproduce the known results for the standard Hubbard model. Upon turning on the the ionic term, the energy level spacing distribution of this model continues to obey the Poissonian distribution. Occupation patterns as extracted from OPDM indicate that quasi-particles are sharpened upon increasing the ionic potential. This is evidenced by a larger jump in the occupation number distribution.