Low-energy physics and finite-size effects in nonuniform parafermion chains

The nonuniform [Formula: see text] symmetric Kitaev chain, comprising alternating topological and normal regions, hosts localized states known as edge-zero modes (EZMs) at its interfaces. These EZMs can pair to form qubits that are resilient to quantum decoherence, a feature expected to extend to higher symmetric chains, i.e., parafermion chains. However, finite-size effects may impact this ideal picture. Diagnosing these effects requires first a thorough understanding of the low-energy physics where EZMs may emerge. Previous studies have largely focused on uniform chains, with nonuniform cases inferred from these results. While recent work [Narozhny, Sci. Rep. 7, 1447 (2017)] provides an insightful analytical solution for a nonuniform [Formula: see text] chain with two topological regions separated by a normal one, its complexity limits its applicability to chains with more regions or higher symmetries. Here, we present a new approach based on decimating the highest-energy terms, facilitating the scalable analysis of [Formula: see text] chains with any number of regions. We provide analytical results for both [Formula: see text] and [Formula: see text] chains, supported by numerical findings, and identify the critical lengths necessary to preserve well-separated EZMs.