The equilibrium geometries, harmonic vibrational frenquencies, and the dissociation energies of the OCH+-Rg (Rg = He, Ne, Ar, Kr, and Xe) complexes were calculated at the DFT, MP2, MP4, CCSD, and CCSD(T) levels of theory. In the lighter OCH+-Rg (Rg = He, Ne, Ar) rare gas complexes, the DFT and MP4 methods tend to produce longer Rg-H+ distance than the CCSD(T) level value, and the CCSD-calculated Rg-H+ bond lengths are slightly shorter. DFT method is not reliable to study weak interaction in the OCH+-He and OCH+-Ne complexes. A qualitative result can be obtained for OCH+-Ar complex by using the DFT method; however, a higher-level method using a larger basis set is required for the quantitative predictions. For heavier atom (Kr, Xe)-containing complexes, only the CCSD method predicted longer Rg-H+ distance than that obtained at the CCSD(T) level. The DFT method can be applied to obtain the semiquantitative results. The relativistic effects are expected to have minor effect on the geometrical parameters, the H+-C stretching mode, and the dissociation energy. However, the dissociation energies are sensitive to the quality of the basis set. The nature of interaction between the OCH+ ion and Rg atoms was also analyzed in terms of the interaction energy components.