Magnetoelectrics in Disordered Topological Insulator Josephson Junctions

Magnetoelectrics in Disordered Topological Insulator Josephson Junctions

Abstract

We theoretically study the coupling of electric charge and spin polarization in an equilibrium and nonequilibrium electric transport across a two dimensional Josephson configuration comprised of disordered surface channels of a three dimensional topological insulator. In the equilibriun state of the system we predict the Edelstein effect, which is much more pronounced than its counterpart in conventional spin orbit coupled materials. Employing a quasiclassical Keldysh technique, we demonstrate that the ground state of system can be experimentally shifted into arbitrary macroscopic superconducting phase differences other than the standard ’’ or ’’, constituting a -junction, solely by modulating a quasiparticle flow injection into the junction. We propose a feasible experiment where the quasiparticles are injected into the topological insulator surface by means of a normal electrode and voltage gradient so that oppositely oriented stationary spin densities can be developed along the interfaces and allow for directly making use of the spin-momentum locking nature of Dirac fermions in the surface channels. The -state is proportional to the voltage difference applied between the injector electrode and superconducting terminals that calibrates the injection rate of particles and, therefore, the shift.

pacs:
74.50.+r, 73.20.-r, 73.63.-b

I introduction

Topological state of matter has been a striking topic during the past decade and has attracted extensive attention both theoretically and experimentally rev1 (); rev2 (). The study of nontrivial electronic structure topology goes back to the prediction of the time-reversal symmetry protected electronic states bound to the contact of two semiconductors with mutually inverted bands pankratov () and exotic particles in liquid volovik (), which now offer promising prospects to practical platforms such as spintronics and quantum computation rev1 (); rev2 (). The extensive and growing research efforts in this context have led to the exploration of topological insulators (TIs) ti1 (); ti2 (); ti3 (); ti4 (); ti5 (), for instance. Surface channels of a three dimensional topological insulator are conductors pankratov () while bulk material itself is insulator with a gap in its band structure. In contrast to numerous theoretical works so far rev1 (); rev2 (), these surface channels may not be fully ballistic and can host unavoidable defects or other quasiparticle scattering resources such as impurities burkov (); Schwab (); zyuzin (). The recent experiments have also revealed this fact that the conduction of the surface states are limited with a finite resistance nowack (); spanton (). This key observation clearly proves the fundamental importance of considering the contributions and influences of disordered motions of moving particles in the surface channels nowack (); spanton (); burkov (); Schwab (); zyuzin ().

The surface channels may be employed as an appealing opportunity that can support spin-momentum locked modes, namely, the orientation of particle’s spin is locked to its momentum direction rev1 (); rev2 (). The motion of low energy quasiparticles in these surface states are governed by a Dirac equation, meaning the particles’ velocity is energy independent in low energy bands. The interplay of the Dirac particles with superconductivity can result in rich and intriguing physics such as particle to antiparticle conversion rel1 (); rel2 (). While the surface states themselves may not be superconducting, they may however adopt superconducting properties by proximity to a superconducting electrode exp_ti1 (); exp_ti2 (); exp_ti3 (); exp_ti4 (); exp_ti5 (); exp_ti6 (). Therefore, such a platform can provide the unique possibility to materialize the interplay of relativity and superconductivity in laboratory and possibly make use of its advantages in functional devices and technological ways phi0 ().

Motivated by the recent experimental progresses in fabricating TI-superconductor (S) heterostructures exp_ti1 (); exp_ti2 (); exp_ti3 (); exp_ti4 (); exp_ti5 (); exp_ti6 (), we here study one of the manifestations of interplay between the superconductivity and spin-momentum locking nature of Dirac fermions on the surface of a 3D TI: The coupling between spin and charge degrees of freedom in a two dimensional Josephson junction made of disordered surface states of a three dimensional TI under equilibrium and nonequilibrium conditions.

