Classification and construction of interacting fractonic higher-order topological phases

The notion of higher-order topological phases can have interesting generalizations to systems with subsystem symmetries that exhibit fractonic dynamics for charged excitations. In this work, we systematically study the higher-order topological phases protected by a combination of subsystem symmetries and ordinary global symmetries in two and three-dimensional interacting boson systems, with some interacting fermionic examples. In particular, we discover some important results, including (i) an Abelian subsystem symmetry group that does not allow for a nontrivial higher-order symmetry-protected topological phase in 2D without global symmetry. (ii) inhomogeneous subsystem symmetries do not allow for that either.