Crossed Andreev reflection in zigzag phosphorene nanoribbon based ferromagnet/superconductor/ferromagnet junctions

We study the crossed Andreev reflection in zigzag phosphorene nanoribbon based ferromagnet/superconductor/ferromagnet junction. Only edge states, which are entirely detached from the bulk gap, involved in the transport processes. The perfect crossed Andreev reflection, with the maximal nonlocal conductance -2e2/h\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-2e^{2}/h$$\end{document}, is addressed by setting the electric potentials of the leads and device properly. At this situation, the local Andreev reflection and the electron tunneling are completely eliminated, the incoming electrons can only be reflected as electrons or transmitted as holes, corresponding to the electron reflection and the crossed Andreev reflection respectively.The nonlocal conductance oscillates periodically with the length and the electric potential of the superconductor. Our study shows that the phosphorene based junction can be used as the quantum device to generate entangled-electrons.

www.nature.com/scientificreports/ especially by the electric field 19 . Furthermore, the edge bands are entirely detached from the bulk bands, which helps them to eliminate the trouble caused by the bulk bands. Such magnificent properties of edge bands make it good for building the electronic device. Note that the Josephson effect in phosphorene based superconductor junctions were investigated 20 . Researcher found the supercurrent is highly anisotropic, due to the anisotropic band structure of the system. Furthermore, the supercurrent is also affected by the electric potential and the length of the superconductor region. However, compared to graphene, according to our investigation, the study of CAR in phosphorene based junction is still a few.In this Letter, we study the crossed Andreev reflection in zigzag phosphorene nanoribbon based ferromagnet/superconductor/-ferromagnet (F/S/F) junction. The ferromagnetism and the superconductivity of phosphorene nanoribbon can be achieved in proximity to EuO 22 and s-wave superconductor respectively. According to the first principle calculation, the exchange energy can reach 184 meV. Normally, the electron tunneling (ET) and local Andreev reflection (LAR) will compete with the CAR in the scattering processes inevitably 23 , making the nonlocal conductance weaken. In our study, the scattering processes in the junction can be consisted of CAR and electron reflection (ER), where the ET and LAR are completely blocked by setting the electric potentials of both leads and device properly. The maximal nonlocal conductance of CAR can reach −2e 2 /h , corresponding to a perfect CAR. Also, we found out that the nonlocal conductance is oscillating periodically with the length and electric potential of the superconductor. Our finding shows that the phosphorene based junction can be used as the quantum device to generate entangled-electrons.
This work is organized as follows. The structure and model of the junction are presented in "Structure and model" section; the calculation results and the corresponding analysis are presented in "Results and discussion" section. In "Conclusion" section, the conclusion of this work is given. Figure 1 is the structure we studied. The F/S/F junction is based on ZPR, with the superconductor connected to two ferromagnetic leads. The ferromagnet and superconductor can be achieved by the proximity effect. The primitive cell is shown in the regular frame with the size X × Y , where X = 3.27 Å and Y = 4.43 Å 18 are the length of the cell in x and y direction respectively. The size of the device region is L × T = N L X × N T Y , where N L and N T denote the number of unit cell in x and y direction respectively. For a convenient description, N L and N T are referred to as the length and the width of the device region respectively, the width of leads is the same as device region. Four atoms in a unit cell are labelled as A,B,C and D. A and B(C and D) belong to the bottom(upper) layer of the phosphorene nanoribbon. As we know, the effective Hamiltonian of the bare phosphorene nanoribbon can be written in a tight binding form where c A i σ ( c B i σ , c C i σ , c D i σ ) denotes annihilating an electron on site A i ( B i , C i , D i ) with spin σ =↑ / ↓ . There are 5 hopping parameters in H 0σ , shown in Fig. 1. The parameters read as t 1 = −1.220 eV , t 2 = 3.665 eV , t 3 = −0.205 eV , t 4 = −0.105 eV , and t 5 = 0.055 eV 24 . Actually, compare with the model in Ref. 24 , t 3 and t 5 are neglected, because they bring no important effect on the transport properties of phosphorene nanoribbon. i, j in H 0σ denotes that an electron or a hole hops from atom j in one layer to atom i in another layer,shown in Fig. 1.

Structure and model
Bogliubov-de Gennas (BdG) equation is used to study the CAR in ZPR superconductor junction. The tightbinding Hamiltonian of the junction can be written as www.nature.com/scientificreports/ H T stands for the couple between superconductor region and ferromagnetic leads, where c † S c F represents a particle hops from the lead to the device region, and F, S denotes the hopping happens between nearest atoms. H S /H F1(F2) denotes the Hamiltonian of superconductor/left(right) ferromagnet, is the Heaviside step function. in H S is the order parameter of the superconductor. ε ex in H F1(F2) represents the exchange energy in the ferromagnetic leads, and the value of ε ex is set to be 0.2 eV in this study. V 1(2) /V s is the electric potential of the left(right) lead/superconductor region, which can be tuned by a top gate voltage.
The calculation in this study is done by Kwant 25 , a software for quantum transport calculation. The tightbinding model of the system can be built by Kwant's package, and the structure can be visualized by setting system parameters. In detail, the lattice built by Kwant is the phosphorene lattice shown in Fig. 1. Note that the lattice of monolayer phosphorene described by C 2h point group. The leads belong to the same lattice and they are semi-infinite. After that, the scattering matrix [S] can be obtained. In our study, the expression of the wave function in the right lead is e S eh (E F )ψ e , where e denotes the injecting electrons in the left lead, h denotes the outgoing holes in the right lead and E F is the Fermi energy. S eh denotes the element of scattering matrix, and the transmission can be written as T eh = |S eh (E F )| 2 . For simplicity, we can simply just consider the CAR process form the left lead to the right lead 11,14,26,27 . According to the Ref. 26 , the nonlocal conductance can be given as where I R can be measured as the current formed in the right lead and V L is the bias applied between the left lead and the device. Kwant uses the efficient method of solving linear equations, making the calculation speed faster than the lattice Green function approach. The wavefunctions and the band structures of the leads can be obtained. For more details please refer to Ref. 25 .

