Tunneling Conductance in Ferromagnetic Metal / Normal Metal / Spin-Singlet s-Wave Ferromagnetic Superconductor Junctions

In the framework of scattering theory, we study the tunneling conductance in a system including two junctions, ferromagnetic metal/normal metal/ferromagnetic superconductor, where ferromagnetic superconductor is in spin-singlet s-wave pairing state. The non-magnetic normal metal is placed in the intermediate layer with the thickness (L) which varies from 1 nm to 10000 nm. The interesting result which we have found is the existence of oscillations in conductance curves. The period of oscillations is independent of FS and FN exchange field while it depends on L. The obtained results can serve as a useful tool to determine the kind of pairing symmetry in ferromagnetic superconductors.


Introduction
In this paper we have studied tunneling conductance in a two-junction system with structure of ferromagnetic metal/normal metal/ferromagnetic superconductor.A constant spin-singlet -wave pairing potential Δ and a constant exchange field ℎ  are our presuppositions for the ferromagnetic superconductor.The main goal of this paper is to show the existence of oscillations in conductance curves with period of oscillations independent of FS and FN exchange field but dependent on .The structure of this paper is as follows.In Section 1 a historical view of the under study subject together with a short report of the results is presented.In Section 2 the model and the theoretical framework of our system are studied.In Section 3 we have presented and discussed the results of the paper.
For a long time, coexistence of superconductivity and ferromagnetism has been a challenging issue in the history of superconductivity.It was believed that superconductivity and ferromagnetism were two opposite phenomena and coexistence between them was impossible.The view toward ferromagnetism and superconductivity changed when ferromagnetic superconductors were discovered in 2000-2001.
The discovery of superconductivity in some ferromagnetic materials such as  2 [1],  2 [2], and ℎ [3] has attracted much attention in the field of superconductivity.
The questions are mostly raised about the mechanism of cooper pairing in ferromagnetic superconductors.At the first glance, it seems that the pairing symmetry of such materials should be spin triplet equal spin pairing.But some research groups have shown that one could not ignore spin-singlet wave opposite spin pairing [4][5][6][7][8][9][10][11].For example, in [7], Singh and Mazin by calculating the electronic structure, as well as the zone center phonons and their coupling with electrons, have shown that  2 has a spin-singlet -wave pairing state.Also, in [8], due to kinetic exchange, the coexistence between ferromagnetism and -wave pairing has been shown.Also, Dahal et al. [9], based on interaction of electrons mediated by localized spins have shown that the pair symmetry in  2 is spin-singlet -wave state.(Thus, in this paper, we consider spin-singlet -wave pairing for ferromagnetic superconductors.)Other pairing symmetries such as spin triplet opposite spin pairing (OSP) state and spin triplet equal spin pairing (ESP) state are suggested for ferromagnetic superconductors [12,13].These two pairing states are different from SWP state, because their order parameters have a sign Advances in Condensed Matter Physics change on Fermi surface, which leads to the formation of a midgap state while there is not midgap state in the case of SWP state.The formation of midgap states affects the tunneling spectra, so they make difference between different symmetries in terms of the kind of pair symmetries [12].
On the other hand, tunneling spectroscopy can provide much information about the excitation energy, the type of pairing symmetry, and other physical properties of superconductors.For the first time, the theory of tunneling spectroscopy by taking into account Andreev reflection [14] was presented by Blonder et al. (BTK theory) in 1982 [15].After the discovery of high temperature superconductors, Tanaka and Kashiwaya have extended BTK theory to unconventional -wave superconductors in 1995 [16,17].The presence of zero-bias conductance peaks in tunneling spectroscopy of normal metal/unconventional superconductor junctions is accepted as a strong indication of -wave symmetry in high temperature superconductors.
In order to determine the type of pairing symmetry in ferromagnetic superconductors, the detailed comparison between theoretical results and experimental data is required.Anyway, as it is mentioned in [18], the study of tunneling conductance in junctions which consist of ferromagnetic superconductors can settle the problem.Till now, the study of tunneling spectroscopy for ferromagnetic superconductors has been restricted to one junction (see [12,13,19,20]).
In this paper, we theoretically study tunneling conductance in ferromagnetic metal/normal metal/ferromagnetic superconductor (FM/NM/FSc) double junction where ferromagnetic superconductor (FS) is in spin-singlet -wave pairing (SWP) state and has an exchange field ℎ  .The nonmagnetic normal metal is placed in an intermediate layer with a zero barrier potential and the thickness of  which varies from 1 nm to 10000 nm.We find that the tunneling conductance oscillates as a function of bias voltage.The period of oscillations is independent of ℎ  and ℎ  (exchange field of ferromagnetic metal), while it is inversely proportional to .From the obtained results, we can obtain much information about the symmetry of pair potential of ferromagnetic superconductors which is a controversial issue in the field of superconductivity.

