In this paper, we study the Moore-Gibson-Thompson equation in R N , which is a third order in time equation that arises in viscous thermally relaxing fluids and also in viscoelastic materials (then under the name of standard linear viscoelastic model). First, we use some Lyapunov functionals in the Fourier space to show that, under certain assumptions on some parameters in the equation, an energy norm related with the solution decays with a rate (1 + t) −N/4 . But this does not give the decay rate of the solution itself. Hence, in the second part of the paper, we show an explicit representation of the solution in the frequency domain by analyzing the eigenvalues of the Fourier image of the solution and writing the solution accordingly. We use this eigenvalues expansion method to give the decay rate of the solution (and its derivatives), which (for the solution) results in (1 + t) 1−N/4 for N = 1, 2 and (1 + t) 1/2−N/4 when N ≥ 3.