Charge dynamics in the Mott insulating phase of the ionic Hubbard model

We extend to charge and bond operators the transformation that maps the ionic Hubbard model at half filling onto an effective spin Hamiltonian. Using these operators we calculate the amplitude of the charge density wave in different dimensions. In one dimension, the charge-charge correlations at large distance d decay as 1/(d^3 ln^{3/2}d), in spite of the presence of a charge gap, as a consequence of remaining charge-spin coupling. Bond-bond correlations decay as (-1)^d 1/(d ln^{3/2}d) as in the usual Hubbard model.