
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39256344
52066
10.1038/s41467-024-52066-2
Article
Robust poor man’s Majorana zero modes using Yu-Shiba-Rusinov states
http://orcid.org/0000-0002-4250-6566
Zatelli Francesco 1
van Driel David 1
http://orcid.org/0000-0002-7158-1469
Xu Di 1
Wang Guanzhong 1
Liu Chun-Xiao 1
Bordin Alberto 1
Roovers Bart 1
http://orcid.org/0000-0002-7611-8115
Mazur Grzegorz P. 1
http://orcid.org/0000-0001-7461-5428
van Loo Nick 1
Wolff Jan C. 1
http://orcid.org/0000-0003-0593-6062
Bozkurt A. Mert 1
Badawy Ghada 2
Gazibegovic Sasa 2
http://orcid.org/0000-0002-8264-6862
Bakkers Erik P. A. M. 2
Wimmer Michael 1
http://orcid.org/0000-0002-7046-8181
Kouwenhoven Leo P. l.p.kouwenhoven@tudelft.nl

1
http://orcid.org/0000-0003-0966-2564
Dvir Tom 1
1 https://ror.org/02e2c7k09 grid.5292.c 0000 0001 2097 4740 QuTech and Kavli Institute of NanoScience, Delft University of Technology, Delft, The Netherlands
2 https://ror.org/02c2kyt77 grid.6852.9 0000 0004 0398 8763 Department of Applied Physics, Eindhoven University of Technology, Eindhoven, The Netherlands
11 9 2024
11 9 2024
2024
15 79332 7 2024
21 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Kitaev chains in quantum dot-superconductor arrays are a promising platform for the realization of topological superconductivity. As recently demonstrated, even a two-site chain can host Majorana zero modes known as “poor man’s Majorana”. Harnessing the potential of these states for quantum information processing, however, requires increasing their robustness to external perturbations. Here, we form a two-site Kitaev chain using Yu-Shiba-Rusinov states in proximitized quantum dots. By deterministically tuning the hybridization between the quantum dots and the superconductor, we observe poor man’s Majorana states with a gap larger than 70 μeV. The sensitivity to charge fluctuations is also greatly reduced compared to Kitaev chains made with non-proximitized dots. The systematic control and improved energy scales of poor man’s Majorana states realized with Yu-Shiba-Rusinov states will benefit the realization of longer Kitaev chains, parity qubits, and the demonstration of non-Abelian physics.

A Kitaev chain formed by two quantum dots coupled via a superconductor support the so-called poor man’s Majorana bound states. Here, the authors form a minimal Kitaev chain using Yu-Shiba-Rusinov states and show that the resulting bound states are more robust than in the case of unproximitized quantum dots.

Subject terms

Condensed-matter physics
Quantum information
501100003246 Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Netherlands Organisation for Scientific Research) Microsoft Corporation Station Qissue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Kitaev chains based on quantum dots (QDs) coupled via a hybrid semiconductor-superconductor heterostructure are a promising avenue for the creation of Majorana bound states1,2. Even a minimal chain, consisting of only two QDs, supports fine-tuned Majorana zero modes known as “poor man’s Majoranas” (PMMs)3. These PMM states do not benefit from topological protection, but already exhibit robustness to local perturbations and quadratic protection from global fluctuations in the chemical potential3. Moreover, PMM states obey non-Abelian exchange statistics, thus providing a favorable platform for braiding and Majorana-based qubit experiments in the near future4–7.

One approach to form PMM states is to couple spin-polarized QDs via both crossed Andreev reflection (CAR) and elastic co-tunneling (ECT)2,3,8,9, as recently demonstrated in InSb/Al nanowires10. Realizing such a system is challenging since the coupling between the QDs and the superconducting reservoir that enables ECT and CAR needs to be balanced. Too small a coupling results in negligible interaction between the QDs, while too large a coupling prevents their charge from being a good quantum number. While ref. 10 provided the necessary proof of concept, the gap separating PMMs from excited states was not much higher than the electron temperature. Furthermore, appreciable noise in the measured spectra indicated sensitivity to charge fluctuations, setting dim prospects for a qubit realization. Finally, a strong coherent coupling between the QDs was observed only in a limited number of devices, posing a significant yield challenge in realizing longer Kitaev chains and coupling them.

In this work, we show that it is possible to overcome these challenges by using spin-polarized Yu-Shiba-Rusinov (YSR) states11–20 instead of non-proximitized QDs as the sites of a Kitaev chain21–30. We systematically control the hybridization between two QDs and a superconductor to form strongly coupled YSR states, which we use to realize a robust two-site Kitaev chain21,31. We demonstrate that we can tune the system to a sweet spot hosting PMMs without relying on microscopic, device-specific details. Importantly, the PMMs do not hybridize despite their significant wavefunction overlap in the superconductor originating from its strong coupling to the QDs21,31. We show that the sensitivity of the resulting PMMs to charge fluctuations affecting the QDs decreases by two orders of magnitude compared to a QD-based Kitaev chain10. In addition, we measure a gap between ground and excited states of Egap = 76 μeV, a threefold increase compared to our previous report10.

We present results obtained using a hybrid semiconducting InSb nanowire partly covered with a superconducting Al film. However, our approach does not rely on the details of the material platform. Complementary results are also reported in a parallel study on a 2D hybrid platform using an InAsSb/Al two-dimensional electron gas32, demonstrating the wider applicability of our approach.

Results

Fabrication and setup

Figure 1a shows a scanning electron microscope image of the reported device. It consists of an InSb nanowire (green) placed on top of a series of bottom gates33. The middle part of the nanowire is covered by a thin Al shell (blue), forming a superconductor-semiconductor hybrid whose electrochemical potential is controlled by a plunger gate (VH). On both sides of the hybrid segment, QDs are formed in the nanowire using three gates each. The electrochemical potential of each QD is controlled by a plunger gate (VLP and VRP for the left and right QDs, respectively), and the couplings between the QDs and the hybrid are controlled by tunnel gates (VLT and VRT for the left and right QDs, respectively). A normal lead is attached to each QD, separated by another gate-defined tunnel barrier. The superconducting lead is kept grounded at all times. In addition, the two normal leads are connected to off-chip multiplexed resonators for fast RF reflectometry measurements34, using the setup described in ref. 35. Each lead is voltage-biased independently with respect to the grounded Al, with voltages VL and VR on the left and right leads, respectively. The currents (IL and IR on the left and right leads), the local conductances (GLL = dIL/dVL, GRR = dIR/dVR), and the non-local conductances (GRL = dIR/dVL, GLR = dIL/dVR) are measured simultaneously. When the full conductance matrix is measured, we correct for line resistance as described in the Methods section. Further fabrication and setup details can be found in our previous publications36–38. The experiment is conducted in a dilution refrigerator with a base temperature of 30 mK. A magnetic field of B = 150 mT is applied along the nanowire axis, inducing a Zeeman splitting of ~200 μeV in the QDs (Supplementary Fig. 1).Fig. 1 YSR states formed by hybridizing QDs with an ABS.

