# Quasiparticle Interference Patterns in a Topological Superconductor

###### Abstract

In light of recent proposals to realize a topological superconductor on the surface of strong topological insulators, we study impurity and vortex scattering in two dimensional topological superconductivity. We develop a theory of quasiparticle interference in a model of the surface of a three dimensional strong topological insulator with a pairing term added. We consider a variety of different scatterers, including magnetic and nonmagnetic impurity as well as a local pairing order parameter suppression associated with the presence of a vortex core. Similar to the case of a surface of a three dimensional topological insulator without pairing, our results for non-magnetic impurity can be explained by the absence of back scattering, as expected for a Dirac cone structure. In the superconducting case, doping away from the Dirac point leads to a doubling of the contours of constant energy. This is in contrast to the unpaired case where the chemical potential simply adds to the bias voltage and shifts the energy. This doubling of contours results in multiplying the number of possible scattering processes in each energy. Interestingly, we find that some processes are dominant in the impurity case while others are dominant in the vortex case. Moreover, the two types of processes lead to a different dependence on the chemical potential.

## I Introduction

The search for a Majorana mode in a condensed matter system is currently a very active area of research. Some very promising experimental results have been reported in one dimensionMourik et al. (2012); Nadj-Perge et al. (2014) while in two dimensions appealing theoretical proposals existFu and Kane (2008); Tanaka et al. (2009). The above systems are based on the combination of momentum spin locking with superconducting pairing. In two dimensions this combination is provided by heterostructures of a superconducting layer in contact with either a three dimensional topological insulator (3DSTI)Fu and Kane (2008) or a spin-orbit coupled semiconductorSau et al. (2010); Alicea (2010). While the superconductor furnishes the pairing, the spin orbit coupled layer provides the momentum-spin lockingKoren et al. (2011); Sacépé et al. (2011); Qu et al. (2012); Williams et al. (2012); Cho et al. (2013); Xu and et al. (2013). Another type of proposal makes use of an innate tendency for developing superconductivity in materials with spin orbit coupling (SOC)Zhao et al. (2014); Farrell and Pereg-Barnea (2013); Farrell and Pereg-Barnea (2014a, b).

Given the above, it is timely to search for unique properties of a system where superconductivity arises from an underlying Dirac-like band structure. In this Paper we develop a theory of quasiparticle interference (QPI) patterns in such systems. The QPI patterns are of direct experimental relevance as they can be measured using the technique of Fourier-transform scanning tunnelling spectroscopy (FT-STS)Hoffman et al. (2002); He et al. (2014); Vershinin et al. (2003). This method measures the local density of states (LDOS) of a sample in the vicinity of a single impurity. Theoretically, the pattern which is observed is dictated by the underlying clean system (i.e without the impurity) and so properties of the clean system can be deduced from such a pattern and readily compared with calculationsGuo and Franz (2010); Plamadeala et al. (2010); Pereg-Barnea and Franz (2008); Pereg-Barnea and MacDonald (2008).

We consider the surface of a strong topological insulator with pairing added. We do this by treating two complementary models. The first is a single Dirac cone in the continuum and should capture universal properties of a Dirac band structure in the presence of pairing. Second, we scrutinize our continuum model results by looking at a more physical lattice model. This model has been developedMarchand and Franz (2012) for the surface of a strong topological insulator. It specifically avoids the doubling theorem by including both surfaces of the material. We calculate the QPI pattern for a variety of scatterers, including a charge defect, a magnetic defect and a defect in the superconducting order parameter.

The QPI for a strong topological insulator in the presence of charge and magnetic impurities has been calculated by one of us in the pastGuo and Franz (2010). With the chemical potential tuned to the Dirac point, we obtain results consistent with this previous study with the exception that QPI pattern is only observable for energies above the superconducting gap , as illustrated in Fig. 1. As noted previously, unlike in normal metals where QPI shows distinctive peaks, patterns in STIs with non-magnetic impurities are non-singular and exhibit only an edge. We find this behavior to persist when the SC order is included. An interesting difference however occurs when we consider disorder in the SC order parameter amplitude which will generically arise in the presence or pair-breaking impurities of vortices. In this case, QPIs display a sharp peak, as illustrated in Fig. 1c.