In the equilibrium state, we predict that the system responses to a dc supercurrent flowing across the junction by developing a stationary spin density oriented along the junction interfaces. This phenomenon can be viewed as a direct magnetoelectric effect i.e. the Edelstein effect. In the conventional intrinsic spin orbit coupled metals, analogous effect was first theoretically predicted in Refs. aronov89, ; edelstein90, and later observed experimentally in Refs. kato04, ; silov04, . It was shown that spin polarization can be produced in the spin orbit coupled metals by externally applied electric field. aronov89 (); edelstein90 (); kato04 (); silov04 (). The magnetoelectric polarizability was also discussed in the normal phase of topological insulators rev1 (); rev2 (); essin09 (). We note that similar magnetoelectric effect was also predicted for bulk superconductors and superconducting heterostructures by means of supercurrent due to the spin orbit interactions edelstein95 (); edelstein05 (); malshukov08 (). Here we develop a theory to this phenomenon in S-TI-S heterostructures and show that a topological insulator can support much more pronounced signatures than those of conventional spin orbit coupled metals.

We also demonstrate an inverse coupling between the charge and spin degrees of freedom in S-TI-S heterostructures. In particular, it is shown that the injection of quasiparticle current at the middle of junction in opposite directions yields spontaneous magnetizations with opposite orientations in each segment. The current-induced magnetization is directly coupled to the phase difference between the superconducting terminals. This phenomenon renders the minimum of junction free energy into a superconducting phase difference between the S terminals , which is generally unequal to or . The is proportional to the voltage difference between the injector electrode and superconducting terminals and inverse of Fermi velocity on the surface states. -junctions were also predicted in ferromagnetic Josephson junctions with spin orbit interaction krive04 (); braude07 (); reinoso08 (); buzdin08 (); kuzmanovski16 () and may be viewed as an inverse magnetoelectronic effect in such hybrids konschelle15 ().

Nearly all of the past theoretical works on the spin orbit coupled Josephson structures involves ferromagnetism or an external magnetic field as key ingredients to establish the magnetoelectric effect. In this paper, however, we demonstrate that an internal Zeeman-like term can be generated and electrically controlled in S-TI-S heterostructures (in the absence of any ferromagnetic elements or an external magnetic field) simply by means of quasiparticles injection. The Zeeman-like term causes an anomalous phase shift in the supercurrent that can be electrically controlled via the quasiparticle flow injection generated by a voltage difference.

So far, the electric control of the superconducting critical temperature bobkova11 (); bobkova15 (); ouassou16 (), the magnitude of the Josephson current heikkila00 (); bobkova12 (), switching between and states volkov95 (); wilhelm98 (); baselmans99 (); huang02 (); crosser08 (); bobkova10 (), and spontaneously accumulated spin currents ald15 () have been discussed in hybrid structures with spin orbit coupling. The electric control of the anomalous phase shift was experimentally realized in superconducting heterostructures on the basis of a quantum wirephi0 (). This kind of controll also was theoretically proposed in S/silicene/S heterostructures kuzmanovski16 (). Nonetheless, in these proposals, the anomalous phase shift occurs only in the presence of an externally applied magnetic field or ferromagnetic elements. Our findings offer the ability of integrating the superconducting topological insulator nanostructures into electronic circuits without the requirement to apply magnetic field or making use of magnetic elements.

The paper is organized as follows. In Sec. II, we summarize the theoretical framework used and basic assumptions made for studying a disordered topological insulator Josephson junction under nonequilibrium conditions. The direct magnetoelectric effect in S-TI-S heterostructures is presented in Sec. III. The electrically controllable inverse magnetoelectric effect is discussed in Sec. IV. We finally summarize concluding remarks in Sec. V.

Ii nonequilibrium Keldysh technique

Numerous systems in laboratory deal with nonequilibrium processes. To theoretically describe a system that experiences nonequilibrium situations, one powerful approach is the Keldysh technique within the Green function framework keldysh (); mahan (). In this section, we start with the Hamiltonian of surface states of a three dimensional topological insulator and generalize the Keldysh technique to a Josephson structure made of the disordered surface states of a three dimensional TI under nonequilibrium. The system under consideration is schematically shown in Fig. 1(a).

The Hamiltonian that describes the Rashba type surface states in the presence of an in-plane exchange field reads:

(1)
(2)

in which , is the Fermi velocity, is a unit vector normal to the surface of TI (see Fig. 1(a)), is the chemical potential, and is a vector of Pauli matrices in the spin space. Note that the Dirac type of Hamiltonian would only change the notations zyuzin (). We assume that the system involves nonmagnetic impurities that can be described by a Gaussian scattering potential:

