Non-local 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 work, 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 wavefunction 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 $ν< d$, we analytically motivate and numerically demonstrate that, in the fully self-consistent model, algebraic edge mode decay with system size persists for all interaction 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.