Moving the chemical potential away from the Dirac point does not result in any angular dependence of the properties of this pattern. However, we do see some interesting behavior in these patterns. We find a contour of singular LDOS features. As the chemical potential is increased (decreased) with respect to the Dirac point, the radius of the contour in the impurity induced QPI patterns is also increased (decreased) and changes linearly with the chemical potential. On the other hand, in the order parameter suppression QPI pattern, the radius of the singular contour is independent of the chemical potential. We trace this difference back to the fact that different quasiparticle scattering processes are favored in these two situations. It can also be argued that this particular dependence on for these two types of scatterers comes from the underlying Dirac band structure. These observations in the continuum model can be understood heuristically via chirality angle arguments and are supported by our lattice model calculations.

Although tuning the chemical potential in these systems may be challenging from an experimental point of view, the underlying result (e.g. the behavior of the QPI singularity) can be observed via other means. As an example, we suggest a simplified scenario where the sample we consider is placed in a capacitor. This adds a bias voltage between the edges of the material. The voltage acts as a chemical potential with differing signs on each edge of the sample. By changing the bias one can observe the distinct behavior the QPI due to an impurity or a vortex. There has recently been an experimental realization of a set-up analogous to this capacitor systemChang et al. (2015). In this work it was shown that the chemical potential on each edge of thin film 3DSTI can be tuned independently of the other edge through the use of dual-gate structures. Our results are then directly applicable to an FT-STS measurements on such a device when the chemical potential on one edge is tuned to be opposite to each other edge.

The rest of this paper is organized as follows. In the next section we give an overview of the methods we have followed to obtain our results. We present our method for calculating the QPI pattern, define each type of perturbation, present our two model systems and discuss numerical details. Section III then presents our results for the QPI patterns. We begin by tuning the chemical potential to the Dirac point, where the QPI pattern for our continuum model can be calculated exactly (the details of this are relegated to the Appendix). This not only helps us compare our results with a previous study, but also allows us to develop some intuition. We then move on to the main result of this paper, the behavior of QPI patterns at finite chemical potential. Here we present numerical data as well as further discussion. We then discuss an alternative to tuning the chemical potential and perform an explicit calculation using our lattice model.

## Ii Methodology

### ii.1 Green’s function for a single impurity and the Born approximation

We begin by discussing a general expression for the total Green’s function. The approach taken here will be similar to previous workGuo and Franz (2010), but appropriately generalized to include superconductivity and to be suitable for application to our lattice model. In either case we take a total Hamiltonian , the first term describes the underlying clean system while the latter term describes the perturbation. The clean system’s Hamiltonian is given by:

(1) |

here we have defined where labels additional degrees of freedom in the model (spin, sub lattice, etc.). The matrix is given by

(2) |

where and are matrices in the space of the degrees of freedom implicit in before, describes the normal band structure of the model and describes the pairing. Next we consider a local perturbation which couples to the electronic degrees of freedom in some general form:

(3) |

where is the number of lattice sites; is the Fourier transform of the perturbation potential. Its matrix structure depends on the impurity we consider and will be discussed in detail below.

Next we consider calculating the Green’s function for this model. We define the Green’s function as follows

(4) |

where is the imaginary time and is the Heisenberg picture version of . Fourier transforming to Matsubara frequency and defining the -matrix we write

where is the Green’s function of the clean system. In the case of a local perturbation, where the potential is a -function in real space, the potential and the -matrix are momentum independent and the self consistency equation for the -matrix can be solved:

(6) |

In principle poles in the -matrix reveal information about bound states in the system. However, in this study we are interested in the QPI patterns in momentum space. Since the -matrix has no momentum dependence, the results of a weak perturbations will not differ from the simpler Born approximation. To this end, we keep only the leading order term in which leaves

Moreover, in the Born approximation, which is appropriate away from impurity bound state energies, the momentum dependence of the potential separates from that of the Greens functions convolution when calculating the Fourier transformed density of statesCapriotti et al. (2003).

### ii.2 Local density of states and impurity types

Here we will describe how we calculate the various LDOS patterns in the work to come. We are considering FT-STS experiments where the tip of the STM is normal (not magnetic or superconducting) and measures the LDOS, theoretically this quantity is given as follows

(8) |

