with \(f^{(0)}\)  as the equilibrium distribution, \(\phi \)  as the out of equilibrium distribution deviation, \(\tau _R^0\)  as the momentum independent part and \(n\)  as a number specifying the power of the scaled energy.
In order to solve Eq.(REF ), here the well known iterative technique of gradient expansion, the Chapman-Enskog (CE) method has been adopted [1]}. With this, the first order corrections in energy
density, particle number density, pressure, particle flux and energy flow are respectively given by,
\(&\delta \epsilon ^{(1)}=c_{\Lambda }^1(\partial \cdot u)~,~~\delta {\rho }^{(1)}=c_{\Gamma }^1(\partial \cdot u)~,~~\delta P^{(1)}=c_{\Omega }^1(\partial \cdot u)~,\nonumber \\&W^{(1)\alpha }=-c_{\Sigma }\hat{h}(\nabla ^{\alpha }T/T-Du^{\alpha })~,~~V^{(1)\alpha }=c_{\Xi }(\nabla ^{\alpha }\tilde{\mu }).\) 
