# Novel Majorana mode and magnetoresistance in ferromagnetic superconducting topological insulator

## Abstract

Among the potential applications of topological insulators, we investigate theoretically the effect of coexistence of proximity-induced ferromagnetism and superconductivity on the surface states of 3-dimensional topological insulator, where the superconducting electron-hole excitations can be significantly affected by the magnetization of ferromagnetic order. We find that, Majorana mode energy, as a verified feature of TI F/S structure, along the interface sensitively depends on the magnitude of magnetization in FS region, while its slope in perpendicular incidence presents steep and no change. Since the superconducting gap is renormalized by a factor , hence Andreev reflection is more or less suppressed, and, in particular, resulting subgap tunneling conductance is more sensitive to the magnitude of magnetizations in FS and F regions. Furthermore, an interesting scenario happens at the antiparallel configuration of magnetizations and resulting in magnetoresistance in N/F/FS junction, which can be controlled and decreased by tuning the magnetization magnitude in FS region.

PACS: 74.45.+c; 85.75.-d; 73.20.-r

Keywords: topological insulator; ferromagnetic superconductivity; Andreev reflection; Majorana mode; tunneling conductance

## 1 Introduction

Topological insulators (TIs) represent new type of material which has emerged in the last few years as one of the most actively research subjects in condensed matter physics. They are characterized by a full insulating gap in the bulk and gapless edge or surface states, which are protected by the time-reversal symmetry [1, 2, 3, 4]. Regarding Bernevig and Hughes prediction [3, 5], TIs have been experimentally observed with such properties that host bound states on their surface, e.g. in 3-dimensional topological insulators (3DTI) , , and alloy, and also in the quantum well heterostructure [6, 7, 8]. These states form a band-gap closing Dirac cone on each surface, and lead to a conducting state with properties unlike any other known electronic systems. In particular, conformity of the conduction and valence bands to each other in and around Dirac points in the first Brillouin zone, possessing an odd number of Dirac points, description of fermionic excitations as massless two-dimensional chiral Dirac fermions, depending chirality on the spin of electron, having the significant electron-phonon scattering on the surface, owning very low room-temperature electron mobility are the peculiar properties of electronic structure of TIs. Interestingly, the charge carriers in the surface states can behave as massive Dirac fermions [9] due to its proximity to a ferromagnetic material, that the vertical component of the magnetic vector potential may be proportional to the effective mass of Dirac fermion. The experimentally observed proximity-induced superconductivity on the surface state is another interesting dynamical feature occuring in 3DTI, see Refs. [10, 11, 12].

More importantly, the coexistence of superconductivity and ferromagnetism as one of potential interests for spintronics and high magnetic field applications has firstly been predicted by Fulde and Ferrel [13], and Larkin and Ovchinnikov [14] as FFLO state. This effect can be in compliance with standard BCS theory for phonon-mediated -wave superconductivity, because the ferromagnetic exchange field is expected to prevent spin-singlet Cooper pairing, (see, Ref. [15] as a prior work). The magnetic polarization of a pair electron caused by a ferromagnetic material can lead to the different momentum of Cooper pair occurring in a ferromagnetic superconducting (FS) segment. It seems to be in contrast to the formation of a typical cooper pair, where two electrons may be in opposite spin direction with the same momentum. However, Bergeret et.al. [16] and Li et.al. [17] have studied the effect of superconductor/ferromagnetic bilayer on the critical Josephson current, where the orientation of ferromagnetic exchange field strongly affects the critical current. Also, the effect of superconductivity in coexistence with ferromagnetism has been studied on the superconducting gap equation for two case of singlet -wave and triplet -wave symmetries [18]. The authors have reconsidered the Clogston-Chandrasekhar limiting [19, 20]. According to the Clogston criterion in the conventional FS mixture, the normal state is regained as soon as the ferromagnetic exchange field exceeds at zero temperature. To be empirically, the [21] has been discovered to be the first ferromagnetic superconductor, which superconductivity is found to occur in a small temperature interval with adjusted ferromagnetic phase. Also, superconductivity is detected in itinerant ferromagnetic in a limited range of pressure and temperature [22].