where acts on particle-hole degrees of freedom, is an infinitesimal and is the double Fourier transform of defined in Eq. (II.1) above. Fourier transforming one finds that the modulation to the local density of states in the Born approximation is

(9) |

where .

Following Ref. [Guo and Franz, 2010], we are interested in a more general interference pattern. Here we imagine that the STM tip can resolve different degrees of freedom (e.g. spin can be resolved by a spin polarized tip). It is therefore possible to insert a general matrix into the trace. This matrix acts to resolve the component of the LDOS pattern of interestGuo and Franz (2010).

We now discuss the various types of perturbations that will be considered. Here we will focus on a physical description of these impurities and leave the formal details to the Appendix. We begin with a simple charge impurity. Here we imagine that the chemical potential at a certain site has been altered. Second, we are interested in magnetic impurities. These alter the local Zeeman splitting on a single site and so couple to the spin of the electrons. In what follows we will combine charge and magnetic impurities into a single heading which we will refer to as impurities. It is then useful to define the QPI patterns for any of these impurities as

(10) |

The matrices are outlined in the appendix. The first label, , denotes the type of STM tip. is a normal (charge) tip while or resolves the component of the electron’s spin along the or direction respectively. Meanwhile, labels the type of impurity we are considering. is a charge impurity while or refers to a magnetic impurity with its spin along the or axis. When we refer to patterns such as it is the patterns above to which we are referring. It should be noted that any physical QPI pattern due to charge/spin scattering can be written as a combination of the above.

The second class of perturbation we are interested in is a local suppression in the superconducting order parameter. We will refer to this as OP suppression. This perturbation can be thought of as a simplified description of a vortex where only the OP suppression at the vortex core is taken into accountPereg-Barnea and Franz (2008). Alternatively, other types of disorder are sometimes accompanied by OP variationsLang et al. (2002). We refer to the LDOS in this case as .

### ii.3 Model Hamiltonians

To discuss the surface of a 3DSTI we will make use of two model Hamiltonians. The first is a continuum Dirac cone model which will be our primary focus. This model is a single Dirac coneGuo and Franz (2010) with a chemical potential and -wave pairing. Symbolically we have

(11) |

with

(12) | |||

where , is the chemical potential and is the OP amplitude. We have adopted units such that . Such a model can be diagonalized and yields the two eigenvalues at each wave vector . Further, finding an exact, closed-form expression for the Green’s function of the above system is possible. We have outlined the details this calculation in the Appendix.

The second model is a lattice model of a 3DSTI where only the surfaces of this system are considered (we will arbitrarily refer to the two surfaces as ’top’ and ’bottom’). This model will be used sparingly and will mostly be employed as a physical consistency check for features we find in the continuum model.

In order to model the single Dirac cone on the surface of a 3DSTI we follow Marchand and FranzMarchand and Franz (2012). Their approach is briefly outlined here. The doubling theorem states that in a time reversal invariant periodic systems Dirac points appear in pairs. Therefore, any two dimensional lattice model can not have an odd number of Dirac points. The 3DSTI is indeed a periodic system. However, it avoids the doubling theorem by placing half of its Dirac points on each surface. To mimic this we employ a two-surface model. This is a minimal way to model a surface with an odd number of Dirac points.

In this approach, one begins with a three dimensional model and integrates out the bulk, leaving only the two edges of the material. Adopting this model gives our clean system’s Hamiltonian:

(13) |

where and and where . The matrix above is an matrix, the sub matrices and act on a space of spin. The basis of the above is where arrows refer to spin, the numbers refer to one of two bands and ’’ and ’’ refer to top and bottom edges. With this basis in mind, we see that couples orbitals within the same surface and couples different surfaces.

We now introduce superconductivity pairing only within the same orbital and the same surface.

(14) |

where . We assume -wave pairing and therefore take . The above has 8 doubly degenerate eigenvalues:

(15) |

whereMarchand and Franz (2012) with and . The clean Green’s function of this system is not analytically tractable and so we rely on numerics for its calculation.

### ii.4 Numerical method

Here we briefly outline the numerical methods we use. We numerically compute the LDOS patterns for particle-like impurities as described in Eq. (10). The idea is to compute the Green’s function in the absence of any impurity from the Hamiltonian using . Then, we produce the convolution with the impurity . Once this is obtained, multiplying it with then taking the trace produces the wanted results.