Formulation
In Figure 1 we have shown the structure of the system which we study in the current paper.Figure 1(a) shows a schematic diagram of FM/NM/FS double junctions.An NM with the thickness of  is sandwiched by FM and FS layers.Figure 1(b) shows a schematic diagram of energy versus wave vector at FM, NM, and FS regions, respectively.The open circles, the filled circles, and the arrows point denote holes, electrons, and the direction of the group velocity, respectively.The incident electron is indicated by 0 and experiences the following processes: (1) Andreev reflection, (2) normal reflection, (3,4) transmission to the NM, (5,6) reflections at the NM/FS interface, and (7,8) transmission to the FS as electron and hole quasiparticles.We study a system consisting of two junctions, ferromagnetic metal/normal metal/ferromagnetic superconductor (FM/NM/FSc); that is, one junction is placed between a ferromagnetic metal and a normal metal, and the other one is between a normal metal and a ferromagnetic superconductor.The junction NM/FSc is located at  = .The nonmagnetic normal metal is placed between two junctions with the thickness of , and ferromagnetic metal (FM) and the ferromagnetic superconductor (FSc) are located at  < 0 and  > , respectively.The exchange field of the system is in the following form: where ℎ  and ℎ  are the exchange fields of FM and FSc, respectively, and Θ() is the step function which is one for  > 0 and zero for  < 0. The effective pair potential for FSc is defined by , where for spin-singlet wave pairing Δ( ⃗ ) is in the following form: Δ( ⃗ ) = Δ  [12].The motion of quasiparticles is described by Bogoliubov-De Gennes (BdG) equations which have the following form: where  0 = −(ℎ 2 /2)∇ 2 −  is the one particle Hamiltonian.
is the energy of quasi-particles which is measured from Fermi energy;   .⃗  (=  x +   ŷ +   ẑ) is the Pauli matrices where x, ŷ, and ẑ are the unit vectors in the -, -, and directions, respectively.⃗ ℎ is the exchange field of system; we assume that it is in the -direction; ⃗ ℎ = ẑℎ() and ℎ() is defined by (1).Now, we are ready to find the eigenfunction and eigenvalues from (2).By replacing (1) in (2) and solving it, the energy spectrum for an electron is obtained from the following formula: and the energy spectrum of a hole is:  = − ± , where  ± is defined in (3).The corresponding wave functions of spin-singlet -wave state are obtained in the following forms [13]: So, the wave function in  >  region whose ferromagnetic superconductor (FS) is placed there is as follows: where is the wave vector of FS layer.

FM
NM FS In intermediate layer (0 <  < ) whose normal metal (NM) is placed, the wave function is obtained by setting ⃗ ℎ = 0 and Δ = 0 in (3).In this case the solution of ( 3) is obtained as follows: where  ± NM =   √1 ± (/  ) is the wave vector of NM region.
In the  < 0 region where ferromagnetic metal (FM) is placed, the wave function is as follows: where  ± FM =   √1 ± (ℎ  /  ) ± (/  ) is the wave vector in FM region and  and  are the Andreev and normal coefficients, respectively.Now, according to BTK theory, we should replace the above wave functions in the following boundary conditions: where  FM ,  NM , and  FS are the wave functions of FM, NM, and FSc; respectively, and   is the derivative of  with respect to ; that is,   = /.By applying the above boundary conditions, we have obtained Andreev and normal coefficients, and as the result, we have obtained the conductance ( ) in the following form: where we have The normalized tunneling conductance is obtained from this formula:  = (∑    )/ N , where  N is the conductance of N/N/N junction.Before we present the results, to examine the validity of our calculation, we compare (10) with previously done researches.In [13], we have calculated zero-bias conductance for NM/FSc junction.By placing  = 0 in results of [13] and  = 0 and in (10), we come up with the same result.Also, in [21], the tunneling conductance is calculated for FM/Sc junction.In the limit of  → 0, the results of [21] are reproduced.Also, in the limit of ℎ  → 0 the results of [22] are reproduced which shows the validity of our calculations.

