Often, exotic phases appear in the phase diagrams between conventional phases. Their elementary excitations are of particular interest. Here, we consider the example of the ionic Hubbard model in one dimension. This model is a band insulator (BI) for
weak interaction and a Mott insulator (MI) for strong interaction. Inbetween, a spontaneously dimerized insulator (SDI) occurs which is governed by energetically low-lying charge and spin degrees of freedom. Applying a systematically controlled version of the continuous unitary transformations (CUTs) we are able to determine the dispersions of the elementary charge and spin excitations and of their most relevant bound states on equal footing. The key idea is to start from an externally dimerized system using the relative weak interdimer coupling as small expansion parameter which finally is set to unity to recover the original model.
A detailed study of the one-dimensional ionic Hubbard model with interaction $U$ is presented. We focus on the band insulating (BI) phase and the spontaneously dimerized insulating (SDI) phase which appears on increasing $U$. By a recently introduced
continuous unitary transformation [Krull et al. Phys. Rev. B {bf 86}, 125113 (2012)] we are able to describe the system even close to the phase transition from BI to SDI although the bare perturbative series diverges before the transition is reached. First, the dispersion of single fermionic quasiparticles is determined in the full Brillouin zone. Second, we describe the binding phenomena between two fermionic quasiparticles leading to an $S=0$ and to an $S=1$ exciton. The latter corresponds to the lowest spin excitation and defines the spin gap which remains finite through the transition from BI to SDI. The former becomes soft at the transition indicating that the SDI corresponds to a condensate of these $S=0$ excitons. This view is confirmed by a BCS mean field theory for the SDI phase.
Effects of truncation in self-similar continuous unitary transformations (S-CUT) are estimated rigorously. We find a formal description via an inhomogeneous flow equation. In this way, we are able to quantify truncation errors within the framework of
the S-CUT and obtain rigorous error bounds for the ground state energy and the highest excited level. These bounds can be lowered exploiting symmetries of the Hamiltonian. We illustrate our approach with results for a toy model of two interacting hard-core bosons and the dimerized S=1/2 Heisenberg chain.