(3)

Here is the mean free time of quasiparticles and is the density of states at the Fermi level of the normal state of TI.

Figure 1: (a) Schematic of the two dimensional S-TI-S hybrid junction. The surface states of a three dimensional topological insulator is sandwiched between two wave superconductors. The junction length along the direction is equal to while its width along is assumed enough large so that the influences of boundaries in this direction are negligible. The junction plane resides in the plane and S-TI interfaces are located at . (b) The spin-momentum locking phenomenon in the conduction band of surface channels.

The advanced (), retarded (), and Keldysh () blocks of Gor’kov Green function in the Keldysh technique are defined as follows:

(4a)
(4b)
(4c)
(4d)
(4e)

where and stand for anticommutator and commutator relations, respectively. Also, the time ordering operator is shown by a step function in time . The other components of Keldysh block i.e. , , and are related to this component Eq. (4e) the same way as those of the retarded and advanced components given above, namely Eqs. (4b)-(4d), to Eq. (4a).

We now introduce the following matrix Green function in the Nambu space:

(5)

where are matrices in the spin space. The elements of these matrices are given above by Eqs. (4a)-(4d). This definition can be extended to the advanced and keldysh blocks and . It is also convenient to introduce a full Green function in the Keldysh space as follows:

(6)

Therefore, by averaging the Green function over the impurity scattering potential in the Born approximation we find the following Gor’kov equation:

(7)

Here, the two time dependent products of operators is equivalent to .

We consider a situation where the chemical potential is the largest energy scale in the system and hence a quasiclassical approximation is the well suited framework to describe the system. We now can introduce the quasiclassical Green function

(8)

where a Fourier transformation of the Green function with respect to the relative space arguments and an integration over are performed. Furthermore, following the standard procedures zyuzin (), we find the Eilenberger equation for :

(9)

in which , , and , are Pauli matrices in the particle-hole space. This is a nonequilibrium generalization to the Eilenberger equation derived in Ref. zyuzin, for the equilibrium situations. Recently, the analogous quasiclassical equation was also discussed for Dirac edge and surface electrons with another type of spin-momentum locking hugdal16 (). In the stationary situations, that we consider throughout the paper, one can perform a Fourier transformation with respect to . Thus, the Eilenberger equation can be expressed by:

(10)

To find a general solution to Eq. (10), the Green function can be expanded by the Pauli matrices as follows;

(11)

where ,, , and depend on and are matrices in the particle-hole and Keldysh space. It can be shown that and are smaller than both and by a factor of order of . Therefore, one can safely neglect these terms in the expansion without missing significant information. The Eilenberger equation in this approximation decomposes into two coupled equations for and :

(12a)
(12b)

where . It is apparent that Eqs. (12a) and (12b) allow for a solution of . This means that the spin structure of full Green function is determined by . Physically, this operator projects the Green function onto the conduction band of TI surface states as schematically depicted in Fig. 1(b). We assume that this band can describe the surface electrons of TI and take the spinless Green function in the rest of our calculations. Now, equations (12a) and (12b) reduce to a single spinless equation so that satisfies an equation of the same form. Also, the normalization condition i.e. is a solution to this equation.

In the diffusive regime, where the quasiparticles’ motion is fully randomized by strong scattering resources so that , the Eilenberger equation reduces to the Usadel equationusadel (). To have a self-contained presentation, we here give a short recap of basic notations and equations generalized for a nonequilibrium situation zyuzin (). In the diffusive regime, the Green function can be expanded through the first two harmonics:

(13)

where the zeroth harmonic (isotropic) is larger than the first harmonic: . Following the derivation steps described in Ref. zyuzin, , one obtains the Usadel equation to :

(14)

where is the diffusion constant and the first harmonic term can be expressed in terms of as follows:

(15)

The operator is defined for the Rashba type band structure of surface states as follows:

(16)

Note that the formalism can be straightforwardly extended to the Dresselhaus type zyuzin (). In the junction configuration systems, the Usadel equation should be supplied by appropriate boundary conditions. Here, we consider low transparent tunnelling processes at the TI-S interfaces shown in Figs. 1(a) and 2. This boundary condition is experimentally relevant and was discussed in Refs. bc1, ; bc2, for conventional metallic junctions. We assume that the wave superconducting terminals with a gap of in their energy spectrums are in the equilibrium state and thus can be described by the bulk solution at the boundaries zyuzin ():

