with nσ (r) ≡ ψ † σ (r) ψσ (r) and σ ≡ -σ. In other words, n ion (r) expresses the density of lo-cal spin-singlet pairs (double-site occupancy). Profiles of n ion (r) are not presented explicitly as they are contained in the equivalent form as the second part of the wave function (5).