Signatures of Majorana bound states and parity effects in two-dimensional chiral p-wave Josephson junctions

We characterize topological features of Josephson junctions formed by coupled mesoscopic chiral p-wave superconducting islands. Through analytic and numerical studies of the low-lying BdG (Bogoliubov-deGennes) spectrum, we identify localized MBS (Majorana bound states) nucleated in Josephson vortices by the application of a perpendicular magnetic field. Additionally, we demonstrate the existence of an extended MBS that is delocalized around the outer perimeter of the coupled islands, which has measurable consequences on the Josephson supercurrent and phase dynamics of the junction. In particular, we predict a change in the critical current diffraction pattern in which the odd integer-flux nodes are lifted in a fermion parity-dependent fashion. We model the competing stochastic effects of thermal noise and macroscopic quantum tunneling within the RCSJ framework and show the emergence of a bimodal critical current distribution. We demonstrate that increasing the parity transition rate suppresses the bimodal nature of the distribution, thus strongly emphasizing the nontrivial parity dependent nature of the many-body ground state. Finally, we consider a trijunction geometry with three islands and discuss possible schemes to braid Majorana bound states by moving the Josephson vortices to which they are bound.