# Phonon bottleneck in graphene-based Josephson junctions at millikelvin temperatures

###### Abstract

We examine the nature of the transitions between the normal and the superconducting branches of superconductor-graphene-superconductor Josephson junctions. We attribute the hysteresis between the switching (superconducting to normal) and retrapping (normal to superconducting) transitions to electron overheating. In particular, we demonstrate that the retrapping current corresponds to the critical current at an elevated temperature, where the heating is caused by the retrapping current itself. The superconducting gap in the leads suppresses the hot electron outflow, allowing us to further study electron thermalization by phonons at low temperatures (K). The relationship between the applied power and the electron temperature was found to be , which we argue is consistent with cooling due to electron-phonon interactions.

###### pacs:

74.45.+c, 72.80.Vp, 63.22.Rc, 65.80.CkThe electron-phonon interaction and the thermal properties of graphene have attracted a lot of attention (for an extensive review, see e.g. Ref. Balandin, .) Most studies focused on relatively high temperatures; however, this year several experiments were performed at temperatures below 10 Kelvin Voutilainen (); Vora (); Yan (); Fong (); Betz (); Price (). Measuring the intrinsic thermal properties of graphene becomes more challenging in this regime. In particular, making normal metal contacts to graphene would provide the dominant thermalization path, effectively shunting the electron cooling by phonons Price (). In order isolate the phonon contribution to thermal transport in graphene below 1K, we have contacted the graphene crystal with electrodes made from lead (Pb), which becomes superconducting below K. Superconducting leads exponentially suppress the thermal transport of hot electrons outside of the sample, allowing one to study thermalization by phonons.

We study several superconductor-graphene-superconductor (SGS) junctions, and focus on the hysteresis between the switching current (from the superconducting to the normal state) and the retrapping current (from the normal to the superconducting state), which is observed despite the fact that our junctions are overdamped tinkham_book (). We demonstrate that the hysteresis is due to electron overheating Pekola_08 (); Albert_11 (). Multiple contacts are made to the same graphene crystal, and resistive heating is applied at one location; simultaneously, we measure the electron temperature by monitoring the change in the critical current through a pair of contacts at a different location on the crystal. By changing the pairs of contacts where heating is applied and where the critical current is measured, we conclude that the graphene crystal is well thermalized, i.e. its electrons maintain a uniform temperature. We evaluate the electron temperature as a function of applied power, which allows us to estimate the electron-phonon scattering rate in graphene at sub-Kelvin temperatures.

Graphene was deposited on the Si/SiO substrate using the standard mechanical exfoliation recipe novoselov_2005 () and verified to be a single layer by Raman spectroscopyRaman (). Six parallel superconducting strips (500 nm wide) were deposited on top of the graphene crystal by thermally evaporating a 4 nm contact layer of palladium (Pd) followed by 120 nm of lead (Pb)Ivan_2011 (). Thereby, a total of 5 superconductor-graphene-superconductor (SGS) junctions were created with lengths of 0.3, 0.7, 1, 1.5 and 2 m, labeled 1 through 5, respectively (Fig. 1a, b). The widths were roughly similar for all the junctions at around m. In order to measure the relatively small superconducting currents, the sample was encased in a copper box thermalized at the millikelvin temperature. The wires connecting the sample to the measurement setup were filtered by cold RC filters and thermocoax cables.

At the base temperature of 35 mK the shorter junctions support a gate-dependent supercurrent (Figure 1c). The switching currents from the superconducting to normal state at the fully open gate voltage of V were nA, nA, and nA for junctions respectively. The two longest junctions (junctions and ) never developed a supercurrent. These values are suppressed compared to the naive estimates of the critical current ( is the normal resistance of the junction) most likely due to disorder, as discussed in our earlier publication Ivan_2011 ().

Nanoscale hybrid Josephson junctions tend to be underdamped due to the large capacitance shunting the junction’s bonding pads to the back gate. (The damping factor of the Josephson junction is related to the critical current , junction capacitance and the shunting resistance as tinkham_book ().) We made the leads connecting the pads to the junctions somewhat resistive (a few hundred ) to partially isolate the junctions from the pad capacitance and to introduce dissipation. As a result, we estimate that our junctions are in the overdamped regime Vion (). Indeed, successive current sweeps produce a very narrow distribution of the switching and retrapping currents (see inset to Figure 2b; also note that the variations of the switching current seen in Figure 1c are reproducible on successive sweeps of .) Therefore, the observed hysteresis between the switching and retrapping current cannot be due to an underdamped dynamics, and we should seek an alternative explanation.

