# Yu-Shiba-Rusinov screening of spins in double quantum dots

A magnetic impurity coupled to a superconductor gives rise to a Yu-Shiba-Rusinov (YSR) state inside the superconducting energy gap YuAPS1965 (); ShibaPTP1968 (); Rusinov1969 (). With increasing exchange coupling the excitation energy of this state eventually crosses zero and the system switches to a YSR groundstate with bound quasiparticles screening the impurity spin by Satori1992 (); Bauer2007 (); ZitkoPRB2011 (); YaoPRB2014 (); Hatter2015 (). Here we explore InAs nanowire double quantum dots tunnel coupled to a superconductor and demonstrate YSR screening of spin-1/2 and spin-1 states. Gating the double dot through 9 different charge states, we show that the honeycomb pattern of zero-bias conductance peaks, archetypal of double dots coupled to normal leads Jeong2001 (); vanderWiel2002 (); Chorley2012 (), is replaced by lines of zero-energy YSR states. These enclose regions of YSR-screened dot spins displaying distinctive spectral features, and their characteristic shape and topology change markedly with tunnel coupling strengths. We find excellent agreement with a simple zero-bandwidth approximation, and with numerical renormalization group calculations for the two-orbital Anderson model.

Yu-Shiba-Rusinov states can be imaged in a direct manner by scanning-tunneling spectrocopy of magnetic adatoms on the surface of a superconductor Heinricharxiv2017 (). Using superconducting tips, high-resolution bias spectroscopy of multiple sub-gap peaks reveals an impressive amount of atomistic details like higher angular momentum scattering channels, crystal-field splitting and magnetic anisotropy Yazdani1997 (); Heinricharxiv2017 (); Ji2008 (); Ruby2016 (); Franke2011 (); Hatter2015 (); Ruby2015 (); Ruby2016 (); ZitkoPRB2011 (). In general, however, it can be an arduous task to model the complex pattern of sub-gap states Ji2008 (); Hatter2015 (); Ruby2016 (), let alone to calculate their precise influence on the conductance Ruby2015 ().

In contrast, the ’atomic physics’ of Coulomb blockaded quantum dots is simple. Changing the back-gate voltage, subsequent levels are filled one-by-one and the different charge states alternate in spin, or Kramers degeneracies for dots with spin-orbit coupling, between singlet and doublet. With normal metal leads, charge states with spin-1/2 exhibit zero-bias Kondo resonances at temperatures below the Kondo temperature, , reflecting a Kondo-screened singlet groundstate. If the leads are superconducting with a large BCS gap, , this resonance is quenched and the groundstate recovers its doublet degeneracy. The system now displays a YSR singlet excitation close to the gap edge, which can be lowered in energy by increasing PilletNP2010 (); Deacon2010a (); Grove-Rasmussen2009 (); Lee2014a (); ChangPRL2013 (); Jellinggaard2016 (). Close to it crosses zero and becomes the YSR-screened singlet groundstate.

YSR states were first discussed in the context of gapless superconductivity arising in the presence of randomly distributed paramagnetic impurities YuAPS1965 (); ShibaPTP1968 (); Rusinov1969 (). However, the ability to assemble spins into dimers, chains and lattices, has recently prompted the exciting idea of engineering YSR molecules YaoPRB2014 (); Kezilebiekearxiv2017 (), YSR sub-gap topological superconductors and spiral magnetic states Heinricharxiv2017 (); Nadj-Perge2014 (); PientkaPRB2013 (); Schecter2016 (). Quantum dots have the advantage of being tunable via electrical gates, and this plays an important role in recent proposals for topological superconductivity in systems of coupled quantum dots FulgaNJPh2013 (); SuNatComm2017 (); ShermanNatNano2017 (). Here we utilize this electrical control to manipulate YSR states in a double quantum dot (DQD) formed in an InAs nanowire. Using multiple finger gates to tune the total DQD spin and the interdot coupling, we demonstrate control of the YSR phase diagram, including electrical tuning between YSR singlets, and a novel YSR doublet arising from the screening of an excited spin triplet.

A scanning electron micrograph of an actual device (Device A) is shown in Fig. 1a, where bottom gates are used to define a normal-DQD-superconductor structure Jellinggaard2016 (). The corresponding schematic is shown in Fig. 1b, where plunger gates labelled g and g control left (QD) and right (QD) quantum dot, respectively, while an auxiliary gate, g, tunes the interdot tunneling barrier. The essential physics of this system can be understood in terms of a simple zero-bandwidth (ZBW) model in which the superconductor is modelled by a single quasiparticle coupled directly to an orbital in QD via and indirectly to QD through . Figure 1c shows the corresponding energy diagram in the regime of dominating on-site charging energies. In the Supplementary Information (Sec. III) we compare this model to NRG calculations to establish its reliability as a quantitative tool.

