Boundary anomaly detection in two-dimensional subsystem symmetry-protected topological phases

We generalize the topological response theory to detect the boundary anomalies of linear subsystem symmetries. This approach allows us to distinguish different subsystem symmetry-protected topological (SSPT) phases and uncover new ones. We focus on the cases where the mixed anomaly exists within the adjacent subsystems. Using numerical simulations, we demonstrate the power of this method by identifying strong and weak ${Z}_{2}^{\ensuremath{\tau}}\ifmmode\times\else\texttimes\fi{}{Z}_{2}^{\ensuremath{\sigma}}$ SSPT phases in a tunable tensor network state. Our analysis reveals an intrinsic ${Z}_{2}$ SSPT phase characterized by its degenerate entanglement spectrum. Furthermore, we extend the anomaly indicator to mixed-state density matrices and show that quantum anomalies of subsystem symmetry can persist under both uniform and alternating disorders. This finding establishes a connection between boundary quantum anomalies in pure and mixed states. Our work provides a numerical method to detect quantum anomalies of subsystem symmetries, offering new insights into the study of topological quantum phases.