a False-colored scanning electron micrograph and measurement circuit of the device. Scale bar is 200 nm. b Illustration of the QD-ABS-QD model. Two QDs are coupled to the same ABS via spin-conserving and spin-flipping tunneling. QD and YSR regimes are defined by the relative energy scales shown below the sketch. c Numerical CSD using zero-bias conductance of a QD coupled to an ABS. “E''/"O” indicate even/odd occupation of the QD. d, e Numerical conductance spectrum of a QD coupled to an ABS as a function of QD chemical potential when ABS is far (d) and close to its energy minimum (e). In e, green/pink arrows indicate avoided crossings due to spin-conserving/spin-flipping tunneling. f, h Same as c–e for the left QD coupled to an ABS, experimentally measured using local conductance. Gray/blue ticks in f mark VH values at which g/h are measured. The right QD is off resonance. In h, we extract a lever arm α ≈ 0.05e by fitting the spectrum with the black dotted line at the zero-energy crossing. This procedure overestimates the actual lever arm since it does not account for capacitance to the normal lead. For comparison, we plot the profile of a typical Coulomb diamond in our device with gray dotted lines (Supplementary Fig. 3a). The VLP range in g is shifted because of a gate jump affecting the left QD. i–k, Same as f–h, but for the right QD. The extracted lever arm is α ≈ 0.04e.

YSR states in quantum dots

We model our system using a three-site model in which the hybrid is considered as a single Andreev bound state (ABS) in the atomic limit15,16 tunnel-coupled to two QDs (see schematics in Fig. 1b)21. The QDs have charging energy U, Zeeman splitting EZ, and chemical potentials μL, μR for the left and right QD, respectively. The ABS has an induced gap Γ, which in the atomic limit15,16 can be identified with its tunnel coupling to the bulk superconductor. Its charging energy is negligible compared to the induced gap Γ due to the screening of the grounded Al film. It also has Zeeman splitting EZH, which is smaller than that of the QDs due to metallization of the ABS39. We ensure EZH < Γ so that the ground state of the ABS is always a BCS singlet. The electron-hole composition of the ABS depends on its chemical potential μH.

In our model, the QDs are coupled to the ABS by spin-conserving and spin-flipping tunneling due to spin-orbit interaction, with amplitudes w and wSO, respectively21. The hybridization between QDs and the ABS becomes significant when w, wSO ~ Γ. As a consequence, the QDs become proximitized and form YSR states20,31, which we distinguish from ABSs because of their large charging energy15–19,40. We refer readers to Supplementary Notes 1 and 2 and the parallel work of ref. 31 for a more detailed discussion of the distinction between ABSs and YSR states, and for theory models to describe the strong coupling regime investigated in this work.

To understand the nature of the YSR states formed in proximitized QDs, we first examine the coupling between a single QD and an ABS in the hybrid segment, following ref. 20. Fig. 1c shows the theoretical zero-bias conductance of a QD-ABS charge stability diagram (CSD), while the second QD is off-resonance. The two vertical features indicate parity transitions of the system, largely corresponding to the consecutive filling of a single orbital of the QD by two electrons. The S-shaped conductance features result from renormalization of the QD energy via hybridization with the ABS. The ABS reaches its energy minimum at charge neutrality, i.e., μH = 0, where its excitation is equal-parts electron and hole. The hybridization with the QD is also maximal here due to their minimal energy separation, evident in the enhanced zero-bias conductance as local Andreev reflection becomes stronger20,31. The QD spectra in panels d and e further reveal its hybridization with the ABS. When the latter is away from charge neutrality, the QD spectrum as a function of μL (Fig. 1d) exhibits straight features reminiscent of Coulomb diamonds. As the ABS approaches μH = 0 (Fig. 1e), such features evolve into an eye-shaped spectrum typical of YSR states19,41–44. The arrows indicate avoided crossings in the excited states of the spectrum produced by spin-conserving (green) and spin-flipping (pink) tunneling between the ABS and the QD levels31.

The hybrid segment of our device features multiple discrete ABSs well-separated from each other in VH (Supplementary Fig. 2). We operate in a VH range containing a single ABS level. In Fig. 1f, we show the zero-bias conductance measured with the left lead as a function of VLP and VH. VRP is fixed to keep the right QD off-resonance. We observe the QD-ABS charge stability diagram features described in panel c, indicating the presence of an ABS reaching its energy minimum at VH ≈ 337 mV. The QD-ABS hybridization is further confirmed by the QD spectrum being in agreement with the model when the ABS is away from (Fig. 1g) and at its energy minimum (Fig. 1h).

The YSR zero-energy excitations in Fig. 1h are our building blocks of a Kitaev chain. Compared to a non-proximitized QD zero-energy crossing, these YSR crossings have noticeably weaker energy dispersion as a function of gate (see dashed lines in panel h). To quantify this observation, we can estimate the lever arm of the YSR excitation at charge degeneracy using its gate-dispersion slope: α ≡ ∂E/∂VLP. In contrast to the above-gap QD lever arm of α ≈ 0.4e in our devices (Supplementary Fig. 3a), this subgap lever arm reduces to α ≈ 0.2e (Supplementary Fig. 3b) when the ABS is detuned from its charge neutrality. Tuning the ABS to its energy minimum further reduces the lever arm to α ≈ 0.05e (Fig. 1h). This signals a strong reduction in the effective charge of the fermionic excitation, attributable to charging energy renormalization, QD-ABS hybridization, and electron-hole superposition45,46, as detailed in Supplementary Fig. 4.

The right QD shows similar behavior when coupled to the same ABS (Fig. 1i–k). However, when the ABS is away from charge neutrality, the subgap conductance is not significantly suppressed (Fig. 1i) and the QD spectrum still shows typical YSR features (Fig. 1h), albeit with a small superconducting coupling. This could be due to residual proximity resulting from direct coupling between the QD and the superconducting film since VH does not affect it appreciably (Supplementary Fig. 5).

Coupled YSR states

Recent theoretical8,21 and experimental9,10,47 works have shown that an ABS can mediate elastic co-tunneling (ECT) and crossed Andreev reflection (CAR) between QDs. These two processes implement the hopping and pairing terms of the original Kitaev chain model1,2. To form PMM states in a two-site chain, the amplitudes of both terms must be equal3. Such control can be achieved by tuning the electrochemical potential of the ABS in the hybrid nanowire8–10,21.

Similar effective ECT and CAR couplings, with respective amplitudes t and Δ, also emerge between YSR states (Fig. 2a) formed by strongly coupling QDs to the same ABS31,48,49. To observe them, we turn to the CSD of two such YSR states, akin to those explored in refs. 50,51. In Fig. 2b, we show the zero-bias conductance matrix measured as a function of VRP and VLP when the ABS is tuned away from its energy minimum. All elements of the conductance matrix show prominent resonances arising from the two charge transitions of each QD. The type of avoided crossings observed in the CSD serves as an indication of the coupled-YSR system’s ground state10,21. Avoided crossings along a negative diagonal, as seen in all four resonances in Fig. 2b, show hybridization between states with the same total charge. This is the ground state of the system when t > Δ. This observation is further confirmed by the negative non-local conductance that is characteristic of ECT10.Fig. 2 ECT and CAR coupling between YSR states.