In Fig. 2a-b, we reproduce the well-known sub-gap state behaviour for a single dot coupled to a superconductor within the ZBW model. The panels show excitation energy as a function of the dimensionless gate voltage (corresponding to the noninteracting average occupation of QD) for weak and strong . As expected the sub-gap excitations cross (do not cross) zero energy for weak (strong) coupling. The ground state of the system for odd occupancy is thus a doublet or a YSR singlet (screened spin) Satori1992 (); Deacon2010a (); Jellinggaard2016 ().

In Figs. 2f-i we extend the ZBW model to a DQD (finite ) and calculate stability diagrams for increasing . For weak coupling the characteristic honeycomb pattern is observed similar to DQD in the normal state vanderWiel2002 (). However, as increases, entirely new types of stability diagrams emerge. In Fig. 2g, the pattern resembles two mirrored arcs, where the lack of zero energy excitations as a function of for even occupation of QD is due to a doublet to singlet transition in QD (see Figs. 2a,b). Moreover, as the coupling increases even further, the ground state in the region becomes a YSR doublet altering the stability diagram to vertically shaped rectangular regions.

To understand this behaviour, we show the states of the system in the (11) region in Fig. 2c. Two electrons in the DQD may form either a singlet or a triplet state with energy splitting . Due to the superconductor, a third state may also exist in the gap. In analogy with the QD-S system, where a doublet state may be screened to a YSR singlet, the triplet state may be screened to form a YSR doublet (called ) ZitkoPRB2011 (). The energy of and versus coupling is plotted in Fig. 2c, and the latter eventually becomes the ground state at strong coupling.

The relevant ground states and corresponding phase diagram ( versus ) of the DQD-S system with two ground state transitions is shown in Fig. 2d-e. The first occurs when the system transitions from a honeycomb pattern to the case where the spin in QD is screened. The latter regime we call partially screened (PS) since only some of the charge states are affected. The second transition happens when the in (11) becomes the ground state. This regime we name fully screened (FS) since all possible screened states are ground states of the system.

With the qualitative behaviour of this system in place, we explore the different regimes experimentally. The honeycomb regime is presented in the Supplementary Information (Sec. II.B, Device B), while below we focus on the stronger coupled regimes. Figure 3a shows linear conductance versus plunger gates for a two-orbital DQD shell (i.e. one spin-degenerate level in each dot). A pattern of two arcs are observed resembling the partially screened regime. To verify that the conductance resonances originate from sub-gap states, gate traces for different fillings of the two dots are measured. Figures 3c-e trace out the filling of electrons in QD along the red arrows in Fig. 3a, keeping the electron number in QD constant. The sub-gap spectroscopy plots c,e for even filling of QD show similar behaviour, differing from d with odd occupation. When fixing (sweeping) the occupation of QD (QD), the qualitative behaviour is switched (Figs. 3f-h corresponding to green arrows in Fig. 3a). For even occupancy in QD (f,h) no zero-bias crossing is observed, while the opposite is true for odd occupancy (g). In particular Figs. 3d,g are interesting, since they involve the (11) charge state region. In contrast to single dot systems, the singlet ground state shows different behaviour whether tuning the electrochemical potential of the dot close to the superconductor or the normal lead, i.e. concave and convex excitation behaviour versus gate voltage in the (11) state. The experimental data clearly confirm that the resonances in the stability diagram originate from sub-gap excitations. The stability diagram generated by our DQD-S model for realistic parameters reproducing the experimental behaviour is shown in Fig. 3b, and corresponding gate traces for fixed occupations are shown in Figs. 3i-n. The qualitative agreement between theory and experiment is striking and even subtleties like the asymmetry of the sub-gap resonance splitting in j (see arrows) are reproduced.

The transition between different YSR states can also be driven by changing the singlet-triplet (-) splitting by tuning (cf. Fig. 2e). In Fig. 4e-h we show calculated diagrams for in a parameter range, where the ground state in the (11) charge state transitions from (h) to (e), i.e. from double-arc, to vertical-lines diagrams. The corresponding measured stability diagrams for the two orbitals analyzed in Fig. 3 are shown in Fig. 4a-d, where the gate voltage between the two dots are tuned to more negative values (decreasing ). The effect of this tuning qualitatively follows the expectation of the model: a transition from to in the fully screened regime where all spin states are YSR screened.

