Pith. sign in

REVIEW 2 major objections 4 minor

All four anticrossing gaps in a silicon double quantum dot encode the full set of tunnel and valley phases.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-13 01:33 UTC pith:TTQ2OG4H

load-bearing objection Solid experimental completion of the Si DQD valley-phase picture: all four anticrossings via high-res DAXS give the six parameters, including the previously missing interdot-intervalley phase, with transparent fitting and quantified method error. the 2 major comments →

arxiv 2607.09638 v2 pith:TTQ2OG4H submitted 2026-07-10 quant-ph cond-mat.mes-hall

Complete measurement of tunnel- and valley-coupling parameters in a silicon double quantum dot

classification quant-ph cond-mat.mes-hall
keywords silicon quantum dotsvalley phasestunnel couplingDAXSintervalley couplingSi/SiGedouble quantum dot
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

In silicon quantum dots the electron can tunnel both within a valley and between the two valleys of the conduction band. Those intervalley processes carry complex phases fixed by the random Ge alloy in the well. Until now those phases were not fully measured. This work uses high-resolution delta-axis spectroscopy to resolve every anticrossing among the ground and first excited valley states of both dots, then inverts the four measured gaps into the complete six-parameter Hamiltonian, including the previously inaccessible interdot-intervalley coupling and its phase. The authors further show that simple dimensionless ratios of the gaps themselves equal the valley-phase difference and the ratio of intervalley to intervalley tunnel strength. Because the phases change when the dots are moved or the barrier is raised, the same method maps how the atomic landscape imprints itself on every quantity that depends on valley phase: spin-orbit coupling, g-factors, and two-qubit exchange. Completing that map is presented as a necessary step toward reliable, scalable silicon spin qubits.

Core claim

High-resolution DAXS spectra that resolve all four anticrossings of a single-electron silicon double quantum dot allow simultaneous extraction of the six Hamiltonian parameters tc, |ΔL|, |ΔR|, |ΔRL|, Δφ and ΔφRL. In the large-valley-splitting limit the measured gap ratios satisfy tan(Δφ/2) ≈ (Ege + Eeg)/(Egg + Eee) and |ΔRL|/tc ≈ √(F−/F+), so the valley phases are directly readable from the spectrum. Both phases evolve with gate voltage as the dots sample different alloy disorder, and |ΔRL| is non-negligible over the full experimental range.

What carries the argument

Delta-axis spectroscopy (DAXS) combined with the four-level valley Hamiltonian. DAXS maps energy versus detuning by square-wave pulsing along the common-mode axis δ; the four measured anticrossing gaps are then inverted through the closed-form relations that connect the ground/excited-valley tunnel couplings to the underlying ±z valley phases.

Load-bearing premise

That a four-level model with one orbital per dot, plus the choice that the parity of the tunnel-coupling signs is odd and that the ordinary tunnel coupling is larger than the interdot-intervalley coupling, is enough to invert every measured spectrum.

What would settle it

Repeat the same high-resolution DAXS protocol on a device whose valley splittings are deliberately made comparable to or smaller than the tunnel couplings; if the extracted phases and |ΔRL|/tc ratios then disagree with independent microwave spectroscopy or with the same four-level inversion, the completeness claim fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Sample-to-sample scatter in g-factors, spin-orbit matrix elements and exchange can now be attributed to measurable local valley phases rather than treated as unexplained device variation.
  • Gate-voltage maps of Δφ and ΔφRL become a practical tool for locating spatial regions of enhanced or suppressed valley splitting.
  • Two-qubit gate design can incorporate the measured |ΔRL| and its phase instead of assuming pure valley-conserving tunneling.
  • Statistical sampling of valley phases across a wafer becomes feasible once DAXS is automated on larger arrays.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the interdot-intervalley phase is fixed by the alloy under the barrier, deliberate barrier-gate shaping could be used to engineer a target ΔφRL for exchange-based gates.
  • The same four-gap inversion applied at finite magnetic field would simultaneously yield the valley-dependent g-tensors, closing the remaining free parameters of the spin-valley Hamiltonian.
  • If the parity-sector ambiguity proves resolvable by weak higher-orbital spectroscopy, the method could become fully model-independent.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript reports a complete experimental extraction of the six single-electron Hamiltonian parameters of a Si/SiGe double quantum dot (tc, |ΔL|, |ΔR|, |ΔRL|, Δφ, ΔφRL), including the previously unmeasured interdot-intervalley coupling and the two relative valley phases. High-resolution delta-axis spectroscopy (DAXS) resolves all four anticrossing gaps between ground and excited valley states; these gaps are fit to a four-level model and inverted via closed-form relations (Eqs. 1–4 and S11–S14). The authors show that gap ratios encode the phases (tan(Δφ/2) ≈ (Ege + Eeg)/(Egg + Eee) and |ΔRL|/tc ≈ √(F−/F+)) and track the evolution of the parameters under screening-gate and barrier-gate voltage, attributing the changes to alloy-disorder sampling as the dots move. Method uncertainties from higher orbitals and δ-dependent confinement deformation are quantified with 128 effective-mass disorder realizations and folded into the reported error bars.