(17)

where is the ratio of resistance of the interface barrier per unit area to the resistivity of the TI surface states and is the unit vector normal to the interface (it points towards the superconductor region). In the tunneling low proximity regime we consider throughout our calculations, the opacity of interfaces should be sufficiently large to satisfy this condition. The superconducting bulk solution can be given by:

(18)

where , is an infinitesimal positive value and the system temperature is denoted by .

In the context of quantum transport, the electric current is one of the most important physical quantities that can explain transport experiments. The electric current density across the surface channels of TI junction we consider can be expressed by:

(19)

Rewriting the current density via the quasiclassical Green functions in the diffusive limit we obtain:

(20)

where is the normal state conductivity.

Iii magnetoelectric effect

In this section, we develop a theory of the dissipationless direct magnetoelectric effect to the disordered topological insulator Josephson junction depicted in Fig. 1(a). As it was mentioned earlier in the introduction, in response to a dc Josephson current flowing across the TI junction, a stationary spin density oriented along the junction interface can be generated. To uncover this phenomenon, we first evaluate the average spin polarization in TI which is:

(21)

In terms of the Green function, the components of spin polarization take the following form:

(22)

Within the quasiclassical diffusive regime we consider for the disordered topological insulator, the spin polarization components can be rewritten as

(23)

Substituting from Eq. (15) into Eq. (23), we obtain:

(24)

Comparing the last expression to Eq. (20), we see that

(25)

That is, the electric current and electron spin polarization orientation are perpendicular while their amplitudes are directly proportional.

When a superconducting phase gradient is applied between the superconducting terminals, a dc Josephson current flows across the junction. In the tuneling and low proximity limit, the supercurrent has a standard sinusoidal relation in which is the critical supercurrent zyuzin (). Therefore, considering the spin polarization relation given above, we see that the spin polarization can be simply controlled by the superconducting phase difference . This is a generic phenomenon and occurs in ballistic systems as well. To derive an expression for the spin density in the ballistic regime, we consider a surface channel with superconductivity where the superconducting phase experiences a coordinate gradient so that . The spin density in momentum representation and in the limit when the superconducting gap is much smaller than the temperature is given by

(26)

where . The Green function in the Matsubara representation reads:

(27)

Following similar steps as Ref. edelstein95, , one obtains the spin density in the real space representation as:

(28)

where Considering one finds that at the charge neutrality point at which the spin density reduces to:

(29)

while within regime one arrives at:

(30)

As seen, in both regimes, the spin polarization and the direction of the current (which is proportional to the phase gradient ) are orthogonal to each other.

We note that similar direct magnetoelectric effect was also predicted for superconductors and superconducting heterostructures with intrinsic spin orbit interactions edelstein95 (); edelstein05 (); malshukov08 (). In these systems, a Rashba spin orbit interaction is considered to be present and a generic relation to the spin polarization obtained . Nonetheless, the coefficient in realistic systems of this type is practically negligible and proportional to , where is the splitting of the conducting bands due to the spin orbit interaction and is the Fermi energy.

This issue however can be resolved in the surface channels of a 3D TI. The surface electrons are fully spin-momentum locked and therefore they have only one conduction band. Hence, the relation between the spin polarization and supercurrent contains no reducing coefficient of order of and the effect should be much stronger than that of the conventional spin orbit coupled materials with discussed so far in the literature.

Iv inverse magnetoelectric effect

In this section, we first show that a nonequilibrium quasiparticle injection on the surface of 3D TI induces a Zeeman-like field which is proportional to the voltage difference between an injector electrode and superconducting terminals that drives the quasiparticle flow. We then calculate the dc Josephson current through a disordered topological insulator junction where a flow of nonequilibrium quasiparticle current is injected into the TI region. We discuss the appearance of states which is controlable by means of the quasiparticle flow injection.

iv.1 Zeeman-like field induced by quasiparticle injection

