The analytic theory kVSC (Eq. 17) predicts that the VSC rate enhancement (Eq. 19) scales as kVSC /k 0 ∝ Ω 2 R as the light-matter coupling increases, then plateau when Ω R becomes large. This is in excellent agreement with the numerically exact HEOM simulations and provides a non-linear relation between the change of the effective free energy barrier and the light-matter coupling strength (Fig. ), which has been observed in the VSC experiments . The theory kVSC (Eq. 17) also predicts changes in both effective activation enthalpy and entropy (as observed in the experiments ), which agrees well with the numerical exact HEOM results within all the parameter regimes we explored. Although a re-scale parameter α = 0.5 is needed to bring the numerical values of the FGR rate constants in consistency with the HEOM exact results, the overall scaling relations with respect to Ω R , ω c , τ c , and T reported in the paper are rather general and impressive, which should not be restricted to the detailed shape of the potential energy surface or environmental spectral density functions. We further generalized the kVSC expression to consider the many mode effects, and the resulting theories (Eq. 33 for 1D FP cavity and Eq. 34 for 2D cavity) predict the normal-incidence resonance condition: the peak of the rate constant enhancement occurs when k ∥ = 0 and ω c = ω 0 . Last but not least, our current theory also predicts that for two chemically similar reactions, if one satisfies k 1 ≪ k 2 , k 3 but the other does not (due to the low reaction barrier), then there will be a VSC effect for the first reaction but not for the second one. This is because for the second reaction, |ν L ⟩ → |ν ′ L ⟩ is no longer rate limiting, and the cavity modification of this process will no longer influence the apparent rate constant. This might be the explanation for the recently observed null effects in VSC experiments .