Each term in the perturbed Green’s function is therefore a convolution of two functions of momentum, with each of these functions being an entry in the clean system’s Green’s function matrix. In order to minimize run time and obtain high resolution LDOS maps we preform the convolution using the fast Fourier transform algorithm (FFT). This amounts to first using FFT to express the two functions in real space, performing a direct product and then using FFT again to get back to momentum space. This lowers the run time from to where is the number of points in the Brillouin zone.

## Iii Quasiparticle Interference Patterns in a Topological Superconductor

### iii.1 Zero Chemical Potential

At the particle-hole symmetric point, i.e., , the spectrum possesses Lorentz invariance. This enables us to find a closed form solution to the continuum model, for all of the perturbations we consider. This is done using the exact Green’s function calculated in the appendix along with a standard Feynman parameterizations trick (see e.g [Peskin and Schroeder, 1995]) which we will not repeat here. For the impurities these patterns are most easily presented by defining the matrix , the variables and and and the functions

(16) | |||||

(17) |

and we have set hereafter for simplicity. With these definitions we find

(18) |

where is an ultraviolet cut-off and . We note that with some small exceptions the above results are almost identical to those for a strong topological insulator without superconductivityGuo and Franz (2010) with the replacement .

The OP suppression can also be calculated in closed form at and yields the results

For simplicity we have presented to above results in Matsubara space, the proper analytic continuation must be employed to obtain the physical QPI patterns.

We have compared the exact results above with those we have obtained from the numerical calculation at using the model in Eq. (13) with superconductivity added. We find that the two patterns agree remarkably well showing similar angular and radial features (the specifics of these features will be discussed below). We have showcased five of these patterns in density plots in the top of Fig. 2. We now move on to discuss three properties of the QPI patterns we have found in the continuum model, keeping in mind the agreement between the lattice and continuum results.

First, let us discuss the radial features of the QPI patterns. Inspecting the above results we notice that the function is singular at a length of while has a kink at this value (its derivative should be singular). To understand this value let us recall that the bulk energy bands for the continuum model are given by . With the two bands become degenerate and contours of constant energy for a given frequency are then given by . Our results for the LDOS therefore show maximum response at twice this wave vector length. Physically, this corresponds to scattering across the diameter of the contours of constant energy as is expected for such a circularly symmetric problemPereg-Barnea and MacDonald (2008). This change in the singular value of with the frequency can be seen in the bottom plots of Fig. 2. Here we plot the LDOS for several perturbations and probes. The singularities (or kink in the case of ) clearly vary with .

Next, we to study the effect of superconductivity on the QPI patterns. In short, the occurrence of in the critical radius of the LDOS described above is a signature of superconductivity. Without superconductivity the peak would occur at and would thus persist all the way down to Guo and Franz (2010). With superconductivity present all features disappear once , i.e. when we probe energy scales within the gap. Thus having the pairing present in the system shifts the radius of this major feature to lower values. We have explored this result in Fig. 2 by plotting several different LDOS patterns for varying . For we see that the singularity/kink in these patterns moves to smaller values of as decreases. For , or the solid curve in this figure, we see no signal at all.

Finally, we discuss the angular features of the QPI patterns we have calculated. Studying Eq. (18), we see that one requires the input to be non-zero to find a pattern that is not circularly symmetric. In the case we have either zero (for ) or a circularly symmetric function function when . Recall, the result represents a non-magnetic impurity probed with a normal STM tip. The other results above would be those that are obtained from other perturbations and a normal STM tip. Thus our underlying chiral system does not show any angular dependence along the circular singularity in the QPI maps. This is the case for a normal STM tip at regardless of the perturbation type. That being said, as in Ref. [Guo and Franz, 2010] we can see angular dependence in the patterns, which corresponds to a spin-filtered STM tip. To look at these angular features we have presented density plots of our calculations on the lattice model in the top of Fig. 2.

### iii.2 General Chemical Potential

