Thermodynamic stability of interacting fermions system with matrix eigenvalue method

A matrix eigenvalue method is applied to analyse the thermodynamic stability of two-component interacting fermions. The non-relativistic and ultra-relativistic d=1, 2,3 dimensions have been discussed in detail, respectively. The corresponding stability region has been given according to the two-body interaction strength and the particle number density ratio.