The First and second rungs of the Jacob's ladder, i.e., LDA and GGA method, are typically employed in the solid-state calculations while third and fourth rungs, e.g., Meta-GGA and hybrid functional are generally used in the quantum chemistry simulations for the prediction of molecular properties . The meta-GGA method is a semi-local approximation which is developed through the extension of GGA with consideration of the second derivative of the density . As it comprises similar physical information and higher stability, generally orbital kinetic energy density is used at the place of the second derivative of density. The mathematical expression of meta-GGA is given as;