The gate-dispersion of sub-gap excitations also shows good overall correspondence between ZBW modeling and experiment. We measured the sub-gap spectra in the fully screened regime along the red and green arrows in c and a. The first case c is almost at the transition, where the singlet and doublet states are degenerate in (11). Figures 4i,j show sub-gap states versus and , respectively, with a zero bias peak at V reflecting a degeneracy at this value of t. The corresponding ZBW modelling in Figs. 4k,l (for meV) places the system just barely in the fully screened regime with a (11) groundstate and a nearby excitation dispersing very much like in the measurement. A triplet state is predicted inside the gap, and should be accessible from the groundstate. As demonstrated in the Supplementary Information (Sec. III.A.), this is confirmed by more accurate NRG calculations, which however reveal a strong suppression of spectral weight on this state, explaining why it may be difficult to observe in experiment. For even lower , case a, Figures 4m-p again show good overall correspondence between experiment and theory, except for the triplet state, which should be weak, and in this case hardly resolved within the linewidth broadening in the data. Future experiments with hard gap superconductors or improved resolution may eventually lead to capability to detect even such low-weight spectral features. A detailed discussion of nonlinear conductance and broadening of the YSR sub-gap spectra is provided in the Supplementary Information (Sec. IV).

## Methods

The devices are made by defining bottom gate Au/Ti (12/5 nm) electrodes (pitch 55 nm) on a silicon substrate capped with 500 nm SiO followed by atomic layer deposition of 3x8 nm HfO. InAs nanowires (70 nm in diameter) appropriately aligned on bottom gate structures are contacted by Au/Ti (90/5 nm) normal and Al/Ti (95/5 nm) superconducting electrodes separated by nm. The superconducting film has a critical field of around 85 mT. More details on device fabrication can be found in Ref. Jellinggaard2016, . The samples are mounted in an Oxford Instruments Triton 200 dilution refrigerator with base temperature of around 30 mK and are measured with standard lockin techniques. For the data (Fig. 3,4) shown in the partially screened regime, the voltages on the gates define a double dot potential with values (V) , , , , , , and . Here the gate numbers correspond to gates from left to right in Fig. 1a. Gates 3, 4, and 5 are thus the left plunger, the tunnel barrier and right plunger gates, which are tuned within some range of the values stated.

## Acknowledgements

We would like to thank M.C.Hels. Research was supported by the Center for Quantum Devices, The Danish National Research Foundation, Carlsberg Foundation, The Independent Research Fund Denmark (Natural Sciences), the FP7 FET-Open SE2ND project and the Slovenian Research Agency (ARRS) under Program P1-0044 and J1-7259.

## Author contributions

K.G-R. and A.J. made the measurements, A.J. fabricated the devices and M.H.M. grew the nanowires. K.G-R., A.J. and J.N. designed the double quantum dot experiments, G.S., K.G-R. and J.P. made the ZBW analysis, G.S. and J.P. the conductance asymmetry calculations, and R.Z. provided the NRG analysis. K.G-R., G.S., J.P., R.Z. and J.N. participated in discussions, analysis and wrote the paper.

## References