Moving to non-zero chemical potential makes calculating the LDOS exactly for the continuum model intractable. Therefore we must rely on numerics for both the continuum and lattice model. Our results for the lattice model are presented in Fig. 3. These results lead us to make three observations. First, in the figure we notice that in addition to one major circular pattern, we can see a very subtle secondary pattern in the case of magnetic impurities. Second, all patterns essentially respect the same angular symmetry as at . Third, the radius of the major pattern in the impurity QPI patterns changes noticeably when compared to the results whereas the OP suppression pattern does not.

Let us put these observations onto some more solid ground. We begin with the observation that the angular symmetry does not change. We have thus established that tuning the chemical potential away from the Dirac point does not lead to angular signatures in either the impurity or OP suppression perturbation. We now dedicate the rest of this subsection to explaining the other two observations above.

To understand the above dependence, let us recall an observation from the last subsection; the major radial features in the LDOS patterns appear at values of corresponding to impurity scattering between states at the same quasiparticle energy. Additionally, this scattering occurs across the diameter of the contour of constant energy. At we have degenerate contours of constant energy. When is nonzero both of our models develop multiple contours of constant energy. In the continuum model, the radii of these contours are given by . With these two contours one can imagine 4 quasiparticle scattering processes across the diameter of these circles. These processes are illustrated in Fig. 4 and labelled through .

With these different processes in mind we make the following two observations: (1) magnetic and non-magnetic impurities favor scattering processes between the same contour (intra-contour), (2) the OP suppression perturbation favors scattering processes between different contours (inter-contour). The QPI patterns of impurity scaterrers show a major singularity at for with a minor singularity at . For the major singularity occurs at with the secondary minor singularity at . Meanwhile, the pattern of an OP suppression shows singularities at when and when

In relation to the above we find interesting mathematical relations for the radii of the singularities in the LDOS patterns. The impurity QPI patterns have their major singularity at while the OP suppression pattern shows a singularity at . Thus the intra-contour process increases linearly with while the inter-contour scattering wave vector length is independent of .

We can understand the above results by heuristic arguments. The two constant energy contours seen in the case are the result of band reflection as depicted in Fig. 5. In this figure the Dirac cone is intermitted by a gap and reflected about the horizontal line due to the introduction of superconductivity. On each of these contours one can define a chirality angle . This chirality angle affects the wave function in two ways. First, even without superconductivity, spin-orbit coupling locks the spin direction to the momentum direction. Second, the pairing inherits this chirality angle, becoming effectively a -wave superconductor. Its order parameter winding appears in the coherence factors. The notion of chirality is depicted by arrows in Fig. 5.

To describe this analytically we note that in our choice of basis the Hamiltonian has the form

(20) |

and the positive energy eigenvectors have the form

(21) |

here , and .

The two branches of energy above (denoted with a or a ) have chirality directions which wind in opposite directions in the Brillouin zone (see Fig. 5). Notably, the chirality direction for a given branch is completely inverted upon sending . This has important implications for our system where (as we have already discussed) scattering across the diameter of contours of constant energy is favored.

Looking at the impurity potentials and considering its action on one can see that it reverses the chirality angle. This means that will have the same chirality angle as . This explains the favoring of the and transitions.

Meanwhile, the OP suppression potential replaces both the spin and particle hole degree of freedom. As a result, the wavefunction is almost in tact (except for the Bogoliubov coherence factors being exchanged). This means will have the same chirality angle as which leads to a suppression of scattering across the same contour as here the initial state and the final state have opposite chirality angles. Furthermore, these results lead to enhancement of inter-contour scattering. This is because there is a place on the other contour with the same chirality angle. For energies above the Dirac point we are cutting each branch only once and so the two contours we are left with have chirality angles winding in opposite directions. This leads to scattering processes. For energies below the Dirac point we cut the same energy branch twice and so we will have two contours with chirality angles winding in the same direction. For this case processes are favored.

We explore the above results in our continuum model in the left of Fig. 6. Here we have plotted and for several different values of the chemical potential. We see very clearly that the impurity scatterers (e.g. the pattern) have singularities that occur at values that increase linearly with . Meanwhile, the pattern for the OP suppression in this figure illustrates clearly that the singularity in this pattern does not depend on .

We now contend that the functional form with respect to of , the radius for the OP suppression, and , the radius for the regular impurity, (i.e. linear and constant respectively) is a result of the underlying Dirac band structure. To argue this we approximate the spectrum by where is the dispersion of the underlying band structure. If we wish to invert this equation to find as a function of , i.e. to find the equation for the contours of constant energy, we obtain .

