where H mol is the molecular Hamiltonian, H sol considers the solvent DoFs, H cav corresponds to the radiation modes inside the cavity, and H cbath describes a photon bath outside the cavity. For convenience, we set ̵ h = 1 throughout this work. For the molecular system, we consider an aggregate with N mol molecules, where each molecule is represented by a single reactive vibrational DoF, and the Hamiltonian is given by