In order to investigate the nature of the hysteresis, the graphene flake was heated locally by sending a current through one of the five junctions, while measuring the critical current of another junction Voutilainen (). For example, junction 1 would be measured while the heating current is sent through junction 5 (Figure 2a). At zero heating current, the values of the switching and retrapping currents and are the largest, and both decrease as the heating power is increased. The difference between and also decreases and the hysteresis disappears after some intermediate value of .

While the switching current is reached when the measured junction is in the superconducting state, the retrapping current is realized when the measured junction is in the normal state. Therefore, at the retrapping transition, additional power is applied to graphene. Taking into account this additional power, we plot and in Figure 2b. (Here, is the power applied at the heater junction note1 ().) and clearly fall on the same curve, demonstrating that at the same total dissipated power, and coincide.

The curves of and versus map onto each other regardless of the junction that was measured (junctions 1 or 2) and for different gate voltages applied to graphene. Therefore, we can conclude that the retrapping current is suppressed compared to the switching current due to the heat dissipated by the retrapping current itself; the hysteresis between the switching and retrapping currents is thus due to the self-heating of the measured junction. Indeed, it has been demonstrated that self-heating is the origin of the hysteretic behavior in overdamped SNS junctions Pekola_08 (); Albert_11 (). This conclusion allows us to characterize the electron cooling in graphene in the remainder of the paper.

Instead of heating the graphene crystal locally by the current , we can heat the entire sample. We measure the switching current in the two shortest junctions versus temperature at several different gate voltages. By comparing the vs. data with the vs. data, we are able to extract the temperature up to which the electrons in the measured junction are heated for a given . This way, a plot of vs. has been extracted for measured junction 1 while heating either junction 3 or 5 (Figure 3a). Similarly, in Fig 3b, we measured junction 2 while heating junction 1. (In the latter case, large enough current was passed through junction 1 to overcome its own supercurrent. Once that happens, graphene heats to about 200 mK, which explains the gap in the data at lower temperatures.)

The three curves in Figure 3 are remarkably close, although the three setups have different distances and configurations of the heater and the measured junctions, implying that the entire region of graphene crystal in the vicinity of the junctions is thermalized to a uniform electron temperature. Empirically, the curves can be fitted by the same power law . This dependence rules out cooling by the hot electron diffusion into the superconducting leads. Indeed, the dependence would be exponentially suppressed in temperature due to the gap in the quasiparticle density of states of the superconductor. Even at 1 K, the dissipated power ( W) is several orders of magnitude larger than that estimated from the Wiedemann-Franz law, properly adjusted for the suppressed density of states in the leads ( W). This situation is different from the regime realized in Ref. Voutilainen, , where a relatively small superconducting gap in Al leads allows them to provide a dominant thermalization path. We conclude that in our sample electrons in the graphene crystal are thermally decoupled from the leads and must be cooled by phonons.

Furthermore, the observed cannot be limited by the Kapitza resistance between graphene and the substrate, which is estimated to be a few orders of magnitude smaller. (As an upper limit for the resistance, we take the values measured at 40 K in Ref. Chen_09, and extend the cooling power according to or .) Hence, graphene lattice must be well thermalized with the substrate. We conclude that the bottleneck for cooling in our system is the weak thermal coupling between electrons and phonons in graphene.

Following the derivation of Ref. DasSarma_08, , we estimate the electron-phonon cooling power in graphene as . Here, is the area of graphene, is the deformation potential constant, is the electron density, is graphene mass density, is the speed of sound, and is the electron temperature. We assume that at low temperatures the wavelength of the emitted phonon , which enters the scattering matrix element, should be replaced by a distance of the order of a hundred nm (e.g. the distance between the leads, or the electron mean free path). We also took into account that Ref. DasSarma_08, calculated the transport time while we are interested in the scattering time, because each phonon emission event typically results in the electron energy loss of the order of .

Plugging in the numbers from Refs. DasSarma_08, ; Kim_10, , we estimate the cooling power at , within an order of magnitude from the measured value. [This is a rather crude estimate, in part because we do not know how large is the thermalized area of graphene crystal, which extends beyond the immediate vicinity of the contacts, see Figure 1a. This area determines , taken to be m in this estimate.]

Recent measurements of Refs. Fong, ; Betz, report cooling power for hot electron , indicating that the emitted phonon wavelength is not cut off, as we assumed above. The difference may stem both from the lower temperatures in our measurement and a more restricted sample geometry resulting in shorter in our case. The cooling power of Ref. Fong, , extrapolated down to 1 K would yield 0.07 W/m, which in fact is close to our result for m.