Results and discussion
The band structures of the ferromagnetic lead had been shown in Fig. 2. The electron and hole bands are drawn together, which are plot in solid and dot lines respectively. Note that, there exist four electron (hole) edge bands in the band gap, considering the spin degeneracy. The electron edge bands are separated into two double degenerate bands by the exchange field, which can be found in Fig. 2a. As ε ex is set as 0.2 eV, the electrons in upper/bottom electron edge bands are spin-up/spin-down electrons. Here, we set V 1 < 0 and V 2 > 0 to get the band structure of left and right lead in Fig. 2b respectively. According to the Hamiltonian of the ferromagnetic lead, when the extra electric potential is set to be V < 0 , the electron bands will be pushed down and the hole bands will be uplifted. By doing so, the LAR and ET can be eliminated in the scattering processes. This can be understood from Fig. 2b, at the Fermi energy E F , there are no conducting modes to reflect holes and transmit electrons; Furthermore, suppose there is no spin flipping in the processes, the bands of spin-up hole and the spin-down electron will not contribute to the scattering processes when the spin-up electron is injecting into the system from left lead. Therefore, only two situations happen when electrons injecting into the system: reflect electrons or transmit holes for forming ER or CAR respectively.
The calculation results are shown in Fig. 3. The number of unit cell of the device region are N L = 40 , N T = 40 . Without losing generality, the electric potential of left (right) lead is set to be V 1 = −0.4 eV ( V 2 = 0.6 eV ), so that there is only spin-up electron injecting into the system around E = 0 eV . The electric potential and the superconductor gap of the superconductor are set to be V s = −164 meV and = 1.5 meV respectively. Note that, the scattering coefficients T eh , R ee , T ee and R eh denote CAR, ER, ET and LAR process respectively. Agree with our analysis above, the scattering coefficients of LAR and ET are zero, i.e., R eh = 0 and T ee = 0 , for there is no conducting modes for these two processes around E = 0 eV . Since there are two electrons injecting into the junction simultaneously and only CAR and ER existed in the scattering processes at the energy range we studied, the relation T eh + R ee = 2 is obtained. The scattering coefficients T eh and R ee in Fig. 3 are asymmetric with respect to the E coordinate. Two peaks of T eh appear at E = −4 and E = 4 in Fig. 3, with the values 0.75 and 2 respectively. Note that, the maximal T eh can achieve 2 at E = 4 , corresponding to the nonlocal conductance −2e 2 /h . Such phenomenon is referred to as the perfect CAR, where the incoming electrons from one lead are entirely converting to outgoing holes in another lead. Figure 4 shows the scattering coefficients of CAR (and ER) as a function of the length N L and the electric potential V s of superconductor. By the way, as the change of N T of the junction will barely affect the band gap and T eh , we defaultly fixed N T = 40 in this study. The calculation results are shown in Fig. 4. In Fig. 4a, T eh as a www.nature.com/scientificreports/ function of N L , with V s and E are fixed at −149 meV and 2.16 respectively. As we see in Fig. 4a, the T eh and R ee are oscillating periodically with N L , with the period N L = 5 . The peak's value of T eh is changing with N L . We also found that there exists a peak at N L = 45 , with T eh = 2 , in Fig. 4a. Then, when we change V s and fix N L = 40 and E = 2.16 , a similar effect happens. The T eh oscillates with a period V s = 15 meV , which means that the CAR is quite sensitive to the change of electric potential of the superconductor. We also found T eh = 2 around V s = −164 meV . The perfect CAR in Fig. 4 is attributed to the resonance effect of the quantum transport, the resonance condition is kL = nπ , k and L are the wave vector and the length of the device region respectively. The mechanism of such periodic osocillation had reported in Ref. 28 . These properties make the CAR in the junction controllable.  www.nature.com/scientificreports/ The nonlocal conductance are presented at Fig. 5a and d. In Fig. 5a, when N L increases from 45 to 47, the CAR peak moves from the higher energy to a lower one, with the N L changing. Such phenomenon also happens in changing V s , which is shown in Fig. 5d. These properties can be easily found in Fig. 5b and e. When the energy is in the range [−7.5 , 7.5 ] , there is a conductance "gap" at N L = 43/48/53 and V s around −138/ − 153/ − 168 meV , the corresponding T eh are suppressed and close to zero. Considering Fig. 5b/e and c/f, the T eh and R ee are complementary, with T eh + R ee = 2.

Conclusion
In this work, we studied the crossed Andreev reflection in the ferromagnet/superconductor/ferromagnet junction based on zigzag phosphorene nanoribbon. By setting the electric potential of both leads differently, the scattering processes consists of electron reflection and crossed Andreev reflection, where the local Andreev reflection and electron tunneling are suppressed. The maximal nonlocal conductance can reach −2e 2 /h , which is referred to as the perfect CAR. We also found that the transmission coefficient of crossed Andreev reflection oscillates periodically with the length and the electric potential of the superconductor. These properties make the CAR controllable. Our finding will be helpful to devise the quantum device for generating entangled-electrons.