# Shot noise of the conductance through a superconducting barrier in graphene

###### Abstract

We investigated the conductance and shot noise properties of quasiparticle-transport through a superconducting barrier in graphene. Based on the Blonder, Tinkham, and Klapwijk (BTK) formulation, the theory to investigate the transport properties in the superconductive graphene is developed. In comparison, we considered the two cases that are the transport in the presence and absence of the specular Andreev reflection. It is shown that the conductance and shot noise exhibit essentially different features in the two cases. It is found that the shot noise is suppressed as a result of more tunneling channels contributing to the transport when the superconducting gate is applied. The dependence of the shot noise behavior on the potential strength and the width of the superconducting barrier differs in the two cases. In the presence of the specular Andreev reflection, the shot noise spectrum is more sensitive to the potential strength and the width of the superconducting barrier. In both of the two cases, total transmission occurs at a certain parameter setting, which contributes greatly to the conductance and suppresses the shot noise at the same time.

###### pacs:

74.20.-z, 72.70.+m, 72.10.-d## I Introduction

In the past decade, there has been a great deal of interest in studying the physical properties of graphene both theoretically and experimentallyBeenakkerPRL2006 (); BeenakkerRMP2008 (); Castro NetoRMP2009 (). In graphene, the low-energy excitations are massless and chiral Dirac fermions with linear energy dispersion near the Dirac point. Because of the unique energy dispersion in graphene, there are some special phenomena such as the unconventional quantum Hall effect, Klein paradox, etcCastro NetoRMP2009 (); BaiPRB2009 (); BaiPRB2007 (); BaiSuperMicros2011 (); KatsnelsonNature2006 (). Particularly relevant to the present work, the presence of the Dirac Fermions in the graphene-based superconductor junctions results in the specular Andreev reflectionBeenakkerPRL2006 (); BeenakkerRMP2008 ().

The Andreev reflectionAndreevJETP1964 () describes the tunneling phenomenon of electron excitation converting into hole excitation by the superconducting pair potential. The interface between a metal and a superconductor can reflect a negatively charged electron incident from the metal into a positively charged hole, while the missing charge enters the superconductor forming a Cooper pair. Usually the hole is reflected back along the path of the incident electron in the conventional materials, which is called “retro-Andreev reflection”. However, the Andreev reflection in undoped graphene is specular—the so-called ”specular Andreev reflection”, in which the reflected angle is invertedBeenakkerPRL2006 (); BeenakkerRMP2008 (); BaiSuperMicros2011 (); BenjaminPRB2008 (). Because the graphene needs to be described by the Dirac-like equation rather than the usual Schrodinger equation, the superconductive graphene needs to be described by the Dirac–Bogoliubov–de Gennes equation rather than the Bogoliubov–de Gennes equation for usual superconductorsBeenakkerPRL2006 (); BeenakkerRMP2008 (); BhattacharjeePRL2006 (); BaiPLA2010 (); BaiPLA2008 (); BaiAPL2008 (). However, graphene is not a natural superconductor. Recent research has shown that superconductivity can be induced in a graphene layer in the presence of a superconducting electrode by means of the proximity effectBeenakkerPRL2006 (); TitovPRB2006 (); BhattacharjeePRL2006 (). So far, the physical properties of the graphene-based superconductor junctions have been extensively studied and many important results have been obtained. Most of them focus on the transport properties of the graphene/superconductive graphene interface and only the conductance is consideredBeenakkerPRL2006 (); BaiSuperMicros2011 (); BhattacharjeePRL2006 (); BaiPLA2010 (); JongPRB1994 (); BeenakkerPRB1992 (); BlonderPRB1982 (), which prompts us to investigate the transport properties through a superconducting barrier in graphene and especially focus on the shot noise spectrum. In this work, based on the Dirac–Bogoliubov–de Gennes equation and the Blonder, Tinkham, and Klapwijk (BTK) formulation, we derive the shot noise formula and provide numerical results of the transport properties in the graphene/superconductive graphene/graphene (G/SG/G) heterostructure. Dependence of the shot noise spectrum on the structure parameters and its physical mechanisms are discussed.

## Ii Theory and Model

We consider transmission through the G/SG/G heterostructure occupying the - plane, the schematic of which is shown in Fig. 1. The growth direction of the graphene is taken to be the -axis. The left G region extends from to and the right G region extends from to . The superconductive region occupies the region. The electron and hole excitations are described by the Dirac–Bogoliubov–de Gennes equationBeenakkerPRL2006 (); BaiSuperMicros2011 (); BhattacharjeePRL2006 (); BaiPLA2010 (); JongPRB1994 (); BeenakkerPRB1992 (); BlonderPRB1982 ():

(1) |

Here, are the four-component wave functions for the electron and hole spinors. The index denotes or for the electrons or holes near the Dirac and points. takes values for . denotes the Fermi energy. and denote the two inequivalent sites in the hexagonal lattice of graphene. The Hamiltonian is given by

