
2
the model that is used for the underlying N-S-F-S-N junc-
tion. In Sec. II we recap the substance of WTDs. The re-
sults are presented in Sec. III, where we elaborate WTDs
for specific types of hybrid quantum conductors, closely
related to the Fu and Kane structure [7]. We study in
details the WTDs for both, ABS and MBS and identify
signatures to distinguish them. Conclusion and outlook
are given in Sec. IV and several technical details are avail-
able in Appendices.
I. MODEL
The heterostructure we consider consists of helical
counter propagating edge states of a Quantum spin hall
insulator, proximity coupled to a spatially restricted s-
wave superconductor and a magnetic field [see Fig. 1
(a)]. Similar setups have been studied in the litera-
ture [6,7,49–51]. In the basis Ψ = (ψ↑, ψ↓, ψ†
↑, ψ†
↓), the
Hamiltonian of the system under consideration is the fol-
lowing
H=H0H∗
SC
HSC −H∗
0,(1)
where H0=HTI +HZ. The Hamiltonian HTI =
vFkxσ1−µis the kinetic term of the edge states, with
Fermi-velocity vF(~= 1), momentum kx=−i∂xand
chemical potential µ, while HZ=σ1∆Z(x) is the Zee-
man field pointing in x-direction. Furthermore, HSC =
−i∆(x)σ2includes the superconductor, which is assumed
to be grounded, and the Pauli matrices σ(τ) act in
spin (Nambu) space. The spatial arrangements of ∆
and ∆Zare chosen to have a ferromagnetic region en-
closed by two superconducting regions, such that ∆(x) =
∆[Θ(x)−Θ(x−L)]+∆eiφ[Θ(x−L−Lm)−Θ(x−2L−Lm)]
and ∆Z(x)=∆Z[Θ(x−L)−Θ(x−L−Lm)]. The pairing
mechanisms open a gap in the spectrum of the edge states
[see schematic in Fig. 1(b)], which is either at the Dirac-
point (DP) (∆Z) or at the Fermi level crossing (∆). In
addition we also allow the possibility to add a magnetic
flux φto the system [see Fig. 1(a)].
We check the characterization of different states ap-
pearing in the heterostructure, by first presenting the
scattering features as a basis to identify distinct signa-
tures in the WTD. Whereby we use the seminal work of
Blonder, Tinkham and Klapwijk (BTK) to calculate the
corresponding scattering coefficients [52]. For our inter-
ests we study two cases, first, no Zeeman field is applied
(trivial phase) and the middle region hosts freely propa-
gating states, and second, a non vanishing Zeeman field
(topological phase above a certain treshold) couples the
two spin species and opens a gap in the middle region, re-
sulting in the emergence of MBSs on the interfaces of the
superconducting and magnetic regions [7]. We assume in
this work a rather long junction (Lm> vF/∆), such that
for ∆Z= 0 multiple bound states are well defined within
the superconducting gap E < ∆ [see Fig. 2(a)]. Cal-
culating the wavefunction with the continuity conditions
on the interfaces at x= 0, L, L +Lmand x= 2L+Lm
results in the necessary transport coefficients. With a
preserved TRS an incoming right moving electron from
the source is (∆Z= 0) protected from backscattering,
such that there are only two non vanishing processes.
Namely, the co-tunneling of an electron (te) form source
to drain and the local Andreev-reflection (rh) on the SC
interface at x= 0. The breaking of TRS in turn then
also allows normal electron reflection (re), including a
spin flip, and the transmission of holes (th). While the
gap closing and reopening is usually used as an indica-
tor of the topological phase transition, we use only the
normal electron reflection signature, accompanied by its
probability at zero energy from zero to one. Details are
presented in App. A. Based on those energy dependent
coefficients, we are highly interested in the corresponding
WTDs, which we introduce in the next section.
II. WAITING TIME DISTRIBUTION
For the sake of clarity, we recap in this section, the
formalism used for the calculation of WTDs [26–29] as
well as some well established results that will serve as
reference. The WTD W(τ) denotes the probability dis-
tribution for the time delay τbetween the detection of
two consecutive charge carriers. They can be electrons
or holes for instance. This quantities gives precise and
transparent informations about correlations in a trans-
port process. In general, it is customary to calculate it
from the Idle Time Probability (ITP) Π(τ), namely the
probability to detect zero particles during a time interval
τ. For stationary processes these two distributions are
connected by the following expression
W(τ) = hτi∂2Π(τ)
∂τ2. (2)
For non-interacting systems, the ITP is given by the de-
terminant formula [28]
Π(τ) = det1−Qτ, (3)
where Qτis a projector over the time window τwhose ex-
pression depends on the detection scheme (measurement
of consecutives electrons, holes, ...) and the scattering
matrix of the system. Explicit formulae are given be-
low for the processes of interest. In addition, the mean
waiting time is given by
hτi=−∂Π(τ)
∂τ
−1
τ=0
.(4)
To continue, we specify a detection procedure for the
WTD. Thus, we calculate the ITP with the projection
of a single particle scattering state on an discrete
voltage window around the Fermi level. In the range
of the applied voltage V, the linear energy spectrum
is discretized into Ncompartments with wave vectors
k=n
N
eV
~vF. A stationary process is then reached by the