Plasma-neutral interaction in kinetic models for the divertor region