In this subsection, we demonstrate that an electric current flowing through the surface states of the 3D TI can induce a Zeeman-like term. We consider a highly simple model which is sufficient to reveal this effect. In order to more simplify our discussions, we consider fully ballistic surface states without any disorder. In this case, the Hamiltonian is given by Eq. (1) where . Since we have employed the quasiclassical approximation where , we only consider the conduction band and assume in Eq. III that

(31)

Moreover, we include the electron-electron interaction at the Hartree-Fock level. Thus, the self-energy reads:

(32)

where

(33a)
(33b)

We also assume that contains the spin independent part of self-energy and focus on the spin part of . The Green function with interactions reads:

(34)

where and

(35)

By now substituting the Green function Eq. (34) into Eqs. (33a) and (33b) and assuming that we find that the spin dependent part of the self-energy is:

(36)

In equilibrium at , the spin dependent part of the self-energy vanishes due to the integration over momentum. However, in the presence of bias voltage the integral results in a nonvanishing finite value. If we assume an electric current in the direction of a ballistic system, we then are able to substitute in Eq. (36). In the linear approximation with respect to the dimensionless interaction constant together with we find:

(37)

where we set and expand Eq. (36) around to reach at a liner response theory. Hence, the only nonzero component is :

(38)

In conclusion, based on a simplified model, we have shown that an electric current can generate spin dependent exchange that results in an effective Zeeman-like term. This induced exchange field can be controlled by the applied voltage to inject quasiparticles. A more detailed study of the particular connection between the induced Zeeman-like term and the injected quasiparticles’ current is an interesting problem, however, is beyond the scope of the present paper and can be addressed elsewhere.

iv.2 Josephson current through a 3D TI under quasiparticle flow injection

The flow injection of quasiparticles into the junction generates a nonequilibrium situation in the TI region. The Josephson current in this nonequilibrium situation is determined not only by the condensate wave functions in the TI region but also by the nonequilibrium distribution function in this segment. The current can be calculated by Eq. (20). It is convenient to calculate it at the interfaces of TI-S by exploiting the boundary conditions Eq. (17):

(39)

where signs refer to the boundary conditions at . The Keldysh component of the full Green function can be expressed via the retarded and advanced components with the distribution function:

(40)

where the distribution function is a diagonal matrix in the particle-hole space:

(41)

There is a general relation (due to the particle-hole symmetry) between and . For our case of spinless fermions it reduces to .

Utilizing the distribution function, the current Eq. (39) can be expressed by:

(42)

Here, we assume that the left and right superconducting terminals possess macroscopic phases so that the phase difference between the superconductors is . The anomalous components of Green function at the interfaces are denoted by and . In order to find the Josephson current to the leading order with respect to the interface conductance we expand the Green function around the bulk solution in the TI region. To obtain Eq. (42) we have already taken into account this regime and set to the first order in .

Figure 2: Schematic of the topological insulator Josephson junction proposed to generate a -junction by means of quasiparticle flow injection. The junction resides in the plane (the same as Fig. 1(a)) and the injector electrode is located at while the superconducting interfaces are assumed at . The quasiparticle flow is controlled by a voltage difference between the superconducting terminals and injector electrode.

The experimental setup we propose with the electronic current injection is schematically shown in Fig. 2. The setup is made of two superconductor terminals connected via the surface channels of a 3D TI at . The junction plane resides in the plane and a normal electrode that injects quasiparticles into the junction is attached at . The quasiparticles flow can be induced by applying a voltage difference between the superconducting terminals and normal injector electrode. We assume that the voltage difference between the two superconducting terminals is zero so that we are able to eliminate the ambiguities and complications concerning nonequilibrium phenomena that arise when the two superconductors possess a finite voltage difference. Therefore, the quasiparticles current has opposite directions within each region of the TI surface separated by the injector electrode at . The direction of quasiparticles current is shown by white arrows in Fig. 2 and we assume that the current is uniform along the direction parallel to the S-TI interfaces. As discussed in passing, the electric current within TI in the direction generates a finite spin polarization along the direction. Therefore, as it was shown in the previous subsection, this electric current injection induces an effective Zeeman-like energy in the TI region where is the voltage difference applied between the normal injector electrode and the superconductors. The resulting exchange fields are shown in Fig. 2. The proposed setup generates two domains with opposite orientations for the effective Zeeman-like term which are locked to the direction of injected electric currents. This effective Zeeman-like domain scenario can be described by:

