Nonlocal Edge Mode Hybridization in the Long-Range Interacting Kitaev Chain

In one-dimensional p-wave superconductors with short-range interactions, topologically protected Majorana modes emerge, whose mass decays exponentially with system size, as first shown by Kitaev. In this Letter, we extend this prototypical model by including power law long-range interactions within a self-consistent framework, leading to the self-consistent long-range Kitaev chain (seco-LRKC). In this model, the gap matrix acquires a rich structure where short-range superconducting correlations coexist with long-range correlations that are exponentially localized at both chain edges simultaneously. As a direct consequence, the topological edge modes hybridize even if their wave function overlap vanishes, and the edge mode mass inherits the asymptotic scaling of the interaction. In contrast to models with imposed power law pairing, where massive Dirac modes emerge for exponents ν0, despite exponential wave function localization. While the edge mode remains massless in the thermodynamic limit, finite-size corrections can be experimentally relevant in mesoscopic systems with effective long-range interactions that decay sufficiently slowly. Our findings have direct implications for quantum simulations of ultracold microwave-shielded dipolar molecules in one-dimensional optical lattices.