Abstract
Through the analytic solutions of the Bogoliubov de Gennes (BdG) equations the effect of a static and homogeneous magnetic field applied parallel to the interface of an NIS (N: Normal metal, S: superconductor and I: Insulator) junction on the differential conductance is calculated. For a d xy - symmetry we obtain zero bias conductance peak that can be split by a magnetic field. The shift of the zero bias conductance peak depends on the spread (beta) of the tunneling electrons in k space, on the magnitude of the applied field H and on the ratio between the Fermi energy of the superconductor and the normal region, E FS/E FN. Finally we estimate the minimum value of the magnetic field, Hmin, that splits the zero bias conductance peak. In general Hmin depends on beta, E FS/E FN, the strength of the insulating barrier Z and on the temperature T.
Zero bias conductance peak; d-wave superconductivity; Magnetic field
SURFACES, INTERFACES, AND THIN FILMS
Effect of a magnetic field on tunneling conductance in normal metal d-wave superconductor interfaces
William J. Herrera; J. Virgilio Niño
Departamento de Física, Universidad Nacional de Colombia, Bogotá- Colombia
ABSTRACT
Through the analytic solutions of the Bogoliubov de Gennes (BdG) equations the effect of a static and homogeneous magnetic field applied parallel to the interface of an NIS (N: Normal metal, S: superconductor and I: Insulator) junction on the differential conductance is calculated. For a dxy - symmetry we obtain zero bias conductance peak that can be split by a magnetic field. The shift of the zero bias conductance peak depends on the spread (b) of the tunneling electrons in k space, on the magnitude of the applied field H and on the ratio between the Fermi energy of the superconductor and the normal region, EFS/EFN. Finally we estimate the minimum value of the magnetic field, Hmin, that splits the zero bias conductance peak. In general Hmin depends on b, EFS/EFN, the strength of the insulating barrier Z and on the temperature T.
Keywords: Zero bias conductance peak; d-wave superconductivity; Magnetic field
I. INTRODUCTION
In high temperature superconductivity different experiments have been interpreted assuming a d-wave symmetry of the pair potential, see for example [1]. A feature of d symmetry is the zero-bias conductance peak (ZBCP) observed in tunneling conductance in NIS junctions [1]-[6]. The ZBCP appears when the angle between the a axis of the superconductor and the vector normal to the interface is nonzero; it is a maximum when this angle is ±p/4, (110) orientation. This peak is due to the formation of the zero energy states (ZES) that are originated by the Andreev reflection at the interface; it undergoes a difference in phase of p due to the anisotropy of the pair potential [7]-[10]. Different experiments have shown that the ZBCP can be split due to application of a magnetic field [2]-[5]. From numerical solutions of the Eilenberger's equations, it has been shown that the effect of the magnetic field is to produce a shift of the ZES proportional to the applied magnetic field [11]-[12], other experiments show [13]-[14] that the effect of the magnetic field is to decrease of the heigth of the ZBCP. In this work we solve the BdG equations for this system and show that the splitting depends on Z, b,c and T, The above mentioned experimental characteristics are explained in this work.
II. THEORY
The quasiparticles in a superconductor are described by the BdG equations. For steady states and anisotropic superconductors these equations are [15]
where He(r1) = ( -iÑ-eA(r1))2/2m+V( r1) -µ is an electronic hamiltonian, with A(r1) the vector potential associated with the magnetic fields present in the system, V(r1) the scalar potential and m the chemical potential, (r1,r2) is the pair potential, u(r1)and v(r1) are the wave function for the electron- and hole-like components of a quasiparticle. The insulating barrier of height V0 and thickness d is located in x > 0, this barrier can be modeled by a delta function, V(x) = U0d(x), where U0 = V0d. We concentrate on cuprate superconductor junctions. It is supposed that the quasiparticle moves on the CuO2 plane with the a and b axes in the x-y plane, the interfaces are normal to the x-axis, see Fig. 1. The Fourier transform of the pair potential is modeled by
where R = r1-r2,r = (r1+r2)/2 , Q(x) the Heaviside function and D() is the pair potential that undergoes a quasiparticle with momentum k, for s-symmetry D() = D0 and for dx-y-symmetry D() = D0sin2qs,where qs is the quasiparticle angle in the superconductor region qs = sin-1(ky/| k| ). The magnetic field H, is applied parallel to the z-axis, therefore the vector potential can be written as A(r) = Ay(x). We consider the situation of a high Tc superconductor where the coherence length of the pair potential x is much smaller than the penetration length l of the magnetic field and approximate Ay = -Hl. As the potentials depends only on x, the solutions of the BdG equations can be written as Y(r) = .
Considering an incoming electron from the normal region, the wave functions (x) and (x) are given taking into account the in the Andreev approximation [15] by
where
The wavenumbers kFN and kFS are determined by the Fermi energy in the normal and superconducting regions respectively, = EFS/EFN = c2. The quasiparticles with and wavenumber move in the pair potential D+ and D- respectively
The effect of the magnetic field is an energy shift that depends on ky and H. One finds a, b, C and D using the boundary conditions in x = 0. The electron-electron and electron-hole reflection coefficients are respectively, Re = |b| 2, Rh = |a| 2.
III. DIFFERENTIAL CONDUCTANCE
Using the model developed by Blonder et al.[16] the differential conductance for an angle qS and for T = 0K is calculated from Re and Rh coefficients as
with
sN is the differential conductance when D = 0 (NIN junction), Z is the strength of the barrier and qN is the quasiparticle angle in the normal region and is determined by the momentum conservation condition in y direction
Firstly the case of c = 1 (qN = qSº q) is analyzed. The relative total differential conductance is found by integration in the k-space as
where
and b is related to the spread of tunneling electron in k- space, it is given by
It is important to note that the insulating barrier is characterized by the parameter Z and b. If in the equation (6) H = 0, our results agree with [7]-[8]. When H ¹ 0 the differential conductance sS(eV,q) has a shift given by VA(q). The Gaussian distribution, , diminishes 99% for a angle given by qc = . The maximum peak energy shift is determined by VA(qmax), where qmax is the maximum angle for which an electron tunnels the insulating barrier. This angle is p/2-e (e® 0) if 0 < b < 4.2 or approximately qc if b > 4.2.
The maximum value of VA is max{D(q)} = D0, therefore the maximum magnetic field is
with Hc the bulk critic magnetic field.
The average shift is dV = á| VA|ñ = elH á| sinq|ñ /m. For b >> 1,dV is approximately
The average shift dV is strongly affected by b, this behavior is due to the fact that as b increases, the tunneling cone diminishes and therefore the shift energy decreases. Fig. 2 shows sR at different values of b, the splitting of ZBCP increases as b decreases. In the inset of the Fig. 2 it is contrated the average shift obtained in Eq. (14) with the numerically value. Fig. 3 shows sR for different values of H with Z =3, ZBCP is split and the shift depends on the magnetic field and b, as shown the inset of Fig. 3. Fig. 4(a) shows sR for different values of H with Z =0.5. Is observed that beyond some value the magnetic field ( Hmin ) the ZBCP is split. This value depends on Z and b as shown the figures 4(a) and 4(b). For an estimation of Hmin one can compare dV with the width G0 of the ZBCP. For an angle q, the width is given by
Averaging over q the width is approximately
As b increases , decreases, this is because the tunneling cone decreases and therefore average gap diminishes. The splitting appears when dV ~ /2, from this relation we obtained Hmin
with p ~ 1. Hmin is independent of b , this is due to the fact that as b increases, the width of ZBCP decreases but also dV decreases and these effects compensate each other. It is important to stand out that although in this case Hmin is independent of b, Hmax depends on b. If H > Hmin the magnetic field induces a splitting of the ZBCP, but if Hmin > Hmax it is not possible to observe this splitting, as is shown in figures 4(a) and 4(b) for b = 5 and b = 20 respectively, and Hmin@ 1.4Hc with p = 0.8. If b = 5, Hmax = 1.2Hc < Hmin and the splitting does not take place. If b = 20, Hmax@ 2.2Hc and ZBCP is split. In both cases if H < Hmin the effect of the magnetic field is to decrease the height of the ZBCP.
Now we study the effect of the Fermi energy difference between the normal and superconductor regions. This difference is quantified by the parameter c = kFS/kFN, the average angle in Eq. (10) is modified by
where qm is the maximum angle in the superconductor, forc < 1 , qm = p/2 and for c > 1, qm = sin-1(1/c) . If bc2 >> 1 from Eq's. (10), (14), (15), dV and 0 are
In this case Hmin and Hmax are
In Figs. 5(a) and 5 (b) sR is plotted against V for c = 0.5 and c = 1.5 respectively. From Eq. (22) Hmin = 0.8Hc and Hmin = 1.6Hc in agreement with the numerical results shown in the figures.
Finally we study the differential conductance at finite temperature, in this case sR is calculated from
with f(E) the Fermi-Dirac distribution function at temperature T: Fig. 6 shows sR for T = 0.1Tc and for different values of Z. The main effect of T is to increase the width of the peak to H = 0 and therefore to increase Hmin. In order to analyze the width of the ZBCP we considerer first the case kBT << 0, we find that width of the ZBCP is
and Hmin is given by
Therefore the value Hmin at finite temperature is larger than Hmin for T = 0. The relative differenceDH/Hmin(0) = (Hmin(T)-Hmin(0))/Hmin(0) is proportional to .gif" align="absmiddle">, and therefore in this case the relative increase of Hmin is always less than one, as shown in Fig, 6(a) where Hmin@ 1.6Hc ; from Eq's. (17) and (22) Hmin(0) = 1.4Hc and Hmin(T) = 1.57Hc as shown in the figures. For kBT >> 0 (T)
In this case the width is determined main by the temperature ~ KBT, Hmin(T) is proportional to T and b1/2, and can be larger than Hmin(0), in contrast to the previous case , see Fig. 6(b).
IV. CONCLUSIONS
Starting from the solutions of the Bogoliubov of Gennes equations in an NIS junctions we have determined the effect of a magnetic field on the differential conductance. For dxy symmetries the differential conductance presents a ZBCP. The magnetic field induces a splitting of the ZBCP when H > Hmin, where Hmin depends on Z, c,and T, in this case the shift of the peak is proportional to H/(b1/2c). When H < Hmin the effect of the magnetic field is to decrease of the height of the ZBCP. If Hmin > Hmax, the splitting does not take place when the magnetic field is applied; the value of Hmax depends on c,and b. This results can be used for the interpretation of the tunneling characteristic in high Tc superconducting junctions in applied magnetic fields.
Acknowledgments
The authors have received support from División de Investigaciones de la Universidad Nacional de Colombia sede Bogotá.
Received on 8 December, 2005
- [1] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72 969 (2000).
- [2] M. Covington, M. Aprili, E. Paraoanu, F. Xu, J. Zhu and C. A. Mirkin, Phys. Rev. Lett. 79, 277(1997).
- [3] M. Aprili, E. Badica and L. H. Greene, Phys. Rev. Lett. 83, 4630 (1999).
- [4] A. Alff, H. Takashima, S. Kashiwaya, N. Terada,H. Ihara, Y. Tanaka, M. Koyanagi, and K. Kajimura, Phys. Rev. B 55 R14757 (1997).
- [5] R. Beck, Y. Dagan, A. Milner, A. Gerber and G. Deutscher, Phys. Rev. B. 69,144506 (2004).
- [6] I. Iguchi, W. Wang, M. Yamazaki, Y. Tanaka and S. Kashiwaya, Phys. Rev. B, 62, R6131 (2000).
- [7] S. Kashiwaya, Y. Tanaka, M. Koyanagi, H. Takashima and K. Kajimura, Phys. Rev. B 51, 1350 (1995).
- [8] S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63,1641 (2000)
- [9] M. B. Walker, P. Pairor, Phys. Rev. B 60, 10395 (1999).
- [10] T. Löfwander, V. S. Shumeiko and G. Wendin, Supercond. Sci. Technol. 14, R53 (2001).
- [11] M. Fogelström, D. Rainer and J. A. Sauls, Phys. Rev. Lett. 79 281 (1997).
- [12] Y. Tanaka, H. Tsuchiura, Y. Tanuma and S. Kashiwaya, J. Phys. Soc. Jpn. 71 271, (2002).
- [13] J. W. Ekin, Yizi Xu, S. Mao, T. Venkatesan, D. W. Face, M. Eddy and S. A. Wolf. Phys. Rev. B 56, 13746 (1997).
- [14] A. Biswas, P. Fournier, M. M. Qazilbash, V. N. Smolyaninova, H. Balci, and R. L. Greene, Phys. Rev. Lett. 88, 207004 (2002).
- [15] C. Bruder, Phys. Rev. B 41, 4017 (1990).
- [16] G. E. Blonder, M. Tinkham and T. M. Klapwijk, Phys. Rev. B 25, 4515 (1982).
Publication Dates
-
Publication in this collection
04 Dec 2006 -
Date of issue
Sept 2006
History
-
Received
08 Dec 2005