Grad's method of moments is employed to develop higher-order Grad moment equations-up to first 26-moments-for dilute granular gases within the framework of the (inelastic) Boltzmann equation. The homogeneous cooling state of a freely cooling granular gas is investigated with the Grad 26-moment equations in a semi-linearized setting and it is shown that the granular temperature in the homogeneous cooling state still decays according to Haff's law while the other higher-order moments decay on a faster time scale. The nonlinear terms of fully contracted fourth moment are also considered and, by exploiting the stability analysis of fixed points, it is shown that these nonlinear terms have negligible effect on Haff's law. Furthermore, an even larger Grad moment system which includes the fully contracted sixth moment is also scrutinized and the stability analysis of fixed points is again exploited to conclude that even the inclusion of scalar sixth order moment into the Grad moment system has negligible effect on Haff's law. The constitutive relations for the stress and heat flux (i.e., the Navier-Stokes and Fourier relations) are derived through the Grad 26-moment equations and compared with those obtained via CE expansion and via computer simulations. The linear stability of the homogeneous cooling state is analyzed through the Grad 26-moment system and various sub-systems by decomposing them into longitudinal and transverse systems. It is found that one eigenmode in both longitudinal and transverse systems in the case of inelastic gases is unstable. By comparing the eigenmodes from various theories, it is established that the 13-moment eigenmode theory predicts that the unstable heat mode of the longitudinal system remains unstable for all wavenumbers below a certain coefficient of restitution while any other higher-order moment theory shows that this mode becomes stable above some critical wavenumber for all values of the coefficient of restitution. In particular, the Grad 26-moment theory leads to a smooth profile for the critical wavenumber in contrast to the other considered theories. Furthermore, the critical system size obtained through the Grad 26-moment and existing theories are also in excellent agreement.