Benchmark calculations with a goal to find dispersion-corrected DFT-D3 methods suitable for a reliable estimation of matrix shifts on the vibrational frequency were carried out on the ThO molecule in three rare gas (Rg = Ne, Ar, and Kr) matrices. The matrices were modelled by the explicit approach, in which a single and a double shell of Rg atoms around ThO was considered. The selection of exchange-correlation functionals was based on test calculations on triatomic ThO⋯Rg models. The B3LYP, PBE0, CAM-B3LYP, and LC-ωPBE functionals were found to be the best suited for the estimation of matrix shifts. The single shell of Rg’s around ThO accounted for a major part of the shifts; the addition of a second Rg shell resulted only in a minor improvement. Continuum solvation models considerably overestimated the effect of Rg matrices both when the whole matrix was treated by the model and when the first shell was treated explicitly and the rest with a continuum solvation model.