Regarding several works in the recent few years concerning with the topological insulator-based junctions [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36], which are related to the Andreev process and resulting subgap conductance we proceed, in this paper, to theoretically study the dynamical properties of Dirac-like charge carriers in the surface states of 3DTI under influence of both superconducting and ferromagnetic orders via the introducing the proper form of corresponding Dirac spinors, which are principally distinct from those given in Ref. [29]. The magnetization induction opens a gap at the Dirac point (no inducing any finite center of mass momentum to the Cooper pair), whereas the superconducting correlations causes an energy gap at the Fermi level in the 3DTI. It will be particularly interesting to investigate the topological insulator superconducting electron-hole excitations in the presence of a exchange field. We assume that the Fermi level is close to the Dirac point, and the ferromagnet has a magnetization . The chirality conservation of charge carriers on the surface states in the presence of magnetization (due to opening the band gap) allows to use a finite magnitude of . In the absence of topological insulator, the spin-splitting caused by magnetization gives rise to limiting the magnitude of in a FS structure. These excitations, therefore, are found to play a crucial role in Andreev reflection (AR) process leading to the tunneling conductance below the renormalized superconducting gap. Particularly, we pay attention to the formation of Majorana bound energy mode, as an interesting feature in topological insulator ferromagnet/superconductor interface, depending on the magnetization of FS hybrid structure. We present, in section 2, the explicit signature of magnetization in low-energy effective Dirac-Bogogliubov-de Gennes (DBdG) Hamiltonian. The electron(hole) quasiparticle dispersion energy is analytically calculated, which seems to exhibit qualitatively distinct behavior in hole excitations () by varying the magnitude of magnetization. By considering the magnetization is ever less than chemical potential in FS region, the superconducting wavevector and corresponding eigenstates are derived analytically. Section 3 is devoted to unveil the above key point of FS energy excitation, Majorana mode energy, Andreev process and resulting tunneling conductance in N/F/FS junction and respective discussions. In the last section, the main characteristics of proposed structure are summarized.

## 2 Theoretical Formalism

### 2.1 Topological insulator FS effective Hamiltonian

In order to investigate how both superconductivity and ferromagnetism induction to the surface state affects the electron-hole excitations in a 3DTI hybrid structure, we consider magnetization contribution to the DBdG equation. Let us focus first on the Hubbard model Hamiltonian [37] that is included the effective exchange field follows from:

(1) |

where denotes the effective attractive interaction between arbitrary electrons, labeled by the integer and with spins and . The matrices are responsible for the hopping between different neighboring sites, and and indicate the second quantized fermion and number operators, respectively. Here, is the vector of Pauli matrix. Using the Hartree-Fock-Gorkov approximation and Bogoliubov-Valatin transformation [38], the Bogoliubov-de Gennes Hamiltonian describing dynamics of Bogoliubov quasiparticles is found. In Nambu basis, that electron(hole) state is given by , the BdG Hamiltonian for a -wave spin singlet superconducting gap in the presence of an exchange splitting can be written as:

(2) |

where denotes the non-superconducting Schrodinger-type part, and is superconducting order parameter. In the simplest model, can be chosen to be real to describe time-reversed states. The effective exchange field by rotating our spin reference frame can be gain as . The four corresponding levels of a singlet superconductor in a spin magnetic field is obtained with , where is the normal state energy for . However, dependence of superconducting order parameter on the exchange energy can be exactly derived from self-consistency condition [18]:

(3) |

where is the conventional order parameter in absence of ferromagnetic effect, and are the Boltzmann constant and temperature, respectively. The exchange splitting dependence of superconducting gap indicates that equation (3) has no functionality of at zero temperature. This takes place under an important condition known as Clogston-Chandrasekhar limiting [19, 20]. According to this condition, if the exchange splitting becomes greater than a critical value , then the normal state has a lower energy than the superconducting state. This means that a phase transition from the superconducting to normal states is possible when the exchange splitting is increased at zero temperature.

We now proceed to treat such a ferromagnetic superconductivity coexistence at the Dirac point of a 3DTI. It should be stressed that the dressed Dirac fermions with an exchange field in topologically conserved surface state have to be in superconducting state. Here, the influence of exchange field interacts in a fundamentally different way comparing to the conventional topologically trivial system, where the exchange field splits the energy bands of the majority and minority spins. A strong TI is a material that the conducting surface states at an odd number of Dirac points in the Brillouin zone close the insulating bulk gap unless time-reversal symmetry is broken. Candidate Dirac-type materials include the semiconducting alloy , as well as HgTe and under uniaxial strain [39]. In the simplest case, there is a single Dirac point in the surface Fermi circle and general effective Hamiltonian is modeled as , where indicates the surface Fermi velocity, and is the chemical potential. Under the influence of a ferromagnetic proximity effect, the Hamiltonian for the two-dimensional surface states of a 3DTI reads as:

