In a previous work we have presented a numerical procedure for the calculation of the internal chemical hardness tensor at the molecular orbital resolution level from standard density functional calculations. In this article we describe an improvement of our method using the thermal extensions of density functional theory. Furthermore, new concepts are introduced in the orbitally resolved theory of chemical reactivity. Traditional molecular orbital theories of chemical reactivity are based only on considerations concerning the highest occupied molecular orbitals (HOMOs) and the lowest unoccupied molecular orbitals (LUIMOs) of molecules, supposed to describe the behavior towards electrophiles, respectively, nucleophiles. By applying our methodology to two test molecular systems, namely water and ferrocene, we show how chemical reactivity can be differentiated against hard and soft electrophiles (acids) and hard and soft nucleophiles (bases). As a by-product of the numerical algorithms being used, a self-consistent method for calculating the molecular chemical potential is also described.