pith. sign in

arxiv: 2601.14181 · v2 · submitted 2026-01-20 · ❄️ cond-mat.mtrl-sci

Faster grain-boundary diffusion with a higher activation enthalpy than bulk diffusion in ionic space-charge layers

Timon F. Kielgas , Roger A. De Souza This is my paper

Pith reviewed 2026-05-16 12:16 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords grain boundary diffusionspace-charge layersactivation enthalpyperovskite oxidescation vacanciesdefect associatesionic conductors
0
0 comments X

The pith

Grain-boundary cation diffusion in acceptor-doped perovskites can be faster than bulk diffusion even when its activation enthalpy is higher.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper demonstrates through continuum simulations that the ratio r of grain-boundary to bulk activation enthalpy for cation diffusion can exceed 1 in acceptor-doped AB O3 perovskites. This outcome arises when bulk diffusion is carried by a small concentration of faster neutral defect associates, while grain-boundary diffusion benefits from the strong accumulation of slower isolated cation vacancies inside negative space-charge layers. The neutral associates remain unaffected by the electrostatic potential, so the enhanced vacancy population at the boundaries produces faster overall transport despite the higher energy barrier per jump. A reader would care because earlier simulations had shown r approaching but not exceeding unity, leaving some experimental reports of r greater than 1 unexplained.

Core claim

Cation diffusion in these perovskites occurs by two related mechanisms: slower charged isolated cation vacancies and faster neutral associates of cation and anion vacancies. In the negative space-charge layers at grain boundaries the isolated vacancies are strongly accumulated while the neutral associates are unaffected by the potential. Continuum simulations that first solve Poisson’s equation for the space-charge profile and then solve the diffusion equation across a two-dimensional bicrystal show that, when bulk diffusion is dominated by the associates and grain-boundary diffusion is dominated by the enhanced isolated vacancies, the effective activation enthalpy for grain-boundary paths,

What carries the argument

The differential electrostatic response of charged isolated cation vacancies (accumulated in negative space-charge layers) versus neutral defect associates (unaffected by potential), implemented by sequential solution of Poisson’s equation and the diffusion equation in a bicrystal geometry.

If this is right

  • The ratio r can exceed unity when bulk transport is associate-mediated and grain-boundary transport is vacancy-mediated within the space-charge region.
  • The magnitude of r greater than 1 depends on the relative concentrations of the two defect species and on the strength of the space-charge potential.
  • Experimental observation of r greater than 1 requires doping levels that keep the associate fraction small in the bulk while still allowing measurable grain-boundary enhancement.
  • Measurement scatter in activation enthalpies can mask the distinction between r less than 1 and r greater than 1.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same two-mechanism picture could explain why r values vary across different acceptor-doped perovskites depending on dopant size and concentration.
  • Changing the dopant level or temperature range might allow experimental tuning between regimes where r is below or above unity.
  • The result suggests that space-charge layers can invert the usual expectation that grain boundaries lower the activation barrier for diffusion.

Load-bearing premise

Bulk diffusion is carried by a small concentration of faster neutral defect associates while grain-boundary diffusion is carried by the space-charge-enhanced concentration of slower isolated cation vacancies, with the associates remaining neutral and unaffected by the potential.

What would settle it

An experimental measurement that finds the concentration of neutral defect associates to be large in the bulk or that shows the associates themselves are repelled or attracted by the space-charge potential would eliminate the conditions needed for r greater than 1.

read the original abstract

Faster diffusion of cations along grain boundaries is reported in the literature for a variety of acceptor-doped $AB\mathrm{O}_{3}$ perovskite-type oxides. The ratio $r$ of the activation enthalpy of grain-boundary diffusion ($\Delta H^\mathrm{gb}$) to the activation enthalpy of bulk diffusion ($\Delta H^\mathrm{b}$) is seen experimentally to lie in the range $0.7 < r = \Delta H^\mathrm{gb} / \Delta H^\mathrm{b} < 1.3$, albeit with substantial errors. In a previous publication [Parras and De Souza, Acta Mater. 195 (2020) 383] it was shown through a set of continuum simulations that cation-vacancy accumulation within negative space-charge layers at grain boundaries in acceptor-doped perovskites will give rise to faster grain-boundary diffusion of cations, with the associated values of $r$ approaching but not exceeding unity. In the present study, we demonstrate by means of continuum simulations that $r > 1$ is possible for faster cation diffusion along grain boundaries in an acceptor-doped perovskite. The specific case we consider is cation diffusion occurring by two related mechanisms, by slower (charged) isolated cation vacancies and by faster (neutral) defect associates of cation and anion vacancies. Within the negative space-charge layers, the isolated cation vacancies are strongly accumulated, whereas the neutral associates are unaffected. We calculate diffusion profiles for a two-dimensional bicrystal geometry by solving, first, Poisson's equation, and subsequently, the diffusion equation. We find that, if a small concentration of faster defect associates is responsible for bulk diffusion, and a hugely enhanced concentration of slower isolated vacancies yields faster diffusion along space-charge layers, $r>1$ is obtained. The conditions under which $r > 1$ may be observed are described, and issues with experimental confirmation are discussed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript uses continuum simulations of a 2D bicrystal to show that the ratio r = ΔH^gb / ΔH^b for cation diffusion in acceptor-doped ABO3 perovskites can exceed unity. Bulk diffusion is assumed to be carried by a low concentration of fast neutral cation-anion vacancy associates, while grain-boundary diffusion is carried by the space-charge-enhanced population of slower isolated cation vacancies; Poisson's equation is solved first for the equilibrium potential, after which the diffusion equation is integrated with neutral associates remaining uniform.

