In a charge-stabilized colloidal solution, the large colloidal
particles are surrounded by microions that are up to four orders
of magnitude smaller than the colloidal particles. Because of this size
asymmetry, it is desirable to obtain an effective one-component
description of the mixture where the colloidal particle plus its
ionic atmosphere is treated as one, dressed particle. The effective
pair potential between these dressed particles is a screened
Coulomb potential. The screening depends, of course, on the
density distribution of the small ions around and between the big
colloidal particles. If the colloidal charge and the concentration
of the ions is not too high, this distribution can be
approximately determined from the linearized Poisson-Boltzmann
equation, and the resulting effective pair potentials are Yukawa
potentials. In concentrated suspensions, however, the full,
non-linear Poisson-Boltzmann equation must be solved to determine
the density distribution of the small ions. In this article, we
suggest a way to obtain effective pair potentials for
this case. We solve the non-linear Poisson-Boltzmann equation
around a colloidal particle that is displaced a certain distance
from the centre of its Wigner-Seitz cell. From the resulting
density profile of the ions, we determine the total force acting
on the shifted particle as a function of the displacement. From
this function one can then estimate the non-linearly
screened pair forces, and, thus, the effective pair potentials.