I searched and found out that the thermal broadening is given as 3bkT/(4a2)3bkT/(4a^2) , but this way need some data to find out aa and bb , as the aa and bb is just a number given by a model. I can't find a clue to find the thermal broadening just with Fermi distribution function and without any experimental data. I am working on this problem by trying to integrate the Fermi distribution function, as it seems to be the only thing I can do with this function, but that's all I could do.