Significance. If the result holds, it supplies a concrete defect-chemistry scenario that can produce the experimentally observed r values up to ~1.3, extending the authors' prior work in which r approached but did not exceed 1. The approach is standard, the governing equations are solved sequentially without circularity, and the outcome is directly traceable to the relative concentrations of the two defect species. The work therefore offers a falsifiable prediction for the conditions under which r > 1 should appear.

minor comments (2)
  1. [Results] The specific numerical values chosen for the bulk associate fraction and the space-charge enhancement factor of isolated vacancies (used to obtain r > 1) should be placed in the context of measured or calculated defect concentrations for typical acceptor-doped perovskites.
  2. [Discussion] The discussion of experimental confirmation would benefit from a short table or paragraph listing candidate materials and measurement techniques that could distinguish the proposed two-mechanism scenario from single-mechanism models.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, the recognition that our continuum simulations provide a falsifiable defect-chemistry scenario for r > 1, and the recommendation of minor revision. We have no major comments to address point by point, as none were specified in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives r > 1 by explicitly solving Poisson's equation for the equilibrium space-charge potential (yielding position-dependent isolated-vacancy concentrations) and then integrating the diffusion equation over a 2D bicrystal geometry. The key assumptions—bulk diffusion carried by a low concentration of fast neutral associates while GB diffusion is carried by the accumulated slower charged vacancies—are stated as inputs and produce the result directly; no parameter is fitted to the target ratio and then renamed a prediction. The citation to the authors' prior work supplies only the contrasting r ≤ 1 background case and is not required to establish the new simulation outcome. The derivation chain is therefore self-contained against the stated electrostatic and transport equations.

Axiom & Free-Parameter Ledger

2 free parameters · 3 axioms · 0 invented entities

The central claim rests on standard continuum electrostatics and diffusion applied to known defect species in acceptor-doped perovskites; the key free parameters are the relative concentrations chosen to produce the r > 1 regime.

free parameters (2)
  • relative concentration of neutral defect associates
    Assumed small enough that bulk diffusion is dominated by the associates rather than isolated vacancies
  • enhancement factor of isolated vacancies in space-charge layer
    Assumed large enough that the accumulated slow vacancies produce net faster grain-boundary diffusion
axioms (3)
  • standard math Poisson's equation governs the electrostatic potential and resulting defect distributions in space-charge layers
    Invoked to calculate accumulation of charged isolated vacancies
  • domain assumption Neutral defect associates carry zero net charge and are therefore unaffected by the space-charge potential
    Fundamental to why only isolated vacancies accumulate while associates remain constant
  • domain assumption Cation diffusion occurs by two parallel mechanisms with different rates and charge states
    Core premise required to obtain r > 1

