REVIEW 3 major objections 5 minor 1 cited by
Resolving Andreev spin qubits in germanium-based Josephson junctions
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A germanium Josephson junction's geometry can suppress its Andreev spin-qubit frequency below the thermal limit, and a narrower channel with higher filling can make it resolvable.
desk verdict A useful, honest theory paper explaining a null result in Ge Andreev spin qubits, with caveats about the assumed pairing symmetry. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the quasi-one-dimensional heavy-hole band structure of a germanium two-dimensional hole gas, confined by widths $L_y$ and $L_z$, with Rashba spin-orbit coupling. In the regime where only the heavy-hole band is filled, the model yields two right- and two left-moving bands with different Fermi velocities $v_1$ and $v_2$; this velocity difference is the origin of the nondegenerate Andreev levels. The paper defines the Andreev bound-state energy $\epsilon_j = \Delta \cos(\phi/2)\big/(1 + (L\Delta/v_{Fj})\sin(\phi/2))$ and the qubit frequency $\omega_q = |\epsilon_1 - \epsilon_2|$. The analytical result is supplemented by a tight-binding Hamiltonian that includes the full spin-3/2 Luttinger-Kohn model, Birs-Pikus strain, Zeeman coupling, and interface transparency, and it is this numerical model that produces the geometry- and filling-dependence of $\omega_q$.
What would settle it
Measure the tunnelling spectrum of a germanium Josephson junction with $L_y \approx 55\,\mathrm{nm}$, $L_z \approx 25\,\mathrm{nm}$, and filling $\bar{\mu} \approx 0.66\Delta_{HL}$ at $\phi = \pi/2$, and look for two spin-split Andreev levels separated by at least 1 GHz; a null result would contradict the model's core prediction.
Extended reading notes
Core claim
The central claim is that the Andreev spin qubit frequency in a germanium Josephson junction is controlled by the junction geometry through the difference in Fermi velocities of the two heavy-hole bands. The authors derive an analytic formula, $\epsilon_j = \Delta \cos(\phi/2)\big/(1 + (L\Delta/v_{Fj})\sin(\phi/2))$, showing that the two Andreev states are split only when the inner and outer heavy-hole bands have different velocities; the qubit frequency is $\omega_q = |\epsilon_1 - \epsilon_2|$. For an experimentally realized $350 \times 150 \times 25\,\mathrm{nm}^3$ junction with $\Delta = 50\,\mu\mathrm{eV}$ and transparency $r = 0.995$, this frequency is about 150 MHz at $\phi = \pi/2$, below the 30 mK thermal limit of about 625 MHz, explaining why spin-resolved Andreev states were not seen. The tight-binding calculation, which includes strain, magnetic field, and imperfect interfaces, predicts that a narrower channel ($L_y = 55\,\mathrm{nm}$) and higher filling ($\bar{\mu} = 0.66 \Delta_{HL}$) raise $\omega_q$ to about 1.3 GHz, making the qubit resolvable.
Load-bearing premise
The calculation assumes the superconductivity induced in germanium is a spin-3/2 singlet that pairs each heavy-hole or light-hole band with its time-reversed partner; if the real proximity pairing has a different spin structure or a heavy/light-hole asymmetry, the predicted qubit frequencies and the explanation of the null result would change.
Editorial extensions
If this is right
- If the model is right, the null result of Ref. [17] is explained by geometry: at $L_x = 350\,\mathrm{nm}$, $L_y = 150\,\mathrm{nm}$, and $\Delta = 50\,\mu\mathrm{eV}$, the qubit frequency falls below the 30 mK thermal limit and cannot be resolved.
- Design rules follow: $L_y$ should sit near the characteristic electric-field length (roughly 7 to 55 nm depending on field), $L_x$ near the coherence length (about 200 nm to 1 $\mu$m), filling near $0.66$ to $0.8\Delta_{HL}$, strain minimized, and interface transparency maximized.
- A resolvable germanium Andreev spin qubit with $\omega_q \approx 1.3\,\mathrm{GHz}$ should be achievable without changing materials, using the experimentally observed gap and transparency of Ref. [17].
- Because the mechanism is generic to two Fermi velocities, the same geometric optimization should apply to planar InAs junctions, with higher confined subbands playing the role of the light-hole band.
- The predicted $g$-factor anisotropy, with $g_y \approx g_z$ near $L_y \approx L_z$ and $g_x$ small, suggests magnetic sweet spots where decoherence from nuclear noise could be minimized.
Reading between the lines
- The paper's mechanism implies that earlier null searches for spin-resolved Andreev states may have been blinded by geometry; re-analyzing existing devices with $L_y \ll L_z$ could reveal the qubit without new materials work.
- The strong geometric dependence of $\omega_q$ also means that the qubit frequency of a germanium Andreev spin qubit could serve as a sensitive local probe of confinement, strain, and filling, useful for device metrology.
- The singlet-pairing assumption is the main uncertainty: if the proximity-induced pairing is not a spin-3/2 singlet, the two-level basis and the reported frequencies could shift, so a microscopic theory of the Ge/Al interface pairing would sharpen or overturn the design guidance.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a theoretical study of Andreev spin qubits (ASQs) in Ge-based Josephson junctions. Starting from the Luttinger-Kohn Hamiltonian with transverse confinement, the authors construct a one-dimensional multiband model and a tight-binding model of a superconductor-normal-superconductor junction, including Rashba spin-orbit coupling, magnetic field, strain, and interface transparency. They derive an effective two-level Andreev Hamiltonian and study the qubit frequency as a function of junction geometry (Ly, Lz, Lx), filling, strain, and magnetic field. They use the model to argue that the failure to resolve spin-resolved Andreev states in a recent Ge experiment (Ref. 17) can be explained by a qubit frequency (~150 MHz) below the thermal energy (~625 MHz), and they propose design changes (Ly~55 nm, filling ~0.66 Delta_HL) that would raise the frequency to ~1.3 GHz.
Significance. If the central model is accepted, the paper provides a concrete and testable explanation for a puzzling null result and gives design guidance for Ge-based ASQs. The tight-binding model is described in sufficient detail to be reproduced, it reproduces expected qualitative features (different Fermi velocities, HH-LH gap, linear dispersion near phi=pi), and the parametric dependence on geometry and filling is computed directly from the model rather than fitted. The main quantitative conclusions, however, rest on a specific assumption about the spin structure of the proximity-induced pairing, which the authors themselves identify as non-unique; this makes the numerical values conditional rather than definitive. The paper is nevertheless valuable as a systematic study of the geometry dependence of ASQ frequencies in a hole-gas platform.
major comments (3)
- [Sec. II.B and Sec. IV] The reduction of the 8x8 problem to four 2x2 blocks and the resulting Andreev energies in Eq. (11) rely on the spin-3/2 singlet pairing H_SC = i tau_y exp(i pi J_y) Delta cos(phi) - i tau_x exp(i pi J_y) Delta sin(phi), which pairs each particle solution Phi^+_{nu,e} with a specific hole solution Phi^-_{nu,h}. This pairing choice is not unique for a hole gas, and the authors concede in Sec. IV that it depends on the microscopic coupling of heavy and light holes to the superconductor and that small HH-LH pairing differences can alter the qubit frequency while large ones can destroy the gap. Since every reported omega_q, including the 150 MHz and 1.3 GHz values in Sec. III.C, is a direct output of this assumption, the central claims are conditional on a pairing structure not established by experiment or first principles. I request either a microscopic estimate of the proximity pairing matrix along the lines of Refs. 25 and 26, or an explicit sensitivity analysis over the HH/LH pairing amplitude ratio.
- [Sec. III.C] The comparison to the experiment of Ref. 17 fixes the interface transparency parameter r=0.995 so that the Andreev level energies match experiment, and no sensitivity scan in r is shown. The conclusion that the null result is explained by omega_q ~ 150 MHz being below 625 MHz is partially robust because the r=1 limit with bulk Delta also gives only ~300 MHz, but the quantitative design recommendation of omega_q ~ 1.3 GHz is a single point in parameter space and needs a sensitivity analysis in both r and the pairing amplitude. Please provide omega_q as a function of r and state how the inferred resolvability threshold depends on this fit parameter.
- [Eq. (13)] In Eq. (13), the expression for B contains four identical sine terms [3 sin(theta_1+theta_2+pi/3) + sin(theta_1+theta_2+pi/3) + sin(theta_1+theta_2+pi/3) + sin(theta_1+theta_2+pi/3)]/4, which reduces to (3/2) sin(theta_1+theta_2+pi/3) and makes the displayed magnetic-field matrix element incorrect as written. This affects the analytic model used for the magnetic-field dependence in Fig. 2(b). Although the numerical tight-binding results in Sec. III do not use Eq. (13), the analytical derivation is presented as supporting intuition and should be corrected.
minor comments (5)
- [Eq. (6)] As printed, Eq. (6) appears to contain a product (k^2/m_s) tilde_J_x^2 tilde_J_y^2 and a bare -omega_z, whereas consistency with Eq. (3) requires separate terms -(k^2/m_s) tilde_J_x^2, -omega_z tilde_J_y^2, and -omega_y tilde_J_z^2.
- [Sec. II.B] There are small grammatical errors: "two and two-by-two Hamiltonians" should be "two-by-two Hamiltonians", and "are are separated" should be "are separated".
- [Fig. 5 and Fig. 7] The symbol E_s is used for strain in Sec. III.A and Fig. 7, but Fig. 5's inset appears to use "Es" for the characteristic electric-field length; these should use distinct notation to avoid confusion.
- [Sec. III.A and Sec. III.C] The coherence length is estimated as xi ~ 200 nm in Sec. III.A, but Sec. III.C recommends Lx ~ xi ~ 1 micron; please reconcile these two estimates.
- [Eq. (15)] The tight-binding superconducting term in Eq. (15) uses a bare J without defining it; it should be explicit that this is the same exp(i pi J_y) factor used in the analytical Hamiltonian, or a separate definition should be given.
Circularity Check
No significant circularity: the central qubit-frequency predictions are direct outputs of a Luttinger-Kohn/tight-binding model; the experimental comparison calibrates one transparency parameter to the measured Andreev energy but does not fit the predicted splitting.
full rationale
The derivation chain is self-contained. The analytical Andreev spectrum uses Eq. (11), taken from the independent Ref. [22], and evaluates it with Fermi velocities obtained from the Luttinger-Kohn Hamiltonian (Eqs. (1)-(4)); the qubit frequency omega_q = |epsilon_1 - epsilon_2| is then a direct output of the band-structure input. The numerical tight-binding model (Eqs. (14)-(17)) is a direct simulation with Luttinger parameters, confinement lengths, and an explicit singlet-pairing assumption, and yields the geometry/filling dependence in Secs. III.A-III.B. In the experimental comparison (Sec. III.C), the single free parameter r = 0.995 is calibrated to match the measured Andreev-state energy (~3 GHz); the reported omega_q ~ 150 MHz is a distinct output of the calibrated model, not the fitted quantity, so the explanation of the null result is a postdiction rather than a reduction to the input. The authors' Sec. IV caveat that the spin-3/2 singlet pairing assumption 'is not unique' is an acknowledged modeling limitation, not a circular step. The only self-citation (Ref. 15) is used for qualitative phase/transparency behavior and nuclear-decoupling expectations, and is not load-bearing for the central geometry/filling predictions.
Assumptions & free parameters
free parameters (1)
- Interface transparency parameter r =
0.995
assumptions (5)
- standard math Bulk Ge is described by the Luttinger-Kohn Hamiltonian with Luttinger parameters gamma1 = 13.35, gamma2 = 4.25, gamma3 = 5.69 (Eq. 1).
- domain assumption Confinement along the y and z axes is an infinite square well of widths Ly and Lz (Eq. 3 and tight-binding setup).
- domain assumption The induced superconductivity is singlet pairing in spin-3/2 space with a step-function phase profile (Eq. 9 and HSC).
- domain assumption Only the heavy-hole band is filled; light holes are neglected in the Andreev derivation.
- domain assumption Magnetic field and barrier transparency are treated perturbatively in the analytical model (Sec. II.B).
Cite this review
Pith. "Pith review of Resolving Andreev spin qubits in germanium-based Josephson junctions." pith.science (2026). https://pith.science/paper/DNJU7SRU
@misc{pith2026250613988,
author = {Pith},
title = {Pith review of: Resolving Andreev spin qubits in germanium-based Josephson junctions},
year = {2026},
howpublished = {\url{https://pith.science/paper/DNJU7SRU}},
note = {Machine review of arXiv:2506.13988}
}
read the original abstract
Andreev spin qubits (ASQs) are a promising platform for quantum information processing which benefit from both the small footprint of semiconducting spin qubits and the long range connectivity of superconducting qubits. While state-of-the-art experiments have developed ASQs in InAs nanowires, these realizations are coherence-time limited by nuclear magnetic noise which cannot be removed by isotopic purification. In Ge-based Josephson junctions, which can be isotopically purified, Andreev states have been experimentally observed but spin-resolved Andreev states remain elusive. Here, we theoretically demonstrate that the geometry of the Josephson junction can limit the qubit frequency to values below typical experimental temperatures and render the ASQ effectively invisible. ASQs could be experimentally resolved by judiciously choosing the geometry of the junction and filling of the underlying Ge. Our comprehensive study of ASQ frequency on in situ and ex situ experimentally controllable parameters provides design guidance of Ge-based Josephson junctions and paves the way towards realization of high-coherence ASQs.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 1 Pith paper
-
A 66-Gb/s/5.5-W RISC-V Many-Core Cluster for 5G+ Software-Defined Radio Uplinks
A 1024-core RISC-V cluster is claimed to process 5G NR uplink (PUSCH) at 66 Gb/s in 1.7 ms at under 6 W, from a placed-and-routed simulation in 12-nm CMOS.
Reference graph
Works this paper leans on
-
[1]
author author G. Scappucci , author C. Kloeffel , author F. A. \ Zwanenburg , author D. Loss , author M. Myronov , author J.-J. \ Zhang , author S. De Franceschi , author G. Katsaros , \ and\ author M. Veldhorst ,\ @noop journal journal Nature Reviews Materials \ volume 6 ,\ pages 926 ( year 2021 ) NoStop
work page 2021
-
[2]
author author H. Watzinger , author J. Kuku c ka , author L. Vuku s i \'c , author F. Gao , author T. Wang , author F. Sch \"a ffler , author J.-J. \ Zhang , \ and\ author G. Katsaros ,\ @noop journal journal Nature Communications \ volume 9 ,\ pages 3902 ( year 2018 ) NoStop
work page 2018
-
[3]
author author N. W. \ Hendrickx , author D. P. \ Franke , author A. Sammak , author G. Scappucci , \ and\ author M. Veldhorst ,\ @noop journal journal Nature \ volume 577 ,\ pages 487 ( year 2020 ) NoStop
work page 2020
-
[4]
van Riggelen-Doelman , author C.-A
author author F. van Riggelen-Doelman , author C.-A. \ Wang , author S. L. \ de Snoo , author W. I. L. \ Lawrie , author N. W. \ Hendrickx , author M. Rimbach-Russ , author A. Sammak , author G. Scappucci , author C. D \'e prez , \ and\ author M. Veldhorst ,\ @noop journal journal Nature Communications \ volume 15 ,\ pages 5716 ( year 2024 ) NoStop
work page 2024
-
[5]
author author K. M. \ Itoh \ and\ author H. Watanabe ,\ @noop journal journal MRS communications \ volume 4 ,\ pages 143 ( year 2014 ) NoStop
work page 2014
-
[6]
author author S. Bosco \ and\ author D. Loss ,\ @noop journal journal Phys. Rev. Lett. \ volume 127 ,\ pages 190501 ( year 2021 ) NoStop
work page 2021
-
[7]
author author N. W. \ Hendrickx , author L. Massai , author M. Mergenthaler , author F. J. \ Schupp , author S. Paredes , author S. W. \ Bedell , author G. Salis , \ and\ author A. Fuhrer ,\ @noop journal journal Nature Materials \ volume 23 ,\ pages 920 ( year 2024 ) NoStop
work page 2024
-
[8]
author author N. W. \ Hendrickx , author D. P. \ Franke , author A. Sammak , author M. Kouwenhoven , author D. Sabbagh , author L. Yeoh , author R. Li , author M. L. V. \ Tagliaferri , author M. Virgilio , author G. Capellini , author G. Scappucci , \ and\ author M. Veldhorst ,\ @noop journal journal Nature Communications \ volume 9 ,\ pages 2835 ( year 2...
work page 2018
Show all 29 references
-
[9]
Tosato , author V
author author A. Tosato , author V. Levajac , author J.-Y. \ Wang , author C. J. \ Boor , author F. Borsoi , author M. Botifoll , author C. N. \ Borja , author S. Mart \' -S \'a nchez , author J. Arbiol , author A. Sammak , author M. Veldhorst , \ and\ author G. Scappucci ,\ @...
2023
-
[10]
Lakic , author W
author author L. Lakic , author W. I. L. \ Lawrie , author D. van Driel , author L. E. A. \ Stehouwer , author Y. Su , author M. Veldhorst , author G. Scappucci , author F. Kuemmeth , \ and\ author A. Chatterjee ,\ @noop journal journal Nature Materials \ ( year 2025 ) NoStop
2025
-
[11]
author author N. M. \ Chtchelkatchev \ and\ author Y. V. \ Nazarov ,\ 10.1103/PhysRevLett.90.226806 journal journal Phys. Rev. Lett. \ volume 90 ,\ pages 226806 ( year 2003 ) NoStop
2003 doi
-
[12]
Padurariu \ and\ author Y
author author C. Padurariu \ and\ author Y. V. \ Nazarov ,\ 10.1103/PhysRevB.81.144519 journal journal Phys. Rev. B \ volume 81 ,\ pages 144519 ( year 2010 ) NoStop
2010 doi
-
[13]
Hays , author V
author author M. Hays , author V. Fatemi , author D. Bouman , author J. Cerrillo , author S. Diamond , author K. Serniak , author T. Connolly , author P. Krogstrup , author J. Nygård , author A. L. \ Yeyati , author A. Geresdi , \ and\ author M. H. \ Devoret ,\ 10.1126/science...
-
[14]
Pita-Vidal , author A
author author M. Pita-Vidal , author A. Bargerbos , author R. Žitko , author L. J. \ Splitthoff , author L. Grünhaupt , author J. J. \ Wesdorp , author Y. Liu , author L. P. \ Kouwenhoven , author R. Aguado , author B. van Heck , author A. Kou , \ and\ author C. K. \ Andersen ...
-
[15]
Hoffman , author M
author author S. Hoffman , author M. Hays , author K. Serniak , author T. Hazard , \ and\ author C. Tahan ,\ 10.1103/PhysRevB.111.045304 journal journal Phys. Rev. B \ volume 111 ,\ pages 045304 ( year 2025 ) NoStop
2025 doi
-
[16]
Pita-Vidal , author J
author author M. Pita-Vidal , author J. J. \ Wesdorp , \ and\ author C. K. \ Andersen ,\ 10.1103/PRXQuantum.6.010308 journal journal PRX Quantum \ volume 6 ,\ pages 010308 ( year 2025 ) NoStop
2025 doi
-
[17]
Hinderling , author S
author author M. Hinderling , author S. C. \ ten Kate , author M. Coraiola , author D. Haxell , author M. Stiefel , author M. Mergenthaler , author S. Paredes , author S. Bedell , author D. Sabonis , \ and\ author F. Nichele ,\ @noop journal journal PRX Quantum \ volume 5 ,\ p...
2024
-
[18]
author author J. M. \ Luttinger ,\ 10.1103/PhysRev.102.1030 journal journal Phys. Rev. \ volume 102 ,\ pages 1030 ( year 1956 ) NoStop
1956 doi
-
[19]
author author B. A. \ Bernevig \ and\ author S.-C. \ Zhang ,\ @noop journal journal Phys. Rev. Lett. \ volume 95 ,\ pages 016801 ( year 2005 ) NoStop
2005
-
[20]
Mao , author M
author author L. Mao , author M. Gong , author E. Dumitrescu , author S. Tewari , \ and\ author C. Zhang ,\ @noop journal journal Phys. Rev. Lett. \ volume 108 ,\ pages 177001 ( year 2012 ) NoStop
2012
-
[21]
note While we have considered a linear spin-orbit interaction, our results are consistent with a cubic spin-orbit interaction with the substitution k k^3 . Stop
-
[22]
Park \ and\ author A
author author S. Park \ and\ author A. L. \ Yeyati ,\ 10.1103/PhysRevB.96.125416 journal journal Phys. Rev. B \ volume 96 ,\ pages 125416 ( year 2017 ) NoStop
2017 doi
-
[23]
note While we have assumed only the HH band is filled, the calculation proceeds analogously for LH bands. Stop
-
[24]
Bosco , author M
author author S. Bosco , author M. Benito , author C. Adelsberger , \ and\ author D. Loss ,\ @noop journal journal Phys. Rev. B \ volume 104 ,\ pages 115425 ( year 2021 ) NoStop
2021
-
[25]
Babkin , author B
author author S. Babkin , author B. Joecker , author K. Flensberg , author M. Serbyn , \ and\ author J. Danon ,\ @noop title Superconducting proximity effect in two-dimensional hole gases , \ ( year 2024 ) NoStop
2024
-
[26]
author author D. M. \ Pino , author R. Seoane-Souto , author M. J. \ Calder \'o n , author R. Aguado , \ and\ author J. C. \ Abadillo-Uriel ,\ @noop title Theory of superconducting proximity effect in hole-based hybrid semiconductor-superconductor devices , \ ( year 2024 ) NoStop
2024
-
[27]
Machnikowski , author K
author author P. Machnikowski , author K. Gawarecki , \ and\ author L. Cywi n \' n ski ,\ @noop journal journal Phys. Rev. B \ volume 100 ,\ pages 085305 ( year 2019 ) NoStop
2019
-
[28]
Testelin , author F
author author C. Testelin , author F. Bernardot , author B. Eble , \ and\ author M. Chamarro ,\ @noop journal journal Phys. Rev. B \ volume 79 ,\ pages 195440 ( year 2009 ) NoStop
2009
-
[29]
Fischer , author W
author author J. Fischer , author W. A. \ Coish , author D. V. \ Bulaev , \ and\ author D. Loss ,\ 10.1103/PhysRevB.78.155329 journal journal Phys. Rev. B \ volume 78 ,\ pages 155329 ( year 2008 ) NoStop
2008 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.