(43)

To find the anomalous Green function within the TI region, we solve the linearized Usadel equation in the surface channels:

(44a)
(44b)

where refer to the retarded and advanced components.

The Usadel equation should be solved together with the boundary conditions Eq. (17) at S-TI interfaces:

(45a)
(45b)

and at the interface between the domains :

(46a)
(46b)
(46c)
(46d)

where . The solution to component takes the following form:

(47)

in which and . To calculate the current across the junction, one needs also to obtain a solution to . The solution to this component however can be simply given as (, ). By substituting the obtained and into the current definition (42), we are able to analyze the transport properties of the system proposed. The current relation (42) can be decomposed to two components

(48)

where

(49)

and is expressed by the remaining parts of Eq. (42) that contains the anomalous Green functions. Note that because of opposite signs at the S-TI interfaces, the contribution of component to the Josephson current vanishes. Physically, this component represents a part of the injected quasiparticle current. Specifically, it can be seen that is equal to zero in the equilibrium simply because . Therefore, component contains the Josephson current and we focus on this component of the current in the following. Since the anomalous components of Green function are of the first order in , we thus only need the distribution function in the zeroth order with respect to this parameter for calculating . In this regime we have

(50a)
(50b)

These expressions imply that the main contribution to the resistance of the system originates from the tunnel process at the TI-S interfaces. The following combination of the distribution functions enables us to even decompose component of current to two parts and :

(51)
(52a)
(52b)

in which and at the right hand side of Eq. (52a) pertinent to the left and right interfaces. As seen, reverses its direction at the left and right S-TI interfaces. Consequently, this component is unable to contribute to the current between the superconductors. This is also a part of the injected electric current and does have no effect on the Josephson current. In contrast, the current component does not change its sign at the superconductor interfaces, and therefore represents the nonvanishing current flow through the junction. This term reveals an electrically controllable -junction that we discuss in detail in the next subsection.

iv.3 Controllable -Josephson junction

The Josephson current flowing through the junction is expressed by Eq. (52b). The ground state of the system corresponds to zero current and it is apparent that this condition is satisfied at a nonzero phase difference between the superconducting leads namely at: . According to the previous section, this anomalous phase shift is proportional to the voltage bias between the superconducting terminals and the additional normal injector electrode. Therefore, the ground state of the Josephson junction can be experimentally switched on and off by simply controlling the quasiparticle flow injection.

Figure 3: Critical supercurrent as a function of voltage difference between the normal injector electrode and superconducting terminals . In panel (a) we set the length of TI region equal to while in panel (b) .

The quasiparticle injection controls not only the anomalous phase shift but also modifies the critical current of the Josephson junction. To illustrate this fact, we plot the corresponding critical current in Fig. 3 vs the voltage difference for two different lengths of the TI region. The panel (a) exhibits the critical supercurrent vs the voltage difference where the junction thickness is while in panel (b) we set . Here and this quantity is a superconducting coherence if the diffusion constant in the superconducting terminals and TI surface are equal. The current is measured in units of . As seen, the nonequilibrium conditions, generated by the quasiparticles injection into the TI region, strongly alter the critical Josephson current flowing through the junction and can result in supercurrent reversals. We see that the critical current in a junction of thickness shows a single sign reversal while larger thicknesses can result in multiple sign reversals. Note that sign reversals of supercurrent due to the quasiparticles injection was also theoretically predicted in conventional diffusive metalsvolkov95 (); wilhelm98 () and observed in experiments baselmans99 (); huang02 (); crosser08 (). However, the novelty of our finding is a transition between and ground states rather than transition. Our proposal offers new venues to explore more about the superconductor-topological insulator heterostructures, identify their topological characteristics, and potentially utilize them for technological functionalities.

V Conclusions

