Numerical simulation of far infrared emission under population inversion of hole sub-bands