Now let us think about a circularly symmetric dispersion for simplicity. Then . Let us assume that we are close enough to the center of the Brillouin zone so that the leading order term in an expansion of is valid. We then take where is some number and is a constant. Then the circles of constant energy have a radius given by

(22) |

We then see that, provided and are of comparable size, the only type of dispersion that yields energy contours that depend on in a linear matter is , or a Dirac-like dispersion. Thus this linear scaling of the major singularity in the QPI patterns is a property of a linearly dispersing band structure. Furthermore, the independence on of scattering from one contour to another is also a consequence of having a linear dispersion. The radius of singularities for such a process will be

(23) |

The above is independent for only.

The above argument of course only holds provided the values we are interested in are suitably close to the -point (or, more generally, wherever the Dirac point occurs) so that our linear approximation is valid. That is to say, in a realistic system will contain other subleading contributions on top of the linear Dirac like term. To explore how well our results hold in such a system we have calculated and for several different values of the chemical potential in our lattice model. These results are presented on the right of Fig. 6. In this figure we see that our observation holds very well provided that is kept small enough. For intermediate values of deviations increase as moves away from the Dirac point.

### iii.3 System in a Capacitor

In the previous subsection we have outlined an interesting dependence of the radii of the major feature in the LDOS on the chemical potential. The chemical potential can be tuned in STIs by chemical doping or, in a thin flake, by electrostatic gating. To illustrate the effect of a non-zero chemical potential in a simple setting we imagine placing a system inside a capacitor which has an effect of biasing the system so that there is a potential at the top and at the bottom. As mentioned in the introduction of this work such a scenario could be realized through the use of dual-gate structuresChang et al. (2015). We describe this by adding

(24) |

to our clean Hamiltonian in Eq. (13).

In our lattice model the two surfaces (top and bottom) are only weakly coupled and so we may think of this as introducing a chemical potential to the two surfaces. This chemical potential is equal in magnitude but opposite in sign on each edge. Let us, without loss of generality, assume the positive bias is introduced on top of the sample (notice that enters the Hamiltonian different by a minus sign from the chemical potential). Thinking of the two edges as isolated Dirac points, our findings thus far in the paper then simply predict the following: The radius of the major peak in the LDOS pattern for an impurity QPI should decrease on the top of the sample (as if ), and decrease on the bottom (as if ). Meanwhile, an OP suppression should look the same on the top and bottom of the sample.

In addition to the qualitative prediction made above, the strength of the bias in this setup is, at least in principle, tuneable via external means.

In Fig 7 we present the LDOS on the top and bottom of a sample for several values of this bias voltage. We see that the exact scenario described above plays out, namely as is increased the singularity in the impurity QPI pattern moves to larger on one edge and smaller on the other. In addition, we see that the OP suppression placed on the top and bottom of the sample give identical results and are essentially independent of this bias voltage.

## Iv Conclusions

We have studied the QPI patterns induced by different local perturbations in both a lattice and a continuum model for a SC surface of a three dimensional topological insulator. Our results for a half filled band (such that the chemical potential is at the Dirac point) are similar to those calculated in the past for a strong topological insulator surface. The existence of superconductivity in the system gaps out the low-lying excitations and shows up qualitatively in the radius of the singularity/kink that occurs in the QPI pattern. The presence of reduces the critical radius of this singularity. For STM bias values below the superconductive gap no signal can be observed. For non-magnetic impurities the QPI response is weak and consists of an edge similar to the normal state result. Remarkably though, we find that when disorder in the order parameter amplitude is included (one that will generically be present) the edge transforms into a peak, which should be more easily observed in experiment. The shape of the singularity – edge vs. peak – can thus be used as an indicator of the dominant source of quasiparticle scattering in the sample.

With a finite chemical potential we find almost no change in the angular features of the QPI pattern. We do however find that the quasiparticle scattering processes contributing to impurity and OP suppression scattering are different. As a result, the singular features of the impurity QPI pattern depend linearly on the chemical potential and that those of the OP suppression are independent of the chemical potential. We argued that this functional dependence is unique to a Dirac-like spectrum and showed that it approximately holds in our more complicated lattice model. We have also proposed and verified with a calculation an alternative method to tune the chemical potential by placing the sample in a capacitor.

