The 1D Kitaev model in the topological phase, with open boundary conditions, hosts strong Majorana zero modes. These are fermion parity-odd operators that almost commute with the Hamiltonian and manifest in long coherence times for edge degrees of freedom. We obtain higher-dimensional counterparts of such Majorana operators by explicitly computing their closed form expressions in models describing 2D and 3D higher-order superconductors. Due to the existence of such strong Majorana zero modes, the coherence time of the infinite temperature autocorrelation function of the corner Majorana operators in these models diverges with the linear system size. In the presence of a certain class of orbital-selective dissipative dynamics, the coherence times of half of the corner Majorana operators is enhanced, while the time correlations corresponding to the remaining corner Majoranas decay much faster as compared with the unitary case. We numerically demonstrate robustness of the coherence times to the presence of disorder.
Read full abstract