Significance. If the extraction is reliable, the work closes a long-standing gap in Si quantum-dot physics: valley phases control tunnel-coupling ratios, spin–orbit matrix elements, valley–orbit mixing, and g-factor variations, yet had never been measured in full. The demonstration that all four anticrossings are required and that |ΔRL|/tc is non-negligible (≈ 0.2–0.25) revises the standard modeling assumption that interdot-intervalley coupling can be ignored. The transparent pipeline (global DAXS → high-resolution 1-D scans → four-level least-squares → closed-form inversion), Monte-Carlo measurement-error propagation, and independent effective-mass validation of method uncertainties constitute a reusable metrology tool for mapping valley disorder across devices. The results are therefore of immediate interest for qubit variability, shuttling, and g-factor engineering in Si/SiGe.

major comments (2)
  1. Supp. S1D and Fig. S1: the discrete parity-sector choice s = −1 (and the convention tc > |ΔRL|) is load-bearing for the reported phases and |ΔRL|. The authors correctly show that s = +1 places the extracted |ΔRL|/|Δi| ratios deep in the tails of the disorder distributions while s = −1 places them near the modes. This statistical argument is persuasive for the present data set, but the manuscript should state more explicitly that the choice remains a modeling assumption whose validity must be re-checked whenever the measured |ΔRL|/tc approaches unity or the valley splittings become comparable to the tunnel couplings. A short sentence in the main text (near the discussion of Fig. 3i,j) would suffice.
  2. Methods and Supp. S4: lever-arm calibration is performed at low interdot tunnel coupling and then rescaled for the DAXS regime by matching electron temperature. The residual uncertainty in this rescaling is absorbed into the method-error budget via the simulated δ-dependent confinement deformation, which is appropriate. However, the main text never quantifies how large the lever-arm change actually is (“typically less than 10 %” is only in Methods). Adding the measured percentage change (or an upper bound) would let readers judge whether the method-error floor is dominated by lever arms or by higher-orbital hybridization.
minor comments (4)
  1. Fig. 2i: the estimated right-dot position under BC is surprising; a one-sentence reminder that the high BC voltage plus alloy disorder can pull the wavefunction under the barrier would help non-specialist readers.
  2. Eqs. (1)–(4) and the subsequent gap-ratio approximations: the sign convention Eμν = ±2 t~μν is mentioned only in the Supplementary Materials. A brief note in the main text (or a footnote) would prevent confusion when readers compare Fig. 2f with the raw gaps.
  3. Abstract and conclusion: the claim that valley phases control Landé g-factors is correct but is not demonstrated in the present data set. Softening the language to “parameters that depend on them, including …” (already used) is fine; avoid implying that g-factor maps were extracted here.
  4. Typographical consistency: “Landég-factors” appears without a space in the abstract; “delta-axis spectroscopy (DAXS)” is introduced cleanly but later sometimes written “DAXS spectroscopy.”

Circularity Check

1 steps flagged

No significant circularity: parameters are inverted from measured anticrossing energies via a self-contained four-level model; minor self-citation of the authors' prior DAXS technique is not load-bearing for the physics claims.

specific steps
  1. self citation load bearing [Methods / Ref. [29]; main text p. 3]
    "We perform DAXS here by applying periodic, square wave voltage pulses along the δ tuning axis... DAXS is a recently developed technique that directly maps the energy states of a DQD as a function of detuning [29]."

    The measurement technique that supplies all raw energy levels is the authors' own prior work. This is a normal self-citation of a method paper and is not load-bearing for the inversion or the valley-phase claims; the physics extraction stands on the new spectra and the self-contained four-level algebra.

full rationale

The derivation chain is experimental extraction followed by algebraic inversion. High-resolution DAXS supplies 24 energy-level positions (4 anticrossings imes 3 detunings imes 2 levels). These are least-squares fitted to the spectrum of the four-level Hamiltonian H' (Eq. S6 / main-text Fig. 1c), yielding the four real tunnel couplings ˜tµ u and the two valley splittings. Closed-form inversion (Eqs. S11–S14, derived in Supp. S1 from the unitary change of basis) then returns tc, |ΔRL|, Δϕ, ΔϕRL. The approximate gap-ratio encodings tan(Δϕ/2) ≈ (Ege+Eeg)/(Egg+Eee) and |ΔRL|/tc ≈ √(F−/F+) are rearrangements of the same Eqs. (1)–(4) in the large-valley-splitting limit; comparing them to the full-fit values (Figs. 2f,h) is therefore a consistency check on the data, not a prediction forced by a fitted input. The discrete parity-sector choice s = −1 (and the convention tc > |ΔRL|) is resolved by comparing the extracted |ΔRL|/|Δi| ratios against independent effective-mass disorder distributions (Supp. S3C, Fig. 3j / S1); those distributions are not fitted to the present data. Method uncertainties arising from higher orbitals and δ-dependent confinement deformation are quantified by 128 independent microscopic simulations (Supp. S4) and folded into the reported error bars; they do not redefine the central quantities. Self-citations (prior DAXS demonstration, alloy-disorder statistics) supply the measurement technique and the theoretical distributions used for post-hoc comparison, but the load-bearing inversion itself is performed on new spectra with equations derived inside the paper. No self-definitional loop, no fitted-input-called-prediction, and no uniqueness theorem imported from the authors appear.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The experimental extraction rests on standard effective-mass and electrostatic modeling of Si/SiGe quantum dots plus a small number of gauge and parity conventions needed to invert the four-level spectrum. No new particles or forces are postulated; free parameters are the six Hamiltonian quantities fitted to each DAXS spectrum and the material parameters of the disorder model used only for uncertainty estimation and statistical comparison.

