Dualities between 2+1d fusion surface models from braided fusion categories

Fusion surface models generalize the concept of anyon chains to 2+1 dimensions, utilizing fusion 2-categories as their input. We investigate bond-algebraic dualities in these systems and show that distinct module tensor categories $\mathcal{M}$ over the same braided fusion category $\mathcal{B}$ give rise to dual lattice models. This extends the 1+1d result that dualities in anyon chains are classified by module categories over fusion categories. We analyze two concrete examples: (i) a $\text{Rep}(S_3)$ model with a constrained Hilbert space, dual to the spin-$\tfrac{1}{2}$ XXZ model on the honeycomb lattice, and (ii) a bilayer Kitaev honeycomb model, dual to a spin-$\tfrac{1}{2}$ model with XXZ and Ising interactions. Unlike regular $\mathcal{M}=\mathcal{B}$ fusion surface models, which conserve only 1-form symmetries, models constructed from $\mathcal{M} \neq \mathcal{B}$ can exhibit both 1-form and 0-form symmetries, including non-invertible ones.