where the ferromagnetic contribution corresponds to an exchange field . It has been shown [9] that transverse components of the magnetization on the surface are responsible to shift the position of the Fermi surface of band dispersion, while its perpendicular component to the surface induces an energy gap between conduction and valence bands.

In what follows, we will employ the relativistic generalization of BdG Hamiltonian, which is interacted by the effective exchange field to obtain the dispersion relation of FS dressed Dirac electrons in a topological insulator:

(4) |

The superconducting order parameter now depends on both spin and momentum symmetry of the Cooper pair, that the gap matrix for spin-singlet can be given as , where is the uniform amplitude of the superconducting gap and phase guarantees the globally broken symmetry. By diagonalizing this Hamiltonian we arrive at an energy-momentum quartic equation. Without lose of essential physics, we suppose the component of magnetization vector along the transport direction to be zero for simplicity. Also, we set , since the analytical calculations become unwieldy otherwise. The dispersion relation resulted from Eq. (4) for electron-hole excitations is found to be of the form:

(5) |

where, the parameter denotes the electron-like and hole-like excitations, while distinguishes the conduction and valence bands. We might expect several anomalous properties from the above superconducting excitations, which is investigated in detail in the next section. Equation (5) is clearly reduced to the standard eigenvalues for superconductor topological insulator in the absence of exchange field as (see Ref. [9]), . The mean-field conditions are satisfied as long as . In this condition, the exact form of superconducting wavevector of charge carriers can be acquired from the eigenstates .

The Hamiltonian Eq. (4) can be solved to obtain the electron (hole) eigenstates for FS topological insulator. The wavefunctions including a contribution of both electron-like and hole-like quasiparticles are analytically found as:

(6) |

where we define

Note that, the solution is allowed as long as the Zeeman field being lower than chemical potential .

### 2.2 FS interplay at the TI interface

We consider Andreev reflection in a hybrid N/F/FS structure formed on the surface of a 3DTI which coexistence between ferromagnet and superconductor is assumed to be induced by means of the proximity effect. The wide topological insulator junction is taken along the -axis with the FS region for , F region for and N region for . The superconducting order parameter vanishes identically in N and F regions, and we can neglect its spatial variation in the FS region close to the interface. The magnetization vectors of both sections is taken, in general, , which can be at the parallel or antiparallel configuration, as shown in Fig. 1. In the scattering process follows from the Blonder-Tinkham-Klapwijk (BTK) formula [40], we find the reflection amplitudes from the boundary condition at the interface. In ferromagnetic case, right- and left- moving electrons (holes) with energy excitation below the superconducting gap, transmitted (normal reflected) from the N region and reflected (Andreev reflected) at the FS interface. Thus, the leftover charge is transferred into the FS region as a Cooper pair at Fermi level. At energy excitation above the normalized superconducting gap resulted from Eq. (5) (see, in particular, Fig. 2) quasiparticle states can directly tunnel into the superconducting section. The reflected hole leading to AR can be actually controlled by the doping level in order to take place possible specular Andreev reflection. Particularly, we have to determine (via the dynamical features of system) the allowed values of Fermi energy in three regions. We set the Fermi energy to zero in F region. The electron(or hole) transmitted to the FS region angle may be accordingly obtained from the fact of conservation of transverse wavevector under quasiparticle scattering at the interface:

(7) |

where and are the chemical potential and incidence angle in N region, respectively. As an important point, the electron(hole) angle of incidence in all regions may be span the range from to around the normal axis. Regarding the Eqs. (7), the angle needs to be meaningful when the chemical potential of FS region takes a magnitude greater than its value in N region (). On the other hand, we previously applied the condition , as an experimentally used manner to calculate the wave functions Eq. (6).

By introducing the normal and Andreev reflection coefficients and the scattering coefficients in F region, the total wave function inside the N and F region can be written as:

(8) |

where the eigenvectors can be found in Appendix A. The probability amplitude of reflections in Eq. (8) are calculated from the continuity of the wavefunctions at the interface. The wave function in FS region is defined as . Finally, we find the following analytical expressions for the reflection coefficients, that the auxiliary quantities is described in Appendix A:

