1. Introduction
Diode is one of the basic devices in semiconductor technology, and is characterized by different amplitudes of currents flowing in opposite directions. It depends on the lack of an inversion symmetry center in a system and is one of the central building blocks required to build powerful instruments, such as a transistor. In many superconductors-based systems, a similar effect named as superconductor or Josephson diode effect (JDE) appears due to broken time-reversal and inversion symmetries [
1,
2,
3,
4,
5,
6]. The JDE can result in direction-dependent (nonreciprocal) critical currents driven by phase difference between superconductors and is fascinating in various applications due to its connection to fundamental properties of diverse superconducting systems [
1,
2,
3,
4,
5,
6,
7,
8]. Generally, time-reversal symmetry may be broken by an external magnetic field, which exerts impacts on specific term in the systems’ Hamiltonian related to inversion symmetry breaking [
9], such as spin-orbit interaction (SOI). Therefore, the JDE can be used to detect SOI strength [
10] or existence of a topological phase in nanowires in proximity-contacted with a superconductor [
5,
10,
11]. The JDE is also promising in design of electronic devices including photodetectors, ac/dc converters, superconducting qubits [
8], and racetrack memory devices [
12] etc.
The JDE was demonstrated as early as 1970s in superconducting quantum interference devices (SQUIDs) based on superconductor bridges [
13] and Josephson junctions [
14]. It was also observed in non-centrosymmetric conventional superconductor thin films in the device geometry [
15]. Since then, many experimental and theoretical platforms for the realization of JDE has been continuously proposed, such as non-centrosymmetric superconductors [
14,
15], stacks of different superconductors with broken inversion symmetries [
16], Andreev molecules [
17], artificial superlattices [
2], topological semimetals and insulators [
5], insulator heterostructure devices [
18], nanowires [
10], and disordered systems [
19]. Another efficient platform proposed to achieve high diode efficiency relies on a Josephson current interferometer, in which conjunct Josephson junctions with nonsinusoidal current-phase relations form a SQUID. In such systems, the Josephson currents are contributed from higher harmonics other than the usual
-periodic ones with a magnetic flux penetrating through the SQUID loop [
6,
8,
14]. It has recently been demonstrated in two-dimensional electron and many three- and four-terminal setups in which the diode efficiencies at equilibrium can reach up to about 30% [
20,
21].
Conventional s-wave superconductors in proximity contacted with semiconductor nanowires, which have strong Rahsba spin-orbit interaction and are subjected to external magnetic fields, have been demonstrated for the realization of topological superconductivity [
6,
8,
10,
11], a topological state that hosts Majorana bound states (MBSs) [
22,
23]. The MBSs are charge-neutral topologically protected quasiparticles that are spatially separated, i.e., character of spatial nonlocality [
24]. They offer an attractive way for constituting Majorana qubits [
25,
26] allowing to store information in a nonlocal manner and immune to decoherence by a local disturbance [
27]. Along with the intensive investigations on preparation and detection of MBSs, exploiting their possible applications is also an active research subject in condensed matter physics. For example, in systems with MBSs side-coupled to quantum dots (QDs) which are connected to external leads, the sign of thermopower which measures the induced bias voltage in response to a temperature difference can be reversed by QDs-MBSs hybridization strength or MBS-MBS overlap amplitude [
28,
29]. Moreover, the magnitudes of both the thermopower and thermoelectric efficiency in such systems can be obviously enhanced [
28,
29,
30]. Since the MBSs often emerge with the help of spin-orbit interaction, strong magnetic field or magnetic materials, they also play an important role in the research field of spintronics [
31]. In recent years, some works have been devoted to the study of Josephson current through a QD connected to two semiconductor nanowires hosting MBSs (MNWs) [
32,
33,
34,
35,
36], which are stimulated by the interesting results found in various system composed of QDs connected to superconductors with normal phase [
37,
38,
39,
40,
41,
42,
43,
44,
45,
46]. It was shown that the Josephson current driven by the topological phase difference is quite stronger than that by normal phase difference, and the bent angle formed by the two MNWs as well as the magnetic fields in the QD will significantly suppress the Josephson current [
32,
33].
Until now, there are two limitations in previous work on Josephson current through MNWs-QD-MNWs systems: one is that only the current’s magnitude and period were studied, whereas less attention has been paid on the control of its direction [
32,
33,
34,
35,
36]; the other is that only single QD was proposed to be inserted between two MNWs, and the interesting quantum interference effect was left untouched [
14,
16,
20,
21,
38,
40,
41,
42,
43,
46,
47]. In view of these, we propose a structure composed of parallel double QDs (DQDs) sandwiched between the left and right MNWs to achieve the JDE based on quantum interference effect arising from the two transport paths through the DQDs, as well as the magnetic flux penetrating through the loop, which is shown in
Figure 1. We emphasize that the present device can be experimentally realized in terms of that reported in Ref. [48], in which DQDs connected to two normal leads couple to each other via both superconductor hosting MBSs and normal tunnel barrier. If the leads in Ref. [48] are replaced by the superconductor connecting the DQDs, their system is exactly the one studied here. In some previous work, DQDs has been proposed to be inserted between conventional superconductor leads to generate spin-correlated electron pairs, and to control the Josephson current and its critical one [
40,
41,
43,
44,
49]. Our studies show that in this MNWs-DQDs-MNWs, the period, magnitude, and the directions of the Josephson current can be effectively adjusted with the help of dots’ energy levels, overlap amplitude between the MBSs, as well as the magnetic flux through the loop. Accordingly, tunable JDE with large value of diode efficiency emerges, and may find real use in design of superconductor-based instruments.
2. Model and Method
The Hamiltonian of the present structure is divided into three parts as
[
32,
33,
39,
40], in which the Hamiltonian of the DQDs and interaction between them is given by
where the creation (annihilation) operator
is for electrons in dot-
i with spin-independent energy level
. In experiments, the dots’ level can be tuned via gate voltages
, and is given by
, where
is the bare energy level in dot-
i. The tunnel coupling strength between the DQDs is
. The Hamiltonian
denotes the left and right MNWs connected to the DQDs, whose explicit expression is as follows [
32,
33,
50],
where
is for direct hybridization strength between the MBSs prepared at the ends of the
-th nanowire. In what follows, we assume
. The creation and annihilation operators of the MBSs satisfies the relationship of
, and
[
32,
51] due to the unique self-conjugate character of the MBSs. Tunnel coupling between the DQDs and the MNWs is described by the Hamiltonian of
, in which [
32,
33]
in which
stands for coupling strength between QD-
i and the MNW-
. Note that there is a phase factor
in
, which arises from the proximity of the MNWs to the superconductor substrates and induces the Josephson current, with
. In the present paper, we set
and
for the sake of clarity [
32]. As usual, we make an unitary transformation to change the Majorana fermion representation to a conventional fermion representation one: [
32,
51]
and
. Then we rewrite the total Hamiltonian in a matrix form to calculate the Green’s functions needed for the Josephson current. In the basis of
, the transformed Hamiltonian
is given by [
32,
33,
34,
35]:
in which the
sub-matrix
diag
with
being the Pauli matrix of the
-component, and the symbol ⊗ denotes the matrix direct product. Similarly, the
sub-matrix
, and the interaction between the DQDs and MBSs are
The dc Josephson current
J tunneling between the two MNWs via the DQDs is calculated in terms of the nonequilibrium Green’s function technique [
32,
37,
38,
39,
40]
where
, and
with the retarded/advanced self-energies due to the MNWs. The self-energies are calculated by
[
19,
26,
29]. The free retarded/advanced Green’s function of the MNW is
. The retarded/advanced Green’s function of the DQDs is then obtained with the help of Dyson’s equation as [
32,
37,
38,
39,
40]
[
32,
39]. In Equation (
6),
is the equilibrium Dirac-Fermi function, where
T and
denote the temperature and Boltzmann constant, respectively.
3. Numerical Results
In numerical calculations, we consider the case of the DQDs couple to the left and right MNWs with equal strengths
, which is set to be the energy unit. In
Figure 2 we present the gate voltage
dependence of the Josephson current
J with different direct hybridization coupling strength between the MBSs
and energy levels
for fixed
and
. First, we diagonalize the Hamiltonian
for the DQDs and transform the QDs’ discrete energy states into a pair of bonding and antibonding states,
in which
is the annihilation operator for the bonding (antibonding) state of the DQDs with energy
where
and
.
For the case of QDs’ energy level configuration of
and
, there is only one current peak associated with the bonding state with energy
. Therefore, the peak position will shift towards lower (higher) energy regime with increasing (decreasing)
. Note in this totally symmetrical structure, i.e., the energy levels of the DQDs are the same (
) and are coupled to the MNWs with equal strengths
, the antibonding state
disappears and only the bonding state contributes to the Josephson current through the DQDs. When the MBSs at the opposite ends of the MNWs are hybridized to each other
, the position of the Josephson current peak remains unchanged as shown in
Figure 2a. With increasing
, the current’s amplitude is reduced which is in consistent with previous results [
32,
33]. Moreover, the width of the peak is narrowed by increasing
. To explain the dependence of the Josephson current on the direct hybridization of the MBSs
, we shows the current-carrying density of states (CCDOS)
varying with respective to energy
. Since at zero temperature, only the states located in negative energy regime contribute to the Josephson current, we only show the behaviors of
for the case of
. When the MBSs at the MNWs are decoupled to each other
, there are three peaks in
, of which two positive and one negative as is indicated by the black solid line in
Figure 2b. When
, a negative peak is induced in
around the Fermi energy of the MNWs
, and its height is increased by increasing
. As a result of it, the current’s amplitude is reduced. The additional current peak originates from the fact that the MBSs in the same MNWs are destroyed by their overlap.[
50,
51]
If the energy levels in the two QDs are different from each other
, there are two current peaks corresponding to the two molecular states
in the DQDs as shown in
Figure 2c,d, in which
are set to be (0,-4) and (-2,4), respectively. The height of the current peak at bonding state
is lower than that at the antibonding state
. When the MBSs in the same MNWs are overlapped (
), the current peaks in
Figure 2c,d are lowered and narrowed as in the case of
Figure 2a. Moreover, the positions of the peaks are slightly changed by
for the cases of
, which is different from the case in
. It should be noted that the reduction of the Josephson current peaks in
Figure 2c,d is mainly induced by charge flowing through the two molecular states
. This is different from the case in
Figure 2a, where the charge flows only through the bonding state
. It means that electrons will transport through the DQDs whenever their energy is in resonant with the bonding and antibonding states, even if
or
is away from the Fermi level.
We now study the influences of the magnetic flux
and phase difference
on
J, critical Josephson currents
, as well as the diode efficiency
[
45] for different dots’ level configurations for fixed value of
. For the case of
, the oscillation period of both the Josephson current
J and its critical counterpart
is
versus either
or
as shown in
Figure 3a,b. The Josephson current
J in
Figure 3a has an abrupt jump from positive to negative value at
. This is identical to the case in systems of a single QD sandwiched between two MNWs.[
32,
33] The current
J in
Figure 3a also has a maximum value at
, and does not change its sign when
is varied. This behavior is also identical to the structure in which the DQDs are connected to conventional non-topological superconductors (S-DQDs-S)[
38,
39]. The positive and negative critical currents as functions of
in
Figure 3b are antisymmetrical with respective to each other, i.e.,
, hence the diode efficiency
as is indicated by the blue dotted line therein. For the cases of dots’ energy levels of (
) and
, the Josephson current
J in
Figure 3c,e are
-period function of
and
, and its sign depends on the values of both
and
These results are similar to those found in S-DQDs-S [
38,
39]. It means that the normal
-period Aharonov-Bohm oscillations such as for
is destroyed and complex periodic interference effects occurs by changing the dots’ levels [
40]. Correspondingly, electron transport processes including current’s amplitude or directions are controllable by adjusting both the magnetic flux and the dots’ energy levels.
The period of the positive and negative critical currents
in
Figure 3d,f is the same as that of
J, which individually corresponds to the cases of (0,-4) and (-2,4). Now
and
are not antisymmetrical to each other, i.e.,
at most value of
, and the phenomenon of JDE emerges accordingly. We find that the diode efficiency is anti-symmetrical with respective to
and exhibits triple-peak configuration, of which two higher ones locate around
and two lower ones at the two sides of
, respectively. It shows that
depends on the energy levels of the QDs and the magnetic flux. For example, the absolute maximum of
for
in
Figure 3d emerged at about
can reach up to
. Whereas for the level configuration of
, the absolute maximum of the diode efficiency is about
at about the same value of
. The above results show that the JDE can be efficiently controlled by combined functions of the QDs’ energy levels and the magnetic flux which induces complex interference effects.
In
Figure 4, we present the impacts of MBS-MBS overlap amplitude
on the diode efficiency
with fixed QDs’ energy levels configuration
. From the figure one can see that the diode efficiency remains as a
periodic function of
in the presence of finite
, and is anti-symmetrical with respective to
. With increasing
, the peaks’ heights around
are lowered with almost unchanged locations. The double-peak configuration around
in the case of
, however, evolves into a triple-peak one for
as is indicated by the red dashed and blue dotted lines. Moreover, the two pairs of peaks, which are shifted individually to magnetic flux values of greater and smaller than
, keep as the same shape as those of
, i.e., one positive and one negative. The heights of the two pairs of the peaks change nonlinearly with respective to
. For the chosen values of
, the maximum of
first increases (by comparing the black dotted line for
with the red dashed line for
), and then decreases (the blue solid line for
). The results displayed in
Figure 4 indicate that by properly adjusting the value of
, one can change either the amplitude or the sign of the diode efficiency. In experiments,
depends on the lengths of the MNWs and superconductor coherence [
24], which can all be used for change the JDE.
Dependence of the Josephsonc current on the gate voltage
and phase difference
for fixed magnetic flux
and different values of bare dots’ energy levels are displayed in
Figure 5a for (
) and 5(c) for (
), respectively. In both of the cases of (
) and (
), the Josephson current varying a function of
is a
periodic function of
, and jumps from positive to negative or vice versa at a particular
. The value of the particular
depends on both the dots’ energy levels through the gate voltage
and the magnetic flux, whereas the amplitude of
J has no obvious change. This result is quite different from that in S-DQDs-S, in which the current is significantly reduced as the dots’ energy levels are tuned away from the Fermi level[
40,
41]. As for the diode efficiency
, it is zero in the case of
regardless of the value of dots’ energy levels
as is shown by the blue solid lines in
Figure 5b,d. The diode efficiency
develops two peaks with opposite signs for the cases of
and
, and their locations and heights are varied when the value of magnetic flux is changed. Here only the results of
are displayed, we have examined that for finite MBS-MBS overlap amplitude
, the behaviors of the diode efficiency are similar to those in
Figure 4 and we do not list them here.