free parameters (3)
  • tc, |ΔL|, |ΔR|, |ΔRL|, Δϕ, ΔϕRL (per gate-voltage setting)
    Six parameters of the four-level Hamiltonian fitted by least-squares to the 24 measured energy levels of each DAXS spectrum; central extracted quantities of the paper.
  • E_Ge = 0.6 eV, n_Ge profile (well 1.7 %, barrier 30 %, w = 1.9 nm)
    Material parameters of the alloy-disorder model used in Supp. S3–S4 for method-uncertainty simulations and |ΔRL|/|Δi| distributions; taken from literature but fixed by hand for the present calculations.
  • ε0, Eoff, κ (spectrum offsets and tilt)
    Nuisance parameters in the four-level fit (Eq. S28) that absorb lever-arm miscalibration and global energy shifts.
axioms (4)
  • domain assumption Single-orbital four-level model (one spatial orbital per dot plus valley) captures the low-energy spectrum near the four anticrossings.
    Stated in Sec. ‘Tunnel couplings and valley couplings’ and Supp. S1; higher orbitals are treated only as a quantified method uncertainty (S4).
  • ad hoc to paper Parity sector s = −1 and the convention tc > |ΔRL| may be chosen without loss of physical content.
    Gauge and inversion ambiguities of Eqs. S11–S14; justified by statistical likelihood against the disorder model (Fig. S1) rather than direct spectral discrimination.
  • domain assumption Alloy disorder is the dominant source of intervalley coupling; the projected disorder fields are circular Gaussian with site-diagonal statistics.
    Standard effective-mass treatment (Supp. S3) used for both uncertainty estimation and the |ΔRL|/|Δi| distributions of Fig. 3j.
  • domain assumption Lever arms can be calibrated from bias triangles and line-shape analysis and remain sufficiently constant under the DAXS pulse amplitudes.
    Methods section; residual tilt is absorbed by the free parameter κ.

pith-pipeline@v1.1.0-grok45 · 36076 in / 2969 out tokens · 30795 ms · 2026-07-13T01:33:17.489131+00:00 · methodology

0 comments
read the original abstract

Tunneling is essential in the initialization, measurement, and control of quantum dot qubits. In silicon, such tunneling connects not only the qubit states but also valley minima in the conduction band on opposite sides of the Brillouin zone, with large consequences for the quantum dot behavior. Here we present a full characterization of the intravalley and intervalley tunnel couplings, including their complex phases -- the valley phases. These phases are shown to control measurable parameters, including the ratios of the gaps at anticrossings between quantum states of a double quantum dot. The valley phases themselves evolve as a function of the quantum dot gate voltages and depend on the underlying atomic structure of the quantum well. Knowledge of the valley phases completes the picture and fills a key gap in our understanding of sample-wide variations of valley couplings and the physical parameters that depend on them, including spin-orbit coupling, valley-orbit mixing, and Land\'e $g$-factors.

Figures

Figures reproduced from arXiv: 2607.09638 by Benjamin D. Woods, Daniel J. King, Jonathan C. Marcks, J. Reily, M. A. Eriksson, Mark Friesen, Minyoung Kim.

Figure 1
Figure 1. Figure 1: b. Here, a single spatial orbital (L/R) is considered in each dot. In addition, each spatial orbital has a valley degree of freedom (±z) due to the two-fold degenerate valley minima in the conduction band of the biaxially strained Si quantum well [36]. The various level cou￾plings are illustrated by arrows in Fig. 1b. These include valley-conserving tunnel coupling tc as well as intradot and interdot valle… view at source ↗
Figure 2
Figure 2. Figure 2: a the anticrossing gaps are labeled, while in Fig. 2b, the locations of the high-resolution DAXS measurements are indicated. Fig. 2c shows a typical high-resolution measurement, where the peaks are fit to derivatives of the Fermi-Dirac distribution. Following the procedure described above, for each of six values of VSG, we fit these peak locations in order to extract the Hamilto￾nian parameters, yielding t… view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.