{"id":"6ced6e77-013a-4d17-a09a-41a1ab799bb0","arxiv_id":"2501.06671","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"In bilayer graphene double quantum dots, the strength and magnetic-field dependence of spin and valley blockade are set by the orbital splitting, the short-range electron-electron interaction, and the difference in valley g-factors between symmetric and antisymmetric two-particle orbital states.","lead":"This paper measures how current is blocked through a pair of graphene quantum dots holding two electrons, and shows that the blockades depend on higher orbital states that were previously ignored in such experiments. Understanding these limits matters for graphene spin and valley qubits, because Pauli blockade is the standard way to read out their quantum state.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative blockade formulas rest on the weak-interdot-coupling assumption; if the actual coupling is large enough to hybridize (1,1) and (0,2) states, the extracted ΔOrb, δ2, and g_v^a are not uniquely determined.","rationale":"I agree with the reader's CONDITIONAL verdict and with the identification of the independent-dot, negligible-coupling premise as the load-bearing assumption. The paper reports a real experimental effect, and the qualitative switch from valley to spin blockade with field is convincingly reproduced by the rate-equation simulation, so this is not a reason to reject. However, the strongest claim is quantitative: the blockade extents are given by the formulas preceding Eq. (6), and those formulas are derived from an uncoupled single-dot energy ladder. The authors state the assumption explicitly and then concede, in the discussion of Figs. 4 and 5, that the interdot coupling may be too large to justify it. If the coupling is finite, it renormalizes the resonance positions and can be absorbed into the fitted ΔOrb, δ2, and valley g-factors, so the same data may be consistent with different parameter sets. The baseline uncertainty noted by the reader is real, but it shifts the origin of the detuning axis uniformly and does not, by itself, change the field slopes that determine g_v^a and g_v^s. The finite-coupling issue can change those slopes and therefore more directly threatens the central quantitative formulas. A concrete finite-coupling simulation test would settle whether the extracted parameters are stable; until that test is run, CONDITIONAL is the appropriate verdict.","tokens_in":23035,"tokens_out":11002,"duration_ms":110198,"concrete_test":"Extend the deposited rate-equation code, or use a two-site Hubbard model with the same single-particle energies, to include an interdot tunnel amplitude t_c as an avoided-crossing term between each (1,1) and (0,2) state, with t_c in the 1-10 μeV range estimated from the interdot transition width. Refit the five energy parameters to the detuning-versus-field maps in Figs. 4-5. If the best-fit ΔOrb, δ2, and g_v^a shift by more than the spread already present between triple points A and B (ΔOrb from 0.575 to 0.7 meV, g_v^a from 18 to 19), the quantitative blockade formulas are not robust to the stated model limitation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims, the expressions for ε_VB and ε_SB given before Eq. (6), are derived from detunings of independent single-dot levels. The derivation assumes the (1,1) state is a product of two uncoupled single-particle states ('we assume a negligible mixing between the two QDs due to a small interdot tunnel coupling'), and the authors later concede that the interdot coupling is 'potentially being too large to justify the approximation.' Under that premise, the slope of the valley-blockade boundary is exactly -g_v^a μB and the spin-blockade boundary slope is (g_v^s - g_s - g_v^a) μB. If the interdot tunnel coupling is not negligible, however, (1,1) and (0,2) configurations hybridize, shifting every resonance by an amount that depends on the coupling, the detuning, and the magnetic field. The five fitted parameters (ΔOrb, δ2, g_v^(1), g_v^s, g_v^a) can then absorb the coupling, so the same experimental maps could be reproduced with different parameter sets. The claimed 'limited only by properties of the antisymmetric orbital states' would then be an artifact of the model rather than a directly measured property. The authors themselves list finite interdot coupling as the first explanation for the observed discrepancies between simulation and data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports magnetotransport measurements of a bilayer graphene double quantum dot at the (1,1)↔(2,0) and (1,1)↔(0,2) charge transitions, supported by rate-equation simulations that include both symmetric and antisymmetric two-particle orbital states. The central experimental observation is a magnetic-field-tunable switch from valley blockade at low perpendicular field to spin blockade at higher field, together with resonances attributed to antisymmetric orbital states. The authors derive limiting detuning expressions, ε_VB = Δ_Orb − Δ_SO − g_v^a μ_B B_⊥ and ε_SB = Δ_Orb − Δ_SO − δ_2 + (g_v^s − g_s − g_v^a) μ_B B_⊥, and fit the underlying parameters (Δ_Orb, δ_2, g_v^(1), g_v^s, g_v^a) to reproduce the data. They conclude that the valley-blockade extent is governed by antisymmetric orbital states, while the spin-blockade extent involves properties of both symmetric and antisymmetric orbitals.","tokens_in":23311,"tokens_out":6342,"duration_ms":54871,"significance":"Understanding the limits of Pauli blockade in bilayer graphene double quantum dots is directly relevant for spin- and valley-qubit readout, and the inclusion of antisymmetric two-particle states goes beyond earlier work that considered only the six lowest symmetric states. The qualitative field-driven switch from valley to spin blockade is clearly present in the data and is captured by the simulation, and the paper provides reproducible simulation code and a data-availability statement. However, the quantitative formulas and the parameter extraction rest on a weak-interdot-coupling assumption that the authors themselves question, and the parameters are chosen per triple point without reported uncertainties, so the quantitative central claim is not yet fully supported.","major_comments":[{"comment":"The derivation of the blockade extents assumes that the (1,1) configuration consists of two independent single-particle states with negligible interdot tunnel coupling. In the discussion of the discrepancies between simulation and data the authors state that the interdot coupling is 'potentially being too large to justify the approximation of completely independent single particle states'. Under non-negligible hybridization the resonance condition in Eq. (2) is no longer simply E_(0,2) − E_(1,1) for uncoupled states, and the fitted values of Δ_Orb, δ_2, g_v^s, and g_v^a can absorb the coupling. The central claim that the valley-blockade extent is limited only by properties of the antisymmetric orbital states therefore depends on an assumption that is acknowledged to be questionable. I ask the authors to quantify the interdot tunnel coupling (for example from the stability diagram or from separate measurements) or to include hybridization in the model and demonstrate that the extracted parameters and Eqs. (5)–(6) are stable under this variation.","section":"§4, Eqs. (5)–(6); discussion of Figs. 4 and 5; Appendix B"},{"comment":"The five energy parameters are chosen separately for the two triple points (Δ_Orb = 0.7 meV, δ_2 = 0.34 meV, g_v^a = 19 for triple point A; Δ_Orb = 0.575 meV, δ_2 = 0.2 meV, g_v^a = 18 for triple point B), with no uncertainty estimates or goodness-of-fit measure reported. The difference |g_v^a − g_v^s|, which is central to the claimed role of the antisymmetric valley g-factor, is 1 in one fit and 0 in the other. Moreover, Fig. 10(b) explicitly notes that the valley blockade nearly vanishes in the experimental data at B_⊥ ≈ 0.22 T while persisting in the simulation. These issues mean the quantitative support for Eqs. (5)–(6) is not established; a sensitivity analysis over the parameter set and a direct test of the dependence on g_v^a − g_v^s are needed.","section":"§4, Figs. 4(b) and 10(b,d); Appendix B"},{"comment":"The zero-detuning baseline from which all blockade extents are measured is estimated with an uncertainty of about 15% and shifts slightly with magnetic field due to electrostatic drift and the field-dependent ground-state-to-ground-state transition. These systematic uncertainties are not propagated into the reported slopes and intercepts of ε_VB(B_⊥) and ε_SB(B_⊥). The qualitative switch is unaffected, but the quantitative limits in Eqs. (5)–(6) should be accompanied by uncertainty estimates or a sensitivity analysis before they can be taken as quantitative predictions.","section":"§4, Eq. (5) and baseline definition"}],"minor_comments":[{"comment":"There are several typographical and grammatical errors: 'perfomed' in the second section, 'inot' in Appendix B, 'paramters' in the Fig. 10(d) caption, 'Einital' in Appendix B, and 'to large' in the discussion of discrepancies; these should be corrected.","section":"Throughout"},{"comment":"The sentence referring to 'Figs. 2(c) and 2(c)' should read 'Figs. 2(c) and 2(d)'.","section":"§2, paragraph on triple point B"},{"comment":"The literal placeholder '[REFS]' appears in the text and should be replaced with the appropriate citations for the spin-triplet-to-spin-singlet ground-state transition.","section":"§4, near Eq. (4)"},{"comment":"The data availability statement still contains 'under XXX' as a placeholder; the Zenodo DOI should be inserted before publication.","section":"Appendix C, Data availability"},{"comment":"The phrase 'extend of the valley blockade' should be 'extent of the valley blockade'.","section":"Fig. 8 caption"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about interdot hybridization is legitimate and is, in fact, acknowledged in the manuscript itself. The qualitative observation of a field-tunable valley-to-spin blockade switch is solid, but the quantitative central claim needs stronger support through an estimate of the interdot coupling and a parameter-sensitivity analysis. If the authors can provide those, the paper would be suitable for publication in this journal; I see no fundamental flaw in the experiment itself."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news: this is the first BLG DQD transport analysis that includes the ten antisymmetric two-particle states, and it shows they are what actually opens the blocked regions. The field-tunable valley-to-spin blockade switch is clearly in the data, and the rate-equation model reproduces it, including the positive-bias subtleties that a ground-state-only picture misses. That is a genuine step beyond the earlier Tong et al. work, and the authors are honest about discrepancies between simulation and data. Credit where due: the simulation code is on Zenodo, the parameter choices are explicit, and they walk through the discrepancies instead of hiding them.\n\nThe soft spot is the quantitative layer. Equations (5) and (6) — the centerpiece claims about what limits the valley and spin blockade — are derived assuming the (1,1) state is two independent single-particle dots with negligible interdot tunnel coupling. Later the authors say the coupling is “potentially being too large to justify the approximation.” If the coupling is non-negligible, (1,1) and (0,2) hybridize, and the five fitted parameters can absorb the coupling; the extracted ΔOrb, δ2, and g_a^v are then not uniquely determined, and the claim that valley blockade is “limited only by” antisymmetric orbital properties is as much a model output as a measurement. The stress-test note has it right on that point.\n\nAdd to that: no error bars on the five per-triple-point parameters, an ε=0 baseline with roughly 15% uncertainty that drifts with field, and a Zenodo data link that is still a placeholder. Those are fixable in revision, but they cap how strongly the paper can claim quantitative limits.\n\nThe qualitative phenomenon, though, is an external anchor, not a fit to noise: the same data show the blockade shrinking and a spin-blocked region appearing at the expected field scale. The model is doing real work, not just absorbing noise.\n\nWho should read it: anyone working on graphene DQD qubits and Pauli blockade readout. The paper deserves a serious referee, but I would send it back requesting (1) a test or explicit bound on the interdot-coupling effect—ideally a two-level (1,1)-(0,2) hybridization model—(2) uncertainties or at least a sensitivity analysis over the fitted parameters, and (3) deposition of the data and fixing the placeholder links. With those, the quantitative claims would be defensible.","headline":"Gets the qualitative blockade switch right, but the quantitative limit formulas rest on a weak-coupling assumption the authors admit may be violated; deserves review with serious revision.","tokens_in":23940,"tokens_out":2632,"would_cite":true,"duration_ms":26402,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In bilayer graphene double quantum dots, the spin and valley blockade extent is set by the orbital splitting, electron-electron interactions, and the difference in valley g-factors between symmetric and antisymmetric orbital states.","keywords":["double quantum dots","bilayer graphene","Pauli blockade","valley blockade","spin blockade","electron-electron interaction","orbital splitting","valley g-factor"],"falsifier":"Measure a bilayer graphene double quantum dot with a clearly larger interdot tunnel coupling and record the detuning extent of the spin-blocked region versus perpendicular field; if it deviates from $\\epsilon_{\\mathrm{SB}} = \\Delta_{\\mathrm{Orb}} - \\Delta_{\\mathrm{SO}} - \\delta_2 + (g^s_v - g_s - g^a_v)\\mu_B B_\\perp$ beyond the estimated baseline uncertainty, or if the valley-to-spin switch occurs at a field inconsistent with the triplet-singlet crossing computed from the same parameters, the independent-dots premise is falsified.","tokens_in":22716,"feed_emoji":"🧲","tokens_out":9919,"duration_ms":157224,"temperature":0.7,"pith_summary":"This paper reports magnetotransport measurements of the (1,1) to (2,0) charge transition in electrostatically defined bilayer graphene double quantum dots and argues that the two-particle spin and valley blockade can only be understood if the ten antisymmetric orbital states are included alongside the six symmetric ones. The central quantitative claim is that the valley-blocked detuning window is limited by the orbital splitting, the Kane-Mele spin-orbit splitting, and the valley g-factor of the antisymmetric orbital, while the spin-blocked window additionally depends on the short-range electron-electron splitting and the difference between symmetric and antisymmetric valley g-factors. A rate-equation simulation with one shared parameter set reproduces the main resonances and the magnetic-field-driven switch from valley to spin blockade in both bias directions. This matters because spin and valley qubit readout in bilayer graphene relies on Pauli blockade, so knowing exactly what limits the blockade identifies which splittings a device must engineer.","feed_headline":"Spin and valley blockade tied to orbital splitting","feed_subtitle":"Magnetotransport through 16 two-particle states ties blockade extent to orbital splitting and valley g-factor difference.","key_machinery":"The central object is the 16-state two-particle spectrum of the (2,0)/(0,2) configuration: six orbitally symmetric states, whose spin-triplet valley-singlet ground state is separated from the spin-singlet valley-triplet excited states by the short-range splitting $\\delta_2$, and ten orbitally antisymmetric states, separated from the symmetric manifold by the orbital splitting $\\Delta_{\\mathrm{Orb}}$ and internally by the Kane-Mele spin-orbit coupling $\\Delta_{\\mathrm{SO}}$. The load-bearing calculation is the resonance condition $\\varepsilon(B_\\perp) = E_{(0,2)}(B_\\perp) - E_{(1,1)}(B_\\perp)$, measured relative to the ground-state-to-ground-state baseline, which turns every transition into a line in the detuning-versus-field plane; the blockade is the detuning region below the first transition that requires a spin or valley flip. The supporting machinery is a Pauli rate-equation master equation over all 36 (0,1), (1,1), and (0,2) configurations, with equal tunnel rates modified by ad-hoc spin and valley flip penalties, which produces the simulated transport maps.","core_discovery":"On the paper's own terms, the discovery is that the blockade extents are governed by the orbital and interaction structure of the two-particle states: the valley-blockade extent is $\\epsilon_{\\mathrm{VB}} = \\Delta_{\\mathrm{Orb}} - \\Delta_{\\mathrm{SO}} - g^a_v \\mu_B B_\\perp$, while the spin-blockade extent is $\\epsilon_{\\mathrm{SB}} = \\Delta_{\\mathrm{Orb}} - \\Delta_{\\mathrm{SO}} - \\delta_2 + (g^s_v - g_s - g^a_v)\\mu_B B_\\perp$. Here $\\Delta_{\\mathrm{Orb}}$ separates the symmetric and antisymmetric two-particle orbital states, $\\Delta_{\\mathrm{SO}}$ is the Kane-Mele spin-orbit coupling, $\\delta_2$ is the short-range interaction splitting inside the symmetric manifold, and $g^s_v$, $g^a_v$, and $g_s$ are the symmetric-orbital valley, antisymmetric-orbital valley, and spin g-factors. The valley-blocked region is therefore limited only by properties of the antisymmetric orbital states, whereas the spin-blocked region depends on both orbital species. The authors identify individual measured resonances with transitions into the antisymmetric states and reproduce the field-dependent maps with a rate-equation simulation, concluding that the blockade switch and its extent are a direct readout of these splittings.","pith_inferences":["A direct test of the parameter set would be to fit the slopes of the individual resonances labeled a-d in the detuning-versus-field maps; the model predicts their crossings, which would pin down $g^a_v$ independently of the $\\varepsilon=0$ baseline uncertainty.","In a device with stronger interdot tunnel coupling, the (1,1) and (0,2) states hybridize, so the blockade extents should deviate from the two formulas; measuring that deviation would quantify when the independent-dots approximation breaks.","The same 16-state machinery likely carries over to three-carrier (1,2) to (0,3) transitions, where antisymmetric states could impose analogous limits on higher-order Pauli blockade; this is a prediction the authors did not test.","The unexplained conductance dips near $B_\\perp \\approx 0.05$-$0.1$ T and the valley blockade surviving past $B_{\\mathrm{TS}}$ in one data set suggest that field-dependent tunnel rates or g-factor renormalization, absent from the model, could be checked by repeating the detuning cuts at different barrier gate voltages."],"forward_implications":["The measured valley-blockade extent is a direct spectroscopic measure of $\\Delta_{\\mathrm{Orb}} - \\Delta_{\\mathrm{SO}} - g^a_v \\mu_B B_\\perp$, so one experiment yields the antisymmetric orbital splitting and its valley g-factor.","The field $B_{\\mathrm{TS}}$ where the valley blockade switches into a spin blockade marks the crossing of the spin-triplet valley-singlet and spin-singlet valley-triplet two-particle ground states, giving an independent handle on $\\delta_2$ and the valley g-factors.","For $g^s_v < g^a_v + g_s$, the spin-blockade window closes at the finite field $B_{\\mathrm{Orb}}$ where the relevant antisymmetric states become degenerate, so spin-based Pauli readout works only for $B_{\\mathrm{TS}} < B_\\perp < B_{\\mathrm{Orb}}$.","The symmetric-only six-state model cannot generate the observed resonances or the blockade extents, so any bilayer graphene qubit readout protocol must account for the antisymmetric orbital states in the (2,0)/(0,2) configuration.","The same parameter set reproduces both current directions at both triple points, indicating that including all 16 (1,1) states and all 16 (0,2) states gives a unified description of the blockade switch."],"supporting_citations":[{"why":"provides the two-electron multiplet spectrum of bilayer graphene quantum dots, including the (0,2) energies used for symmetric and antisymmetric orbital states.","marker":"[57]"},{"why":"supplies the tunneling theory and rate-equation treatment of single- and two-electron quantum dot states on which the simulation is built.","marker":"[58]"},{"why":"gives the general method for solving rate equations for electron tunneling through discrete quantum states.","marker":"[69]"},{"why":"provides the transition-rate matrix formalism used to find the stationary state of the master equation.","marker":"[70]"},{"why":"reports the earlier magnetic-field-tunable spin and valley Pauli blockade in bilayer graphene double quantum dots that this work extends.","marker":"[29]"},{"why":"established the particle-hole symmetry protected spin-valley blockade picture that fixes the interpretation of the blockade switch.","marker":"[25]"},{"why":"supplies the Kane-Mele spin-orbit coupling value used in the two-particle state energies.","marker":"[62]"},{"why":"provides the valley lifetime estimates used to argue that valley relaxation, not spin relaxation, limits the valley blockade.","marker":"[64]"},{"why":"catalogues Pauli blockade transitions in bilayer graphene double quantum dots and anchors the resonance assignment extended here.","marker":"[31]"}],"fun_headline_variants":["Orbital splitting sets spin and valley blockade extent in graphene dots","Blockade in graphene double dots limited by orbital structure","Antisymmetric orbitals govern blockade in two-electron graphene dots","Valley g-factor difference caps spin and valley blockade","Spin and valley blockade mapped to orbital and interaction splittings"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central formulas and resonance assignments hold only if the two dots can be treated as independent single-particle dots with negligible interdot tunnel coupling and equal tunnel probabilities for all states, so that each measured resonance is matched to one specific two-particle transition.","fun_headline_variants_meta":{"raw":{"variants":["Orbital splitting sets spin and valley blockade extent in graphene dots","Blockade in graphene double dots limited by orbital structure","Antisymmetric orbitals govern blockade in two-electron graphene dots","Valley g-factor difference caps spin and valley blockade","Spin and valley blockade mapped to orbital and interaction splittings"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1535,"prompt_tokens":975,"completion_tokens":560,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":479}},"tokens_in":591,"tokens_out":560,"duration_ms":5952,"temperature":1.0,"reasoning_tokens":479,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:56:15.661319+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a bilayer graphene double quantum dot with a clearly larger interdot tunnel coupling and record the detuning extent of the spin-blocked region versus perpendicular field; if it deviates from $\\epsilon_{\\mathrm{SB}} = \\Delta_{\\mathrm{Orb}} - \\Delta_{\\mathrm{SO}} - \\delta_2 + (g^s_v - g_s - g^a_v)\\mu_B B_\\perp$ beyond the estimated baseline uncertainty, or if the valley-to-spin switch occurs at a field inconsistent with the triplet-singlet crossing computed from the same parameters, the independent-dots premise is falsified.","supporting_citations":[{"cited_title":"M¨ oller, L","cited_arxiv_id":null,"evidence_quote":"provides the two-electron multiplet spectrum of bilayer graphene quantum dots, including the (0,2) energies used for symmetric and antisymmetric orbital states."},{"cited_title":"Knothe, L","cited_arxiv_id":null,"evidence_quote":"supplies the tunneling theory and rate-equation treatment of single- and two-electron quantum dot states on which the simulation is built."},{"cited_title":"Bonet, M","cited_arxiv_id":null,"evidence_quote":"gives the general method for solving rate equations for electron tunneling through discrete quantum states."},{"cited_title":"Timm, Random transition-rate matrices for the mas- ter equation, Phys","cited_arxiv_id":null,"evidence_quote":"provides the transition-rate matrix formalism used to find the stationary state of the master equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"reports the earlier magnetic-field-tunable spin and valley Pauli blockade in bilayer graphene double quantum dots that this work extends."},{"cited_title":"Banszerus, S","cited_arxiv_id":null,"evidence_quote":"established the particle-hole symmetry protected spin-valley blockade picture that fixes the interpretation of the blockade switch."},{"cited_title":"Kurzmann, Y","cited_arxiv_id":null,"evidence_quote":"supplies the Kane-Mele spin-orbit coupling value used in the two-particle state energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"catalogues Pauli blockade transitions in bilayer graphene double quantum dots and anchors the resonance assignment extended here."}],"review_version":1}