Motivated by recent experiments on hybrid structures of topological insulator-supercondutor exp_ti1 (); exp_ti2 (); exp_ti3 (); exp_ti4 (); exp_ti5 (); exp_ti6 (), we have generalized the quasiclassical approach using the Keldysh technique to disordered TI-S heterostructures under nonequilibrium situations. We study the coupling between electric charge current and spin polarization on the surface of a three dimensional TI. Utilizing the generalized approach, we propose a setup to experimentally achieve a controllable -Josephson junction solely by means of a quasiparticle flow injection without resorting to an external magnetic field or ferromagnetic element. The quasiparticles should be injected into the TI surface from the middle of junction by applying a voltage difference between a normal injector electrode and the superconducting terminals. The quasiparticle currents in opposite directions generate a Zeeman-like domain with opposite orientations. We show that the Zeeman-like term is directly coupled to the Josephson phase difference and shifts the junction ground state into arbitrary values by means of and uncover the influences of on critical supercurrent. We also discuss direct magnetoelectric effect that appears in response to a dc Josephson current on the surface channels and show that the induced electron spin polarization is much more pronounced than that of a conventional spin orbit coupled material. Our work demonstrates the great potential of topological insulators in magnetoelectrics by superconducting hybrid configurations due to the strong spin-momentum locking property of TI surface channels.

Acknowledgements.
A.M.B and I.V.B were supported by Grant of the Russian Scientific Foundation No. 14-12-01290.