a Illustration of the effective two-site system of YSR states coupled via ECT (t) and CAR (Δ). b–d Conductance matrices of CSDs measured at different VH. By tuning the electrochemical potential of the ABS in the hybrid, it is possible to continuously vary t and Δ. In b, avoided crossings along the antidiagonal are observed, indicating t > Δ. The opposite regime is shown in d, with all the avoided crossings along the diagonal, indicating t < Δ. The crossover between these two regimes is shown in c, where two crossings indicate t ≈ Δ.

Increasing the value of VH leads to a change in the CSD, shown in Fig. 2c. Here, the bottom-left and the top-right resonances show avoided crossings along a positive diagonal associated with Δ > t. In the top-left and bottom-right, the resonance lines cross each other, indicating a PMM sweet spot, with t ≈ Δ. Finally, bringing the ABS close to its energy minimum by a further increase of VH tunes all of the quadrants of the CSD to the Δ > t regime (Fig. 2d). The theoretical model reproduces the observed evolution of the CSDs (Supplementary Fig. 6a–c).

Gate control of CAR and ECT

ABS-mediated ECT and CAR couplings between QDs are controlled by the chemical potential of ABS8,9. At the energy minimum of the ABS, its excitation is equally electron- and hole-like, both parts interfering constructively to enhance CAR and destructively to quench ECT. Finite values of μH lead to an imbalance between the electron and hole parts of the ABS, decreasing the value of the CAR coupling while enhancing ECT. This control over the ECT and CAR amplitudes guarantees a PMM sweet spot when the two QDs are coupled via a single ABS8,9,21. To demonstrate that this description can be extended to YSR states coupled via an ABS, we study how the couplings t and Δ vary as the electrochemical potential of the ABS changes.

The magnitudes of t and Δ can be extracted by measuring the excitation spectrum of the system at μL = μR = 0, when spectral splitting is determined by the couplings alone. We limit the discussion here to the bottom-right crossing of Fig. 2, where the QDs’ spins are opposite. Other crossings exhibit qualitatively similar behavior (Supplementary Fig. 7). At each fixed value of VH, we measure a CSD (Fig. 3a) and set the QD gates to the center of an avoided crossing. There, the subgap spectrum exhibits two sets of electron-hole symmetric peaks at energies ∣t − Δ∣ and t + Δ (Fig. 3b), as detailed in Methods. By fitting the measured spectrum with two pairs of Gaussians symmetric around VL,R = 0, we extract the energy of the excited states and calculate t and Δ. We repeat this procedure for different values of VH and collect the spectra in Fig. 3c, where each line was measured as described above.Fig. 3 ECT and CAR coupling as a function of the ABS chemical potential.

a CSD measured using RF reflectometry. The center of the avoided crossing, marked with a gray cross, is identified as μL = μR = 0. b Spectrum measured on both sides of the device at the center of the avoided crossing shown in a. The inner peaks correspond to excited states at ∣t − Δ∣, while the outer ones correspond to excited states at t + Δ. By summing and subtracting these energies, it is possible to extract t and Δ for each VH value. c Spectrum as a function of VH while μL = μR = 0. Each line of the spectrum was measured as described above. The blue tick indicates the VH value at which a and b were measured. d ECT and CAR couplings extracted from c. e Zoom-in on the zero-energy crossings highlighted in c. The white dotted lines are fits to extract the slope of the linear splitting. The data reported here was collected for the bottom-right charge transition shown in Fig. 2. More details about data processing and data for all charge degeneracies are reported in Supplementary Fig. 7.

Figure 3 d shows the extracted values of ECT and CAR amplitudes. At VH ≈ 322 mV and VH ≈ 355 mV, the outer peaks almost merge with the inner ones, indicating ∣t − Δ∣ ≈ t + Δ. Inspecting the corresponding CSDs (see Fig. 2b for an example, all the CSDs are available in the online repository), we can see that ECT is dominant and, therefore, t ≫ Δ. Upon varying VH, the two peaks split into four well-separated peaks, signaling that CAR is increasing although ECT still prevails. At VH ≈ 333 mV and 345 mV, the two inner peaks merge into a single zero-bias peak. These are two sweet spots where t ≈ Δ. In between them, CAR dominates, as the CSDs can confirm (Fig. 2d). Around VH ≈ 338 mV, the outer peaks merge again, this time indicating Δ ≫ t. This feature is the peak in Δ and dip in t seen in Fig. 3d, because of the aforementioned interference effects8,9. Finally, Supplementary Fig. 6d shows that the same spectral features are reproduced by the theoretical model.

The zero-energy crossings in Fig. 3e allow us to characterize the robustness of the PMM sweet spots to charge fluctuations affecting the ABS and, consequently, t and Δ, causing a splitting of the zero-energy states. Indeed, we observe that each zero-energy state splits with a linear dependence on VH as predicted by theory, indicating the lack of protection against deviations from the condition t = Δ known for a two-site Kitaev chain3. We extract a slope of ∂E∂VH≈7μeV/mV and 3 μeV/mV for the two crossings, comparable to what we measured in ref. 10.

Majorana sweet spot

Finally, we measure the spectrum and stability of the PMM states against perturbation of the QDs. We tune our device to the sweet spot at VH = 333.5 mV, where t ≈ Δ. In Fig. 4a, we show the full conductance matrix of the CSD at the sweet spot. As expected, the resonance lines cross each other and non-local conductance alternates between positive and negative values. The spectrum measured at the sweet spot in the center of the CSD (Fig. 4b) shows, on both sides of the device, a zero-bias conductance peak clearly separated from the excited states. As a result of the much stronger coupling between YSR states, the first excited states reside at energies of Egap ≈ 76 μeV, three times larger than in our previous report10 and significantly above the electron temperature.Fig. 4 Poor man’s Majorana sweet spot for coupled YSR states.

a Conductance matrix of a crossing in the CSD, when t ≈ Δ. b Spectrum of both sides measured at the sweet spot in the center of the crossing of panel a. c Conductance matrix of the spectrum as a function of VRP (along the gray dotted line in a). The right QD is detuned across the sweet spot, while the left one is kept on resonance. d Conductance matrix of the spectrum measured detuning both QDs simultaneously along the antidiagonal along the pink dashed line in a. e, f GRR as a function of the simultaneous detunings of both sites, δVLP and δVRP, away from the sweet spot, using QDs (e, replotted from ref. 10) and YSR states (f, from the same dataset presented in d). For comparison, we have plotted the two gate dispersions with the same scale in both plots. The green dotted lines correspond to the expected energy splitting of a PMM realized with QDs, the blue dashed lines to that of a PMM realized with YSR states, as detailed in the Methods section.

