The IR vibrational properties and the corresponding reflectance spectra of the six most common members of the garnet family (pyrope Mg3Al2Si3O12, almandine Fe3Al2Si3O12, spessartine Mn3Al2Si3O12, grossular Ca3Al2Si3O12, uvarovite Ca3Cr2Si3O12, and andradite Ca3Fe2Si3O12) were simulated at the ab initio level with the CRYSTAL09 code by using a large all-electron Gaussian-type basis set and the B3LYP hybrid functional. The 17 IR active F1u transverse optical (TO) and longitudinal optical (LO) frequencies, the oscillator strengths, the high frequency and static dielectric constants, and the reflectance spectrum were computed. The agreement with experiments for the TO and LO peaks is always excellent, the mean absolute difference for the whole set of data (overall 178 peaks) being 5 cm−1. Oscillator strengths, calculated from the mass-weighted effective Born charges, are found in semi-quantitative agreement with the experimental data. The reflectance spectra, simulated through the classical dispersion relation, reproduce the experimental curves extremely well. The availability of the full set of simulated frequencies and intensities, obtained by using uniform computational tools (computer code, variational basis sets, density functional), permits the establishment of correlations between IR wavenumbers and structural features suggested, but only partially documented, in the past.