 Solution

The Hückel orbital charge density and Hückel pi bond order giveEach site receives $1/2$ electron from the lowest orbital and one from the occupied nonbonding subspace, so $n_i=3/2$ at every carbon, summing to six. Relative to one neutral pi electron per carbon this is a pi charge of $-e/2$ at each site.

For adjacent sites, only the lowest orbital contributes, giving $p_{12}=p_{23}=p_{34}=p_{41}=1/2$. For opposite sites the nonbonding contribution reverses the sign: $p_{13}=p_{24}=-1/2$. These non-neighbour negative density-matrix entries are not extra sigma bonds or extra resonance integrals in the Hamiltonian. If only skeletal nearest-neighbour pi bond orders are requested, the four values $1/2$ are the relevant ones.