Probing superconductivity with tunneling spectroscopy in rhombohedral graphene

Motivated by experiments on rhombohedral tetralayer graphene showing signs of superconductivity emerging from a valley-polarized normal state, we here analyze theoretically how scanning tunneling spectroscopy can be used to probe the superconducting order parameter of the system. To describe different pairing scenarios on equal footing, we develop a microscopic tunneling approach that can capture arbitrary, including finite-momentum, superconducting order parameters and low-symmetry normal-state Hamiltonians. Our analysis shows that the broken time-reversal symmetry in a single valley leads to unique features in the weak-tunneling regime that are different for commensurate and incommensurate Cooper pair momenta. We further uncover an unconventional spatial dependence of the Andreev conductance, allowing to distinguish between three topologically distinct classes of single-$\mathbf{q}$ pairing states in the system, and compute the signatures of a competing translational-symmetry breaking three-$\mathbf{q}$ ''moir\'e superconductor''.