## V Acknowledgements

The authors are thankful for useful discussions with J. Hoffman Financial support for this work was provided by the NSERC and FQRNT (TPB), the Vanier Canada Graduate Scholarship (AF) and CIfAR. Some of the numerical calculations for this work were performed using CLUMEQ/McGill HPC supercomputing resources.

## References

- Mourik et al. (2012) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012).
- Nadj-Perge et al. (2014) S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yazdani, Science (2014).
- Fu and Kane (2008) L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008).
- Tanaka et al. (2009) Y. Tanaka, T. Yokoyama, and N. Nagaosa, Phys. Rev. Lett. 103, 107002 (2009).
- Sau et al. (2010) J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. Rev. Lett. 104, 040502 (2010).
- Alicea (2010) J. Alicea, Phys. Rev. B 81, 125318 (2010).
- Koren et al. (2011) G. Koren, T. Kirzhner, E. Lahoud, K. B. Chashka, and A. Kanigel, Phys. Rev. B 84, 224521 (2011), URL http://link.aps.org/doi/10.1103/PhysRevB.84.224521.
- Sacépé et al. (2011) B. Sacépé, J. B. Oostinga, J. Li, A. Ubaldini, N. J. G. Couto, E. Giannini, and A. F. Morpurgo, Nat. Comm. 2, 575 (2011).
- Qu et al. (2012) F. Qu, F. Yang, J. Shen, Y. Ding, J. Chen, Z. Ji, G. Liu, J. Fan, X. Jing, C. Yang, et al., Scientific Reports 2, 339 (2012).
- Williams et al. (2012) J. R. Williams, A. J. Bestwick, P. Gallagher, S. S. Hong, Y. Cui, A. S. Bleich, J. G. Analytis, I. R. Fisher, and D. Goldhaber-Gordon, Phys. Rev. Lett. 109, 056803 (2012), URL http://link.aps.org/doi/10.1103/PhysRevLett.109.056803.
- Cho et al. (2013) S. Cho, B. Dellabetta, A. Yang, J. Schneeloch, Z. Xu, T. Valla, G. Gu, M. J. Gilbert, and N. Mason, Nat. Comm. 4, 1689 (2013).
- Xu and et al. (2013) J.-P. Xu and et al., unpublished ArXiv/1312.7110 (2013).
- Zhao et al. (2014) L. Zhao, H. Deng, I. Korzhovska, J. Secor, M. Begliarbekov, Z. Chen, E. Andrade, E. Rosenthal, A. Pasupathy, V. Oganesyan, et al., unpublished arXiv:1408.1046 (2014).
- Farrell and Pereg-Barnea (2013) A. Farrell and T. Pereg-Barnea, Phys. Rev. B 87, 214517 (2013).
- Farrell and Pereg-Barnea (2014a) A. Farrell and T. Pereg-Barnea, Phys. Rev. B 89, 035112 (2014a), URL http://link.aps.org/doi/10.1103/PhysRevB.89.035112.
- Farrell and Pereg-Barnea (2014b) A. Farrell and T. Pereg-Barnea, Phys. Rev. B 90, 144518 (2014b), URL http://link.aps.org/doi/10.1103/PhysRevB.90.144518.
- Hoffman et al. (2002) J. Hoffman, K. McElroy, D.-H. Lee, K. M. Lang, H. Eisaki, S. Uchida, and J. C. Davis, Science 297, 1148 (2002).
- He et al. (2014) Y. He, Y. Yin, M. Zech, A. Soumyanarayanan, M. M. Yee, T. Williams, M. C. Boyer, K. Chatterjee, W. D. Wise, I. Zeljkovic, et al., Science 344, 608 (2014), eprint http://www.sciencemag.org/content/344/6184/608.full.pdf, URL http://www.sciencemag.org/content/344/6184/608.abstract.
- Vershinin et al. (2003) M. Vershinin, S. Misra, Y. Abe, Y. Ando, and Y. A., Bull. Am. Phys. Soc. 48, 372 (2003).
- Guo and Franz (2010) H.-M. Guo and M. Franz, Phys. Rev. B 81, 041102 (2010).
- Plamadeala et al. (2010) E. Plamadeala, T. Pereg-Barnea, and G. Refael, Phys. Rev. B 81, 134513 (2010).
- Pereg-Barnea and Franz (2008) T. Pereg-Barnea and M. Franz, Phys. Rev. B 78, 020509 (2008).
- Pereg-Barnea and MacDonald (2008) T. Pereg-Barnea and A. H. MacDonald, Phys. Rev. B 78, 014201 (2008).
- Marchand and Franz (2012) D. J. J. Marchand and M. Franz, Phys. Rev. B 86, 155146 (2012).
- Chang et al. (2015) C.-Z. Chang, Z. Zhang, K. Li, X. Feng, J. Zhang, M. Guo, Y. Feng, J. Wang, L.-L. Wang, X.-C. Ma, et al., Nano Letters 15, 1090 (2015).
- Capriotti et al. (2003) L. Capriotti, D. J. Scalapino, and R. D. Sedgewick, Phys. Rev. B 68, 014508 (2003), URL http://link.aps.org/doi/10.1103/PhysRevB.68.014508.
- Lang et al. (2002) K. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, H. Eisaki, S. Uchida, and J. C. Davis, Nature 415, 412 (2002).
- Peskin and Schroeder (1995) M. Peskin and D. Schroeder, An Introduction to Quantum Field Theory, Advanced book classics (Addison-Wesley Publishing Company, 1995), ISBN 9780201503975.