(2) |

where denotes the Fermi velocity of the quasiparticles in graphene and takes the values of for .

For , the superconducting electrode on the top of the graphene layer induces a nonzero pair potential via the proximity effect. We model the pair potential as

(3) |

where and are the amplitude and the phase of the induced superconductive order parameter. The electrostatic potential in the G and SG regions can be tuned independently by a gate voltage or by doping. We take

(4) |

Eq. (1) can be solved straightforwardly to yield the wave functions in the G and SG regions, respectively. In the G region, for the electrons and holes traveling in the directions with a transverse momentum and energy , the wave functions are given by

(5) |

with

(6) |

is the incident angle of the electron and is the reflection angle of the hole. Note that for the Andreev process to take place, the maximum incident angle of the electron is given by

(7) |

In the SG region, the quasiparticles are mixtures of the electrons and holes. The wave functions of these quasiparticles moving along the -direction with the transverse momentum and energy has the form of

(8) |

The parameters , , , and are defined by

(9) |

Taking into account both the Andreev and normal reflection processes, the wave functions in the left G, SG, and right G regions can be written as

(10) |

respectively. Here, and are the amplitudes of the normal and Andreev reflections, respectively; and are the amplitudes of the normal and Andreev transmissions, respectively. , , , and are the amplitudes of electronlike and holeslike quasiparticles in the SG region. All the amplitudes in Eq. (10) can be determined by demanding wave function continuity at the interfaces. These boundary conditions are given by

(11) |

The scattering amplitudes can be obtained by numerically solving these continuity equations.

The electron and hole operators of the outgoing states are related to the electron and hole operators of the incoming states via the scattering matrixBlanterRep2000 (); MuzykantskiiPRB1994 (),

(12) |

where the element gives the outgoing electron current amplitude in response to an incoming electron current amplitude and gives the outgoing hole current amplitude in response to an incoming electron current amplitude. The generalized current operator for the electrons and holes in the left electrode can be written asBlanterRep2000 (); MuzykantskiiPRB1994 (); AnantramPRB1996 ()

(13) |

Substituting the scattering matrix in Eq. (12), we can obtain

(14) |

where

(15) |

The general expression for the current fluctuations between contacts and is

(16) |

with its Fourier transform

(17) |

We restrict our consideration to coherent tunneling and neglect the Coulomb interaction. In the zero-frequency limit with , providing all the information above, we can express the noise power as

(18) |

By introducing the distribution function for the electrons and that for the holes , the zero-temperature conductance can be obtained from the current operator and from the usual quantum statistical assumptions for the averages and correlations of the electron and hole operators in the normal reservoirs as

(19) |

The shot noise power can be obtained as

(20) |

where is the ballistic conductance of metallic graphene, is the bias voltage, and denotes the number of available channels for a graphene sample of width .

## Iii Numerical results and discussion

Now we present numerical results for the Andreev reflection coefficients and the tunneling conductance for the junction with . In this condition, there is a large mismatch of Fermi surfaces on the G and SG sides. Such a mismatch is well known to act as an effective barrierBaiSuperMicros2011 (); BhattacharjeePRL2006 (); BaiPLA2010 (). The transport properties of the two cases of () and are significantly different. In the case of , i.e., the incident electron and the reflected hole both lie in the conduction band, which results in the “retro-Andreev reflection”, while in the case of , only the “specular Andreev reflection” takes place, since the incident electron in the conduction band is converted into the reflected hole in the valence bandBeenakkerPRL2006 (); BaiSuperMicros2011 (); BhattacharjeePRL2006 (); BaiPLA2010 (). In general, it is very difficult to reach the regime in experimentBeenakkerPRL2006 (). So we only consider the case of and the condition of comparable and , in the latter of which both normal Andreev reflection and specular Andreev reflection play roles and we can see that the retroreflection crosses over to specular Andreev reflection.

Firstly, we consider dependence of the transport properties on the SG-barrier thickness. The tunneling conductance through the SG-barrier as a function of the thickness for different potential strengths is shown in Fig. 2. Dimensionless thickness is used. The solid and dashed lines correspond to and 10, respectively. When the bias voltage is small, the conductance exhibits oscillation features. Their oscillation amplitudes decay with increase of the thickness of the SG layer. The oscillation period of the shot noise power is the same as that of the conductance. The shot noise characterises correlations of the current. It can be seen from the curves that the growing trend of the shot noise is opposite to that of the conductance. When the conductance reaches the maximum value, the shot noise approaches the minimum value. This can be interpreted by the relation between the shot noise and the properties of the scatterer. A coherent conductor with all the transmission channels open (The open channel means that the transmission probability is close to one.) has minimal shot noise with the Fano factor approaching 0. A coherent conductor with all closed channels (The closed channel means that the transmission probability is close to zero.) has maximal shot noise with the Fano factor approaching 1. The strength of the shot noise is in the middle of the two limits when a conductor has open and closed channels coexistent. Therefore large transmission probabilities enhance the conductance and suppress the shot noise. It can also be seen in the panels (b) and (d) that for the case of the conductance and shot noise decrease with the increase of the thickness of the SG layer. This is because that the parameter is proportional to for identical and it is in the exponential form of in the wave function as shown in Eq. (8). It is larger for the case of than for the case of , which results in the quick decrease of the shot noise and the conductance. On the other hand, we can see from Fig. 2 that with the increase of the potential strength, the tunneling conductance decreased, which illustrates that the Fermi surface mismatch between the normal and superconducting regions suppresses transmission. We also considered the dependence of the tunneling conductance on the value of for small bias voltages. As expected, we found that the oscillation amplitude decreases monotonically with the increase of in the case of and finally approaches a constant value.

