We also explored the reaction at a larger light-matter interaction strength, η c = 0.005 a.u. . A detailed analysis of the reaction dynamics is presented in the SI. A larger η c entails considering more highly excited photonic states (e.g. N cav v > 6 for η c = 0.005 a.u. ) to obtain accurate results. The reaction rates exhibit a significant increase across the entire frequency range, as depicted in Fig. , which displays the ratio k/k 0 as a function of ω c . Additionally, we observe that the relative heights of two lower-frequency peaks are changed. The most prominent peak now appears at 825cm -1 with a shoulder around 750cm -1 . This shift is attributed to the stronger coupling to the cavity photon mode and successively to the cavity bath, resulting in a further broadening of molecular vibrational levels. As a consequence, a cavity photon with a frequency 825cm -1 can induce both the vibrational transitions v mol 0/1 ↔ v mol 3/2 as well as v mol 2/3 ↔ v mol 5/4 , and thereby maximizes the rate enhancement. In other words, a stronger light-matter interaction can amalgamate nearby peaks corresponding to different vibrational transitions in an anharmonic system and intensify the rate modification. These explanations can be applied to elucidate a rate profile with multiple broad peaks observed in the previous theoretical works. As a side note, it is important to note that not all cavityinduced cotunneling pathways necessarily manifest as peaks in the cavity frequency-dependent rate profile. If the tunneling component of a cavity-induced cotunneling pathway is considerably slower than the reaction rate outside the cavity, a corresponding feature may not be observable in the rate profile. An illustrative example is provided in Appendix B. Multiple peaks in the rate profile have not been observed in the experiments. However, we predict that they would be measurable, e.g. in a strongly anharmonic system where the selection rules for the harmonic oscillator break down.