Mirror-protected Majorana zero modes in f-wave multilayer graphene superconductors

Inspired by recent experimental discoveries of superconductivity in multilayer graphene, we study models of f-wave superconductivity on the honeycomb lattice with arbitrary numbers of layers. For odd numbers of layers, these systems are topologically nontrivial, characterized by a mirror-projected winding number ν_± = ±1. Along each mirror-preserving edge in armchair nanoribbons, there are two protected Majorana zero modes. These modes are present even if the sample is finite in both directions, such as in rectangular and hexagonal flakes. Crucially, zero modes can also be confined to vortex cores. Finally, we apply these models to twisted bilayer and trilayer systems, which also feature boundary-projected and vortex-confined zero modes. Since vortices are experimentally accessible by local scanning probes, our study suggests that superconducting multilayer graphene systems are promising platforms to create and manipulate Majorana zero modes.