(9) |

The reflection amplitudes measurements under the BTK formalism enables us to capture the tunneling conductance through the junction:

(10) |

where the critical angle of incidence is determined depending on the doping of F region. The quantity is a renormalization factor corresponding to the ballistic conductance of normal metallic junction.

## 3 Results and Discussion

### 3.1 Energy excitation and Majorana mode

In this section, we proceed to analyze in detail the dynamical features of Dirac-like charge carriers in 3DTI with ferromagnetic and superconducting orders deposited on top of it. We assume that the Fermi level controlled by the chemical potential is close to the Dirac point. In this case, it is expected the signature of to be significant. In Fig. 2, we demonstrate the FS 3DTI electron-hole excitations. A net superconducting gap is obtained in Dirac points (for , where is Fermi wavevector) when we set . Increasing up to its possible maximum value results in three outcomes: i) the superconducting excitations, which is renormalized by a factor , disappear in hole branch (). It means that for the greater magnetizations, if we consent the superconductivity in FS 3DTI still exists, there is almost vanishing quantum state for reflected hole by Andreev process in the valence band, ii) Dirac point is shifted towards smaller FS quasiparticle electron-hole wavevectors, iii) the superconducting gap decreases slowly, where the variation of net gap is very low . The Andreev process, therefore, is believed to inconsiderably supress. The signature of these valence band excitations can be clearly shown in AR, where the Majorana mode may also be formed at the 3DTI F/FS interface [23, 25].

As a verified result, considering the topological insulator interface between the ferromagnetic insulator and conventional superconductor leads to the appearance of the chiral Majorana mode as an Andreev bound state. In other words, the Majorana mode and Andreev reflection are strongly related to each other. The latter can be realized by the fact of looking for bound energies produced by the perfect AR, which yields the following solution:

where we define

We have checked numerically that sign of is changed by . Thus, the sign of Andreev resonance states may be changed by reversing the direction of , and it corresponds to the chirality of Majorana mode energies. As shown in Fig. 3, the slope of the energy curves of around become steep and show no change with the increase of for fixed , while it exhibits significantly decreasing behavior with the increase of for fixed . The dispersion of Majorana modes along the interface () decreases with the increase of both magnetizations of FS and F regions. Note that, due to the presence of it needs to consider the Fermi level mismatch between normal and FS sections, i.e. . Then, the above contributions can be considerable in Andreev process and resulting subgap tunneling conductance.

### 3.2 Tunneling conductance

From the angle-resolved Andreev and normal reflection probabilities using Eq. (9), we see from Fig. 4(a) the main contribution of AR belongs to the angle of incidence in zero bias. It, therefore, is expected to achieve the lower zero bias conductance, as shown in Fig. 4(b) and (c). Furthermore, varying has no significant influence on AR in zero bias owing to the very small decrease of the renormalized superconducting gap with the increase of the , while the increasing results in more suppression of AR. The latter can be understood by the increase of band gap in Dirac point in F region. The resulting normalized angle-averaged tunneling conductance curves are reported in Figs. 4(b) and (c) for two parallel and antiparallel configurations of magnetizations in FS and F regions. Zero bias conductance peak disappears with the decrease of the , and instead of it a high conductance peak appears in bias . This result should be compared to that is obtained in Ref. [9, 24]. Interestingly, by increasing the the magnitude of subgap bias , for which the new peak takes place, is limited, as seen from Fig. 4(c). Thus, parameter can be considered as a “bias-limitation coefficient”. These features have been obtained when the direction of magnetizations in F and FS regions are at the parallel configuration. The fundamentally distinct scenarios we find for the case of antiparallel configuration of magnetizations. In this case, first, the tunneling subgap conductance is enhanced, secondly, the zero bias conductance peak presented in parallel case is replaced by a deep, see, in detail, Fig. 4(d). Dynamically description, when the direction of is inverted we, indeed, meet with an inverse energy gap in Dirac point of 3DTI giving rise to enhancing the conductance peak respective bias energy in low values of . One can express that the zero bias conductance originates from the chiral Majorana mode, which significantly depends on the . The chirality actually corresponds to the sign of , while the magnitudes of zero bias conductance at the parallel and antiparallel configurations are the same. Hence, the both deep and peak of conductance curves in antiparallel case are significance being influenced by the inverted gap caused by the .