Results
Figure 2 shows the conductance curves for  = 1000 nm and ℎ  /Δ = 0.2, 0.4, 0.6.As it can be seen in Figure 2(a), the conductance shows a periodic behavior especially in subgap region ( < Δ).The origin of this periodic behavior comes back to the interference of electron waves in the intermediate layer [23,24].We can interpret this periodic behavior by using the energy diagram of Figure 1(b).In fact the incident electron propagates through the NM layer toward the NM/FS interface as an electron (3 in Figure 1(b)) and is diffused into a hole (6 in Figure 1(b)) by the FS order parameter.The FS order parameter can couple an excited electron with an electron inside the Fermi sea, and it leaves a hole excitation.This hole is scattered back into the NM layer, but it does not interfere with the primary electron.In order to take place an interference at the NM layer, the hole should be reflected at the NM/FS interface, be propagated into the electron state (5 in Figure 1(b)) by the FS order parameter, and then propagate into the NM layer as an electron.So, the interference at the NM layer can occur with the primary electron (3 in Figure 1(b)), and it produces an oscillation in the conductance curves at Figure 2(a).The effect of midgap states on conductance curves is an interesting subject.For example, in the case of SWP state there is not midgap state at the interface.So, we have not zero-bias conductance peak (ZBCP), while as it is mentioned in [12,13] ZBCP exists for the case of OSP and ESP states.Also, we can see with the increase in exchange field, the magnitude of conductance decreases.This is because with the increase in exchange field, the symmetry of spin-singlet wave pairing breaks and so conductance decreases.In Figure 2(a), we can see that in Δ ≪ , normalized conductance tends to 1.It is because when |Δ| is small in comparison with the energy of bias voltage, the FS tends to NM, and so normalized conductance which is defined by ( = (∑    )/  =  F/N/FS / F/N/N ) tends to 1.
In Figures 2(b) and 2(c), we study the period of conductance oscillation.Analytically by replacing (−  ) =  in (9), we find the following equation for the period of oscillations: which is obtained for   = 10 nm −1 .We note that (11) may change when we have a two-dimensional junction with other symmetries.For example, in the case of a -wave symmetry, the interferences take place at two directions, and the period of oscillations may be obtained from ( −   ) cos  = .Therefore, (11) has a different form in the case of other symmetries.We examine (11) numerically.We measure the period of oscillations from Figure 2(a) and we find a good agreement between numerical data and analytic result (see (11)).Solid line in Figure 2(b) is plotted from (11) and circles in Figure 2(b) are obtained from numerical data.As it can be seen from this figure there is good correspondence between numerical and analytical results.
Figure 2(c) shows that the period of oscillations does not depend on the FS exchange field.We have plotted different - curves for different exchange fields and we have found that the period of oscillations is independent of exchange field.This is because the period of oscillations just depends on  and   while it does not depend on exchange fields.A point that should be noted is that Figure 2(c) is plotted for the values of ℎ  < 0.7Δ.This is because at zero temperature for the values of ℎ  > 0.7Δ ferromagnetic superconducting state converts to normal metal state.In fact, as it is mentioned in [13,25], by calculating the difference in thermodynamic potential between superconducting state and normal state, Ω  − Ω  , we see that for ℎ  > 0.7Δ the thermodynamic potential of normal metal is lower than the thermodynamic potential of superconducting state, Ω  < Ω  .So in a first order phase transition process, the superconducting state falls into normal state, and in the presence of an exchange field, superconducting state exists just for ℎ  < 0.7Δ.
Beside spin-singlet -wave pairing (SWP) state, the other symmetries such as spin triplet opposite spin pairing (OSP) state and spin triplet equal spin pairing (ESP) state are suggested for ferromagnetic superconductors [12,13].As it is mentioned at [12,13], one of the characteristic features of OSP and ESP states is the existence of zero-bias conductance peaks in conductance curves.This comes back to the type of their pair potentials.In OSP and ESP states the sign change of pair potential at Fermi surface leads to the formation of midgap Andreev resonant states, and they are responsible for the formation of zero bias conductance peaks.There is not conductance peak at zero-bias for SWP state, because its pair potential does not change at the Fermi surface.In fact the kind of pair potential strongly affects the tunneling conductance and so our results can help distinguish the type of symmetry in ferromagnetic superconductors which is still unknown for physicists.

Figure 1 :
Figure 1: (a) Schematic diagram of a FM/NM/FS double junctions.(b) Schematic diagram of energy versus wave vector of a FM/NM/FS double junctions for a spin-up incident electron.

Figure 2 :
Figure 2: (a) Normalized tunneling conductance versus bias voltage for different exchange fields, (b) period of conductance oscillations versus thickness of intermediate layer, and (c) period of conductance oscillations versus exchange field of ferromagnetic superconductor.