In Figs. 3 and 4, we provide numerical results of the tunneling conductance as a function of the bias voltage . Similar results to Beenakker’sBeenakkerPRL2006 () are obtained. In the limit of , all the conductances have the same value of . For the case of , a sharp change in the conductance occurs at and all the conductances vanish at , which is shown in Fig. 3. This is because of that no Andreev reflection occurs for all the incident angles (the critical angle of incidence ) when . For the small bias voltages before the turning point , the conductance decreases with the increase of the SG energy gap. The conductance curves exhibit oscillatory behavior in the region of , which is different from the condition considered by Beenakker. By analyzing the transmission coefficients of the system, we found that in the condition of the transmission spectrum demonstrates oscillatory behavior and the oscillating period increases with the increase of the SG energy gap. The oscillatory behavior originates from the effect of the quantum-mechanical interference between the electron-like and hole-like quasiparticlesMcMillanandPRL1966 (); PasanaiPhysica.C2014; DongPRB2003 () in the SG barrier. This effect gives rise to oscillations in the reflection and transmission probabilities for the incident energies larger than the gap energy.

In Fig. 4, numerical results of the conductance for different SG-barrier thickness are provided. It can be seen that the conductance is sensitive to the SG-barrier thickness, especially in the condition of . In this condition, the oscillation period increases and the oscillation amplitude decreases sharply with the increase of the thickness. This is also a result of the exponential term in the wave function. It can be interpreted by the limit. In this limit, the components of the wave function with the exponential term are nonphysical, therefore only two of the components of wave function in the SG-region are physical. The model reduces to that of the the G/SG-junction, in other words, the model proposed by BeenakkerBeenakkerPRL2006 () is obtained.

We also considered the shot noise properties of the conductance through the G/SG/G structure. Numerical results of the shot noise and the Fano factor are provided in Fig. 5. The parameters of panels (a) and (c) are the same with Fig. 4 and teh parameters of (b) and (b) are the same with Fig. 3. We can see from the curves that the shot noise oscillates greatly in the region of , while it increases monotonically in the region . These behaviors are similar to the tunneling conductance. In tunneling through the G/SG/G-structure, transmission is enhanced by the active hole channels. As a result, the values of the shot noise are small in comparison with the Poisson value, which is corresponding to the uncorrelated transport. When the hole channels in addition to the electron channels contribute to the transport, the interference effect is strong and the shot noise is significantly suppressed. In the region of , the impact of the proximity effect in the SG-region is strong giving rise to strong conductance and the values of and approach 0. In the case of , the amplitude of oscillation decreases with the increase of ; in the case of , the amplitude of oscillation increases with the increase of , which originates from the same reason as the conductance.

## Iv Conclusions

Based on the Dirac–Bogoliubov–de Gennes equation and the scattering theory, we investigated the transport properties of the relativistic electrons and holes through the G/SG/G junction. We have deduced the analytical formulas of the tunneling conductance and the shot noise. Numerical results of the tunneling conductance and shot noise in the system are provided. We compared the two cases, one of which is in the presence of the specular Andreev reflection and the other of which is in the absence of the specular Andreev reflection. The physical results can be summarized as follows. Firstly, the conductance increases with the increase of the thickness of the SG-layer and the shot noise is suppressed by the conductance in the case of . Secondly, the potential strength significantly affects the transport properties. It suppresses the conductance and enhances the shot noise. Thirdly, we obtained similar results with the model of the G/SG junction in the condition . In the limit of , the results of Beenakker for a G/SG junction can be reproduced. Fourthly, the thickness of the SG-layer affects the conductance and the shot noise more prominently in the condition , causing the decrease of the oscillation amplitudes and the characteristic features of the specular Andreev reflection. In conclusion, the conductance is a combined result of the Andreev reflection and the specular Andreev reflection, which can be tuned by the system parameters; the shot noise is suppressed by the SG-barrier because of the contribution of the hole channels in addition to the electron channels.

## V Acknowledgements

This project was supported by the National Natural Science Foundation of China (No. 11004063), and the Fundamental Research Funds for the Central Universities, SCUT (No. 2014ZG0044).

## References