{"id":"8b78267e-90ac-4bfb-88b1-2018403a495e","arxiv_id":"2607.13657","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Calcium doping of BaCeO3 increases proton uptake and CO2/SO2 resistance because the smaller, more electronegative Ca weakens A-O bonding and contracts the lattice, while slightly slowing proton rotation.","lead":"Using density-functional-theory calculations, the authors explain why replacing some barium with calcium in BaCeO3 improves both proton conduction and chemical stability: calcium's smaller size and higher electronegativity lower the energy needed to create oxygen vacancies, and they also make carbon dioxide and sulfur dioxide bind less strongly. This gives materials designers a simple rule — small-radius, high-electronegativity A-site dopants — for tuning proton-conducting fu","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Electronegativity-induced basicity mechanism is internally inconsistent: lower O Bader charges near Ca imply less basic oxide ions, yet the paper attributes more exothermic hydration to increased basicity.","rationale":"The reader's weakest assumption concerned the reliability of Bader charges (O3 anomaly, zero-flux artifacts). My concern is more fundamental: even if the Bader charges are accurate, they show the opposite trend required by the paper's basicity argument. Section C explicitly claims that Ca becomes more electron-rich and 'more basic,' enhancing proton uptake, but the proton acceptor in the hydration reaction is an oxide ion. Fig. 5a shows those oxide ions have reduced negative charge near Ca, i.e., lower basicity. The paper never reconciles this. This is an internal inconsistency in the central mechanism, not a numerical uncertainty. The more exothermic hydration in BCCO is a solid computational observation, but the causal attribution to Ca electronegativity is unsupported and, on standard chemical grounds, backwards. The title and abstract promote 'high-electronegativity A-site dopants' as the key design principle, so this flaw strikes at the headline claim. A revised paper that attributes the hydration enhancement to lattice contraction or hydrogen-bonding effects, and removes the electronegativity-basicity claim, might be salvageable, but the current manuscript's central mechanistic argument should not be accepted.","tokens_in":18880,"tokens_out":12588,"duration_ms":118797,"concrete_test":"Compute a direct oxide-basicity descriptor for the proton-accepting oxygens used in Section C: the O 2p band center (or proton affinity) at O1/O2 in BCCO and the equivalent O1/O3 in BCO from the same DFT calculations. If the O 2p center in BCCO is lower (more stable) or the proton affinity is smaller than in BCO, oxide basicity is decreased, directly contradicting the paper's mechanism. Alternatively, compute Bader charges on the protonated oxygens in the hydrated cells; if these charges are less negative in BCCO, the more exothermic hydration cannot be attributed to increased oxide basicity.","verdict_should_be":"REJECT","load_bearing_attack":"The central mechanistic claim (Abstract, Section C) is that Ca's higher electronegativity weakens A-O ionic bonding and reduces charge transfer, making the A-site more electron-rich and thereby increasing the basicity that drives more exothermic hydration (ΔE_hyd = -2.41 eV → -3.03 eV). But the paper's own Bader analysis (Fig. 5a) shows that oxygen atoms near Ca have lower (less negative) Bader charges than those near Ba, meaning the oxide ions are less electron-rich. Protons bind to oxide ions, not to A-site cations; the relevant basicity is that of oxygen. A less negative oxide ion is a weaker base for protonation, not a stronger one. The paper's appeal to the Ca Bader charge (8.55 vs 8.33, Fig. A4) conflates A-site electron richness with oxide basicity and does not resolve the contradiction. The more exothermic hydration energy is a valid DFT result, but the proposed electronegativity→basicity mechanism cannot explain it; the effect likely arises from lattice contraction and hydrogen-bonding geometry, not from Ca electronegativity. Since this mechanism is the paper's central explanatory claim, the argument fails at the conceptual level even if the computed energetics are correct.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses DFT (PBE+U, CI-NEB, AIMD) to investigate how A-site Ca substitution in BaCeO3 affects proton incorporation, diffusion, and chemical/thermal stability. It reports that Ca doping lowers oxygen-vacancy formation energies, makes water incorporation more exothermic, reduces CO2/SO2 adsorption strength, induces lattice contraction, and promotes Ca surface segregation. The authors attribute these effects to Ca's higher electronegativity (weakening A-O ionic bonding and increasing basicity) and smaller ionic radius (lattice contraction). They conclude that A-site Ca doping provides a dual strategy for improving both proton conductivity and stability in Ba-based electrolytes.","tokens_in":19056,"tokens_out":3973,"duration_ms":44027,"significance":"If the proposed mechanism were correct, this work would offer a clear design rule for A-site doping in proton-conducting perovskites: small-radius, high-electronegativity dopants simultaneously improve proton uptake and stability. The paper is computationally solid in its core energies and barriers, with standard settings, and it provides quantitative comparisons (e.g., ΔE_vac decreases, ΔE_hyd changes from -2.41 to -3.03 eV, CO2 adsorption shifts by 0.19 eV). These numbers are useful as numerical observations. The central explanatory claim, however, contains an internal inconsistency between the Bader-charge evidence and the proposed basicity mechanism, which is load-bearing for the abstract and summary. The thermal-stability conclusion also rests on a single 10 ps AIMD trajectory without statistical quantification. The paper's value is partly preserved by the raw energetics, but the mechanistic narrative needs substantial revision.","major_comments":[{"comment":"The claimed electronegativity→basicity→hydration mechanism is internally inconsistent with the paper's own Bader analysis. Fig. 5(a) shows that oxygen atoms near Ca (O1-O8) have lower Bader charges (less negative) than the same oxygens near Ba. Protons bind to oxide ions, so a less negative oxygen is a weaker base, not a stronger one. The appeal to the Ca Bader charge (8.55 vs 8.33) in Fig. A4 concerns the A-site cation, not the proton-accepting oxygen, and therefore does not explain the more exothermic ΔE_hyd. The -3.03 eV vs -2.41 eV hydration energy is a valid DFT result, but the paper's proposed basicity mechanism cannot explain it; lattice contraction or hydrogen-bond geometry may be the actual driver. This is a central claim of the abstract and Section C, so it must be corrected or substantially reframed with a mechanism consistent with the calculated charge redistribution.","section":"Section C; Eq. (2); Fig. 5(a); Fig. A4"},{"comment":"The conclusion that 'Ca doping further enhances the thermal stability' is based on a single 10 ps AIMD trajectory for each system. The paper reports that the amplitude of energy fluctuation in BCCO is 'significantly smaller' than in BCO, but no statistical analysis is provided: no standard deviations, no multiple independent trajectories, no structural order parameters, and no assessment of finite-size or thermostat effects. The energy-fluctuation criterion is not a robust measure of thermodynamic stability on this timescale. Either quantitative error bars should be provided, or the claim should be softened to a qualitative statement about lattice dynamics/phonon hardening, which is already supported by the phonon analysis.","section":"Section E; Fig. 11"},{"comment":"The central bond-weakening inference relies on Bader charges, but the analysis is not conclusive and contains an unexplained outlier. The O3 site shows the opposite Bader-charge trend, and the paper dismisses it as a 'zero-flux surface partitioning artifact' with no supporting evidence; reference [60] is a general Bader method paper, not a demonstration of this artifact in BaCeO3. Meanwhile, the ICOHP in Fig. 3(b) shows the Ca-O covalent overlap slightly increasing (0.12 to 0.14), which does not directly support 'weaker A-O ionic bonding.' The Bader charges indicate less negative oxygen charge, which is a statement about charge transfer, not directly about bond strength. More direct evidence—e.g., charge-density differences, projected COHP, or a full Bader analysis for all oxygen sites including O9-O11—is needed to establish the proposed electronegativity mechanism.","section":"Section B; Fig. 5(a); Appendix A4"}],"minor_comments":[{"comment":"The title uses 'A-Site' with a capital S; for consistency with the text it should be 'A-site'. Also, the abstract's phrase 'increasing the basicity' is ambiguous about whether the A-site or the oxide ion is meant.","section":"Title/Abstract"},{"comment":"The notation E_perfect is misspelled ('perf ect' in the text). Also, Eq. (2) uses Etot(VO(1)) and Etot(VO(2)) without specifying whether these are the same vacancy configuration or two distinct configurations; the definition is clear from the text but should be stated explicitly with the equation.","section":"Computational Method, Eq. (1)"},{"comment":"The figure caption and text refer to a 'CO3−2' group; this should be CO3^2− (or CO3^{2−}) to correctly indicate the carbonate ion.","section":"Fig. 10"},{"comment":"The sentence 'The reduced bond thereby decreases the E_vac' appears to have a missing word; it should probably read 'The weakened bond' or 'The reduced ICOHP magnitude' rather than 'reduced bond.'","section":"Section B, Fig. 5(b)"},{"comment":"A few references have incomplete metadata (e.g., ref. [59] lists only the journal and year without an article title), and ref. [91] is a data repository link that is acceptable but should be cited consistently.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a solid computational core and useful quantitative results, but the central mechanistic narrative has a logical inconsistency between the Bader-charge data and the proposed basicity argument. This inconsistency is fixable by reframing the mechanism, but it is a substantive issue rather than a cosmetic one. The thermal-stability claim also needs a quantitative basis. Given the paper's otherwise sound methodology and the practical importance of A-site doping in proton conductors, major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a qualified pass for me. The DFT machinery is standard and the energy landscape — vacancy formation, hydration, rotation/transfer barriers, CO2/SO2 adsorption, Ca segregation — is a useful systematic dataset for Ca-doped BaCeO3. The full ten-step migration pathway with separate rotation and transfer barriers is the most valuable piece; it convincingly shows Ca raises the rate-limiting rotation barrier and that the conductivity gain must come from higher proton concentration, not faster diffusion. That is an honest, non-obvious result.\n\nThe problem is the central mechanism. The abstract says higher electronegativity of Ca weakens A-O ionic bonding and enhances proton uptake by increasing basicity. But their own Bader analysis (Fig. 5a) shows the oxygen atoms nearest Ca have lower (less negative) charges than oxygen near Ba. That means the oxide ions are less electron-rich, i.e., weaker proton acceptors, not stronger. The paper tries to rescue this by pointing to the Ca Bader charge (8.55 vs 8.33) and the charge-density isosurface in Fig. A4, but protons bind to oxide ions, not to A-site cations. The hydration energy may indeed be more exothermic (-3.03 vs -2.41 eV), but the explanation cannot be the claimed increase in oxide basicity. The stress-test note gets this right. The authors need to remove or rewrite that causal chain — likely the effect comes from lattice contraction and hydrogen-bond geometry, which they already discuss elsewhere.\n\nOther soft spots are less severe but real: the thermal-stability claim rests on 10 ps AIMD energy-fluctuation amplitudes without statistical analysis; the grain-boundary resistance conclusion is inferred from a slab segregation energy, not from a grain-boundary calculation; and the modeled 12.5% Ca lies well above the experimental 1–5% optimum. None of these sabotage the core energetics, but they should be reported as extrapolations.\n\nBottom line: a competent, useful computational study with one load-bearing interpretive mistake. It deserves a serious referee if the authors can reconcile the Bader charges with a chemically sensible basicity argument or drop it. I would not desk-reject.","headline":"Solid DFT dataset on Ca-doped BaCeO3, but the electronegativity–basicity mechanism contradicts their own Bader charges and needs rethinking.","tokens_in":19637,"tokens_out":6642,"would_cite":true,"duration_ms":70002,"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":"This paper claims that substituting Ca for Ba on the A site of BaCeO3 simultaneously raises proton uptake and chemical stability, because Ca's higher electronegativity weakens A-O ionic bonding while its smaller radius contracts the lattice","keywords":["BaCeO3 perovskite","proton-conducting electrolyte","A-site doping","oxygen vacancy formation","hydration energy","proton diffusion barrier","chemical stability","first-principles calculations"],"falsifier":"Replace Ca with Sr in the same supercell and recompute oxygen-vacancy formation and hydration energies: Sr has a similar electronegativity shift but a smaller radius change than Ca, so the paper's mechanism predicts a weaker lattice-contraction effect but a comparable electronegativity effect. If the vacancy and hydration energies do not split accordingly, the dual-role attribution fails.","tokens_in":18659,"feed_emoji":"🔋","tokens_out":6597,"duration_ms":64969,"temperature":0.7,"pith_summary":"BaCeO3-based electrolytes conduct protons well but react with CO2, SO2, and moisture, and common B-site doping fixes one problem at the expense of the other. The paper argues that putting a small amount of calcium on the A site, replacing barium, gives both improvements at once. Calcium's higher electronegativity weakens the ionic bond between the A-site cation and oxygen, lowering the energy to create oxygen vacancies and making the material more basic; both effects make it easier for water to split and protons to enter the lattice. Calcium's smaller radius contracts the lattice, which stiffens the structure, improves thermal stability, lowers the proton-transfer barrier in transfer-limited systems, and helps keep positively charged oxygen vacancies away from grain boundaries. The calculated vacancy-formation, hydration, adsorption, and migration-barrier changes are all consistent with experiments showing better conductivity and stability in Ca-doped BaCeO3.","feed_headline":"Calcium doping boosts proton uptake and acid resistance","feed_subtitle":"A tiny Ca-for-Ba swap raises proton uptake and acid resistance — a design rule for fuel-cell ceramics.","key_machinery":"The carrying device is the A-site dopant itself, treated through two structural descriptors: electronegativity difference (Ca vs Ba) and ionic radius difference. Charge-density partitioning and orbital-overlap analysis are used to claim that Ca-O ionic bonding is weaker than Ba-O, making oxygen easier to remove and the A-site cation more electron-rich and basic; the smaller Ca radius lowers the tolerance factor, increasing octahedral tilting, which shortens the key hydrogen bond H···Of, lengthens H–Ce distances, and hardens phonon modes. Those structural descriptors then connect the same dopant to proton uptake, proton diffusion, grain-boundary resistance, and chemical and thermal stability.","core_discovery":"On the paper's own terms, the discovery is: in Ba0.875Ca0.125CeO3, every inequivalent oxygen site has a lower oxygen-vacancy formation energy than in BaCeO3 (the decrease growing with proximity to Ca), hydration becomes markedly more exothermic (−2.41 → −3.03 eV, and −2.83 → −3.53 eV when Y acceptor doping is included), and CO2/SO2 adsorption on the BaO-terminated (001) surface weakens by about 0.19 eV. The rotation barrier that limits long-range proton diffusion in BaCeO3 rises under Ca doping, while intra- and inter-octahedral transfer barriers fall; the author reasons that the lattice contraction shortens hydrogen bonds and elongates H–Ce distances. These numbers are interpreted through a","pith_inferences":["Editorial inference: A natural next calculation is a Sr-substituted series: Sr has a comparable electronegativity shift but a much smaller lattice-contraction effect than Ca, so it would separate the two proposed roles and map the design space.","Editorial inference: Because rotation is the bottleneck in BaCeO3, the stability gain from Ca could be combined with co-doping or strain that suppresses the rotation barrier, potentially recovering bulk diffusion without giving up acid resistance.","Editorial inference: The predicted 0.19 eV weakening of CO2 adsorption is directly testable by comparing CO2 temperature-programmed desorption from doped and undoped powders; the doped surface should desorb at lower temperature."],"forward_implications":["The same dual-role mechanism should operate for other small, high-electronegativity A-site cations in Ba-based perovskites, making A-site doping a general design route rather than a Ca-specific fix.","In perovskites where proton transfer, not rotation, limits diffusion, the same doping lowers the rate-limiting barrier, so conductivity gains would come from both concentration and mobility.","Surface segregation of Ca suppresses oxygen-vacancy accumulation at grain boundaries, predicting lower grain-boundary resistance and better performance in polycrystalline electrolytes.","Weakened carbonate and sulfate formation predicts longer chemical stability under CO2- and SO2-containing atmospheres, directly addressing the failure mode that blocks BaCeO3 use.","The lattice-contraction effect on thermal stability suggests that Ca doping can broaden the safe operating temperature window of proton-conducting ceramics."],"fun_headline_variants":["Ca doping lifts proton uptake, cuts acid adsorption","Ca swap: faster protons, weaker acid grip","Tiny Ca dopant strengthens perovskite fuel-cell ceramics","Electronegativity rule: Ca improves proton transport and stability"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The causal chain from electronegativity to weakened A-O bonding rests on a charge-density partitioning that the paper itself notes gives the opposite trend for one oxygen site; if that partitioning is unreliable for these mixed ionic/covalent bonds, the microscopic argument loses its evidence even though the calculated energies remain as numerical observations.","fun_headline_variants_meta":{"raw":{"variants":["Ca doping lifts proton uptake, cuts acid adsorption","Ca swap: faster protons, weaker acid grip","Tiny Ca dopant strengthens perovskite fuel-cell ceramics","Electronegativity rule: Ca improves proton transport and stability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000947,"raw_usage":{"total_tokens":3913,"prompt_tokens":809,"completion_tokens":3104,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":553,"completion_tokens_details":{"reasoning_tokens":3041}},"tokens_in":553,"tokens_out":3104,"duration_ms":26591,"temperature":1.0,"reasoning_tokens":3041,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T04:29:52.599417+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Replace Ca with Sr in the same supercell and recompute oxygen-vacancy formation and hydration energies: Sr has a similar electronegativity shift but a smaller radius change than Ca, so the paper's mechanism predicts a weaker lattice-contraction effect but a comparable electronegativity effect. If the vacancy and hydration energies do not split accordingly, the dual-role attribution fails.","supporting_citations":[],"review_version":1}