Remarkably, the importance of above findings can be featured by the capture of magnetoresistance (MR) of the topological insulator junction. The magnetization (specialy in FS region) dependence of MR is presented in Fig. 5, where we observe a considerable MR peak for extra values of (e.g. in figure). Importantly, increasing the weakens the MR peak, since, regarding the superconducting excitations in Fig. 2, the AR is more or less suppressed in the presence of and Fermi wavevector mismatch also causes to decrease the -bias conductance peak at the antiparallel configuration. According to the conductance curves, there is no MR in zero bias.

## 4 Conclusion

In summary, we have investigated the influence of ferromagnetic superconducting orders coexistence in the surface state of topological insulator. The topological insulator superconducting electron-hole excitations in the presence of magnetization have led to achieve qualitatively distinct transport properties in tunneling N/F/FS junction. One of key findings of the present work is that the resulting subgap conductance has been found to be strongly sensitive to the parallel or antiparallel configuration of magnetization directions in FS and F regions. Thus, this feature has actually led to present the magnetoresistance peak for bias energy close to the renormalized superconducting gap , which the bias limitation coefficient includes the magnetization of FS region . Particularly, we have found the presence of Majorana mode at the F/FS interface to be controlled by the tuning of magnetizations magnitude. However, these results have been obtained in the case of and , which is relevant to the experimental regime.

## APPENDIX A: Normal and Andreev reflection amplitudes

To complete calculation of probability of reflections in N/F/FS junction, we write down right and left moving electron and hole spinors in F and N region:

where we define . By matching boundary conditions on and at and and at , the reflection amplitudes are obtained. We introduce auxiliary quantities in Eq. (9) as:

with

h.goudarzi@urmia.ac.ir ; goudarzia@phys.msu.ru

m.khezerlou@urmia.ac.ir

Figure captions

Figure 1 (color online) Sketch of the topological insulator-based N/F/FS junction. The magnetization vectors in F and FS regions can be at the parallel or antiparallel configuration.