One of the hallmarks of PMM states is their stability against local perturbations of the QDs’ electrochemical potential. In Fig. 4c, we measure the spectrum of the QDs varying VRP while keeping the left QD on resonance. The observed zero-bias conductance peak persists until the next charge degeneracies of the detuned QD. The range of this region is determined by its charging energy and level spacing. For our device, this corresponds to a gate range of  ~ 10 mV (Fig. 2c). Tuning VLP while keeping the right QD on resonance (Supplementary Fig. 8b) shows the same qualitative features. The same behavior can be reproducibly observed with other QD and ABS orbitals, as shown in Supplementary Fig. 9.

If the electrochemical potentials of both QDs are detuned from the sweet spot, PMM states are expected to split quadratically. This is verified in Fig. 4d, where we measure the spectrum while detuning both QDs along the antidiagonal path shown in Fig. 4a. The spectrum taken along the diagonal path is shown in Supplementary Fig. 8c. Numerical simulations reproducing these measurements are reported in Supplementary Fig. 6e–i. For comparison, we plot the quadratic dispersion measurements of the QD-based vs YSR-based PMMs side by side in panels e and f. The curvature of the energy-gate dispersion close to the sweet spot is directly proportional to the dephasing rate resulting from charge noise affecting both QDs52–55 and is therefore a measure of the Majorana states’ robustness against it. Comparing the overlaid curves in Fig. 4e, f, we find the gate dispersion curvature reported in this work to be a factor of  ~ 150 lower than the non-proximitized case (details in Methods). This striking reduction can be fully explained by the decreased lever arm of the YSR states, α, and the increased Egap, since the curvature expected from the theory model is α2/Egap (see “Methods” section).

To illustrate the effect of reduced charge dispersion on the coherence of a potential poor man’s Majorana qubit, we calculate the dephasing rates using realistic charge noise estimations and data presented above (see Supplementary Note 3). While perturbation of QDs’ potentials is expected to be the dominant mechanism of energy splitting and thus dephasing for non-proximitized PMMs, the drastic reduction in gate dispersion of YSR-based PMMs makes it negligible compared to dephasing caused by deviations from the t = Δ condition. The estimated T2*~10ns is now limited by noise affecting the couplings t and Δ, nearly an order of magnitude higher than the charge-noise-limited T2* extracted from data in our previous report10. Importantly, the expected dephasing time is now also much longer than the adiabatic limit ℏ/Egap ~ 10 ps, which sets an upper bound on how fast Majorana states can be manipulated without populating the excited states5,52.

Discussion

In conclusion, we have realized a two-site Kitaev chain hosting PMM states by coupling YSR states via ECT and CAR. We have formed YSR states by coupling QDs to a common ABS. We have shown how controlling their hybridization with the ABS allows us to deterministically achieve the strong coupling regime between the YSR states. By tuning the electrochemical potential of the ABS, we can systematically control the ECT and CAR couplings in accordance with a single-ABS theoretical model. We have demonstrated that, despite the strong hybridization between QDs and ABS, it is possible to tune the system to a regime consistent with PMM states. We emphasize that the strong QD-ABS coupling causes the PMMs’ wavefunction to substantially overlap in the hybrid. Nevertheless, such overlap does not cause hybridization between the PMMs themselves21,31. The resulting PMM states have two significant improvements over those in a non-proximitized QD chain. First, the stronger coupling between the YSR states triples the gap between the PMM and excited states, protecting the former from finite temperature excitations53 and enabling faster adiabatic operations5,52. Second, the reduced charge dispersion of the YSR states enhances the robustness of the PMMs against charge noise affecting both QDs by two orders of magnitude. Thanks to these, even a two-site Kitaev chain realized using YSR states should suffice for a prototypical Majorana qubit and verification of non-Abelian properties with fusion and braiding experiments. Despite the present lack of protection against tunnel-coupling noise, the expected coherence of Majorana qubits made from PMM states is increased by close to an order of magnitude compared to the first report10. In the future, increasing the number of sites can mitigate noise affecting the tunnel-coupling rates. Estimations using parameters of the YSR-based PMMs suggest that a Majorana qubit realized with Kitaev chains as short as 3 to 5 sites could already achieve dephasing times comparable to those predicted for continuous nanowires2,53,56. Finally, the systematic control of QD-ABS and subsequently YSR-YSR hybridization demonstrated in this work will be useful for Andreev spin qubits57,58, long-distance coupling of spin qubits via a superconductor59,60, and analog quantum simulations of Fermi-Hubbard systems61–63 with superconductivity.

Methods

Device fabrication

The InSb/Al hybrid nanowire device presented in this work was fabricated using the shadow-wall lithography technique36,37.

A substrate is patterned with Ti/Pd gates. 10 nm of AlOx and 10nm of HfOx are deposited by ALD as gate dielectric. HSQ shadow walls are then patterned. Nanowires are deposited and pushed next to the shadow walls using an optical micro-manipulator. An 8nm Al shell is deposited at alternating angles of 15∘ and 45∘ with respect to the substrate, followed by a capping layer of 20 nm of AlOx. Normal leads in ohmic contact with the nanowire are fabricated by Ar milling and evaporation of Cr/Au.

Additional details about the substrate fabrication and the Al deposition are described in ref. 38.

Transport measurements and data processing

The measurements are done in a dilution refrigerator with a base temperature of 30 mK. A magnetic field of 150 mT is applied along the nanowire axis except in Supplementary Fig. 1c, d.

The three-terminal setup used to measure the device is illustrated in Fig. 1a. The SC lead is always kept grounded. The two normal leads can be voltage-biased independently. When a bias is applied to one side, the other one is kept grounded. The currents IL and IR are measured separately. Digital multimeters and lock-in amplifiers are used to read the voltage outputs of the current meters. AC excitations of 5 μV RMS are applied on each side with different frequencies (39 Hz on the left and 29 Hz on the right), except for Supplementary Fig. 2 and Supplementary Fig. 5d–f, where excitations of 10 μV RMS were used.

Off-chip multiplexed resonators34 connected to the two normal leads are used for fast RF reflectometry measurements. This measurement scheme was employed to speed up the tune-up of the device, for Fig. 3, and for Supplementary Fig. 7, as explained in Fig. 3a. Additional details about the reflectometry setup are described in ref. 35.

Spinless PMM spectrum

The low-energy spectrum of the coupled YSR system arising from the model described in Supplementary Note 1 can be reduced to that of a spinless PMM model3,6,8,21,31.

In this simplified model, the non negative energy eigenvalues are31 E±=μ+2+μ−2+t2+Δ2±2(μ+2+Δ2)(μ−2+t2),

where μ±=μL±μR2, μL and μR are the chemical potentials of the two sites, and t, Δ are the couplings. The spectrum is symmetric around zero energy because of particle-hole symmetry. If μL = μR = 0, the eigenvalues reduce to E+ = t + Δ and E− = ∣t − Δ∣. These are the excitation energies that we use to extract t and Δ from the measured spectra in Fig. 3.