## Appendix A Impurity Potentials

Here we will describe the mathematical expressions we have used to describe the impurity potential. Starting with the charge impurity we define where acts on spin and any other degree of freedom in the system and acts on Nambu space. The operator acts, in a suitable manner, on any other degrees of freedom that may be present (i.e. not spin or particle-hole). For example, in our lattice model acts on the top-bottom surface degree of freedom while in the continuum model there are no other degrees of freedom left. In the lattice model we consider an impurity localized on a particular edge and so is diagonal with entries of on the impurity edge and on the other edge.

Secondly, we are interested in magnetic impurities. These alter the local Zeeman splitting on a single site and so couple to the spin of the electrons. We consider three separate impurities (one for each cartesian direction) while noting that a general magnetic impurity can be written as a linear combination of these potentials. As such we define (in Nambu space) where .

Finally, we consider our treatment of the OP suppression perturbation. This is modelled by taking where we remind the reader that act on spin degrees of freedom, on Nambu space and is as defined above. We calculate the OP QPI pattern by using this potential in Eq. (9).

Let us close this subsection with a brief description of the role of the probe potential, . In the main text we have imagined an experimental set-up where the physical STM tip is capable of resolving the component of the electron spin along a particular projection. In this case we are interested in only certain components of the change in the Greens function matrix. To find these relevant components and how they contribute to an interference pattern we must place an operator in the trace in Eq. (9) which acts to find the proper contribution. For example, if we are interested in an STM tip which resolves the -component of the electron spin we would add a to Eq. (9). In general, we call this matrix where is defined in the same way as above. For we are considering only a normal STM tip, while for we consider an STM tip capable of resolving the spin in the , , or direction respectively. In either case we additionally consider the physically relevant case of the tip only resolving degrees of freedom on one edge of the system; the tip can only make contact with one edge or the other.

## Appendix B Exact Green’s Function for Continuum Model

Here we will outline our calculation of the exact, clean Green’s function for the simple Dirac continuum model. To find the Green function of we begin by diagonalizing . This is accomplished by making the transformation where

(25) |

where is the phase of . In the above the operators ( ) annihilate electrons from the upper (lower) Dirac cone. Using this transformation we can rewrite where and

(26) |

In terms of this transformation the Green’s function can be written

(27) |

where is the Green function in the basis. After making this change of basis we can write as

(28) |

where

(29) |

where . We see from the above matrix that the two bands are completely decoupled from each other and we have -wave pairing on each Dirac cone. We have the two independent systems

(30) | |||

Defining the two Green’s function where it is straightforward to show

(31) |

where . Our Full Green function is then

(32) |

where

(33) | |||

Now we apply the matrices to get back to the Green’s function in a spin basis. It reads

(34) |

where and . Performing the matrix multiplication one can show that

(35) | |||

and

(36) | |||

where and .