Redox behaviour of Fe impurities in BaTiO$_3$ based on many-body calculations

Hamza Zerdoumi; Hao Wang; Hongbin Zhang; Ruiwen Xie; Zhiyuan Li

arxiv: 2605.24747 · v1 · pith:YKYQI36Qnew · submitted 2026-05-23 · ❄️ cond-mat.mtrl-sci

Redox behaviour of Fe impurities in BaTiO₃ based on many-body calculations

Zhiyuan Li , Hamza Zerdoumi , Hao Wang , Ruiwen Xie , Hongbin Zhang This is my paper

Pith reviewed 2026-06-30 12:46 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci

keywords Fe impuritiesBaTiO3mixed valenceoxygen vacanciescharge transferredox behaviorDFT many-bodyXPS

0 comments

The pith

Fe impurities in BaTiO3 exist mostly as mixed Fe2+ and Fe3+ states whose ratio shifts with oxygen vacancies through d-orbital charge transfer.

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

The paper uses DFT-based many-body calculations to map the electronic states of iron atoms placed inside barium titanate. It establishes that these impurities spend most of their time in a combination of Fe2+ (3d6) and Fe3+ (3d5) electron configurations rather than a single fixed valence. Oxygen vacancies move the balance toward Fe2+ by altering the energy of specific orbitals. The driving process is charge transfer triggered when the d3z2-r2 orbitals shift relative to the rest of the 3d manifold. Calculated x-ray photoemission spectra from the resulting local Hamiltonian line up with measured data, indicating the valence mixture is observable.

Core claim

Iron impurities in BaTiO3 exhibit a mixed valence nature, comprising mostly Fe2+ (3d6) and Fe3+ (3d5) configurations, and such configurations can be tuned via oxygen vacancies which favor Fe2+. The origin of such a redox behavior can be attributed to the charge transfer caused by shifting of the d3z2-r2 orbitals. X-ray photoemission spectroscopy evaluated by solving the Wannier function-derived local atomic Hamiltonian using the crystal field multiplet approach shows good agreement with recent experimental measurements.

What carries the argument

Shifting of the d3z2-r2 orbitals that produces charge transfer between the Fe 3d states and the surrounding lattice.

If this is right

Oxygen vacancies can be introduced to adjust the average valence of Fe impurities toward the divalent state.
The mixed valence produces distinct spectroscopic signatures that can be calculated from the local Hamiltonian.
The same orbital-shift mechanism governs the redox response when Fe is substituted into the perovskite lattice.
Similar vacancy-induced valence changes are expected in other transition-metal-doped oxide hosts where crystal-field splitting is comparable.

Where Pith is reading between the lines

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

Controlling oxygen content during growth or annealing could serve as a practical lever to set the Fe valence fraction for desired electrical or magnetic response.
The orbital-specific charge-transfer channel identified here may appear in other d5/d6 impurity systems once the local symmetry and vacancy positions are comparable.
Time-resolved or in-operando spectroscopy under varying oxygen partial pressure would provide a direct test of the predicted vacancy-driven shift.

Load-bearing premise

The DFT-based many-body methods plus crystal-field multiplet treatment correctly capture the relative energies of the Fe configurations and the orbital shifts without large errors from functional choice or basis construction.

What would settle it

Direct experimental measurement of Fe 2p core-level spectra or valence-band features in BaTiO3 samples with systematically varied oxygen-vacancy concentrations that show no increase in Fe2+ weight when vacancies are added would falsify the predicted tuning mechanism.

read the original abstract

Based on detailed electronic structure and spectroscopy obtained using DFT-based many-body techniques, the redox behavior of Fe impurities in BaTiO$_3$ is investigated. It is observed that Fe impurities exhibit a mixed valence nature, comprising mostly Fe$^{2+}$ ($3d^6$) and Fe$^{3+}$ ($3d^5$) configurations, and such configurations can be tuned via oxygen vacancies which favor Fe$^{2+}$. The origin of such a redox behavior can be attributed to the charge transfer caused by shifting of the $d_{3z^2-r^2}$ orbitals. Furthermore, x-ray photoemission spectroscopy is evaluated by solving the Wannier function-derived local atomic Hamiltonian using the crystal field multiplet approach, with good agreements with recent experimental measurements.

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.

Desk Editor's Note private letter to a colleague

The paper pins Fe redox in BaTiO3 on d_{3z^2-r^2} shifts producing mixed 2+/3+ valence tunable by oxygen vacancies, but the result sits on untested DFT choices.

read the letter

The main point is that Fe impurities in BaTiO3 sit in a mixed Fe^{2+}(3d^6)/Fe^{3+}(3d^5) state that oxygen vacancies push toward Fe^{2+}, with the shift traced to the d_{3z^2-r^2} orbital moving relative to the oxygen bands. The authors compute this with DFT-based many-body methods, extract a local Hamiltonian via Wannier functions, and feed it into a crystal-field multiplet calculation that reproduces recent XPS data.

