Topological Phase Diagram of Generalized SSH Models with Interactions

We investigate interacting Su-Schrieffer-Heeger (SSH) chains with two- and three-site unit cells using density matrix renormalization group (DMRG) simulations. By selecting appropriate filling fractions and sweeping across interaction strength $J_z$ and dimerization $\delta$, we map out their phase diagrams and identify transition lines via entanglement entropy and magnetization measurements. In the two-site model, we observe the emergence of an interaction-induced antiferromagnetic intermediate phase between the topologically trivial and non-trivial regimes, as well as a critical region at negative $J_z$ with suppressed magnetization and finite-size scaling of entanglement entropy. In contrast, the three-site model lacks an intermediate phase and exhibits asymmetric edge localization and antiferromagnetic ordering in both positive and negative $J_z$ regimes. We further examine the response of edge states to Ising perturbations. In the two-site model, zero-energy edge modes are topologically protected and remain robust up to a finite interaction strength. However, in the three-site model, where the edge states reside at finite energy, this protection breaks down. Despite this, the edge-localized nature of these states survives in the form of polarized modes whose spatial profiles reflect the non-interacting limit.