This paper describes the present state of the development of the computational model of the ITER edge plasma. Neutral–neutral collisions and molecular dynamics are introduced into the self-consistent scheme. First results of ITER modelling including these effects indicate that the operational window for the ITER divertor shifts towards higher neutral pressure in the private-flux region, retaining the operational flexibility determined in the previous analyses.