1. Introduction
Superconductivity [
1] is at the heart of many quantum technological applications [
2]. For example, SQUIDS [
1] and nano-SQUIDS [
3] represent extremely useful sensors for magnetic fields, hence finding application in the most diverse settings. Additionally, QuBits based on superconductors [
4] are prominent in quantum computation, a field gaining more and more relevance as quantum supremacy appears to be a concrete goal. Superconductive correlations are crucial for the generation of Majorana zero modes [
5,
6,
7] and parafermions [
8,
9], that could alleviate the demanding needs of error correction: Indeed, such topological quasiparticles are characterized by non-abelian exchange statistics – which is expected to open the way to topologically protected quantum computation protocols – although it is fair to remark that such non-abelian statistics has not been experimentally detected so far. A further renowned application of superconductivity is superconducting spintronics [
10]. The aim of this sub-field of spintronics is to merge the advantages of spintronics on standard electronics, with the properties of Cooper pairs (CPs). Furthermore, superconductors find room in caloritronics [
11]. The list of applications is however continuously updating.
Recently, superconducting structures with broken time-reversal and inversion symmetry [
12,
13,
14] have been shown to allow for two connected and promising behaviors: the
Josephson effect [
15,
16,
17,
18,
19,
20,
21,
22] and the superconducting diode effect [
23,
24,
25]. The first consists in the presence of a finite Josephson current, called anomalous Josephson current, in the absence of a phase difference between the two superconductors and has recently been experimentally observed [
26]. It is remarkable because it can be used to design phase batteries [
27,
28], and to drive superconducting circuits [
29] and superconducting memories [
30]. The second, experimentally observed as well [
31,
32,
33], inspires great perspectives, given the analogy to its semiconducting counterpart that basically opened the way to modern electronics [
34]. While the breaking of inversion symmetry can be provided by the geometry of the nanostructure [
35,
36] or by the microscopic lattice [
37,
38], the breaking of time-reversal symmetry is usually, although not always [
39,
40,
41,
42], driven externally by means of applied magnetic fields.
Many of these striking advances in the functionalities of superconducting structures are related to the superconducting proximity effect [
1], which consists in the induction of a superconducting pairing in non-superconducting materials put in close contact to a superconductor. This has remarkable significance, because CPs in the proximitized material do not need to have the same properties as in the proximitizing superconductor. As an implication, despite most known superconductors host zero-momentum spin-singlet CPs, very exotic superconducting states can be implemented in proximitized materials. In these regards, materials with strong spin-orbit coupling play a special role. Spin-orbit coupled quantum wires, for example, enabled the realization of spinless
p-wave superconductivity [
6] with the purpose of realizing Majorana fermions. A second class of materials with strong spin-orbit coupling is that of topological insulators (TIs), for which the effects of proximization with superconductors are countless and extremely promising [
43]. Within this vast literature, a major role is played by two-dimensional topological insulators (2DTIs), where a gapped two-dimensional bulk coexists with symmetry protected metallic edge channels at the boundaries of the structure. 2DTIs have been successfully proximitized, and intriguing effects, such as missing Shapiro steps, have been observed [
44].
In this Article, we further analyze the properties of proximitized 2DTIs. Inspired by recent developments in the field of non-proximitized TIs, we propose a novel type of Josephson junction (JJ). Its main building block is represented by a narrow constriction between the edge states characterizing the 2DTI. Such edge states are reconstructed, meaning that the metallic states characterized by different spin polarization are spatially separated. Notably, we analytically show that the system hosts a Josephson effect in the absence of external magnetic fields. Moreover we show that, surprisingly, the effect is reinforced by a temperature increase in a range of parameters. Indeed, for low temperature and weak tunneling between the topological edge channels, the anomalous current has a small zero temperature contribution and a quadratic one. We interpret such temperature dependence on the basis of a simple perturbative argument.
The rest of the Article is structured as follows: In
Section 2, we introduce the main ideas of our proposal by qualitatively describing the setup, its symmetries, and its working mechanisms. Subsequently, we provide a theoretical model for the system and its properties. In
Section 3, we provide a quantitative analysis of the emerging
effect. In
Section 4, we discuss and interpret our findings, even on the basis of perturbative expansions and simple arguments. Finally, in
Section 5, we summarize and draw our conclusions.
2. Methods
The setup we consider is shown in
Figure 1 and consists of a JJ with a 2DTI as a normal part [
44,
45].
The chemical potential of the TI is tuned inside the bulk gap, such that the relevant degrees of freedom are the ones related to the one-dimensional edge channels [
46,
47,
48]. Crucially, such channels are helical,
i.e. on each edge electrons with opposite momentum have opposite spin projection. With reference to the picture, on the upper edge, right (left) moving electrons have spin up (down) projection. On the lower edge, the situation is reversed,
i.e. electrons with spin down (up) move right (left). At this stage, inversion symmetry is already broken at the level of the single edges, but time-reversal symmetry is not [
49]. To break the latter, we follow a road provided by the recent advancements in the creation of constrictions between topological edge channels [
50]. The mechanism, which has three conceptual steps, is the following:
By means of etching, the topological edge channels are brought at a distance that is comparable with their localization length. Electrons can then tunnel between the upper and lower edge. This can happen in two ways [
51,
52,
53]: via spin-preserving backscattering and via spin-flipping forward scattering.
The constriction is made long with respect to the inverse Fermi momentum. This step is, to some extent, unavoidable: Currently, constrictions between helical edge states have only been performed in thick HgTe quantum wells, where the Dirac point is hidden in the bulk valence band [
50]. The spin-preserving backscattering is hence irrelevant, and the only process remaining is the spin-flipping forward scattering [
50], denoted by
(with
) in
Figure 1. Even alone, such process is in any case still time-reversal invariant.
The third step is to induce edge reconstruction: As proposed in Ref. [
54], if the potential confining the edge channels is soft, the spin up and spin down channels can separate in real space, thus creating an unbalance in the tunneling rates for spin up and spin down electrons (see
Figure 1). Time-reversal symmetry is hence also broken, as well as inversion symmetry at the level of the full structure.
Thanks to these three considerations, the JJ with reconstructed helical edges as a link is hence expected to exhibit a Josephson effect in the absence of external fields.
Such anomalous effect is in a way similar to the one previously discussed in [
55], where, however, a magnetic field is required. Moreover, in that case, the effect is mostly present when the weak link is a single edge, while it tends to vanish when both edges are considered.
We now provide a quantitative analysis of what we have just described. To begin, we introduce the fermionic operators
that annihilate an electron at position
x propagating in the
-direction channel of the
edge [
56,
57]. We set
for the right (left) direction of motion and
for the upper (lower) edge. Due to the spin-momentum locking, these two indices completely define the edge states, since the spin polarization is
. For notational convenience, we introduce the spinor
. The Hamiltonian can hence be recast in the quadratic form
Here,
describes the 2DTI and
the coupling to the superconductors, that are considered of
s-wave type. Explicitly, we have
where
L is the separation between the superconductors and
is the Hamiltonian density of the edge channels, made up of two contributions:
, that are the kinetic term and the forward scattering term, respectively. In detail,
In the equations above,
and
(with
) represent the identity and the three Pauli matrices acting in the right/left-movers and the upper/lower edge space, respectively, and
is the Fermi velocity. The couplings
and
parameterize the tunneling between the edges: when
the edges are reconstructed and time-reversal symmetry is broken. Our model does not include magnetic impurities, that would introduce intra-edge backscattering, due to the fact that long and ballistic edges can now routinely be synthetized. Furthermore, we neglect the possible presence of charge puddles in the bulk of the 2DTI since there is little room for them in a narrow sample. To summarize, the non-superconducting part just discussed represents a pair of one-dimensional Dirac cones, in which the two right-moving branches are shifted with respect to the left-moving ones proportionally to
. Moreover, branches with fixed chirality do not have fixed spin projection anymore, due to the spin non-conserving tunneling.
The lateral proximitizing superconductors, assumed to have a large superconducting gap with respect to all the other energy scales, influence the system via the perturbative Hamiltonian
. Besides the large gap hypothesis, its derivation requires a weak coupling between the edges and the superconductors, and the specific model of the latters. In our case, it is the standard BCS one, since in experiments on 2DTIs niobium is the prime example of employed superconductor [
44]. What one obtains is that the coupling to the right and left superconductors (
) is given by [
58,
59]
In Eq. (
4),
labels the right/left superconductor, and
are collective indices:
. Moreover, we have introduced
, where
is the coherence length in the edges, which will be the short distance cutoff of our system,
, and
is the superconducting gap. This splitting makes processes with tunneling of spin-triplet CPs into or out from the same edge possible without violating the Pauli principle. It is here evident the general appeal of proximity-induced superconductivity: the symmetry of the superconducting pairing in the proximitized system can be different from the original one in the proximitizing superconductor. Finally, the last summation in Eq. (
4) has been antisymmetrized for each
term, so that
for
, and
for
. Following this requirement, the last sum runs over only 10 terms, whose coefficients read as [
58]
Here,
, where
t parametrizes the magnitude of tunneling across the superconductor/TI interface and
is the normal density of states per spin in the superconductors at the Fermi energy, is the typical amplitude of the CP tunneling processes,
are the phases of the superconductors (
),
is the Fermi momentum,
and
are coefficients related to the occurrence of spin flips [
60] and tunneling into different edges (crossed Andreev reflection) [
61], respectively. Of particular interest for the following is
, since it allows for the propagation of triplet CPs on a single edge and hence implements processes lacking inversion symmetry. A finite
, expected due to the strong spin-orbit coupling characterizing HgTe quantum wells, is hence essential for the occurrence of the
Josephson effect in the system. Although expected as well,
is not crucial for the following. As an only comment about such term, it is here worth to mention that
is expected to become significant if a perpendicular magnetic field is piercing the junction. Indeed, the presence of
affects the periodicity of the critical supercurrent as a function of the magnetic flux [
59], since it implements a physical behavior similar to the one of the so called nano-SQUIDs [
3]. The main classes of superconducting tunnelings are schematically shown in
Figure 2.
Having set the Hamiltonian of the system up, we can now move to the main object of our investigation: the DC Josephson current. To perform the calculation we rely on a standard perturbative Kubo-like approach that takes into account the Hamiltonian of the edge channels as the unperturbed system and the tunneling from the superconductors up to order
, the lowest non-zero contribution (see [
59] for all the details of the method). The calculation, lengthy and cumbersome but straightforward in essence, leads to the Josephson current
The function
has the following structure
and, from now on, an overbar denotes complex conjugation. Here, the integral evaluates the amplitude associated to the transport of electrons in the initial state
and in the final state
through the junction, hence encoding all the interesting information. More quantitatively,
arises as a combination of products of two Green functions (see [
59] for the details) calculated over the Hamiltonian in the absence of superconductors. For convenience, we have rescaled all the parameters via the introduction of the adimensional quantities
with
T the temperature and
the Boltzmann constant. In Eq. (
6), the natural scale is set by
, and the phase
is the difference between the bare phases of the two superconductors.
3. Results
We now proceed to the analysis of the most striking effect related to Eq. (
6), that is the anomalous Josephson current flowing in the system when the phases
are set to be equal. To better characterize such an effect, we first recast the current in the form
where
has been introduced in Eq. (
8), and we have defined
In Eq. (
10), we have hence singled out the parameter
, that is the one describing the anomalous Josephson effect. When
mod
there is no anomalous Josephson effect. The effect is present otherwise. For completeness, we mention that a notable value in this context is
, where one has the so called
-junction, that has less symmetry requirements with respect to the generic case. In this work we however do not put our focus on such special case.
To give a quantitative description of the
effect we predict, we now inspect the current for a vanishing phase difference between the superconductors
. What we explicitly have is then
It is clear from this relation that processes contributing to
symmetrically under the exchange
, do not produce anomalous current (
). Moreover, we find that the processes responsible for the anomalous effect are
and independent of
(more details and our interpretation can be found below). They hence correspond to injections as depicted in the upper right panel of
Figure 2. Quantitatively, we have
where the coefficients
are listed in the table that follows.
Coefficient
|
Corresponding
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
It can be shown that
and that, more generally, they can all be written in an unified form. The explicit form of
and other computational details are provided in
Appendix A for clarity. All in all, the anomalous current characterizing the system reads as
with the function
given by
4. Discussion
The anomalous current is shown in
Figure 3, as a function of temperature and tunneling amplitudes.
First of all, we note that the anomalous current does not have a well defined sign. Indeed, as expected, it is odd under the exchange
, which essentially realizes the inversion of the structure. Accordingly, we find
for
, and in that case Eq. (
9) recovers the standard Josephson effect. This behavior could be expected since for
time-reversal symmetry is preserved. An oscillatory behavior is also present as a function of
and
. This fact is not surprising since oscillations in the amplitude of the DC and AC (non-anomalous) Josephson currents are also present for
[
59] (
Figure 3(a)). As a second interesting fact, thanks to such oscillations,
increases rapidly as soon as an asymmetry between the couplings, although small, is introduced. This is visible in
Figure 3(b), where we set
: the anomalous current assumes a finite value as
deviates, even slightly, from 0.2. Finally, it is here worth noticing a peculiar effect: the anomalous current increases as the temperature is increased, keeping however
T smaller than the induced gap (
Figure 3(c)). This effect can be understood with a simple “perturbative” argument: the scattering between the edges, reported in Eq. (3b), mixes states with the same kinetic energy according to Eq. (
3a). Hence, in a Fermi’s golden rule approach, its effects are strongly suppressed at low temperature due to the Fermi functions. Such a temperature dependence is hence a clear cut signature of the constriction. Lastly, we point out that, going beyond the low-transparency assumption, and hence including higher order contributions of
in the calculation of the current, would result in a superconducting diode effect [
62].
In order to capture the essential physics underlying the temperature behavior, some expansions of Eq. (
13) are particularly helpful. First of all, the anomalous current has a zero-temperature contribution, although hardly visible in
Figure 3(c):
Restoring finite temperature, and in the long junction limit
, we get
On the other hand, without assumptions on
but for small
and
, we can further expand Eq. (
13) as
that is the temperature scaling visible in
Figure 3(c).
We are now in the position of interpreting the temperature dependence of the anomalous current. At zero temperature (see Eq. (
14)), for small
, the inter-edge tunneling events are suppressed by phase space arguments. The coupling between the edges, for
, hence only splits the energies of electrons with different chiralities. In terms of the spin degree of freedom, it is thus analogous to the application of two magnetic fields in the
z direction opposite to each other on the two edges. It hence generates a
effect analogous of the one reported in [
55]. However, thanks to the fact that the effect is now opposite for the two edges, the anomalous Josephson current does not vanish when both edges are taken into account. The energy splitting just commented is independent of the length
of the system and, consistently, the associated current decays as the distance between the superconductors is increased (see Eq. (
15)), as usual for Josephson currents. When the temperature is finite (see Eq. (
16)), a different mechanism for the
effect is added to the one just discussed: The inter-edge tunneling events. Just as the zero temperature contribution, the finite temperature one scales as
. It thus involve two inter-edge tunneling events. The probability of each tunneling event scales as
, so that the finite temperature term scales, with respect to the zero temperature one, with a factor
more. Moreover, each tunneling corresponds to a product
, with
the Fermi function. Under integration, such product contributes with a factor
, and since we are in the presence of a double tunneling, the leading order becomes
. To summarize, the probability of the events contributing to the temperature activated anomalous current increases with both junction length and temperature. Our qualitative interpretation is hence able to capture the behavior of the
effect at finite temperature, its non-zero value at zero temperature, and the fact that the temperature activated processes have a better scaling with the length of the junction.