- (1) Yu, L. Bound State in Superconductors with Paramagnetic Impurities. Acta Phys. Sin. 21, 75 (1965).
- (2) Shiba, H. Classical Spins in Superconductors. Prog. Theor. Phys. 40, 435–451 (1968).
- (3) Rusinov, A. I. Superconductivity near a paramagnetic impurity. JETP Lett. 9, 85–87 (1969). [Zh. Eksp. Teor. Fiz. 9, 146 (1968)].
- (4) Satori, K., Shiba, H., Sakai, O. & Shimizu, Y. Numerical Renormalization Group Study of Magnetic Impurities in Superconductors. J. Phys. Soc. Japan 61, 3239–3254 (1992).
- (5) Bauer, J., Oguri, A. & Hewson, A. C. Spectral properties of locally correlated electrons in a Bardeen-Cooper-Schrieffer superconductor. J. Phys.: Cond. Mat. 19, 486211 (2007).
- (6) Žitko, R., Bodensiek, O. & Pruschke, T. Effects of magnetic anisotropy on the subgap excitations induced by quantum impurities in a superconducting host. Phys. Rev. B 83, 054512 (2011).
- (7) Yao, N. Y. et al. Phase diagram and excitations of a Shiba molecule. Phys. Rev. B 90, 241108 (2014).
- (8) Hatter, N., Heinrich, B. W., Ruby, M., Pascual, J. I. & Franke, K. J. Magnetic anisotropy in Shiba bound states across a quantum phase transition. Nat. Commun. 6, 8988 (2015).
- (9) Jeong, H., Chang, A. M. & Melloch, M. R. The Kondo Effect in an Artificial Quantum Dot Molecule. Science 293, 2221–2223 (2001).
- (10) van der Wiel, W. G. et al. Electron transport through double quantum dots. Rev. Mod. Phys. 75, 1–22 (2002).
- (11) Chorley, S. J. et al. Tunable Kondo physics in a carbon nanotube double quantum dot. Phys. Rev. Lett. 109, 156804 (2012).
- (12) Heinrich, B. W., Pascual, J. I. & Franke, K. J. Single magnetic adsorbates on s-wave superconductors (2017). eprint ArXiv1705.03672.
- (13) Yazdani, A., Jones, B. A., Lutz, C. P., Crommie, M. F. & Eigler, D. M. Probing the local effects of magnetic impurities on superconductivity. Science 275, 1767–1770 (1997).
- (14) Ji, S.-H. et al. High-resolution scanning tunneling spectroscopy of magnetic impurity induced bound states in the superconducting gap of Pb thin films. Phys. Rev. Lett. 100, 226801 (2008).
- (15) Ruby, M., Peng, Y., von Oppen, F., Heinrich, B. W. & Franke, K. J. Orbital Picture of Yu-Shiba-Rusinov Multiplets. Phys. Rev. Lett. 117, 186801 (2016).
- (16) Franke, K. J., Schulze, G. & Pascual, J. I. Competition of Superconducting Phenomena and Kondo Screening at the Nanoscale. Science 332, 940 (2011).
- (17) Ruby, M. et al. Tunneling processes into localized subgap states in superconductors. Phys. Rev. Lett. 115, 087001 (2015).
- (18) Pillet, J.-D. et al. Andreev bound states in supercurrent-carrying carbon nanotubes revealed. Nature Physics 6, 965–969 (2010).
- (19) Deacon, R. S. et al. Tunneling Spectroscopy of Andreev Energy Levels in a Quantum Dot Coupled to a Superconductor. Phys. Rev. Lett. 104, 076805 (2010).
- (20) Grove-Rasmussen, K. et al. Superconductivity-enhanced bias spectroscopy in carbon nanotube quantum dots. Phys. Rev. B 79, 134518 (2009).
- (21) Lee, E. J. H. et al. Spin-resolved Andreev levels and parity crossings in hybrid superconductor-semiconductor nanostructures. Nat. Nanotechnol. 9, 79–84 (2014).
- (22) Chang, W., Manucharyan, V. E., Jespersen, T. S., Nygård, J. & Marcus, C. M. Tunneling Spectroscopy of Quasiparticle Bound States in a Spinful Josephson Junction. Phys. Rev. Lett. 110, 217005 (2013).
- (23) Jellinggaard, A., Grove-Rasmussen, K., Madsen, M. H. & Nygård, J. Tuning Yu-Shiba-Rusinov states in a quantum dot. Phys. Rev. B 94, 064520 (2016).
- (24) Kezilebieke, S., Dvorak, M., Ojanen, T. & Liljeroth, P. Coupled Yu-Shiba-Rusinov states in molecular dimers on NbSe (2017). eprint ArXiv1701.03288.
- (25) Nadj-Perge, S. et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602–607 (2014).
- (26) Pientka, F., Glazman, L. I. & von Oppen, F. Topological superconducting phase in helical shiba chains. Phys. Rev. B 88, 155420 (2013).
- (27) Schecter, M., Flensberg, K., Christensen, M. H., Andersen, B. M. & Paaske, J. Self-organized topological superconductivity in a Yu-Shiba-Rusinov chain. Phys. Rev. B 93, 140503 (2016).
- (28) Fulga, I. C., Haim, A., Akhmerov, A. R. & Oreg, Y. Adaptive tuning of Majorana fermions in a quantum dot chain. New Journal of Physics 15, 045020 (2013).
- (29) Su, Z. et al. Andreev molecules in semiconductor nanowire double quantum dots. Nature Communications 8, 585 (2017).
- (30) Sherman, D. et al. Normal, superconducting and topological regimes of hybrid double quantum dots. Nat Nano 12, 212–217 (2017).