2
II. MICROSCOPIC HAMILTONIAN
The Josephson junction consists of a conventional single
band s-wave superconducting region, a normal insulator re-
gion, and the superconducting SRO region. A schematic of
the regions is shown in Fig. 1(a), and they are denoted A,
B, and C, respectively, in the equations below. Within re-
gions Aand Bwe use the single band kinetic Hamiltonian,
written here for region A,
HA=X
k,iy,δy,σ
ξA(k, δy)c†
A,k,iy,σcA,k,iy+δy,σ (1)
where ξA(k, δy) is the electron dispersion in region Abe-
tween slabs at positions iyand iy+δy, and c†
A,k,iy,σ cre-
ates an electron in slab iyof region Awith momentum
k= (kx, kz) and spin σ. The form of all dispersion terms
as well as the values of the hopping parameters are given
in Appendix A. The superconductivity in region Ais de-
scribed by,
HSC
A=X
k,iy
eiφA∆Ac†
A,k,iy,↑c†
A,−k,iy,↓+ h.c., (2)
where ∆Ais the s-wave order parameter, and φAis the
superconducting phase. Hopping between regions Aand
Bis taken to have the same parameters as hopping within
either of the regions.
The normal state Hamiltonian of region Cis,
HC=X
k,iy,a,σ
ξa
C(k, δy)ca†
C,iy,k,σca
C,iy+δy,k,σ
+X
k,iy,a6=b,σ
ξab
C(k, δy)ca†
C,k,iy,σcb
C,k,iy+δy,σ + h.c. + HSOC,
(3)
which includes intraorbital and interorbital dispersions,
ξa
C(k, δy) and ξab
C(k, δy), respectively, where aand bare
the orbital indices representing the yz,xz, and xy orbitals,
as well as SOC terms. Finally, the Hamiltonian describ-
ing hopping between regions Band C, where the interface
occurs between iy= 1 and 2, has the form,
Hint =X
k,a,σ
ξa
int(k)ca†
C,k,2,σcB,k,1,σ + h.c., (4)
which features orbital dependence in the SRO region, as
denoted by ξa
int(k).
We also consider the effects of mirror symmetry break-
ing in the z-direction (out of plane direction). This effect
is largest near the surface normal to the z-direction, but
imperfections of the interface leading to a broken mirror
plane in the z-direction have previously been proposed to
occur [37], where it was shown that the experimental re-
sults may be explained by a dxz +idyz-wave state if the
tunneling directions tilt out of the xy plane. Additionally,
the growth of Au0.5In0.5directly onto SRO to create the
junction may cause strain in SRO. Any deformations of
the lattice that this leads to may further contribute to bro-
ken mirror symmetry throughout the sample. Dislocations
may also contribute to this mirror symmetry breaking, and
have been found to occur near interfaces [38]. The lack of
mirror symmetry in the z-direction means that hopping be-
tween the xy and xz(yz) orbitals is allowed to be finite in
the y(x) direction. An example of such hopping is shown
in Fig. 1(b) for the xy to xz interorbital hopping. These
hoppings have the form:
hISB
k=−α[2isin kx(cyz†
k,iy,σcxy
k,iyσ−cxy†
k,iy,σcyz
k,iy,σ)
−(cxz†
C,k,iy,σcxy
C,k,iy+1,σ −cxy†
C,k,iy,σcxz
C,k,iy+1,σ + h.c.)],(5)
where αrepresents the hopping integral, which depends on
the strength of the mirror symmetry breaking, and the use
of δy= 1 here represents nearest-neighbor hopping between
slabs.
In the next section, we describe the setup for the Joseph-
son tunneling calculations and apply it with conventional
s- and p-wave superconducting states in the SRO region.
Then, we consider interorbital superconductivity and show
how the current-phase relation (CPR) is affected by the
mirror symmetry breaking terms introduced here, showing
that they may behave like the conventional s-wave state,
or potentially more like the p-wave state, depending on the
nodal structure as well as the strength of the mirror sym-
metry breaking.
III. JOSEPHSON CALCULATIONS
To calculate the Josephson CPR, we use the lattice
Green’s function method presented in Ref. 34, which con-
siders only p-wave pairing to explain experimental results.
Semi-infinite Green’s functions are obtained for the s-wave
and SRO regions using the recursive Green’s function ap-
proach [39]. A single layer of the normal insulator is added
on the surfaces of both of these regions by the Dyson equa-
tions,
ˆ
GB
0(k, iωl)=(iωl−ˆu0(k)−ˆ
t0,−1ˆ
GA
−1(k, iωl)ˆ
t−1,0)−1,(6)
ˆ
GB
1(k, iωl)=(iωl−ˆu1(k)−ˆ
t1,2ˆ
GC
2(k, iωl)ˆ
t2,1)−1,(7)
where the interface is in the xz plane. Here, ˆ
Gm
n(k, iωl) is
the Green’s function of layer nin region m, ˆun(k) is the
part of the Hamiltonian of layer n, and ˆ
tn,n+1 is the part
of the Hamiltonian featuring hopping the the y-direction,
representing hopping between layers nand n+ 1. The left
and right systems are combined using the two equations,
ˆ
G00(k, iωl) = {[ˆ
GB
0(k, iωl)]−1−ˆ
t01 ˆ
GB
1(k, iωl)ˆ
t10}−1,(8)
ˆ
G11(k, iωl) = {[ˆ
GB
1(k, iωl)]−1−ˆ
t10 ˆ
GB
0(k, iωl)ˆ
t01}−1.(9)
These are then used to obtain the nonlocal Green’s func-
tions
ˆ
G01(k, iωl) = ˆ
GB
0(k, iωl)ˆ
t01 ˆ
G11(k, iωl),(10)
ˆ
G10(k, iωl) = ˆ
GB
1(k, iωl)ˆ
t10 ˆ
G00(k, iωl).(11)