What is actually new is the specific orbital attribution plus the quantitative populations under vacancy conditions. Earlier studies on Fe-doped BaTiO3 exist, but this level of detail on which d orbital drives the charge transfer and how vacancies alter the balance does not appear in the cited literature.

The work applies standard tools to a practical material and shows decent spectral agreement, which is useful for anyone modeling defects in perovskites. The multiplet step for XPS is a straightforward way to connect the electronic-structure output to experiment.

The soft spot is the one flagged in the stress test. The d-manifold alignment that produces the reported charge transfer is known to move by 1-2 eV with functional choice, U value, or Wannier downfolding. The abstract gives no sign of testing alternatives, so it is unclear whether the dominant configurations or the vacancy preference would survive a different setup. If the full manuscript contains those checks, the concern shrinks; otherwise it remains load-bearing.

This is for people working on impurity states or defect chemistry in titanates. A reader focused on BaTiO3 or on similar many-body-plus-multiplet workflows will find concrete numbers to compare against. It is worth sending to peer review because it supplies a detailed, falsifiable picture for a common system even if the robustness needs verification.

Referee Report

2 major / 1 minor

Summary. The manuscript uses DFT-based many-body techniques to study Fe impurities in BaTiO3. It reports that Fe exhibits mixed valence dominated by Fe^{2+}(3d^6) and Fe^{3+}(3d^5) configurations that can be tuned by oxygen vacancies (favoring Fe^{2+}), with the redox behavior originating from charge transfer driven by shifts in the d_{3z^2-r^2} orbital. XPS spectra are computed by solving the Wannier-derived local atomic Hamiltonian in the crystal-field multiplet approach and show good agreement with recent experiments.

Significance. If robust, the work supplies a microscopic picture of impurity valence and defect tuning in a prototypical perovskite, relevant to dielectric and ferroelectric applications. The forward-computational approach (no post-hoc fitting to the compared XPS data) and the direct link from many-body electronic structure to multiplet spectra are methodological strengths that increase the result's value if the charge-transfer energetics prove insensitive to technical choices.

major comments (2)

[electronic-structure and many-body methods sections] The central attribution of redox behavior to d_{3z^2-r^2} orbital shift and resulting charge transfer (abstract and results on valence configurations) rests on the alignment of Fe d states relative to O p bands. No tests of sensitivity to DFT functional (LDA/PBE/hybrid), Hubbard U value, or Wannier downfolding procedure are reported; these choices are known to move the d-manifold position by 1-2 eV and can therefore alter the dominant 3d^5 vs. 3d^6 weights or the vacancy preference.
[spectroscopy evaluation section] The crystal-field multiplet XPS calculation (section on spectroscopy) inherits the local Hamiltonian from the preceding DFT+Wannier step. Without reported variations in the Wannier projection window or alternative orbital bases, it is unclear whether the reported spectral agreement is robust or specific to one construction.

minor comments (1)

Notation for the orbital labels and configuration weights should be made consistent between the abstract, main text, and any tables/figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below.

read point-by-point responses

Referee: [electronic-structure and many-body methods sections] The central attribution of redox behavior to d_{3z^2-r^2} orbital shift and resulting charge transfer (abstract and results on valence configurations) rests on the alignment of Fe d states relative to O p bands. No tests of sensitivity to DFT functional (LDA/PBE/hybrid), Hubbard U value, or Wannier downfolding procedure are reported; these choices are known to move the d-manifold position by 1-2 eV and can therefore alter the dominant 3d^5 vs. 3d^6 weights or the vacancy preference.

Authors: We agree that the alignment of the Fe d states relative to the O p bands is central to the reported charge-transfer mechanism and that the results could in principle be sensitive to the DFT functional, U value, and Wannier construction. Our calculations employed the PBE+U approach with U = 4 eV on Fe, a standard parametrization for transition-metal perovskites. To address the concern, the revised manuscript will include additional calculations with U = 3 eV and U = 5 eV as well as a limited set of hybrid-functional results. These will be presented in a new subsection demonstrating that the mixed 3d^5/3d^6 character and the oxygen-vacancy preference for Fe^{2+} remain qualitatively unchanged. We will also report Wannier projections performed with two different energy windows to confirm the robustness of the d_{3z^2-r^2} orbital shift. revision: yes
Referee: [spectroscopy evaluation section] The crystal-field multiplet XPS calculation (section on spectroscopy) inherits the local Hamiltonian from the preceding DFT+Wannier step. Without reported variations in the Wannier projection window or alternative orbital bases, it is unclear whether the reported spectral agreement is robust or specific to one construction.

