Self-consistent determination of the electron velocity distribution function in molecular plasmas