A Monte Carlo simulation has been developed for electron thermalization in a liquid hydrocarbon taking n-hexane as a paradigm. For the first time, both the effect of the geminate cation and that of an external electric field have been included in the predicted electron thermalization distance distribution (ETDD). Also for the first time, both energy losing (direct) and energy gaining (inverse) collisions of the epithermal electron are included in a manner that ensures the equality of their rates, averaged over the Maxwellian distribution, in the long-time limit. Parametric values have been rationalized from the experimental data. Comparison with a free-ion yield experiment in n-hexane at 300 K indicates that about one-third of total epithermal collisions involve the energy exchange of quanta related to intermolecular motion. These energy exchange collisions partition between energy losing to energy gaining processes approximately as 3:2 in the initial stage of thermalization. The important features of the ETDD are displacement from the origin and an exponential tail for which rationalization has been provided. Distortion of the ETDD by an external field does not affect the free-ion yield in the first order; in the second order, an effect quadratic in the field may be seen.