p-Type NbFeSb-based half-Heusler compounds possess good thermoelectric (TE) properties with the benefits of abundantly accessible constituent elements, a high Seebeck coefficient, and a large power factor (PF) but a relatively high κ l . Reducing the κ l is crucial to increase the figure of merit (zT). This is generally done with isovalent doping, alloying, nanostructuring, etc., but these adversely affect carrier mobility. We have performed transport calculations on NbFeSb and NbFeBi to see the same-group substitution effects on thermal conductivity. A usual way to study the phonon transport properties is to solve the Boltzmann transport equation (BTE) within the firstprinciples approach, which has proven experimental reproducibility. The doped or nanostructured system is computationally demanding with the above approach and appeals for more convenient access, such as the deformation potential (DP) method. We study the κ l of pristine NbFeSb and NbFeBi with the first-principles-based BTE and compare it with the DP model. Bi substitution at the Sb site in NbFeSb reduces the κ l (from 4.5 to 2 W/mK at 1000 K) without any complication of nanostructuring/doping and has the same high band degeneracy at the band edges. We also benchmark the DP model prediction so that further doping with aliovalent atoms can be studied with the DP model to obtain the high zT value.