Ionization and recombination of plasma impurities govern their charge state kinetics, which is imposed upon the dynamics of ions that implies a superposition of the appropriate probabilities and causes an impurity charge state dynamics. The latter is considered in terms of a vector field of conditional probabilities and presented by a vector charge state distribution function with coupled equations of the Kolmogorov type. Analytical solutions of a diffusion problem are derived with the basic spatial and temporal dimensionless parameters. Analysis shows that the empirical scaling DA∝ne−1 [K. Krieger, G. Fussmann, and the ASDEX Upgrade Team, Nucl. Fusion 30, 2392 (1990)] can be explained by the ratio of the diffusive and kinetic terms, DA∕(nea2), being used instead of diffusivity, DA. The derived time scales of charge state dynamics are given by a sum of the diffusive and kinetic times. Detailed simulations of charge state dynamics are performed for argon impurity and compared with the reference modeling.