{"id":"b21bb836-12ab-4d26-921f-46180931a933","arxiv_id":"2507.12249","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In tensile-strained Ge/GeSn light-hole qubits, hyperfine coupling is dominated by the Fermi contact term via s-type conduction-band admixtures, and it grows with Sn barrier concentration.","lead":"A theoretical study of hole spins in GeSn/Ge/GeSn quantum dots calculates how strongly nuclear spins disturb the qubit, using atomistic tight-binding and DFT. The key finding is that the hyperfine coupling grows with tin content in the barrier and is driven by tiny electron-like s-orbital admixtures to the hole wavefunction.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The contact-dominance claim rests on bulk DFT hyperfine parameters applied to strained-alloy atoms and TB s* orbitals without a transferability test; an error there would rescale the Sn-dependent Overhauser effect.","rationale":"The reader's verdict is CONDITIONAL, and this stress-test does not move it. The central quantitative claim is genuinely supported by a parameter-free decomposition (Table II) and a monotonic trend (Fig. 3), and the hyperfine parameters are independently derived from DFT rather than fitted, so the work is not circular. The soft spot is the transferability of the bulk DFT contact parameters to the strained QW/alloy environment and to the TB s* basis function. This is load-bearing because Table II shows that removing the contact term collapses the Sn dependence of the perpendicular Overhauser fluctuation, so the contact term is the entire source of the claimed Sn sensitivity. A 10–20% error in |R_S(0)|^2 would shift the quantitative values, while a larger systematic error—especially if s* admixtures are parameterized with the ground-state s radial function—could change the qualitative conclusion. The proposed all-electron/strained-supercell check would settle this. The reader's weakest_assumption identifies the same issue, and the recommended CONDITIONAL verdict with a revision request for qualified wording and an explicit transferability test remains appropriate.","tokens_in":12714,"tokens_out":6366,"duration_ms":85957,"concrete_test":"Recompute the hyperfine parameters used in Table II with all-electron DFT in the actual environments: obtain |R_S(0)|^2 and M_alpha_beta for Ge in a biaxially strained supercell matching the QW (strain about 0.9–1.8%) and for Sn in an ordered Ge0.875Sn0.125 supercell, and separately compute the radial function of the excited s-like state rather than reusing the ground-state s radial function for the s* contact term. Also report the s versus s* decomposition of the total s-type admixture. If |R_S(0)|^2 changes by less than 5% and s* contributes less than 20% of the s-type admixture, the conclusion is robust; if either condition fails, the contact-dominance and Sn-dependence statements must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result—that the Fermi-contact term dominates the Overhauser fluctuations and produces the strong Sn-content dependence (Fig. 3b, Table II)—is carried by the product of the tight-binding s/s* admixture coefficients and the atomic contact parameter |R_S(0)|^2 from Table I. These parameters are obtained from DFT for bulk, unstrained Ge and Sn crystals (Sec. II D), then used unchanged for Ge atoms under biaxial tensile strain in the QW and for Sn atoms in the Ge1-xSnx alloy barrier. No uncertainty is reported, and transferability to the strained/alloy environment is asserted, not tested. In addition, the hyperfine matrix-element formula in Appendix A3 assigns the same |R_S(0)|^2 to every l=0 orbital, including the TB s* state; s* is a fictitious second s orbital, not the physical 4s/5s state whose radial function was extracted from DFT. If s* contributes substantially to the reported 1.4–6.4% s-type admixture, the contact term could be systematically overestimated. This is load-bearing: Table II shows that removing the contact term nearly eliminates the Sn dependence of the perpendicular Overhauser fluctuation (1.58 neV at both 10% and 15% Sn), so the contact term is the entire source of the claimed Sn sensitivity, not a small correction. A quantitative error in |R_S(0)|^2 or an incorrect treatment of s* would therefore directly undermine the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents atomistic tight-binding and DFT calculations of the lowest hole doublet in gate-defined GeSn/Ge/GeSn quantum wells. Strain is obtained with a valence-force-field model, the single-particle states are computed with an sp3d5s* tight-binding Hamiltonian plus a parabolic lateral confinement, and the hyperfine interaction is parameterized using all-electron wave functions reconstructed from PAW-DFT for bulk Ge and Sn. The authors evaluate Overhauser-field fluctuations and find that s-type orbital admixtures, induced by conduction-valence band mixing, allow the Fermi contact term to become the dominant hyperfine channel, with a strong increase of the fluctuations as the Sn content in the barrier grows. Dependences on quantum-well width and lateral confinement strength are also reported.","tokens_in":12963,"tokens_out":6860,"duration_ms":78527,"significance":"If the contact-mediated Sn dependence holds, the paper provides a valuable quantitative prediction for hole-spin qubits in GeSn, where the light-hole ground state and the strong spin-orbit tunability have generated recent interest. The strengths of the work include a fully atomistic treatment that goes beyond k.p models, a clear decomposition of hyperfine channels in Table II, and the use of DFT-derived radial functions near the nucleus instead of hydrogenic approximations. The predicted dependence of Overhauser fluctuations on barrier composition is falsifiable by element- or isotope-sensitive experiments, and the channel decomposition gives design guidance. The calculation is not circular: the tight-binding and DFT parameters are fixed externally and the Overhauser values are outputs. However, the quantitative reliability of the central claim rests on a few load-bearing approximations that need additional sensitivity tests.","major_comments":[{"comment":"The hyperfine matrix element in Appendix A3 assigns the same DFT-derived |R_S(0)|^2 to every l=0 orbital, including the s* orbital. Since s* is an empirical second s-like orbital in the sp3d5s* basis rather than the physical 4s/5s radial function, using the DFT value for s* is an uncontrolled approximation. The paper reports the s-type admixture as the sum of s and s* contributions, and Table II shows that switching off the contact term nearly eliminates the Sn-content dependence in the perpendicular Overhauser fluctuation. An error in the s* contact coupling would therefore directly bias the central claim. Please report the s and s* weights separately and test the sensitivity of the Overhauser fluctuations to how s* is treated, for example by setting its contact parameter to zero or using a separately computed radial overlap.","section":"Appendix A3"},{"comment":"The parameters |R_S(0)|^2 and M_alpha_beta are obtained from bulk, unstrained Ge and Sn crystals and then used without modification for Ge atoms under biaxial tensile strain in the quantum well and for Sn atoms in the GeSn alloy barrier. The central Sn-content trend is governed by the product of the tight-binding s-admixture coefficients and these atomic parameters, especially those of Sn. No uncertainty estimates or transferability tests are provided, and the statement that the PAW reconstruction 'should, in principle' give the correct wave function near the core is an assertion rather than a demonstrated validation. Please quantify how the hyperfine parameters change for Ge under representative tensile strain and for Sn in a GeSn supercell, and estimate the resulting uncertainty in the Overhauser values.","section":"Sec. II D and Table I"},{"comment":"The claim that the hyperfine interaction is 'dominated to a large extent' by the Fermi contact term is not uniformly supported by the data. At 10% Sn the full-model rms values are 2.06 neV (transverse) and 1.44 neV (z), while the no-contact values are 1.58 neV and 0.978 neV; the contact term is clearly dominant only at higher Sn content, where the 15% Sn values are 3.28/2.91 neV versus 1.58/1.32 neV. Please qualify the conclusions to the high-Sn regime, or provide a decomposition that isolates the contact contribution more directly, for example by switching the contact term on and off for s and s* separately.","section":"Sec. III B, Table II, and Conclusions"}],"minor_comments":[{"comment":"The units of |R_S(0)|^2 are given as Å^{-1/3}, but this quantity has dimensions of inverse volume (Å^{-3}). Please correct the typo and indicate whether atomic units are used.","section":"Table I"},{"comment":"The hybrid functional is not fully identified; the text states only that the Hartree-Fock exchange fraction was set to 0.19. Please specify which hybrid functional was used (e.g., PBE0 or HSE with a modified mixing parameter) and whether the fraction was tuned solely to the Ge band gap.","section":"Sec. II D"},{"comment":"The GeSn lattice constant row lists two values, '6.121*' and '6.071', with a footnote about the expanded value used only in the strain VFF simulation. Please clarify which value is used in the VFF calculation and which is used in the tight-binding Hamiltonian, and state the bowing relation used for the alloy lattice constant.","section":"Appendix A1, Table III"},{"comment":"Equation (2) is written without an explicit statistical average over the unpolarized nuclear-spin bath, although the text describes the result as a root-mean-square fluctuation. Please add the angle-bracket notation or otherwise clarify how the thermal average over nuclear spin configurations enters the expression.","section":"Eq. (2)"},{"comment":"The manuscript does not include a data/code availability statement, and many input parameters are deferred to Refs. [20,21]. Given the number of numerical settings (VFF parameters, tight-binding parameters, DFT settings), a self-contained table of the key input parameters would improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central mechanism is plausible. The main concerns are addressable with additional calculations: a sensitivity analysis for the s* contact term, a transferability test of the DFT hyperfine parameters to strained/alloyed environments, and a more qualified statement of contact dominance at low Sn content. I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this paper actually does something new. It computes the hyperfine interaction for a light-hole qubit in a Ge/GeSn quantum dot with an atomistic tight-binding model plus DFT-derived hyperfine parameters, while all prior work on this platform used k·p and ignored hyperfine. The central trend—Overhauser field fluctuations grow strongly with Sn content in the barrier—is internally consistent across the figures and table.\n\nThe paper does well in several respects. The decomposition in Table II is honest: it isolates the contact, d-shell, and Sn contributions. The finding that s-type admixtures (from conduction-valence mixing) open a Fermi contact channel for holes is physically interesting and likely correct in direction. The calculation is not circular—the hyperfine parameters come from independent DFT, not fitted to the Overhauser values.\n\nWhere it gets soft: the 'dominant channel' claim is too strong for low Sn content. At 10% Sn, removing the contact term reduces the perpendicular fluctuation from 2.06 to 1.58 neV—that is about a 23% reduction, not dominance. At 15% Sn it is about half. So the abstract and conclusion overstate the case for the lower end of the composition range. That is a wording issue, easily fixed.\n\nMore substantive is the transferability of the DFT-derived parameters. Table I lists |R_S(0)|^2 for bulk Ge and Sn, with no uncertainty and no test of whether those values hold for Ge under biaxial tensile strain or for Sn in the alloy. The stress-test note also raises the s* orbital: the same |R_S(0)|^2 is used for both s and s*, but s* is a fictitious TB basis state, not the physical 4s/5s radial function. If s* contributes significantly to the reported s-admixture, the contact term could be systematically overestimated. That is a legitimate concern, though not clearly fatal—the qualitative Sn dependence would likely survive, since the s-admixture itself grows with Sn. But the numbers are not as certain as the paper suggests.\n\nBottom line: this is a creditable computational study that deserves a serious referee. I would send it to peer review and ask for a revision that (a) qualifies the dominance claim per composition, (b) reports the s vs s* decomposition, and (c) adds at least a sensitivity check on the hyperfine parameters. For a reader working on hole-spin qubits, it is worth reading.","headline":"Solid first atomistic estimate of hyperfine in GeSn light-hole QDs, with an overstated dominance claim at low Sn and a real but fixable DFT-transferability caveat.","tokens_in":13577,"tokens_out":3641,"would_cite":true,"duration_ms":40021,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.21.La","71.70.Gm","71.15.Mb"],"model":"deepseek-v4-flash","headline":"In light-hole GeSn/Ge quantum dots, hyperfine decoherence is dominated by the Fermi contact interaction, with strength set by the tin fraction in the barriers.","keywords":["light-hole qubit","hyperfine interaction","Overhauser field","GeSn quantum well","tight-binding model","density functional theory","Fermi contact interaction","quantum dot"],"falsifier":"Recompute the Overhauser-field fluctuations using hyperfine parameters obtained from an independent all-electron method (a different DFT code, different functional, or quantum-chemistry calculation) with proper uncertainty estimates; if the contact term no longer dominates, the paper's central claim is falsified. Alternatively, a spin-echo or coherence-time measurement on a GeSn/Ge/GeSn gate-defined dot that tracks the noise floor as the barrier Sn content is varied from 10% to 20% would directly test the predicted growth of the coupling.","tokens_in":12461,"feed_emoji":"⚛️","tokens_out":5995,"duration_ms":68710,"temperature":0.7,"pith_summary":"This paper aims to establish what limits the spin quality of a new type of hole qubit: a gate-defined quantum dot in a tensile-strained Ge well sandwiched between GeSn barriers, where the ground state is a light hole. The authors calculate the nuclear-spin decoherence channel and conclude that it is governed not by the usual dipolar mechanisms but by the Fermi contact interaction, which becomes active only because the hole wave function acquires small s-type atomic orbital admixtures through conduction-valence band mixing. They quantify this by computing Overhauser field fluctuations with a fully atomistic sp3d5s* tight-binding model combined with first-principles hyperfine parameters, and find the coupling grows strongly with the Sn content in the barrier. This matters because it identifies tin composition as the practical knob controlling how much nuclear-spin noise this platform experiences.","feed_headline":"Tin content in barriers sets hole-spin noise in GeSn dots","feed_subtitle":"Atomistic simulations show a Fermi contact term, fed by s-orbital admixtures, dominates the nuclear-spin coupling as Sn rises.","key_machinery":"The machinery is the hyperfine interaction Hamiltonian evaluated in the basis of the two lowest hole states, with the Overhauser-field fluctuations computed from the rms of its components. The key object is the Fermi contact term $\\frac{8\\pi}{3}\\delta(\\boldsymbol{r})\\boldsymbol{S}$ embedded in the hyperfine operator; it is normally negligible for holes because valence orbitals are p-like, but here it becomes the dominant channel through small s-orbital (and s*-orbital) weights in the hole wave function. The paper isolates its contribution by turning off the contact term and the d-shell contributions in controlled model comparisons, and it derives the s-admixture weights from a realistic sp3d5s* tight-binding Hamiltonian with strain, spin-orbit coupling, and the parabolic gate potential.","core_discovery":"The central claim is that for the lowest hole doublet in an electrically defined GeSn/Ge/GeSn quantum dot, the hyperfine interaction is dominated to a large extent by the Fermi contact interaction mediated by s-type orbital admixtures, which arise from conduction-valence band mixing. This is demonstrated by computing Overhauser field fluctuations in a realistic sp3d5s* tight-binding model, with hyperfine matrix elements parameterized from DFT-calculated radial functions instead of hydrogen-like orbitals. The s-type admixture in the light-hole ground state rises from about 1.4% at 10% Sn to 6.4% at 20% Sn in the barrier, and the Overhauser fluctuations increase correspondingly. Systematically removing the contact term from the calculation reduces the fluctuations markedly, especially at 15% Sn, confirming that this channel is the leading source of nuclear-spin coupling.","pith_inferences":["If the contact term truly dominates, isotopic purification of 73Ge alone will not suppress the dominant channel, because the contact coupling acts through s-admixtures spread over the whole dot, including Sn sites; tuning the barrier composition would be the more direct lever.","The same DFT-parameterized contact mechanism should be checked in other tensile-strained group-IV systems, such as GeSn-on-Si quantum wells, where an even smaller band gap should produce larger s-admixtures and possibly a stronger hyperfine coupling.","There may be a design trade-off: higher Sn content helps confine the hole and moves Ge toward a direct band gap (useful for optical interfaces), but it also raises the nuclear-spin noise floor, and this paper's numbers allow that trade to be quantified."],"forward_implications":["Raising the Sn content in the barriers from 10% to 15% increases the transverse Overhauser-field fluctuation by about 60% and the z-component by roughly a factor of 2, so barrier composition is a direct lever on the nuclear-spin noise floor.","The dominant effect of higher Sn content is not the extra Sn nuclei themselves (they contribute only a few percent to the total) but the increased s-type admixture and stronger confinement that come with a deeper well.","The deviation of the transverse-to-z Overhauser ratio from the ideal light-hole value of 2 is largely a fingerprint of the isotropic contact term, which tends to equalize the two components.","Models that ignore conduction-valence band mixing will miss the leading hyperfine coupling channel in this system, so any quantitative prediction of decoherence times in GeSn devices must include that mixing.","The computed growth of Overhauser fluctuations with tin content provides a concrete, testable prediction for future spin-echo or coherence-time measurements on GeSn/Ge/GeSn dots."],"supporting_citations":[{"why":"Supplies the general hole hyperfine-interaction framework and the previous k·p treatment that the tight-binding calculation adapts and contrasts with.","marker":"[15]"},{"why":"Provides the atomistic tight-binding implementation of the hyperfine Hamiltonian and the Overhauser-field fluctuation procedure that this paper follows.","marker":"[20]"},{"why":"Supplies the sp3d5s* tight-binding model parameters and implementation details for GeSn nanostructures used throughout the simulations.","marker":"[21]"},{"why":"Gives the expression for Overhauser-field fluctuations in an unpolarized nuclear spin bath that the authors adopt for their rms calculation.","marker":"[12]"},{"why":"Supplies the nuclear spin quantum numbers, magnetic moments, and natural isotope abundances for Ge and Sn used in the hyperfine sums.","marker":"[13]"},{"why":"Describes the projector augmented-wave method that the paper uses to reconstruct all-electron wave functions from VASP pseudo-wave functions near the atomic core.","marker":"[37]"},{"why":"Establishes the role of d-orbital contributions to hole hyperfine coupling, which the authors separate out and compare against in their analysis.","marker":"[30]"}],"fun_headline_variants":["Sn content amplifies Fermi-contact spin noise in GeSn dots","s-orbital mixing drives hyperfine noise in GeSn hole qubits","Higher Sn content boosts hole-spin noise via contact term","Sn tunes nuclear coupling in strained GeSn light-hole dots","Fermi contact term dominates hole-spin decoherence in GeSn"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The DFT-derived hyperfine parameters, especially the contact density |R_S(0)|^2 for tin, are assumed to be quantitatively accurate and transferable from bulk crystals to the strained GeSn barrier and Ge well; if that number changes appreciably, the dominance of the contact term could disappear.","fun_headline_variants_meta":{"raw":{"variants":["Sn content amplifies Fermi-contact spin noise in GeSn dots","s-orbital mixing drives hyperfine noise in GeSn hole qubits","Higher Sn content boosts hole-spin noise via contact term","Sn tunes nuclear coupling in strained GeSn light-hole dots","Fermi contact term dominates hole-spin decoherence in GeSn"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000389,"raw_usage":{"total_tokens":2035,"prompt_tokens":917,"completion_tokens":1118,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":1038}},"tokens_in":533,"tokens_out":1118,"duration_ms":11501,"temperature":1.0,"reasoning_tokens":1038,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:50:13.972488+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Overhauser-field fluctuations using hyperfine parameters obtained from an independent all-electron method (a different DFT code, different functional, or quantum-chemistry calculation) with proper uncertainty estimates; if the contact term no longer dominates, the paper's central claim is falsified. Alternatively, a spin-echo or coherence-time measurement on a GeSn/Ge/GeSn gate-defined dot that tracks the noise floor as the barrier Sn content is varied from 10% to 20% would directly test the predicted growth of the coupling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the general hole hyperfine-interaction framework and the previous k·p treatment that the tight-binding calculation adapts and contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the atomistic tight-binding implementation of the hyperfine Hamiltonian and the Overhauser-field fluctuation procedure that this paper follows."},{"cited_title":"Del Vecchio and O","cited_arxiv_id":null,"evidence_quote":"Gives the expression for Overhauser-field fluctuations in an unpolarized nuclear spin bath that the authors adopt for their rms calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the nuclear spin quantum numbers, magnetic moments, and natural isotope abundances for Ge and Sn used in the hyperfine sums."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the projector augmented-wave method that the paper uses to reconstruct all-electron wave functions from VASP pseudo-wave functions near the atomic core."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the role of d-orbital contributions to hole hyperfine coupling, which the authors separate out and compare against in their analysis."}],"review_version":1}