Figure 2 (color online) The ferromagnetic superconducting excitation spectra on the surface state of 3DTI for several values of , calculated from Eq. (5). We set the net value of superconducting gap (this value of pair potential is taken only to more clarify the behavior of spectra in Dirac point, although it does not further need to use it in our calculations, since is supposed.

Figure 3 (color online) The dispersion of Majorana modes as a function of the electron incident angle for several values of magnetizations in FS and F regions. The solid lines correspond to and the dashed lines to .

Figure 4(a), (b), (c), (d) (color online) (a) Probability of the normal and Andreev reflections as a function of electron incidence angle at the interface in zero bias with and . The plots show the results for different values of . (b) Normalized tunneling conductance versus bias voltage and magnetization of F region. We set (c) Normalized tunneling conductance versus bias voltage and magnetization of FS region. We set (d) The tunneling conductance as a function of bias voltage for two signs of , corresponding to the parallel and antiparallel configurations in F and FS regions. The solid lines correspond to and marker dashed lines correspond to . We set .

Figure 5 (color online) The magnetoresistance spectra as function of bias voltage, where the influence of and is indicated, separately. We have set in the resulting conductance and magnetoresistance spectra.

### References

- C.L. Kane, E.J. Mele, Phys. Rev. Lett. 95 (2005) 146802.
- B.A. Bernevig, S.C. Zhang,Phys. Rev. Lett. 96 (2006) 106802.
- M.Z. Hasan, C.L. Kane, Rev. Mod. Phys. 82 (2010) 3045.
- X.L. Qi, S.C. Zhang, Rev. Mod. Phys. 83 (2011) 1057.
- B.A. Bernevig, T.L. Hughes, S.C. Zhang, Science 314 (2006) 1757.
- Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y.S. Hor, R.J. Cava, M.Z. Hasan, Nature Physics 5 (2009) 398.
- Y.L. Chen, J.G. Analytis, J.H. Chu, Z.K. Liu, S.K. Mo, X.L. Qi, H.J. Zhang, D.H. Lu, X. Dai, Z. Fang, S.C. Zhang, I.R. Fisher, Z. Hussain, Z.X. Shen, Science 325 (2009) 178.
- M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, L.W. Molenkamp, X.L. Qi, S.C. Zhang Science 318 (2007) 766.
- J. Linder, Y. Tanaka, T. Yokoyama, A. Sudbo, N. Nagaosa. Phys. Rev. B 81 (2010) 184525.
- J. Wang, C.Z. Chang, H. Li, K. He, D. Zhang, M. Singh, X.C. Ma, N. Samarth, M. Xie, Q.K. Xue, M.H.W. Chan, Phys. Rev. B 85 (2012) 045415.
- G. Koren, T. Kirzhner, E. Lahoud, K.B. Chashka, A. Kanigel, Phys. Rev. B 84 (2011) 224521.
- P. Zareapour, A. Hayat, S.Y.F. Zhao, M. Kreshchuk, A. Jain, D.C. Kwok, N. Lee, S.W. Cheong, Zh. Xu, A. Yang, G.D. Gu, Sh. Jia, R.J. Cava, K.S. Burch, Nature communications 3 (2012) 1056.
- P. Fulde, R.A. Ferrell, Phys. Rev. 135 (1964) A550.
- A.I. Larkin, Y.N. Ovchinnikov, Sov. Phys. JETP. 20 (1965) 762.
- N.F. Berk, J.R. Schrieffer, Phys. Rev. Lett. 17 (1966) 433.
- F.S. Bergeret, A.F. Volkov, K.B. Efetov Phys. Rev. Lett. 86 (2001) 3140.
- X. Li, Zh. Zheng, D.Y. Xing, G. Sun, Zh. Dong, Phys. Rev. B 65 (2002) 134507.
- B.J. Powell, J.F. Annett, B.L. Gyorffy, J. Phys. A: Math. Gen. 36 (2003) 9289.
- A.M. Clogston, Phys. Rev. Lett. 9 (1962) 266.
- B.S. Chandrasekhar, Appl. Phys. Lett. 1 (1962) 7.
- S.K. Sinha, et al., Superconductivity in Magnetic and Exotic Materials, Springer Berlin Heidelberg, (1984).
- S.S. Saxena, P. Agarwal, K. Ahilan, F.M. Grosche, R.K.W. Haselwimmer, M.J. Steiner, E. Pugh, I.R. Walker, S.R. Julian, P. Monthoux, G.G. Lonzarich, A. Huxley, I. Sheikin, D. Braithwaite, J. Flouquet, Nature 406 (2000) 587.
- L. Fu, C.L. Kane, Phys. Rev. Lett. 100 (2008) 096407.
- J. Linder, Y. Tanaka, T. Yokoyama, A. Sudbo, N. Nagaosa, Phys. Rev. Lett. 104 (2010) 067001.
- Y. Tanaka, T. Yokoyama, N. Nagaosa, Phys. Rev. Lett. 103 (2009) 107002.
- K.T. Law, P.A. Lee, T.K. Ng, Phys. Rev. Lett. 103 (2009) 237001.
- M. Snelder, M. Veldhorst, A.A. Golubov, A. Brinkman, Phys. Rev. B 87 (2013) 104507.
- J. Nussbaum, T.L. Schmidt, Ch. Bruder, R.P. Tiwari, Phys. Rev. B 90 (2014) 045413.
- P. Burset, B. Lu, G. Tkachov, Y. Tanaka, E.M. Hankiewicz, B. Trauzettel, Phys. Rev. B 92 (2015) 205424.
- T. Choudhari, N. Deo, Physica E 85 (2017) 238-247.
- G. Gupta, H. Lin, A. Bansil, M. Bin Abdul Jalil, G. Liang, Physica E 74 (2015) 10-19.
- M. Snelder, A.A. Golubov, Y. Asano, A. Brinkman, Journal of Physics: Condensed Matter 27 (2015) 315701.
- R. Vali, H.F. Khouzestani, Physica E 68 (2015) 107-111.
- M. Khezerlou, H. Goudarzi, Physica C 508 (2015) 6.
- T. Yokoyama, Sh. Murakami, Physica E 55 (2014) 1-8.
- H. Goudarzi, M. Khezerlou, J. Alilou, J. Super. and Novel Mag. 26 (2013) 3355.
- J. Hubbard, Proc. Roy. Soc. (London) A 276 (1963) 238.
- P.G. de Gennes, Superconductivity of Metals and Alloys, W. A. Benjamin, New York, (1966).
- L. Fu, C.L. Kane, Phys. Rev. B 76 (2007) 045302.
- G.E. Blonder, M. Tinkham, T.M. Klapwijk, Phys. Rev. B 25 (1982) 4515.