References

  1. M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
  2. X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
  3. O. A. Pankratov, S. V. Pakhomov, and B. A. Volkov, Solid State Comm. 61, 93 (1987).
  4. G. E. Volovik, The Universe in a Helium Droplet (Oxford University Press, Oxford, 2003).
  5. B. A. Bernevig, T. L. Hughes, and S-C. Zhang, Science 314, 1757 (2006).
  6. M. Knig, S. Wiedmann, C. Brne, A. Roth, H. Buhmann, L. W. Molenkamp, X-L. Qi, and S-C. Zhang, Science 318, 766 (2007).
  7. D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava and M. Z. Hasan, Nature 452, 970 (2008).
  8. Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal1, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava and M. Z. Hasan, Nat. Phys. 5, 398 (2009).
  9. H. Zhang, C-X. Liu, X-L. Qi, X. Dai, Z. Fang, and S-C. Zhang, Nat. Phys. 5, 438, (2009).
  10. A. A. Burkov and D. G. Hawthorn, Phys. Rev. Lett. 105, 066802 (2010).
  11. P. Schwab, R. Raimondi and C. Gorini, Europhys. Lett, 93, 67004 (2011).
  12. A. A. Zyuzin, M. Alidoust, and D. Loss, Phys. Rev. B 93, 214502 (2016).
  13. K. C. Nowack, E. M. Spanton, M. Baenninger, M. Konig, J. R. Kirtley, B. Kalisky, C. Ames, P. Leubner, C. Brune, H. Buhmann, L. W. Molenkamp, D. Goldhaber-Gordon, and K. A. Moler, Nat. Mat. 12, 787 (2013).
  14. E. M. Spanton, K. C. Nowack, L. Du, G. Sullivan, R.-R. Du, and K. A. Moler, Phys. Rev. Lett. 113, 026804 (2014)
  15. K. Capelle and E. K. U. Gross, Phys. Lett. A 261, 198 (1995).
  16. K. Capelle and E. K. U. Gross, Phys. Rev. B 59, 7140 (1999).
  17. I. Sochnikov, L. Maier, C. A. Watson, J. R. Kirtley, C. Gould, G. Tkachov, E. M. Hankiewicz, C. Brüne, H. Buhmann, Phys. Rev. Lett. 114, 066801 (2015).
  18. 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).
  19. M. Veldhorst, M. Snelder, M. Hoek, T. Gang, V. K. Guduru, X. L. Wang, U. Zeitler, W. G. van der Wiel, A. A. Golubov, H. Hilgenkamp and A. Brinkman, Nat. Mat. 11 417 (2012).
  20. J. B. Oostinga, L. Maier, P. Schüffelgen, D. Knott, C. Ames, C. Brune, G. Tkachov, H. Buhmann, and L. W. Molenkamp, Phys. Rev. X 3, 021007 (2013).
  21. S. Hart, H. Ren, T. Wagner, P. Leubner, M. Mühlbauer, C. Brüne, H. Buhmann, L. W. Molenkamp and A. Yacoby, Nat. Phys. 10, 638 (2014).
  22. S. Lee, X. Zhang, Y. Liang, S. Fackler, J. Yong, X. Wang, J. Paglione, R. L. Greene, I. Takeuchi, arXiv:1604.07455.
  23. D. B. Szombati, S. Nadj-Perge, D. Car, S. R. Plissard, E. P. A. M. Bakkers and L. P. Kouwenhoven, Nat. Phys. 2, 568 (2016).
  24. A. Aronov and Y. Lyanda-Geller, JETP Lett. 50, 431 (1989).
  25. V. Edelstein, Solid State Comm. 73, 233 (1990).
  26. Y.K. Kato, R.C. Myers, A.C. Gossard, and D.D. Awschalom, Phys. Rev. Lett. 93, 176601 (2004).
  27. A.Y. Silov, P.A. Blajnov, J.H. Wolter, R. Hey, K.H. Ploog, and N.S. Averkiev, Applied Phys. Lett. 85, 5929 (2004).
  28. A. M. Essin, J. E. Moore, and D. Vanderbilt, Phys. Rev. Lett. 102, 146805 (2009).
  29. V.M. Edelstein, Phys. Rev. Lett. 75, 2004 (1995).
  30. V.M. Edelstein, Phys. Rev. B 72, 172501 (2005).
  31. A.G. Mal’shukov and C.S. Chu, Phys. Rev. B 78, 104503 (2008).
  32. I.V. Krive, L.Y. Gorelik, R.I. Shekhter, and M. Jonson, Phys. Nizk. Temp. 30, 535 (2004).
  33. V. Braude and Yu.V. Nazarov, Phys. Rev. Lett. 98, 077003 (2007).
  34. A.A. Reynoso, G. Usaj, C.A. Balseiro, D. Feinberg, and M. Avignon, Phys. Rev. Lett. 101, 107001 (2008).
  35. A.I. Buzdin, Phys. Rev. Lett. 101, 107005 (2008).
  36. D. Kuzmanovski, J. Linder, A. Black-Schaffer, arXiv:1605.03197.
  37. F. Konschelle, I.V. Tokatly, and F.S. Bergeret, Phys. Rev. B 92, 125443 (2015).
  38. I.V. Bobkova and A.M. Bobkov, Phys. Rev. B 84, 140508 (2011).
  39. I.V. Bobkova and A.M. Bobkov, JETP Lett. 101, 407 (2015).
  40. J. A. Ouassou, A. Di Bernardo, J. W. A. Robinson, J. Linder, arXiv:1601.07176.
  41. T.T. Heikkil, F.K. Wilhelm, and G. Schn, Europhys. Lett. 51, 434 (2000).
  42. I.V. Bobkova and A.M. Bobkov, Phys. Rev. Lett. 108, 197002 (2012).
  43. A.F. Volkov, Phys. Rev. Lett. 74, 4730 (1995).
  44. F.K. Wilhelm, G. Schn, and A.D. Zaikin, Phys. Rev. Lett. 81, 1682 (1998).
  45. J.J.A. Baselmans, A.F. Morpurgo, B.J. van Wees, and T.M. Klapwijk, Nature (London) 397, 43 (1999).
  46. J. Huang, F. Pierre, T.T. Heikkil, F.K. Wilhelm, and N.O. Birge, Phys. Rev. B 66, 020507 (2002).
  47. M.S. Crosser, J. Huang, F. Pierre, P. Virtanen, T.T. Heikkil, F.K. Wilhelm, and N.O. Birge, Phys. Rev. B 77, 014528 (2008).
  48. I.V. Bobkova and A.M. Bobkov, Phys. Rev. B 82, 024515 (2010).
  49. M. Alidoust and K. Halterman, New J. Phys. 17, 033001 (2015).
  50. L.V. Keldysh, Sov. Phys. JETP 20, 1018 (1965).
  51. G.D. Mahan, in Many-Particle Physics, (Plenum Press, 1990).
  52. H.G. Hugdal, J. Linder, S.H. Jacobsen, arXiv:1606.01249.
  53. K.D. Usadel, Phys. Rev. Lett. 25, 507 (1977).
  54. A. V. Zaitsev, Sov. Phys. JETP 59, 1015 (1984).
  55. M. Y. Kuprianov and V. F. Lukichev, Sov. Phys. JETP 67, 1163 (1988).