The gate dispersion lines plotted in Fig. 4e,f are calculated using the lowest excitation energy. When t = Δ we obtain2 E−=α2δV2+2t2−2t2(t2+α2δV2),

where α is the lever arm converting voltage to chemical potential and δV is the simultaneous detuning of each gate away from the sweet spot. For the PMM realized with QDs, we have used α = 0.33e and t = 12 μeV10. For the PMM realized with YSR states we have used α = 0.05e (Supplementary Fig. 1) and t = 38 μeV (Fig. 4b). Finally, the quadratic splitting of PMM states can be derived by expanding the expression above for small δV, which gives3 E−≈α22tδV2.

Series resistance correction

The effect of series resistance in the fridge line and other parts of the circuit on transport measurements of a three-terminal device of our type is described in ref. 64, as well as how to correct it. In our setup, the resistance of the voltage source is 100 Ω and that of the current meter is 200 Ω. Additional series resistance comes from the fridge lines and the ohmic contacts. Correcting for the voltage fall over these resistors in series to the device requires knowledge of the exact resistance values, which we presently cannot obtain before the next sample exchange. Therefore, we make use of the fact that the bulk Al superconducting gap in local conductance measurements should stay constant across all gate values to arrive at an estimation of the total series resistance in each fridge line.

In Supplementary Fig. 10a, c, it is possible to observe that the superconducting gap appears enlarged when there is finite subgap conductance, indicating unaccounted-for series resistance. To estimate it, we correct the GLL measurement shown in Supplementary Fig. 10a for different trial values of the series resistance. We extract the energies of the coherence peaks and calculate their variance. Finally, we choose the value that minimizes the variance as the optimal series resistance. We find an optimal value of 3.65 kΩ, in addition to the resistance of the voltage source and current meter. The corrected measurement is shown in Supplementary Fig. 10b. All the data processing steps are available in the online repository. In Supplementary Fig. 10c, d, we show that the same series resistance value also corrects an analogous measurement on the right side. In the regimes relevant to this work, the voltage divider effect caused by the series resistance does not strongly affect local conductance. On the other hand, it can be more pronounced for non-local conductance, as discussed in Supplementary Fig. 10e, f. Using the extracted resistance value, we apply the correction method described in ref. 64 to the data presented in Figs. 2–4 and Supplementary Figs. 2, 5d, 7, 8, 9d–k, and 10b, d, f, where the full conductance matrix is measured.

Device tune-up

The device is controlled with seven bottom gates (Fig. 1a). To perform tunnel spectroscopy of the hybrid segment in Supplementary Figs. 2 and 5, we use VLT and VRT to form a tunnel barrier and apply a large positive voltage to VLP, VRP, and the outer gates. VH is used to control the electrochemical potential in the hybrid segment, where discrete states are confined because of the tunnel barriers. QDs are formed by reducing the voltage on the outer gates. This creates additional tunnel barriers next to the normal leads, thus shaping a confining potential. VLP and VRP control the electrochemical potential of the QDs. Finally, the coupling between the QDs and the ABS can be controlled using VLT and VRT. By increasing their voltage, we can achieve the strong coupling regime between the QDs and the ABS in the hybrid as explained in Fig. 1. For all the measurements presented in the main text, the gates defining the tunnel barriers were maintained at the same value.

Supplementary information

Supplementary Information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52066-2.

Acknowledgements

This work has been supported by the Dutch Organization for Scientific Research (NWO) and Microsoft Corporation Station Q. We thank John M. Hornibrook and David J. Reilly for providing the frequency multiplexing chips. We thank Martin Leijnse, Ville Maisi, Pasquale Scarlino, Athanasios Tsintzis, Rubén Seoane Souto, Karsten Flensberg, Srijit Goswami, Sebastiaan L. D. ten Haaf, Qingzhen Wang, Ivan Kulesh, and Yining Zhang for helpful discussions.

Author contributions

D.v.D., F.Z., A.B., G.P.M., N.v.L., and J.C.W. fabricated the device. F.Z., D.v.D., D.X., B.R., and T.D. performed the electrical measurements. T.D. and G.W. designed the experiment. F.Z., D.v.D., G.W., and T.D. analyzed the data. F.Z., D.v.D., T.D., G.W., and L.P.K. prepared the manuscript with input from all authors. T.D. and L.P.K. supervised the project. C.X.L. and A.M.B. developed the theoretical model with input from M.W.; G.B., S.G., and E.P.A.M.B. performed InSb nanowire growth.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The raw data presented in this work, together with the code used to process and plot it, is available at https://zenodo.org/records/10013728.

Code availability

The code used for theory calculations is available at https://zenodo.org/records/10013728.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Francesco Zatelli, David van Driel, Di Xu, Guanzhong Wang.
==== Refs
References

