Topological phases of an interacting Majorana Benalcazar-Bernevig-Hughes model

We study the effects of Coulomb repulsive interactions on a Majorana Benalcazar-Bernevig-Hughes (MBBH) model. The MBBH model belongs to the class of second–order topological superconductors ($HOTS{C}_2$), featuring robust Majorana corner modes. We consider an interacting strip of four chains of length $L$ and perform a DMRG numerical analysis based on a tensor–network approach. Study of the nonlocal fermionic correlations and of the degenerate entanglement spectrum indicates that the topological phases are robust in the presence of interactions, even in the strongly interacting regime.