A selfconsistent calculation of the rigid neutral atom density according to the auxiliary neutral atom model