Infrared spectra of [Re(X)(CO)(3)(alpha-diimine)] (alpha-diimine = 2,2'-bipyridine, X = Cl, NCS, or pyridylimidazo [1,5-a]pyridine, X = Cl) in the ground and the lowest triplet electronic states were calculated by a global hybrid density functional going beyond the harmonic level by means of second-order vibrational perturbation theory (VPT2) and including bulk solvent effects by the polarizable continuum model (PCM). The full-dimensionality (FD) VPT2 is compared with the reduced-dimensionality (RD) model, where only selected vibrational modes are calculated anharmonically. The simulated difference IR spectra (excited state minus ground state) in the nu(CO) region closely match experimental time-resolved infrared (TRIR) spectra. Very good agreement was also obtained for ground-state spectra in the fingerprint region. In comparison with the harmonic simulated spectra, the calculated anharmonic frequencies are closer to experimental values and do not require scaling when the B3LYP functional is used. Several spectral features due to combination bands have been identified by VPT2 simulations in the nu(CO) spectral region, which are of importance for a correct interpretation of TRIR experiments.