pith-pipeline@v0.9.0 · 5653 in / 1556 out tokens · 77225 ms · 2026-05-16T12:16:36.313217+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages · 1 internal anchor

  1. [1]

    Kielgas, a Roger A

    - R X U Q D O 1 D P H Faster grain-boundary diffusion with a higher activation enthalpy than bulk diffusion in ionic space-charge layers Timon F. Kielgas, a Roger A. De Souza ∗a Faster diffusion of cations along grain boundaries is reported in the literature for a variety of acceptor- dopedABO 3 perovskite-type oxides. The ratiorof the activation enthalpy...

    work page 2020

  2. [2]

    Such values are physically inconsistent with the standard picture, and thus sug- gest that there is a phenomenon in oxides, but not in metals, that generatesD gb ≫D b. The one key difference between met- + P V S O B M / B N F < Z F B S > < W P M > 1–9 | 1 arXiv:2601.14181v1 [cond-mat.mtrl-sci] 20 Jan 2026 als and oxides, as pointed out by Kingery ,25 is t...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  3. [3]

    Note that for all but (d) Data for different∆Sa are plotted but overlay each other

    (d) bulk diffusion coefficientDb Sr from error function fit. Note that for all but (d) Data for different∆Sa are plotted but overlay each other. ues of∆S a and∆H a for the formation of (vSrvO)× associates; using defect diffusivities that show higher values for the associate than for the isolated vacancy; using a standard thermodynamic model of space-charg...

    work page 1976

  4. [4]

    4 R. Waser,J. Am. Ceram. Soc., 1989,72, 2234–2240. 5 I. Denk, F. Noll and J. Maier,J. Am. Ceram. Soc., 1997,80, 279–285. 6 R. A. De Souza, J. Fleig, J. Maier, O. Kienzle, Z. Zhang, W. Sigle and M. Rühle,J. Am. Ceram. Soc., 2003,86, 922–

    work page 1989

  5. [5]

    7 R. A. De Souza, C. Voisin, H. Schraknepper, M. Teusner, M. Kessel, P. Dufour, C. Tenailleau and S. Guillemet-Fritsch, Phys. Chem. Chem. Phys., 2014,16, 2568–2575. 8 R. A. De Souza,Phys. Chem. Chem. Phys., 2009,11, 9939–

    work page 2014

  6. [6]

    9 A. L. Usler, F. Ketter and R. A. De Souza,Phys. Chem. Chem. Phys., 2024,26, 8287–8298. 10 P. C. McIntyre,J. Am. Ceram. Soc., 2000,83, 1129–1136. 11 R. Waser,Solid State Ionics, 1995,75, 89–99. 12 R. Waser and R. Hagenbeck,Acta Mater ., 2000,48, 797–825. 13 A. L. Usler and R. A. De Souza,J. Electrochem. Soc., 2021, 168, 056504. 14 G. Garcia-Belmonte and ...

    work page 2024

  7. [7]

    Wærnhus, N

    16 I. Wærnhus, N. Sakai, H. Yokokawa, T. Grande, M.-A. Einarsrud and K. Wiik,Solid State Ionics, 2007,178, 907–914. 17 R. Sažinas, I. Sakaguchi, I. Hasle, J. M. Polfus, R. Haugsrud, M.-A. Einarsrud and T. Grande,Phys. Chem. Chem. Phys., 2017,19, 21878–21886. 18 R. Sažinas, I. Sakaguchi, M.-A. Einarsrud and T. Grande,AIP Adv., 2017,7, 115024. 19 M. Palcut,...

    work page 2007

  8. [8]

    33 R. A. De Souza and J. Maier,Phys. Chem. Chem. Phys., 2003, 5, 740–748. 34 R. A. De Souza, M. S. Islam and E. Ivers-Tiffée,J. Mater . Chem., 1999,9, 1621–1627. 35 A. Walsh, C. R. A. Catlow, A. G. H. Smith, A. A. Sokol and S. M. Woodley ,Phys. Rev. B, 2011,83, 220301. 36 H. J. Heelweg and R. A. De Souza,Phys. Rev. Mater ., 2021,5, 103804. 37 U. N. Gries,...

    work page 2003

  9. [9]

    45 F. A. Kröger and H. J. Vink,Solid State Physics, Academic Press, 1956, vol. 3, pp. 307–435. 46 T. Norby ,J. Korean Ceram. Soc., 2010,47, 19–25. 47 R. De Souza and G. Harrington,Nat. Mater ., 2023,22, 794–

    work page 1956

  10. [10]

    48 A. L. Usler,massaction, 2024,https://www.zenodo.org/. 49 K. Gömann, G. Borchardt, M. Schulz, A. Gömann, W. Maus- Friedrichs, B. Lesage, O. Kaïtasov, S. Hoffmann-Eifert and T. Schneller,Phys. Chem. Chem. Phys., 2005,7, 2053–2060. 50 R. Meyer, R. Waser, J. Helmbold and G. Borchardt,Phys. Rev. Lett., 2003,90, 105901. 51 J. Maier, G. Schwitzgebel and H.-J....

    work page 2024

  11. [11]

    Chung and B

    54 Y.-C. Chung and B. J. Wuensch,J. Appl. Phys., 1996,79, 8323–8329. 55 T. Bak, J. Nowotny and C. C. Sorrel,J. Phys. Chem. Solids, 2004,65, 1229–1241. 56 S. Koerfer, R. A. De Souza, H.-I. Yoo and M. Martin,Solid State Sci., 2008,10, 725–734. 57 P. Erhart and K. Albe,J. Appl. Phys., 2007,102, 084111– 084118. 58 S. Koerfer, B. Dißmann, N. Schier, H.-I. Yoo,...

    work page 1996

  12. [12]

    Bonkowski, M

    59 A. Bonkowski, M. J. Wolf, J. Wu, S. C. Parker, A. Klein and R. A. De Souza,J. Am. Chem. Soc., 2024,146, 23012–23021. 60 O. Schulz, M. Martin, C. Argirusis and G. Borchardt,Phys. Chem. Chem. Phys., 2003,5, 2308–2313. 61 I. V. Belova, G. E. Murch, D. Samuelis and M. Martin,Defect Diffus. Forum, 2007,263, 81–86. + P V S O B M / B N F < Z F B S > < W P M > 1–9 | 9

    work page 2024