In this brief, a new method for extracting carrier mobility from the Monte
Carlo device simulation is presented. Unlike the conventional method which
is based on taking the time-average of the rate of change of momentum befor
e and after the scattering, this method utilizes the spherical harmonic exp
ansion of the carrier distribution function and an analytical expression fo
r the relaxation time. The new method produces a much smoother result than
the time-average method, particularly in low electric field regions.