{"id":"b2c0f75f-d29b-4d02-9936-24dde9ee1d95","arxiv_id":"2607.05948","paper_version":1,"verdict":"CONDITIONAL","confidence":"UNKNOWN","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"The mean inner potential of alloyed and strained semiconductors is modeled as a linear interpolation of endpoint values rescaled by volume, matching DFT calculations to within ~2%.","lead":"The paper shows that the mean inner potential (MIP) of alloyed and strained semiconductors can be modeled by combining linear interpolation between endpoint MIPs with a simple volume rescaling. This matters for electron holography measurements of semiconductor devices, where MIP differences encode composition and strain information.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Model's success likely stems from muffin-tin sphere dominance (~94% of MIP for GaP), not the unverified MASA-for-forward-scattering derivation; this dominance is only checked for one material.","rationale":"The reader correctly identifies the MASA gap and the charge redistribution assumption as concerns. However, I partially disagree that these are the most load-bearing issues. The MASA derivation is unverified but also unnecessary—the paper's own explanation (muffin-tin dominance) provides a simpler and more direct justification that doesn't require MASAs. The DFT validation across three alloy systems with different bonding characters (ionic AlGaAs/InGaP and covalent GeSi) provides genuine independent support for the central claim, and the model clearly outperforms linear interpolation for systems with volume mismatch (Fig. 2b,c). Deviations are within experimental detection limits (0.1–0.4 V). The most important unverified aspect is not the MASA computation but the muffin-tin dominance explanation, which is only quantified for GaP. If the interstitial contribution is significantly larger for Ge or InP, the model's theoretical basis needs revision, even though the empirical fit might still hold. The ambiguity in the combined formula (Eqs. 3+4) for varying-volume alloys is a presentation gap that could affect reproducibility but doesn't undermine the claim if the correct formula was used (which the DFT agreement suggests). The verdict of CONDITIONAL is appropriate—the model is empirically validated but the theoretical justification has gaps. The condition should be: verify muffin-tin dominance across all endpoint materials and explicitly state the combined formula.","tokens_in":6724,"tokens_out":10033,"duration_ms":649261,"concrete_test":"Decompose the DFT-computed MIP into muffin-tin sphere and interstitial contributions for all six pure endpoint materials (AlAs, GaAs, InP, GaP, Ge, Si). If the interstitial fraction exceeds ~10% for any material (particularly Ge or InP), the muffin-tin transferability assumption underlying the model weakens, and the model's success for GeSi or InGaP would require an alternative explanation. Additionally, verify from the Zenodo data [17] that the red curves in Fig. 2 use the formula V = [x·Ω_A·V_A + (1-x)·Ω_B·V_B]/Ω(x) rather than a naive (Ω_ref/Ω(x))·[x·V_A + (1-x)·V_B]; if the latter was used and still matches DFT, the volume mismatch effect is smaller than expected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper derives Eqs. 3 and 4 from the MASA concept for forward scattering, stating twice 'Assuming that we would have computed MASAs for forward scattering' without actually doing so. This is a real theoretical gap, as the reader notes. However, the more load-bearing issue is that the paper's own explanation for why the model works points to a different mechanism: 'We attribute this to the significantly lower contribution of the interstitial regions to the MIP compared to the muffin-tin sphere regions (about 6% for GaP).' If ~94% of the MIP comes from muffin-tin spheres where the potential is atomic-like, linearity in concentration (Eq. 3) and 1/Ω scaling (Eq. 4) follow directly without invoking MASAs at all—the integral within each muffin-tin sphere is approximately transferable from pure materials. The MASA derivation is thus an unnecessary and unverified detour. The real question is whether muffin-tin dominance holds across all tested systems. It is only quantified for GaP (~6% interstitial). For Ge, which is more covalent and has different charge distribution characteristics, the interstitial fraction could be larger. If it exceeds ~10–15%, the assumption that muffin-tin potentials are transferable under alloying weakens, and the model's success for GeSi would need an alternative explanation. The DFT validation (Fig. 2) does show good agreement for GeSi, which is reassuring, but without the muffin-tin decomposition, we cannot confirm the mechanism. Additionally, the paper never explicitly writes the combined formula for alloys with varying volume. The correct derivation gives V_alloy(Ω(x)) = [x·Ω_A·V_A + (1-x)·Ω_B·V_B] / Ω(x), which requires first rescaling pure MIPs to a common volume before applying Eq. 3. A naive sequential application of Eqs. 3 then 4 with an arbitrary reference volume would give a different (incorrect) result for systems with significant lattice mismatch like GeSi (~13% volume difference). The good DFT agreement suggests the correct公式 is被","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This manuscript proposes a simple two-equation model for the mean inner potential (MIP) of alloyed and strained materials: Eq. (3) gives a linear combination of endpoint MIPs for alloying at constant volume, and Eq. (4) rescales the MIP by the volume ratio when the cell volume changes. The model requires no fitted parameters—only the MIPs and lattice parameters of the pure endpoint materials. It is validated against full DFT (WIEN2k, APW+LO) calculations for strained GaP (uniaxial, biaxial, hydrostatic), and for AlGaAs, InGaP, and GeSi alloys, with reported deviations of 1–2%, below typical electron-holography detection limits. The model is also tested against literature data for amorphous carbon density dependence and for biaxially strained GeSi.","tokens_in":6976,"tokens_out":2182,"duration_ms":173649,"significance":"The paper addresses a practical problem in electron microscopy: estimating MIPs of alloys and strained materials without performing full DFT calculations for every composition and strain state. The model is parameter-free and falsifiable, and the authors provide DFT validation across multiple material systems (III-V and group-IV semiconductors) and strain conditions. Data is deposited on Zenodo. The 1–2% agreement with full DFT is a useful result for the electron-holography community, where MIP values are needed for phase interpretation. The contribution is primarily empirical and practical rather than theoretically novel, but it is well-scoped for the intended audience.","major_comments":[{"comment":"The theoretical derivation of Eqs. (3) and (4) is stated to rely on the MASA (modified atomic scattering amplitudes) concept for forward scattering, but MASAs for forward scattering are never actually computed. The phrase 'Assuming that we would have computed MASAs for forward scattering' appears twice (items 1 and 2 in the key-idea list). This is a genuine gap in the theoretical justification. However, the authors' own explanation in the final paragraph—that linearity arises because 'the largest contribution to the MIP arises from the large Coulomb potential within the muffin-tin spheres and thus, becomes linear by definition'—is a more direct and self-contained justification that does not require MASAs at all. The authors should either (a) compute the forward-scattering MASAs to close the loop on the stated derivation, or (b) reframe the theoretical justification around the muffin-tin/","section":null},{"comment":"The muffin-tin dominance argument, which the authors identify as the actual reason the model works (final paragraph: 'about 6% for GaP'), is quantified only for GaP. The DFT validation covers AlGaAs, InGaP, and GeSi, but the interstitial fraction of the MIP is not reported for any of these systems. For GeSi in particular, which is more covalent, the interstitial contribution could differ. Since the muffin-tin transferability assumption is load-bearing for explaining why the model succeeds across all tested systems, the authors should provide the muffin-tin/interstitial decomposition for at least one additional system (ideally GeSi) to confirm that the mechanism is consistent. If the decomposition is not readily available, this limitation should be explicitly acknowledged.","section":null}],"minor_comments":[{"comment":"Eq. (2): the notation uses a superscript 'iso' on V_0 that is not defined in the surrounding text; please clarify that this refers to the isolated-atom approximation.","section":null},{"comment":"Fig. 1 caption: 'Linear relation are found' should read 'Linear relations are found'. Also, 'can be models reasonably' should be 'can be modeled reasonably'.","section":null},{"comment":"The text states maximum deviations of '1.1% (0.17 V), 0.8% (0.14 V) and 2.2% (0.28 V)' for the three strain types in GaP. It would help to state whether these are deviations of the model (Eq. 4) from the linear fit, or from the raw DFT data.","section":null},{"comment":"Table I lists elastic moduli only for Si and Ge, with dashes for AlAs, GaP, and InP. If these were not needed for the calculations presented, a footnote explaining this would avoid confusion.","section":null},{"comment":"The amorphous carbon comparison (final section before summary) compares the model's slope of 5.43 V/(g/cm³) to a refit value of 5.47 V/(g/cm³). The original fit value (5.20 V/(g/cm³)) is also given but the comparison would be clearer if all three values were tabulated.","section":null},{"comment":"Reference [13] and [14] appear to be conference proceedings (BIO Web Conf.); the authors may wish to confirm these are the most appropriate citations.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The reader's circularity concern (that Eqs. 3 and 4 follow 'by construction' from the MIP definition) is partially valid but overstated. While the functional forms are indeed simple consequences of volume-averaging, the non-trivial empirical content is that charge redistribution effects—which Kim et al. [3] showed can change the MIP by ~1.7 V—do not break these relations in practice. The DFT validation across multiple systems is what carries the paper, not the MASA derivation. I agree with the skeptic that the MASA framing is an unnecessary detour; the muffin-tin argument is more honest and direct. The paper is suitable for a methods/applied-physics venue if the theoretical framing is cleaned up."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful reading and constructive comments. Both major points are well-taken. We agree that the MASA framing is incomplete as presented and that the muffin-tin argument is the more self-contained justification. We will reframe accordingly. We also agree that the muffin-tin/interstitial decomposition should be provided for at least one additional system beyond GaP, and we will compute and report it for GeSi.","responses":[{"response":"The referee is correct. The MASA concept is invoked as motivation but is not carried through to an actual computation for forward scattering, and the conditional phrasing ('Assuming that we would have computed MASAs...') makes this gap explicit rather than hiding it. We agree that option (b) is the more appropriate fix: the muffin-tin dominance argument given in the final paragraph is self-contained and does not require the MASA framework at all. The key observation is that the MIP, defined as the volume average of the Coulomb potential in Eq. (1), receives its dominant contribution from the large Coulomb potential inside the muffin-tin spheres around each atomic site. Within these spheres, the potential is approximately atomic-like, so the integral over each sphere is approximately transferable between the pure material and the alloy. The interstitial contribution, which is where charge redistribution due to bonding would manifest, is small (about 6% for GaP). This directly explains why Eq. (3) (linear interpolation at constant volume) and Eq. (4) (volume rescaling) work: the total integral of the Coulomb potential is approximately conserved under alloying at fixed volume, and the MIP changes mainly through the 1/Ω prefactor when the volume changes. We will rewrite the theoretical justification section to lead with this argument and remove the conditional MASA framing from the key-idea list. We will retain a brief reference to the MASA concept as historical context for the idea that muffin-tin and interstitial contributions can be separated, but we will no longer present it as a load-bearing step in the derivation.","revision_made":"yes","referee_comment":"The theoretical derivation of Eqs. (3) and (4) is stated to rely on the MASA concept for forward scattering, but MASAs for forward scattering are never actually computed. The phrase 'Assuming that we would have computed MASAs for forward scattering' appears twice. The authors should either (a) compute the forward-scattering MASAs to close the loop, or (b) reframe the theoretical justification around the muffin-tin argument."},{"response":"This is a fair point. The transferability of the muffin-tin dominance argument across all tested systems is an assumption that we currently support with data from GaP only. GeSi is the most important test case because it is the most covalent system in our validation set, and one would expect a larger interstitial contribution there. We will compute the muffin-tin/interstitial decomposition of the MIP for both pure Ge and pure Si (and, if feasible, for an intermediate GeSi composition) using the same WIEN2k APW+LO framework. We expect the interstitial fraction to be somewhat larger than 6% for GeSi but still a minority contribution, which would be consistent with the model working well for this system. If the decomposition turns out to show a substantially larger interstitial fraction than expected, we will report this honestly and discuss its implications for the model's transferability. In any case, we will add a sentence to the revised manuscript explicitly stating the limitation that the decomposition has been checked for a limited set of systems.","revision_made":"yes","referee_comment":"The muffin-tin dominance argument is quantified only for GaP (~6% interstitial). The DFT validation covers AlGaAs, InGaP, and GeSi, but the interstitial fraction is not reported for these. For GeSi in particular, which is more covalent, the interstitial contribution could differ. The authors should provide the muffin-tin/interstitial decomposition for at least one additional system (ideally GeSi), or explicitly acknowledge the limitation."}],"tokens_in":6490,"tokens_out":1130,"duration_ms":70527,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The bottom line: this paper combines two known results — linear MIP interpolation in composition (Eq. 3) and 1/Ω volume rescaling (Eq. 4) — into a simple, parameter-free model for MIPs in alloyed and strained semiconductors. The model is validated against DFT with 1–2% deviations, below experimental detection limits. It is a practical shortcut for the electron holography community and deserves a serious referee, but the theoretical derivation has a real gap that should be addressed before publication. The DFT data and code are deposited on Zenodo, which is good practice and lends credibility to the validation. The model is genuinely useful: it lets practitioners estimate MIPs from endpoint values and lattice parameters without per-composition DFT. The strained GeSi biaxially stressed to Si is a nice additional test case. The amorphous carbon density-dependence check, while not the main point, is a reasonable sanity check. The authors deserve credit for shipping reproducible data and for testing across ionic (InGaP, AlGaAs) and covalent (GeSi) systems. Now the soft spots. The derivation invokes the MASA concept for forward scattering twice with the phrase 'Assuming that we would have computed MASAs for forward scattering' — but they never actually compute these. This is a genuine theoretical gap. However, I think the stress-test note correctly identifies that the MASA derivation is largely an unnecessary detour. The paper's own explanation for why the model works points to muffin-tin sphere dominance (~94% of MIP for GaP), and if the integral within each muffin-tin sphere is approximately transferable from pure materials, both Eq. 3 and Eq. 4 follow directly without invoking MASAs at all. The real question is whether muffin-tin dominance holds across all tested systems — it is only quantified for GaP. For GeSi, which is more covalent, the interstitial fraction could be larger. The DFT agreement for GeSi is reassuring but does not confirm the mechanism. This is a minor-to-moderate concern: the model works empirically, the explanation is incomplete, and the paper should either compute the forward-scattering MASAs or drop the MASA framing and instead argue from muffin-tin transferability directly. One more thing: the paper never explicitly writes the combined formula for alloys with varying volume. The correct derivation requires first rescaling pure MIPs to a common volume before applying Eq. 3, giving V_alloy = [x·Ω_A·V_A + (1−x)·Ω_B·V_B] / Ω(x). A naive sequential application of Eqs. 3 then 4 with an arbitrary reference volume would give a different result for systems with significant lattice mismatch like GeSi. The good DFT agreement suggests the authors did the right thing computationally, but the paper should state the combined formula explicitly. This paper is for the electron holography and semiconductor characterization community. It is a modeling convenience within an established framework, not a new theoretical advance, but it is a useful one. Recommend conditional acceptance pending revision of the theoretical justification — either compute the MASAs or reframe the argument around muffin-tin transferability, and write the combined formula explicitly.","headline":"Practical MIP model for alloys/strain: useful, validated, but theoretical justification is muddled","tokens_in":7626,"tokens_out":743,"would_cite":true,"duration_ms":146313,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Two-rule model predicts mean inner potential of alloys and strained crystals","keywords":["mean inner potential","alloy modeling","strain","volume rescaling","linear interpolation","density functional theory","electron holography","muffin-tin approximation"],"falsifier":"A material system where bonding-induced charge redistribution contributes comparably to or more than the atomic-core potential — e.g., a highly ionic alloy with large charge transfer — would show MIPs that deviate from the two-rule prediction by more than the ~2.2% maximum deviation seen here, falsifying the model's claim of universal applicability.","tokens_in":6836,"feed_emoji":"","tokens_out":1518,"duration_ms":70038,"temperature":0.7,"pith_summary":"The mean inner potential (MIP) — the average electrostatic potential inside a crystal, relevant to electron holography — can be predicted for alloyed and strained materials using only the MIPs and lattice parameters of the pure endpoint materials. The authors argue that two effects dominate: (1) at constant volume, the MIP of an alloy A_xB_{1-x} is a simple linear combination of the endpoint MIPs, and (2) when volume changes (due to lattice mismatch, strain, or density variation), the MIP rescales proportionally as the ratio of old to new volume. Combining these two rules reproduces DFT-computed MIPs for AlGaAs, InGaP, and GeSi alloys, for biaxially strained GeSi, for uniaxially/biaxially/hydrostatically strained GaP, and even for the density dependence of amorphous carbon, with deviations below ~2.2% — below the detection limit of current electron holography. The physical reason the model works is that the largest contribution to the MIP comes from the Coulomb potential deep within atomic cores (muffin-tin spheres), which is insensitive to bonding environment, while the integral of the potential stays approximately constant when the cell deforms, so the MIP changes mainly because the dividing volume changes.","feed_headline":"Two-rule model predicts mean inner potential of alloys and strained crystals","feed_subtitle":"Linear interpolation plus volume rescaling reproduces DFT results for three alloy systems and strained GaP, needing only endpoint data — no昂","key_machinery":"Mean inner potential (MIP): volume-averaged Coulomb potential of a crystal. Eq. 3: V_0(A_xB_{1-x}) = x*V_0(A) + (1-x)*V_0(B), linear interpolation at constant volume. Eq. 4: V'_0(Ω') = (Ω/Ω')*V_0(Ω), volume rescaling. Modified atomic scattering amplitudes (MASAs): DFT-derived scattering factors that account for charge redistribution, whose forward-scattering extension underpins Eq. 3. Muffin-tin spheres: atomic-core regions whose Coulomb potential dominates the MIP and is insensitive to bonding.","core_discovery":"The paper establishes that the MIP of a mixed or strained crystal factorizes into a composition-dependent linear interpolation (Eq. 3) and a volume-dependent rescaling (Eq. 4), and that this two-equation model captures the bowing in MIP-vs-concentration curves that pure linear interpolation misses. The central object is the MIP itself, defined as the volume average of the crystal Coulomb potential (Eq. 1), and the central mechanism is that this average is dominated by atomic-core contributions that are nearly invariant under alloying, leaving volume change as the primary correction. The authors validate this against full DFT calculations for three alloy systems (one with negligible lattice-m","pith_inferences":["If the model holds for ternary or quaternary alloys (e.g., InGaAsP), the MIP would be a weighted sum of endpoint MIPs with a single volume rescaling — a direct extension the paper does not test.","The model's assumption that charge redistribution is negligible relative to the isolated-atom baseline could break down for materials with strong charge transfer (e.g., highly ionic compounds or metals with delocalized electrons), where the muffin-tin-sphere dominance may not hold.","The ~0.76 V y-intercept found in prior amorphous-carbon data, which the model forces to zero, may indicate a systematic offset from surface dipoles or reference-potential conventions rather than a true density-independent contribution — worth separating experimentally.","For nanoscale objects where surface-to-volume ratios are large, the surface-potential contribution (noted as ~0.2 V by Kim et al.) could become comparable to volume-rescaling effects, setting a size scale below which the model needs correction."],"forward_implications":["Researchers using electron holography on alloyed or strained semiconductor devices can estimate MIPs without running full DFT supercell calculations, using only tabulated endpoint values and lattice parameters.","The model provides a consistency check: if measured MIPs deviate significantly from the two-rule prediction, the discrepancy may signal charge redistribution effects large enough to matter — a diagnostic rather than a nuisance.","The volume-rescaling rule (Eq. 4) is material-agnostic and could be tested on any system where density or pressure changes the MIP, including high-pressure phases or porous materials.","The finding that covalent and ionic alloys behave similarly under the model suggests that bonding character is secondary to atomic-core dominance for MIP, which constrains future charge-redistribution corrections."],"fun_headline_variants":["Alloy and strain effects on mean inner potential captured by two-rule model","Volume rescaling plus linear mixing reproduces alloyed and strained MIP curves","MIP of strained and alloyed crystals splits into composition and volume terms","Two-equation model captures MIP bowing across three alloy systems","Mean inner potential of alloys factors into linear mixing and volume scaling"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The model assumes that charge redistribution due to alloying, straining, or surface relaxation is small enough to neglect, so that the MIP is dominated by atomic-core potentials that are insensitive to bonding environment. If charge redistribution is not negligible, both the linear interpolation and the volume rescaling lose accuracy.","fun_headline_variants_meta":{"raw":{"variants":["Alloy and strain effects on mean inner potential captured by two-rule model","Volume rescaling plus linear mixing reproduces alloyed and strained MIP curves","MIP of strained and alloyed crystals splits into composition and volume terms","Two-equation model captures MIP bowing across three alloy systems","Mean inner potential of alloys factors into linear mixing and volume scaling","DFT validates two-rule MIP model for strained and alloyed crystals","Linear interpolation misses MIP bowing; volume rescaling fixes it","Composition term plus volume correction predicts alloyed and strained MIP"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":713,"prompt_tokens":558,"completion_tokens":155,"prompt_tokens_details":null},"tokens_in":558,"tokens_out":155,"duration_ms":8824,"temperature":1.0,"reasoning_tokens":42,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T20:02:04.931103+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"A material system where bonding-induced charge redistribution contributes comparably to or more than the atomic-core potential — e.g., a highly ionic alloy with large charge transfer — would show MIPs that deviate from the two-rule prediction by more than the ~2.2% maximum deviation seen here, falsifying the model's claim of universal applicability.","supporting_citations":[],"review_version":1}