Authors: The XPS spectra are obtained by solving the local Hamiltonian extracted from the Wannier step, so variations in the projection details could affect the multiplet parameters. The revised manuscript will add a supplementary note that recomputes the XPS spectra using an alternative energy window for the Wannier projection and an extended orbital basis that includes neighboring Ti d states. The principal spectral features, including the satellites associated with the mixed valence, remain in agreement with experiment under these variations, indicating that the comparison is not an artifact of a single construction. revision: yes

Circularity Check

0 steps flagged

No significant circularity; forward calculations compared to external experiment

full rationale

The derivation proceeds from DFT-based many-body electronic structure to Wannier-derived local Hamiltonian, crystal-field multiplet XPS spectra, and attribution of mixed Fe valence and oxygen-vacancy tuning. These steps are presented as predictions from the chosen methodology and are compared to independent recent experimental measurements. No quoted equations or self-citations reduce any load-bearing claim to a fit or definition of its own inputs. The paper is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of standard DFT+many-body approximations for correlated d-states in oxides; no explicit free parameters, new entities, or ad-hoc axioms are stated in the abstract.

axioms (1)

domain assumption DFT-based many-body techniques provide a sufficiently accurate description of the electronic structure and charge transfer for Fe impurities in BaTiO3
Invoked to interpret the mixed valence and orbital-shift mechanism.

pith-pipeline@v0.9.1-grok · 5677 in / 1259 out tokens · 38855 ms · 2026-06-30T12:46:38.267328+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

5 extracted references

[1]

Local-density functional and on-site correlations: The electronic structure of la 2 cuo 4 and lacuo 3.Physical Review B, 49(20):14211, 1994

MT Czyżyk and GA Sawatzky. Local-density functional and on-site correlations: The electronic structure of la 2 cuo 4 and lacuo 3.Physical Review B, 49(20):14211, 1994

1994
[2]

Kristjan Haule, Chuck-Hou Yee, and Kyoo Kim. Dynamical mean-field theory within the full-potential methods: Electronic structure of ceirin 5, cecoin 5, and cerhin 5.Physical Review B—Condensed Matter and Materials Physics, 81(19):195107, 2010

2010
[3]

Covalency in transition-metal oxides within all-electron dynamical mean-field theory.Physical Review B, 90(7):075136, 2014

Kristjan Haule, Turan Birol, and Gabriel Kotliar. Covalency in transition-metal oxides within all-electron dynamical mean-field theory.Physical Review B, 90(7):075136, 2014

2014
[4]

Exact double counting in combining the dynamical mean field theory and the density functional theory.Physical review letters, 115(19):196403, 2015

Kristjan Haule. Exact double counting in combining the dynamical mean field theory and the density functional theory.Physical review letters, 115(19):196403, 2015

2015
[5]

Successive redox- mediated visible-light ferrophotovoltaics.Nature Communications, 11(1):966, 2020

Yuji Noguchi, Yuki Taniguchi, Ryotaro Inoue, and Masaru Miyayama. Successive redox- mediated visible-light ferrophotovoltaics.Nature Communications, 11(1):966, 2020. 7

2020

[1] [1]

Local-density functional and on-site correlations: The electronic structure of la 2 cuo 4 and lacuo 3.Physical Review B, 49(20):14211, 1994

MT Czyżyk and GA Sawatzky. Local-density functional and on-site correlations: The electronic structure of la 2 cuo 4 and lacuo 3.Physical Review B, 49(20):14211, 1994

1994

[2] [2]

Kristjan Haule, Chuck-Hou Yee, and Kyoo Kim. Dynamical mean-field theory within the full-potential methods: Electronic structure of ceirin 5, cecoin 5, and cerhin 5.Physical Review B—Condensed Matter and Materials Physics, 81(19):195107, 2010

2010

[3] [3]

Covalency in transition-metal oxides within all-electron dynamical mean-field theory.Physical Review B, 90(7):075136, 2014

Kristjan Haule, Turan Birol, and Gabriel Kotliar. Covalency in transition-metal oxides within all-electron dynamical mean-field theory.Physical Review B, 90(7):075136, 2014

2014

[4] [4]

Exact double counting in combining the dynamical mean field theory and the density functional theory.Physical review letters, 115(19):196403, 2015

Kristjan Haule. Exact double counting in combining the dynamical mean field theory and the density functional theory.Physical review letters, 115(19):196403, 2015

2015

[5] [5]

Successive redox- mediated visible-light ferrophotovoltaics.Nature Communications, 11(1):966, 2020

Yuji Noguchi, Yuki Taniguchi, Ryotaro Inoue, and Masaru Miyayama. Successive redox- mediated visible-light ferrophotovoltaics.Nature Communications, 11(1):966, 2020. 7

2020