{"id":"c8e05fa8-5d30-4987-a561-c3a8509b9944","arxiv_id":"2412.17750","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A unified EDL model for solid electrolytes with explicit core and space charge layers shows surface defect-energy changes substantially affect potential, capacitance, and parallel ionic conductivity.","lead":"This paper builds a computer model of the electrical double layer in solid electrolytes that treats both the atomically thin core layer and the surrounding space charge region. It shows that changes in defect formation energy at the surface can dominate potential drop, capacitance, and ionic conductivity, especially when the surface stabilizes defects.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sign inconsistency between Eq. (3) and Eq. (7) reverses the meaning of B; the central threshold B<−0.2 for substantial core conductivity is not reproducible as written.","rationale":"The reader correctly identified the constant-migration-energy assumption as a limitation and also noted the Eq. (3)/Eq. (7) sign inconsistency, but treated the migration assumption as the weakest point. I agree the migration assumption matters, yet it is explicitly flagged by the authors, scoped to small EDL potentials, and the paper itself labels the core conductivity contribution as a high estimate. The sign inconsistency is more load-bearing because it sits at the definition of the central control parameter B. The entire quantitative claim—core contribution 'as high as 50% especially when B < −0.2 eV'—depends on which sign of B lowers the surface DFE. If Eq. (7) is taken literally, then B < −0.2 eV corresponds to a higher surface DFE, and the observed enrichment/conductivity behavior would be reversed; the reported threshold would then be an artifact of labeling. The qualitative physics—low surface DFE enriches the core and boosts parallel conductivity—is plausible and independently supported by first-principles DFE trends, so the paper should not be rejected outright. But it cannot be accepted as stated until the sign convention is corrected in Eq. (7), the code is checked against the corrected equations, and the plotted B axis/Fig. 5 results are confirmed to match. This is a concrete, checkable reproducibility issue, and conditional acceptance with a required sign reconciliation is the appropriate verdict. My agreement with the reader is partial: the reader noticed the inconsistency but did not make it the central condition for acceptance.","tokens_in":10398,"tokens_out":4508,"duration_ms":47369,"concrete_test":"Inspect the MOOSE input and source files (Zenodo record 10.5281/zenodo.14538687) to determine the sign multiplying B in the implemented defect chemical potential, then rerun the Fig. 5 scenario at f = 0.5 eV with B = −0.4 eV and B = +0.4 eV. If the implementation follows Eq. (7) as printed, the simulation with B = +0.4 eV should show surface enrichment and a large core contribution r, while B = −0.4 eV should show depletion—directly contradicting the paper's B < −0.2 eV threshold. If the implementation instead follows Eq. (3), then Eq. (7) must be corrected and all B labels in the paper should be verified against the code before the quantitative claim is accepted.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that the core layer contributes substantially to the parallel ionic conductivity when the surface DFE is much lower than the bulk DFE, quantified as 'as high as 50% especially when B < −0.2 eV'. This claim is controlled by the sign convention for B, and the manuscript contradicts itself on that sign. Equation (3) defines the spatially varying standard chemical potential as µ0(x) = µb0 + B exp(−x/λc), and Fig. 1 (caption and plot) uses DFE = A + B exp(−x/λc), so B < 0 means a lower surface DFE. However, the governing equations (7a) and (7b) contain −B+ exp(−x/λc+) and −B− exp(−x/λc−), i.e., µ0(x) = µb0 − B exp(−x/λc). Under Eq. (7), B < 0 produces a higher surface DFE, the opposite of the stated physical setup. Since all simulations, including the conductivity maps in Fig. 5 and the r ≥ 50% statement, are presented as functions of the sign of B, a reader implementing Eq. (7) as written would obtain enrichment for B > 0 and depletion for B < 0, exactly reversing the reported threshold. The surrounding discussion makes it likely the intended convention is Eq. (3)/Fig. 1 and Eq. (7) contains a typographical sign error, but the displayed equations are the formal specification of the model. Until the sign is reconciled and the figures are regenerated or relabeled against the actual code, the quantitative boundary B < −0.2 eV is not a well-defined, reproducible result. The constant-migration-energy assumption, by contrast, is explicitly acknowledged and scoped by the authors as a high estimate, so it is a weaker objection than this internal inconsistency in the parameter that carries the paper's headline claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified one-dimensional continuum model of the electrical double layer (EDL) at solid electrolyte interfaces, explicitly treating both the core layer and the space charge layer. The core layer is represented by an exponentially decaying defect formation energy (DFE) profile, and the model includes site restriction and parameterized defect-defect interactions. The coupled Poisson and algebraic electrochemical potential equations are solved numerically (MOOSE) for two oppositely charged defects. The authors report potential and concentration profiles, surface charge, capacitance, and parallel ionic conductivity, with the central claim that the core layer can contribute up to ~50% of the parallel ionic conductivity, especially when the surface DFE is much lower than the bulk (B < -0.2 eV), regardless of defect interaction strength.","tokens_in":10812,"tokens_out":5250,"duration_ms":45715,"significance":"If the results are correct, they establish that space-charge-only analyses of interfacial ionic conductivity in solid electrolytes are incomplete, and that core-layer parameters (B and lambda_c) are actionable design levers. The paper has several strengths: the exponential DFE profile is motivated by the author's DFT calculations (refs 21, 22), the model includes defect-defect interactions and a concentrated-regime treatment, the numerical implementation is openly available on Zenodo, and the parametric sensitivity study over B and f is clearly framed. However, the central quantitative claim is compromised by an internal sign inconsistency in the governing equations, and the derivation of the interaction potential appears to contain a scaling error. These issues must be resolved before the paper's conclusions can be accepted.","major_comments":[{"comment":"The sign of B is inconsistent between the definition and the governing equations. Eq. (3) defines mu0(x) = mu_b0 + B exp(-lambda_c x), with the text and Fig. 1 stating that B < 0 corresponds to a lower surface DFE. However, Eqs. (7a) and (7b) contain -B+ exp(-x/lambda_c+) and -B- exp(-x/lambda_c-), which is equivalent to mu0 = mu_b0 - B exp(-x/lambda_c). Under this second convention, B < 0 produces a higher surface DFE, the opposite of the stated physical setup. The same -B convention appears in the Supporting Information (Eqs. 1a-1b). The conductivity maps in Fig. 5 and the reported threshold B < -0.2 eV for substantial core contribution are presented as functions of the sign of B, so a reader implementing Eqs. (7) as written would obtain enrichment for B > 0 and depletion for B < 0, reversing the main result. The authors must reconcile this sign (likely a typo in Eq. (7) and the SI), regenerate or relabel all figures accordingly, and verify the GitHub code against the corrected equations.","section":"Sec. 2, Eqs. (3) and (7)"},{"comment":"The derivation of the interaction chemical potential from finite-size scaling is not correct as stated. In standard finite-size scaling of point defects (refs 26, 27), the defect concentration in a supercell of side L scales as c ~ L^{-3}, not L^{-1}. Even if one accepts the paper's claim that c ~ L^{-1}, the free energy f_int ~ L^{-1} + L^{-3} would become f_int ~ c + c^3, so mu_int = df_int/dc would scale as 1 + 3c^2, which does not vanish as c -> 0 and is not proportional to c^2. The proposed form mu_int = f c^2/(1-c)^2 is therefore not a consequence of the cited first-principles calculations. Since the concentrated-regime predictions and the claim that the model is 'consistent with previous first-principles simulations' rest on this functional form, the authors should either correct the derivation, provide a different justification for the proposed form, or explicitly characterize it as a phenomenological choice.","section":"Sec. 2, Eq. (6)"},{"comment":"The central conductivity claim is based on the assumption that migration energies are constant and independent of position, as stated in the text: 'we assume that the migration energies are constant and independent of position in order to calculate the ionic conductivity parallel to the interface.' The authors acknowledge that this assumption is valid only when the EDL potential is small and cite refs 21 and 32, which show significant near-interface migration barrier variations. Nevertheless, the model is applied at potentials up to 1 V and at B values as low as -0.4 eV, and the abstract states without qualification that 'the core contributes substantially to the conductivity when the surface DFE is much lower than the bulk DFE.' Given that the conductivity prediction is the paper's principal new result, the authors should either restrict the claims to the regime where the constant-migration-energy assumption holds, include a sensitivity analysis with position-dependent migration energies, or clearly label the quantitative conductivity results as upper bounds in the abstract and conclusions.","section":"Sec. 3, ionic conductivity parallel to the interface"}],"minor_comments":[{"comment":"The sentence 'A key difference in observed in our EDL model' contains an extraneous 'in'; it should read 'A key difference observed in our EDL model'.","section":"Sec. 3, paragraph before Fig. 4"},{"comment":"The caption lists 'c+ c- c+ c-' in a way that is not immediately clear; labeling each panel with explicit titles or using separate axis labels would improve readability.","section":"Fig. 2 caption"},{"comment":"The standard theory name appears as 'Guoy-Chapman-Stern'; the conventional spelling is 'Gouy-Chapman-Stern'.","section":"Introduction, second paragraph"},{"comment":"The definition of r as the fraction of defects in the core layer integrates only c+ over the core region; the text should explicitly state that this is the fraction for the majority positive carrier and that the identification of r with the conductivity contribution assumes constant migration energies.","section":"Sec. 3, Eq. (9)"},{"comment":"The manuscript motivates the exponential DFE profile using the author's own DFT results (refs 21, 22) but does not show a direct comparison between the model profile and the DFT data for any material; adding such a comparison or reporting fitted B and lambda_c values would strengthen the 'consistent with first-principles' claim.","section":"Sec. 2, Fig. 1 and text"}],"recommendation":"major_revision","confidential_remarks":"The sign inconsistency between Eq. (3) and Eqs. (7) is likely a typographical error, but it is load-bearing for the main result and must be fixed before publication; the associated code should be checked against the corrected equations. The derivation of Eq. (6) also needs correction or reformulation. The constant-migration-energy issue is explicitly acknowledged by the authors, but the abstract overstates the conductivity result without the qualification. If these points are addressed, the manuscript could be a solid contribution to solid-electrolyte interface modeling."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Zeeshan Ahmad's paper is a real step forward for modeling electrical double layers in solid electrolytes. The genuinely new ingredient is treating the core layer as a continuous region with an exponential decay in defect formation energy, rather than as a single monolayer or just ignoring it. That is a physically motivated choice, backed by the author's own DFT results, and it changes the results in a way that matters: the core can carry a large share of the parallel ionic conductivity when the surface DFE is much lower than the bulk. If that holds up, it means space-charge-only analyses are missing a major contribution, and core-layer engineering becomes a design lever.\n\nThe paper also does some things right. The coupled Poisson + equilibrium constraint system is solved with a standard, well-posed numerical method (MOOSE), and the code and input files are on Zenodo. That makes the study reproducible, which is more than most modeling papers offer. The treatment of defect-defect interactions via a concentration-dependent energy, rather than a saturation constraint, is a reasonable choice and yields similar qualitative physics.\n\nThe soft spots are real. The biggest is the sign inconsistency for B. Equation (3) defines mu0(x) = mu_b0 + B exp(-x/lambda_c), so B < 0 means a lower surface DFE. But the governing equations, Eq. (7), contain -B exp(-x/lambda_c), which flips the meaning. A reader implementing Eq. (7) as written would get the opposite of the reported enrichment/depletion behavior, and the quantitative threshold 'B < -0.2 eV' for 50% core conductivity would not be reproducible. I think it's a typographical slip, not a conceptual error — the text and Fig. 1 are consistent with the plus sign — but it is a load-bearing typo that sits in the parameter driving the paper's main claim. It has to be reconciled, and the figures regenerated or relabeled to match the actual code.\n\nThe constant-migration-energy assumption is a limitation, but the authors flag it themselves and call the core contribution a high estimate. That is honest and scoped. The interaction form is a modeling choice with plausible physical justification, but it is not derived from first principles. That is fine for a parameter-sensitivity study.\n\nOverall, the paper is useful, mostly clear, and reproducible once the sign issue is resolved. I would not cite it in its current form, but I would send it to peer review with a request to fix the sign and double-check the conductivity maps. It deserves a careful referee, not a desk reject.","headline":"A genuinely new EDL treatment for the core layer in solid electrolytes, with a reproducible implementation — but a sign inconsistency in the central parameter B guts the headline claim until it is fixed.","tokens_in":11337,"tokens_out":3640,"would_cite":false,"duration_ms":32403,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A unified electrical double layer model shows the core layer can contribute up to half of ionic conductivity in solid electrolytes.","keywords":["electrical double layer","solid electrolyte","space charge layer","core layer","defect formation energy","ionic conductivity","defect-defect interactions","interfaces"],"falsifier":"Recompute the conductivity integral with first-principles migration-energy profiles (which the paper cites as showing significant near-interface variation) or with discrete lattice sites; if the core's share of the integrated carrier density drops below a few percent for $B < -0.2$ eV, the central claim would be falsified. Alternatively, measure the capacitance peak shift and core contribution in a system with engineered surface DFE and compare with the predicted $B$-dependence.","tokens_in":10174,"feed_emoji":"🔋","tokens_out":6928,"duration_ms":54817,"temperature":0.7,"pith_summary":"This paper argues that the standard way of modeling electrical double layers at solid electrolyte interfaces is incomplete: it treats only the space charge layer and ignores the core layer, the atomic-scale region where defect formation energies are altered by the interface. The author builds a model that includes both layers, with the core's defect formation energy varying exponentially with distance, and defect–defect interactions for concentrated conditions. The central finding is that when the surface defect formation energy is much lower than the bulk value, the core layer can account for as much as half of the ionic conductivity parallel to the interface. If this is right, space-charge-only analyses of grain boundary and composite-electrolyte conductivity are missing a major contribution, and engineering the core layer becomes a practical design lever for solid-state batteries and similar devices.","feed_headline":"Solid electrolyte core layer can carry up to 50% of conductivity","feed_subtitle":"Space-charge-only models miss this, so core-layer engineering becomes a real lever.","key_machinery":"The central object is a modified electrochemical potential for each charged defect, $\\mu(x) = \\mu_0(x) + ze\\phi + \\mu_{\\rm int}(c) + kT\\ln[c/(1-c)]$, where the standard chemical potential (the defect formation energy) varies in the core as $\\mu_0(x) = \\mu_b^0 + B\\exp(-x/\\lambda_c)$, and $\\mu_{\\rm int}(c) = f c^2/(1-c)^2$ accounts for defect–defect interactions. Setting $\\mu=0$ at equilibrium (the local chemical potential formalism) and coupling to Poisson's equation gives the potential and concentration profiles; the conductivity parallel to the interface is then taken proportional to the integral of the mobile defect concentration. This machinery lets the author switch from dilute to concentrated regimes and separate core from space charge contributions.","core_discovery":"The central claim is that the electrical double layer in a solid electrolyte must be treated as two coupled regions — a core layer, where the defect formation energy changes smoothly and often substantially from its bulk value, and a space charge layer, where defect concentrations are perturbed — and that doing so changes predictions of capacitance and conductivity. Specifically, the model predicts that when the defect formation energy at the surface is lower than in the bulk ($B < -0.2$ eV), the core layer contributes up to 50% of the total ionic conductivity parallel to the interface, regardless of defect interaction strength. Consequently, conclusions about interface conductivity drawn from space charge layer analysis alone are inaccurate for such interfaces. The core layer also dominates the potential drop: for $B = -0.4$ eV, more than 60% of the potential drop across the EDL occurs within the core.","pith_inferences":["If position-dependent migration barriers were included, the core's conductivity share could shrink or grow, so the 50% figure should be read as an upper bound in the presence of barrier variations.","The same framework could be extended to perpendicular transport and to heterogeneous interfaces where different defect types (kinks, adions) appear in the core, replacing the exponential DFE form with first-principles data.","The model's prediction that core effects activate at dilute bulk concentrations suggests that nominally dilute grain boundaries may behave as concentrated systems, potentially reconciling some discrepancies in measured capacitance.","A direct experimental test would be to engineer the surface DFE (for example by doping or coating) and check whether the conductivity and capacitance respond as predicted with $B$."],"forward_implications":["Space-charge-only models underestimate interface conductivity whenever the surface defect formation energy is lower than the bulk value; the core must be included.","Core-layer properties ($B$ and $\\lambda_c$) are design parameters: lowering the surface DFE and tuning the core length can raise conductivity parallel to the interface.","The model predicts capacitance–voltage curves whose maximum shifts to lower EDL potential when $B<0$, offering a signature that can be compared with experiments.","Defect–defect interactions reduce surface charge and capacitance at high EDL potentials but do not change the potential at which capacitance peaks.","High core defect concentrations can occur even at dilute bulk concentrations, so concentrated-solution effects matter at interfaces even in nominally dilute electrolytes."],"supporting_citations":[{"why":"Supplies the first-principles defect-formation-energy profiles used for the core layer.","marker":"[21]"},{"why":"Extends the same DFE variation data to metal halide perovskites, supporting the exponential form.","marker":"[22]"},{"why":"Provides the local thermodynamic formalism (zero electrochemical potential at equilibrium) that the model is built on.","marker":"[23]"},{"why":"The classical space-charge conductivity framework and constant-migration-energy assumption that this work extends to include the core.","marker":"[10]"},{"why":"The existing concentrated-EDL model for solid interfaces whose treatment of defect interactions this paper contrasts with.","marker":"[6]"},{"why":"Finite-size scaling results motivating the defect-interaction energy form.","marker":"[26]"},{"why":"First-principles treatment of point-defect interactions used to justify the concentration-dependent interaction term.","marker":"[27]"},{"why":"Prior space-charge conductivity analysis whose scope (space charge only) this model generalizes.","marker":"[24]"}],"fun_headline_variants":["Core layer carries up to half of solid electrolyte conductivity","Space-charge-only models miss core layer's 50% conductivity boost","Core layer dominates potential drop in solid electrolytes","Unified EDL model: core layer shapes conductivity and capacitance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole performance claim rests on treating the ion's migration barrier as position-independent, so that conductivity is proportional to the integrated defect concentration, an assumption the author concedes is valid only for small potentials.","fun_headline_variants_meta":{"raw":{"variants":["Core layer carries up to half of solid electrolyte conductivity","Space-charge-only models miss core layer's 50% conductivity boost","Core layer dominates potential drop in solid electrolytes","Unified EDL model: core layer shapes conductivity and capacitance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000525,"raw_usage":{"total_tokens":2523,"prompt_tokens":917,"completion_tokens":1606,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":1539}},"tokens_in":533,"tokens_out":1606,"duration_ms":9543,"temperature":1.0,"reasoning_tokens":1539,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:13:27.827111+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the conductivity integral with first-principles migration-energy profiles (which the paper cites as showing significant near-interface variation) or with discrete lattice sites; if the core's share of the integrated carrier density drops below a few percent for $B < -0.2$ eV, the central claim would be falsified. Alternatively, measure the capacitance peak shift and core contribution in a system with engineered surface DFE and compare with the predicted $B$-dependence.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the first-principles defect-formation-energy profiles used for the core layer."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the same DFE variation data to metal halide perovskites, supporting the exponential form."},{"cited_title":"Local Thermodynamic Formalism for Space Charge in Ionic Crystals","cited_arxiv_id":null,"evidence_quote":"Provides the local thermodynamic formalism (zero electrochemical potential at equilibrium) that the model is built on."},{"cited_title":"Ionic Conduction in Space Charge Regions","cited_arxiv_id":null,"evidence_quote":"The classical space-charge conductivity framework and constant-migration-energy assumption that this work extends to include the core."},{"cited_title":"W.; Swift, J","cited_arxiv_id":null,"evidence_quote":"The existing concentrated-EDL model for solid interfaces whose treatment of defect interactions this paper contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Finite-size scaling results motivating the defect-interaction energy form."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First-principles treatment of point-defect interactions used to justify the concentration-dependent interaction term."},{"cited_title":"D.; Cai, W.; Chueh, W","cited_arxiv_id":null,"evidence_quote":"Prior space-charge conductivity analysis whose scope (space charge only) this model generalizes."}],"review_version":1}