Covariant Interacting Fractons

We show that there exists a natural analogue of the Yang-Mills equations using the Fr\"olicher-Nijenhuis bracket between vector-valued differential forms. The gauge field is a rank-two tensor, and when one constrains it to be symmetric, then the system exhibits fractonic behaviours. In the linearised limit, the constrained equations of motion reduce to those of the covariant fracton model [Bertolini-Maggiore 2022].