We explore the stellar initial mass function (IMF) of a sample of 49 massive quiescent galaxies (MQGs) at 0.9$<$z$<$1.5. We base our analysis on intermediate resolution spectro-photometric data in the GOODS-N field taken in the near-infrared and optical with the HST/WFC3 G141 grism and the Survey for High-z Absorption Red and Dead Sources (SHARDS). To constrain the slope of the IMF, we have measured the TiO$_2$ spectral feature, whose strength depends strongly on the content of low-mass stars, as well as on stellar age. Using ultraviolet to near-infrared individual and stacked spectral energy distributions, we have independently estimated the stellar ages of our galaxies. Knowing the age of the stellar population, we interpret the strong differences in the TiO$_2$ feature as an IMF variation. In particular, for the heaviest z$\sim$1 MQGs (M$>$10$^{11}$Msun) we find an average age of 1.7$\pm$0.3 Gyr and a bottom-heavy IMF ($\Gamma_b$=3.2$\pm$0.2). Lighter MQGs (2$\times$10$^{10}$$<$M$<$10$^{11}$ Msun) at the same redshift are younger on average (1.0$\pm$0.2 Gyr) and present a shallower IMF slope ($\Gamma_b=2.7^{+0.3}_{-0.4}$). Our results are in good agreement with the findings about the IMF slope in early-type galaxies of similar mass in the present-day Universe. This suggests that the IMF, a key characteristic of the stellar populations in galaxies, is bottom-heavier for more massive galaxies and has remained unchanged in the last $\sim$8 Gyr.

(Abridged) We present the Survey for High-z Absorption Red and Dead Sources (SHARDS), an ESO/GTC Large Program carried out with GTC/OSIRIS. SHARDS is an ultra-deep optical spectro-photometric survey of the GOODS-N field (130 arcmin^2) at wavelengths 500 to 950 nm and using 24 contiguous medium-band filters (spectral resolution R 50). The data reach 26.5 mag (>3-sigma level) with sub-arcsec seeing in all bands. SHARDS main goal is obtaining accurate physical properties of interm- and high-z galaxies using well-sampled optical SEDs with sufficient spectral resolution to measure absorption and emission features. Among the different populations of high-z galaxies, SHARDS principal targets are massive quiescent galaxies at z>1. In this paper, we outline the observational strategy and include a detailed discussion of the special reduction and calibration procedures applied to the GTC/OSIRIS data. We present science demonstration results about the detection and study of emission-line galaxies (star-forming and AGN) at z=0-5. We also analyze the SEDs for a sample of 27 quiescent massive galaxies at 1.0<z<1.4. We discuss on the improvements introduced by the SHARDS dataset in the analysis of their SFH and stellar properties. We discuss the systematics arising from the use of different stellar population libraries. We find that the UV-to-MIR SEDs of the massive quiescent galaxies at z=1.0-1.5 are well described by an exponential decaying SFH with scale tau=100-200 Myr, age 1.5-2.0 Gyr, solar or slightly sub-solar metallicity, and moderate extinction, A(V)~0.5 mag. We also find that galaxies with masses above M* are typically older than lighter galaxies, as expected in a downsizing scenario of galaxy formation. This trend is, however, model dependent, i.e., it is significantly more evident in the results obtained with some stellar population synthesis libraries and almost absent in others.