Finally, dependence due to supercollision cooling Levitov () has been just observed in Refs. Betz2, ; Graham, . This regime is realized if the wavevector of the emitted phonon, , exceeds the electron Fermi wavevector ; the emission process is enabled by the disorder. We work in the opposite regime, , so that the theory of Ref. Levitov, cannot be directly applied.

In conclusion, we have shown that the difference between the switching and the retrapping currents in our graphene Josephson junctions is caused by electron overheating in the normal state. The superconducting contacts thermally isolate the graphene crystal from the leads, allowing us to measure the electron temperature rise for a given dissipation power, and hence the electron-phonon energy transfer rate. The observed power law dependence is consistent with theory of electron-phonon interactions in graphene and with other measurements. Due to small electron heat capacitance and their decoupling from phonons, this type of sample may be useful for detector applications Prober ().

## References

- (1) A. A. Balandin, Nature Materials, 19, 569 (2011).
- (2) J. Voutilainen, A. Fay, P. Häkkinen, J. K. Viljas, T. T. Heikkilä, and P. J. Hakonen, Phys. Rev. B 84, 045419 (2011).
- (3) A. C. Betz, F. Vialla, D. Brunel, C. Voisin, M. Picher, A. Cavanna, A. Madouri, G. F‘eve, J.-M. Berroir, B. Plaçais, and E. Pallecchi, Phys. Rev. Lett. 109, 056805 (2012).
- (4) K. Fong and K. Schwab, Phys. Rev. X 2, 4 (2012).
- (5) H. Vora, P. Kumaravadivel, B. Nielsen, and X. Du, Appl. Phys. Lett. 100, 153507 (2012).
- (6) J. Yan, M.-H. Kim, J. A. Elle, A. B. Sushkov, G. S. Jenkins, H. M. Milchberg, M. S. Fuhrer, and H. D. Drew, Nature Nanotech. 7, 472 (2012).
- (7) A. S. Price, S. M. Hornett, A. V. Shytov, E. Hendry, and D. W. Horsell, Phys. Rev. B 85, 161411(2012).
- (8) M. Tinkham, Introduction To Superconductivity (McGraw-Hill, 1996).
- (9) H. Courtois, M. Meschke, J. T. Peltonen, and J. P. Pekola, Phys. Rev. Lett. 101, 067002 (2008).
- (10) P. Li, P. M. Wu, Y. V. Bomze, I. V. Borzenets, G. Finkelstein, A. M. Chang, Phys. Rev. Lett. 107 137004 (2011) and Phys. Rev. B 84, 184508 (2011).
- (11) K. S. Novoselov, D. Jiang, F. Schedin, T.J. Booth, V. V. Khotkevich, S. V. Morozov, and A. K. Geim, PNAS 102, 10451 (2005).
- (12) A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. Novoselov, S. Roth, and A. K. Geim, Phys. Rev. Lett. 97 187401, (2006).
- (13) I. V. Borzenets, U. C. Coskun, S. J. Jones, and G. Finkelstein, Phys. Rev. Lett. 107, 137005 (2011) and IEEE Trans. Appl. Supercond. 22, 1800104 (2012).
- (14) D. Vion, M. Götz, P. Joyez, D. Esteve, and M. H. Devoret, Phys. Rev. Lett., 77,16, 3435, (1996).
- (15) We used as a fitting parameter to ensure the best match between the and curves. The resulting values of are close to the measured two-probe resistances of the heater junctions. The difference between the two-probe resistance and is explained by a more sophisticated set-up used in the heating measurement, where a fraction of the heating current can flow to several of the grounded contacts other than the designated heater ground. These parallel current paths reduce the effective resistance of the heater.
- (16) Z. Chen, W. Jang, W. Bao, C. N. Lau, and C. Dames, Appl. Phys. Lett. 95, 161910 (2009).
- (17) E. H. Hwang and S. Das Sarma, Phys. Rev. B 77, 115449 (2008).
- (18) D. K. Efetov and P. Kim, Phys. Rev. Lett 105, 256805 (2010).
- (19) J. C. W. Song, M. Y. Reizer, and L. S. Levitov, Phys. Rev. Lett. 109, 106602 (2012).
- (20) Supercollision cooling in undoped graphene
- (21) Photocurrent measurements of supercollision cooling in graphene
- (22) C.B. McKitterick, D.E. Prober, and B.S. Karasik, arXiv:1210.5495v3 (2012).