We show the existence of thermal rectification in the graded mass quantum chain of harmonic oscillators with self-consistent reservoirs. Our analytical study allows us to identify the ingredients leading to the effect. The presence of rectification in this effective, simple model (representing graded mass materials, systems that may be constructed in practice) indicates that rectification in graded mass quantum systems may be an ubiquitous phenomenon. Moreover, as the classical version of this model does not present rectification, our results show that, here, rectification is a direct result of the quantum statistics. A fundamental challenge in statistical physics is the derivation of macroscopic phenomenological laws of thermodynamic transport from the underlying microscopic Hamiltonian systems. However, after decades, a first-principle derivation of the Fourier's law of heat conduction, for instance, is still a puzzle [1]. Many works have been devoted to the theme [2], most of them by means of computer simulations. But, in these problems where the central question involves the convergence or divergence of thermal conductivity, sometimes there is a significant difficulty to arrive at precise conclusions from numerical results [3]. Thus, the necessity of analytical studies together with the huge complexity of the associated nonlinear dynamical systems led to several works considering approximative schemes or simplified models [4] -since the pioneering work of Debye, the microscopic models to describe heat conduction are mainly given by systems of anharmonic oscillators, leading to problems without analytical solutions. Many works involve, e.g., the use of approximations such as Boltzmann equations, master equations for effective models, the analysis of Green-Kubo formula, etc. An example of simple (effective) model that can be analytically studied is the classical or quantum harmonic chain of oscillators with self-consistent reservoirs. It has been proposed a while ago [5,6], but it is always revisited [7,8]. In such model, each site of the chain is coupled to a reservoir; the first and last sites are coupled to "real" thermal baths, while the inner reservoirs only mimic the absent anharmonic interactions. The self-consistent condition means that there is no heat flow between an inner reservoir and its site in the steady state: the inner baths act only as a phonon scattering mechanism, such as 1