1. Kitaev AY Unpaired Majorana fermions in quantum wires Phys.-Uspekhi 2001 44 131 10.1070/1063-7869/44/10S/S29
Kitaev, A. Y. Unpaired Majorana fermions in quantum wires. Phys.-Uspekhi 44, 131 (2001).10.1070/1063-7869/44/10S/S29
2. Sau JD Sarma SD Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array Nat. Commun. 2012 3 964 10.1038/ncomms1966 22805571
Sau, J. D. & Sarma, S. D. Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array. Nat. Commun. 3, 964 (2012).22805571 10.1038/ncomms1966
3. Leijnse M Flensberg K Parity qubits and poor man’s Majorana bound states in double quantum dots Phys. Rev. B 2012 86 134528 10.1103/PhysRevB.86.134528
Leijnse, M. & Flensberg, K. Parity qubits and poor man’s Majorana bound states in double quantum dots. Phys. Rev. B 86, 134528 (2012).10.1103/PhysRevB.86.134528
4. Liu C-X Pan H Setiawan F Wimmer M Sau JD Fusion protocol for Majorana modes in coupled quantum dots Phys. Rev. B 2023 108 085437 10.1103/PhysRevB.108.085437
Liu, C.-X., Pan, H., Setiawan, F., Wimmer, M. & Sau, J. D. Fusion protocol for Majorana modes in coupled quantum dots. Phys. Rev. B 108, 085437 (2023).10.1103/PhysRevB.108.085437
5. Boross P Páyi A Braiding-based quantum control of a Majorana qubit built from quantum dots Phys. Rev. B 2024 109 125410 10.1103/PhysRevB.109.125410
Boross, P. & Páyi, A. Braiding-based quantum control of a Majorana qubit built from quantum dots. Phys. Rev. B 109, 125410 (2024).10.1103/PhysRevB.109.125410
6. Tsintzis A Souto RS Flensberg K Danon J Leijnse M Majorana qubits and non-abelian physics in quantum dot–based minimal Kitaev chains PRX Quantum 2024 5 010323 10.1103/PRXQuantum.5.010323
Tsintzis, A., Souto, R. S., Flensberg, K., Danon, J. & Leijnse, M. Majorana qubits and non-abelian physics in quantum dot–based minimal Kitaev chains. PRX Quantum 5, 010323 (2024).10.1103/PRXQuantum.5.010323
7. Pino DM Souto RS Aguado R Minimal Kitaev-transmon qubit based on double quantum dots Phys. Rev. B 2024 109 075101 10.1103/PhysRevB.109.075101
Pino, D. M., Souto, R. S. & Aguado, R. Minimal Kitaev-transmon qubit based on double quantum dots. Phys. Rev. B 109, 075101 (2024).10.1103/PhysRevB.109.075101
8. Liu C-X Wang G Dvir T Wimmer M Tunable superconducting coupling of quantum dots via Andreev bound states in semiconductor-superconductor nanowires Phys. Rev. Lett. 2022 129 267701 10.1103/PhysRevLett.129.267701 36608192
Liu, C.-X., Wang, G., Dvir, T. & Wimmer, M. Tunable superconducting coupling of quantum dots via Andreev bound states in semiconductor-superconductor nanowires. Phys. Rev. Lett. 129, 267701 (2022).36608192 10.1103/PhysRevLett.129.267701
9. Bordin A Tunable crossed Andreev reflection and elastic cotunneling in hybrid nanowires Phys. Rev. X 2023 13 031031
Bordin, A. et al. Tunable crossed Andreev reflection and elastic cotunneling in hybrid nanowires. Phys. Rev. X 13, 031031 (2023).
10. Dvir T Realization of a minimal Kitaev chain in coupled quantum dots Nature 2023 614 445 450 10.1038/s41586-022-05585-1 36792741
Dvir, T. et al. Realization of a minimal Kitaev chain in coupled quantum dots. Nature 614, 445–450 (2023).36792741 10.1038/s41586-022-05585-1
11. Yu L Bound state in superconductors with paramagnetic impurities Acta Phys. Sin. 1965 21 75 91 10.7498/aps.21.75
Yu, L. Bound state in superconductors with paramagnetic impurities. Acta Phys. Sin. 21, 75–91 (1965).10.7498/aps.21.75
12. Shiba H Classical spins in superconductors Prog. Theor. Phys. 1968 40 435 451 10.1143/PTP.40.435
Shiba, H. Classical spins in superconductors. Prog. Theor. Phys. 40, 435–451 (1968).10.1143/PTP.40.435
13. Rusinov, A. Superconductivity near a paramagnetic impurity. JETP Lett. 9, 85–87 (1969).
14. Buitelaar MR Nussbaumer T Schönenberger C Quantum dot in the Kondo regime coupled to superconductors Phys. Rev. Lett. 2002 89 256801 10.1103/PhysRevLett.89.256801 12484909
Buitelaar, M. R., Nussbaumer, T. & Schönenberger, C. Quantum dot in the Kondo regime coupled to superconductors. Phys. Rev. Lett. 89, 256801 (2002).12484909 10.1103/PhysRevLett.89.256801
15. Bauer J Oguri A Hewson AC Spectral properties of locay correlated electrons in a Bardeen-Cooper-Schrieffer superconductor J. Phys.: Condens. Matter 2007 19 486211
Bauer, J., Oguri, A. & Hewson, A. C. Spectral properties of locay correlated electrons in a Bardeen-Cooper-Schrieffer superconductor. J. Phys.: Condens. Matter 19, 486211 (2007).
16. Meng T Florens S Simon P Self-consistent description of Andreev bound states in Josephson quantum dot devices Phys. Rev. B 2009 79 224521 10.1103/PhysRevB.79.224521
Meng, T., Florens, S. & Simon, P. Self-consistent description of Andreev bound states in Josephson quantum dot devices. Phys. Rev. B 79, 224521 (2009).10.1103/PhysRevB.79.224521
17. Grove-Rasmussen K Superconductivity-enhanced bias spectroscopy in carbon nanotube quantum dots Phys. Rev. B 2009 79 134518 10.1103/PhysRevB.79.134518
Grove-Rasmussen, K. et al. Superconductivity-enhanced bias spectroscopy in carbon nanotube quantum dots. Phys. Rev. B 79, 134518 (2009).10.1103/PhysRevB.79.134518
18. Chang W Manucharyan VE Jespersen TS Nygård J Marcus CM Tunneling spectroscopy of quasiparticle bound states in a spinful Josephson junction Phys. Rev. Lett. 2013 110 217005 10.1103/PhysRevLett.110.217005 23745916
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).23745916 10.1103/PhysRevLett.110.217005
19. Jeinggaard A Grove-Rasmussen K Madsen MH Nygård J Tuning Yu-Shiba-Rusinov states in a quantum dot Phys. Rev. B 2016 94 064520 10.1103/PhysRevB.94.064520
Jeinggaard, 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).10.1103/PhysRevB.94.064520
20. Grove-Rasmussen K Yu-Shiba-Rusinov screening of spins in double quantum dots Nat. Commun. 2018 9 2376 10.1038/s41467-018-04683-x 29915280
Grove-Rasmussen, K. et al. Yu-Shiba-Rusinov screening of spins in double quantum dots. Nat. Commun. 9, 2376 (2018).29915280 10.1038/s41467-018-04683-x
21. Tsintzis A Souto RS Leijnse M Creating and detecting poor man’s Majorana bound states in interacting quantum dots Phys. Rev. B 2022 106 201404 10.1103/PhysRevB.106.L201404
Tsintzis, A., Souto, R. S. & Leijnse, M. Creating and detecting poor man’s Majorana bound states in interacting quantum dots. Phys. Rev. B 106, 201404 (2022).10.1103/PhysRevB.106.L201404
22. Fulga IC Haim A Akhmerov AR Oreg Y Adaptive tuning of Majorana fermions in a quantum dot chain N. J. Phys. 2013 15 045020 10.1088/1367-2630/15/4/045020
Fulga, I. C., Haim, A., Akhmerov, A. R. & Oreg, Y. Adaptive tuning of Majorana fermions in a quantum dot chain. N. J. Phys. 15, 045020 (2013).10.1088/1367-2630/15/4/045020
23. Nadj-Perge S Drozdov IK Bernevig BA Yazdani A Proposal for realizing Majorana fermions in chains of magnetic atoms on a superconductor Phys. Rev. B 2013 88 020407 10.1103/PhysRevB.88.020407
Nadj-Perge, S., Drozdov, I. K., Bernevig, B. A. & Yazdani, A. Proposal for realizing Majorana fermions in chains of magnetic atoms on a superconductor. Phys. Rev. B 88, 020407 (2013).10.1103/PhysRevB.88.020407
24. Pientka F Peng Y Glazman L von Oppen F Topological superconducting phase and Majorana bound states in Shiba chains Phys. Scr. 2015 2015 014008 10.1088/0031-8949/2015/T164/014008
Pientka, F., Peng, Y., Glazman, L. & von Oppen, F. Topological superconducting phase and Majorana bound states in Shiba chains. Phys. Scr. 2015, 014008 (2015).10.1088/0031-8949/2015/T164/014008
25. Schecter M Flensberg K Christensen MH Andersen BM Paaske J Self-organized topological superconductivity in a Yu-Shiba-Rusinov chain Phys. Rev. B 2016 93 140503 10.1103/PhysRevB.93.140503
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).10.1103/PhysRevB.93.140503
26. Nadj-Perge S Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor Science 2014 346 602 607 10.1126/science.1259327 25278507
Nadj-Perge, S. et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602–607 (2014).25278507 10.1126/science.1259327
27. Ruby M End states and subgap structure in proximity-coupled chains of magnetic adatoms Phys. Rev. Lett. 2015 115 197204 10.1103/PhysRevLett.115.197204 26588411
Ruby, M. et al. End states and subgap structure in proximity-coupled chains of magnetic adatoms. Phys. Rev. Lett. 115, 197204 (2015).26588411 10.1103/PhysRevLett.115.197204
28. Jeon S Distinguishing a Majorana zero mode using spin-resolved measurements Science 2017 358 772 776 10.1126/science.aan3670 29025997
Jeon, S. et al. Distinguishing a Majorana zero mode using spin-resolved measurements. Science 358, 772–776 (2017).29025997 10.1126/science.aan3670
29. Schneider L Topological Shiba bands in artificial spin chains on superconductors Nat. Phys. 2021 17 943 948 10.1038/s41567-021-01234-y
Schneider, L. et al. Topological Shiba bands in artificial spin chains on superconductors. Nat. Phys. 17, 943–948 (2021).10.1038/s41567-021-01234-y
30. Schneider L Precursors of Majorana modes and their length-dependent energy osciations probed at both ends of atomic Shiba chains Nat. Nanotechnol. 2022 17 384 389 10.1038/s41565-022-01078-4 35256768
Schneider, L. et al. Precursors of Majorana modes and their length-dependent energy osciations probed at both ends of atomic Shiba chains. Nat. Nanotechnol. 17, 384–389 (2022).35256768 10.1038/s41565-022-01078-4
31. Liu C-X Enhancing the excitation gap of a quantum-dot-based Kitaev chain Commun. Phys. 2024 7 1 7 10.1038/s42005-024-01715-5
Liu, C.-X. et al. Enhancing the excitation gap of a quantum-dot-based Kitaev chain. Commun. Phys. 7, 1–7 (2024).10.1038/s42005-024-01715-5
32. ten Haaf SLD A two-site Kitaev chain in a two-dimensional electron gas Nature 2024 630 329 334 10.1038/s41586-024-07434-9 38867129
ten Haaf, S. L. D. et al. A two-site Kitaev chain in a two-dimensional electron gas. Nature 630, 329–334 (2024).38867129 10.1038/s41586-024-07434-9
33. Badawy G High mobility stemless InSb nanowires Nano Lett. 2019 19 3575 3582 10.1021/acs.nanolett.9b00545 31094527
Badawy, G. et al. High mobility stemless InSb nanowires. Nano Lett. 19, 3575–3582 (2019).31094527 10.1021/acs.nanolett.9b00545
34. Hornibrook JM Frequency multiplexing for readout of spin qubits Appl. Phys. Lett. 2014 104 103108 10.1063/1.4868107
Hornibrook, J. M. et al. Frequency multiplexing for readout of spin qubits. Appl. Phys. Lett. 104, 103108 (2014).10.1063/1.4868107
35. Wang J-Y Parametric exploration of zero-energy modes in three-terminal InSb-Al nanowire devices Phys. Rev. B 2022 106 075306 10.1103/PhysRevB.106.075306
Wang, J.-Y. et al. Parametric exploration of zero-energy modes in three-terminal InSb-Al nanowire devices. Phys. Rev. B 106, 075306 (2022).10.1103/PhysRevB.106.075306
36. Heedt S Shadow-wa lithography of baistic superconductor-semiconductor quantum devices Nat. Commun. 2021 12 4914 10.1038/s41467-021-25100-w 34389705
Heedt, S. et al. Shadow-wa lithography of baistic superconductor-semiconductor quantum devices. Nat. Commun. 12, 4914 (2021).34389705 10.1038/s41467-021-25100-w
37. Borsoi F Single-shot fabrication of semiconducting-superconducting nanowire devices Adv. Funct. Mater. 2021 31 2102388 10.1002/adfm.202102388
Borsoi, F. et al. Single-shot fabrication of semiconducting-superconducting nanowire devices. Adv. Funct. Mater. 31, 2102388 (2021).10.1002/adfm.202102388
38. Mazur GP Spin-mixing enhanced proximity effect in aluminum-based superconductor-semiconductor hybrids Adv. Mater. 2022 34 2202034 10.1002/adma.202202034
Mazur, G. P. et al. Spin-mixing enhanced proximity effect in aluminum-based superconductor-semiconductor hybrids. Adv. Mater. 34, 2202034 (2022).10.1002/adma.202202034
39. Reeg C Loss D Klinovaja J Metaization of a Rashba wire by a superconducting layer in the strong-proximity regime Phys. Rev. B 2018 97 165425 10.1103/PhysRevB.97.165425
Reeg, C., Loss, D. & Klinovaja, J. Metaization of a Rashba wire by a superconducting layer in the strong-proximity regime. Phys. Rev. B 97, 165425 (2018).10.1103/PhysRevB.97.165425
40. Kiršanskas G Goldstein M Flensberg K Glazman LI Paaske J Yu-Shiba-Rusinov states in phase-biased superconductor–quantum dot–superconductor junctions Phys. Rev. B 2015 92 235422 10.1103/PhysRevB.92.235422
Kiršanskas, G., Goldstein, M., Flensberg, K., Glazman, L. I. & Paaske, J. Yu-Shiba-Rusinov states in phase-biased superconductor–quantum dot–superconductor junctions. Phys. Rev. B 92, 235422 (2015).10.1103/PhysRevB.92.235422
41. Piet J-D Andreev bound states in supercurrent-carrying carbon nanotubes revealed Nat. Phys. 2010 6 965 969 10.1038/nphys1811
Piet, J.-D. et al. Andreev bound states in supercurrent-carrying carbon nanotubes revealed. Nat. Phys. 6, 965–969 (2010).10.1038/nphys1811
42. Lee EJH Spin-resolved Andreev levels and parity crossings in hybrid superconductor-semiconductor nanostructures Nat. Nanotechnol. 2014 9 79 84 10.1038/nnano.2013.267 24336403
Lee, E. J. H. et al. Spin-resolved Andreev levels and parity crossings in hybrid superconductor-semiconductor nanostructures. Nat. Nanotechnol. 9, 79–84 (2014).24336403 10.1038/nnano.2013.267
43. Scherübl, Z. et al. Large spatial extension of the zero-energy Yu-Shiba-Rusinov state in a magnetic field. Nat. Commun.10.1038/s41467-020-15322-9 (2020).
44. Scherübl, Z. et al. From cooper pair splitting to nonlocal spectroscopy of a Shiba state. Phys. Rev. Res.10.1103/PhysRevResearch.4.023143 (2022).
45. Schindele J Baumgartner A Maurand R Weiss M Schönenberger C Nonlocal spectroscopy of Andreev bound states Phys. Rev. B 2014 89 045422 10.1103/PhysRevB.89.045422
Schindele, J., Baumgartner, A., Maurand, R., Weiss, M. & Schönenberger, C. Nonlocal spectroscopy of Andreev bound states. Phys. Rev. B 89, 045422 (2014).10.1103/PhysRevB.89.045422
46. Gramich J Baumgartner A Schönenberger C Andreev bound states probed in three-terminal quantum dots Phys. Rev. B 2017 96 1 14 10.1103/PhysRevB.96.195418
Gramich, J., Baumgartner, A. & Schönenberger, C. Andreev bound states probed in three-terminal quantum dots. Phys. Rev. B 96, 1–14 (2017).10.1103/PhysRevB.96.195418
47. Wang G Singlet and triplet Cooper pair splitting in hybrid superconducting nanowires Nature 2022 612 448 453 10.1038/s41586-022-05352-2 36418399
Wang, G. et al. Singlet and triplet Cooper pair splitting in hybrid superconducting nanowires. Nature 612, 448–453 (2022).36418399 10.1038/s41586-022-05352-2
48. Trocha, P. & Weymann, I. Spin-resolved Andreev transport through double-quantum-dot cooper pair splitters. Phys. Rev. B10.1103/PhysRevB.91.235424 (2015).
49. Scherübl Z Páyi A Csonka S Transport signatures of an Andreev molecule in a quantum dot–superconductor–quantum dot setup Beilstein J. Nanotechnol. 2019 10 363 378 10.3762/bjnano.10.36 30800576
Scherübl, Z., Páyi, A. & Csonka, S. Transport signatures of an Andreev molecule in a quantum dot–superconductor–quantum dot setup. Beilstein J. Nanotechnol. 10, 363–378 (2019).30800576 10.3762/bjnano.10.36
50. Saldaña, J.C.E. et al. Two-impurity Yu-Shiba-Rusinov states in coupled quantum dots. Phys. Rev. B10.1103/PhysRevB.102.195143 (2020).
51. Kürtössy O Andreev molecule in parael InAs nanowires Nano Lett. 2021 21 7929 7937 10.1021/acs.nanolett.1c01956 34538054
Kürtössy, O. et al. Andreev molecule in parael InAs nanowires. Nano Lett. 21, 7929–7937 (2021).34538054 10.1021/acs.nanolett.1c01956
52. Aasen D Milestones toward Majorana-based quantum computing Phys. Rev. X 2016 6 031016
Aasen, D. et al. Milestones toward Majorana-based quantum computing. Phys. Rev. X 6, 031016 (2016).
53. Knapp C Karzig T Lutchyn RM Nayak C Dephasing of Majorana-based qubits Phys. Rev. B 2018 97 125404 10.1103/PhysRevB.97.125404
Knapp, C., Karzig, T., Lutchyn, R. M. & Nayak, C. Dephasing of Majorana-based qubits. Phys. Rev. B 97, 125404 (2018).10.1103/PhysRevB.97.125404
54. Mishmash RV Bauer B von Oppen F Alicea J Dephasing and leakage dynamics of noisy Majorana-based qubits: Topological versus Andreev Phys. Rev. B 2020 101 075404 10.1103/PhysRevB.101.075404
Mishmash, R. V., Bauer, B., von Oppen, F. & Alicea, J. Dephasing and leakage dynamics of noisy Majorana-based qubits: Topological versus Andreev. Phys. Rev. B 101, 075404 (2020).10.1103/PhysRevB.101.075404
55. Boross P Páyi A Dephasing of Majorana qubits due to quasistatic disorder Phys. Rev. B 2022 105 035413 10.1103/PhysRevB.105.035413
Boross, P. & Páyi, A. Dephasing of Majorana qubits due to quasistatic disorder. Phys. Rev. B 105, 035413 (2022).10.1103/PhysRevB.105.035413
56. Bordin, A. et al. Signatures of Majorana protection in a three-site Kitaev chain. arXiv10.48550/arXiv.2402.19382 (2024).
57. Hays M Coherent manipulation of an Andreev spin qubit Science 2021 373 430 433 10.1126/science.abf0345 34437115
Hays, M. et al. Coherent manipulation of an Andreev spin qubit. Science 373, 430–433 (2021).34437115 10.1126/science.abf0345
58. Pita-Vidal M Direct manipulation of a superconducting spin qubit strongly coupled to a transmon qubit Nat. Phys. 2023 19 1110 1115 10.1038/s41567-023-02071-x
Pita-Vidal, M. et al. Direct manipulation of a superconducting spin qubit strongly coupled to a transmon qubit. Nat. Phys. 19, 1110–1115 (2023).10.1038/s41567-023-02071-x
59. Choi M-S Bruder C Loss D Spin-dependent Josephson current through double quantum dots and measurement of entangled electron states Phys. Rev. B 2000 62 13569 13572 10.1103/PhysRevB.62.13569
Choi, M.-S., Bruder, C. & Loss, D. Spin-dependent Josephson current through double quantum dots and measurement of entangled electron states. Phys. Rev. B 62, 13569–13572 (2000).10.1103/PhysRevB.62.13569
60. Leijnse M Flensberg K Coupling spin qubits via superconductors Phys. Rev. Lett. 2013 111 060501 10.1103/PhysRevLett.111.060501 23971543
Leijnse, M. & Flensberg, K. Coupling spin qubits via superconductors. Phys. Rev. Lett. 111, 060501 (2013).23971543 10.1103/PhysRevLett.111.060501
61. Hensgens T Quantum simulation of a Fermi-Hubbard model using a semiconductor quantum dot array Nature 2017 548 70 73 10.1038/nature23022 28770852
Hensgens, T. et al. Quantum simulation of a Fermi-Hubbard model using a semiconductor quantum dot array. Nature 548, 70–73 (2017).28770852 10.1038/nature23022
62. Dehoain JP Nagaoka ferromagnetism observed in a quantum dot plaquette Nature 2020 579 528 533 10.1038/s41586-020-2051-0 32123352
Dehoain, J. P. et al. Nagaoka ferromagnetism observed in a quantum dot plaquette. Nature 579, 528–533 (2020).32123352 10.1038/s41586-020-2051-0
63. Hsiao T-K Exciton transport in a germanium quantum dot ladder Phys. Rev. X 2024 14 011048
Hsiao, T.-K. et al. Exciton transport in a germanium quantum dot ladder. Phys. Rev. X 14, 011048 (2024).
64. Martinez, E. A. et al. Measurement circuit effects in three-terminal electrical transport measurements. arXiv10.48550/arXiv.2104.02671 (2021).
