We should also emphasize that when multiple reaction pathways exist, neglecting the inference between these paths cannot correctly predict the reaction rate. In a recent work, a variety of mixed quantum-classical approaches are employed to investigate the resonance behavior of the cavity frequencydependent rate profile. Specifically, the reaction coordinate is treated quantum mechanically, while other DoFs are treated classically. It is found that although these mixed quantumclassical methods can capture the shape resonance of the cavity-modified reaction rate, both the reaction rate k 0 outside the cavity and the ratio k/k 0 inside the cavity are significantly overestimated. For an anharmonic system, the transition energy varies between different vibrational levels. In principle, we can observe multiple peaks in the cavity frequency-dependent rate profile, as long as multiple vibrational transitions meet the aforementioned energetic and kinetic criteria. This is the case for Model II with M = 2 a.u. . N cav v = 2 have already captured the key feature in the rate profile, which exhibits three sharp peaks, centered at around 745cm -1 , 900cm -1 , and 1645cm -1 , respectively.