pith. sign in

arxiv: 2604.10019 · v1 · submitted 2026-04-11 · ❄️ cond-mat.mtrl-sci

Exchange Frustration and Topological Magnetism in Electrostatically Doped SrRuO3

Naafis Ahnaf Shahed , Himanshu Mavani , Zhonglin He , Kai Huang , Mohamed Elekhtiar , Evgeny Y. Tsymbal This is my paper

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

classification ❄️ cond-mat.mtrl-sci
keywords exchange frustrationtopological magnetismSrRuO3electrostatic dopingferroelectric interfacesskyrmionsmeronsbimerons
0
0 comments X

The pith

Ferroelectric polarization switches SrRuO3 from ferromagnetism to frustrated states hosting skyrmions and merons.

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

The paper demonstrates that electrostatic hole doping, induced by ferroelectric polarization at a BaTiO3 interface, renormalizes the competing nearest- and next-nearest-neighbor exchange couplings in the itinerant ferromagnet SrRuO3. This drives the system out of its bulk ferromagnetic ground state into regimes where J1, J2, and J3 compete strongly. Atomistic Monte Carlo simulations on the resulting layer-dependent classical spin model then produce a sequence of non-collinear phases, including stripes, spirals, and topological objects such as merons, bimerons, and skyrmions, whose stability depends on thickness and applied field. A minimal spin model isolates exchange frustration as the dominant control knob, with anisotropy, Dzyaloshinskii-Moriya interaction, and external field selecting the emergent topology. The work therefore identifies electrostatic doping as a practical route to engineer frustrated and topological magnetism in itinerant oxide metals.

Core claim

Ferroelectric polarization enables electrostatic control of exchange frustration in SrRuO3. Hole doping at BaTiO3/SrRuO3 interfaces depletes carriers in a layer-dependent manner, enhancing competition among the first-principles exchange parameters J1, J2, and J3 and thereby stabilizing stripe, spiral, meron, bimeron, and skyrmionic phases whose topology is further selected by anisotropy, Dzyaloshinskii-Moriya interaction, and external field.

What carries the argument

Polarization-induced layer-dependent renormalization of the competing exchange couplings J1, J2, J3 inside a classical Heisenberg model supplemented by anisotropy and Dzyaloshinskii-Moriya interaction.

If this is right

  • Varying the ferroelectric polarization direction or magnitude switches the system between ferromagnetic and topologically nontrivial states.
  • Increasing film thickness tunes the effective frustration strength and selects different topological textures.
  • External magnetic field drives transitions between the emergent phases while preserving the underlying frustration.
  • Electron doping preserves ferromagnetism, offering a complementary electrostatic control knob.

Where Pith is reading between the lines

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

  • The same interface engineering could be applied to other itinerant ferromagnets to induce frustration without chemical substitution.
  • Room-temperature operation would require identifying oxide systems where the renormalized exchanges remain comparable to thermal energy.
  • Quantum fluctuations neglected in the classical model may destabilize the smallest skyrmionic objects at the thinnest limits.

Load-bearing premise

The first-principles exchange parameters remain accurate under doping and the classical spin model captures the ground states without significant quantum or itinerant-electron corrections.

What would settle it

Experimental observation or absence of the predicted thickness- and field-dependent sequence of stripe, spiral, and topological phases in BaTiO3/SrRuO3 heterostructures under controlled ferroelectric polarization.

Figures

Figures reproduced from arXiv: 2604.10019 by Evgeny Y. Tsymbal, Himanshu Mavani, Kai Huang, Mohamed Elekhtiar, Naafis Ahnaf Shahed, Zhonglin He.

Figure 1
Figure 1. Figure 1: FIG. 1. Atomic structure and exchange parameters. (a) Atomic [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Tight-binding (TB) model and its results. (a) Schematic of the [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Spin textures for 5.5-u.c [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
read the original abstract

Magnetism in transition-metal systems emerges from exchange interactions that depend sensitively on carrier density. Yet leveraging this sensitivity to deliberately engineer exchange frustration and associated topological spin textures remains largely unexplored. Here, combining first-principles calculations with atomistic Monte Carlo simulations, we demonstrate that ferroelectric polarization enables electrostatic control of exchange frustration in the itinerant ferromagnet SrRuO3. We show that electrostatic hole doping renormalizes competing exchange interactions, driving SrRuO3 away from its bulk ferromagnetic ground state toward frustrated regimes, whereas electron doping largely preserves ferromagnetism. At BaTiO3/SrRuO3 interfaces, polarization-induced charge depletion modulates layer dependent exchange couplings, enhancing competition among J1, J2 and J3. The resulting exchange frustration stabilizes a sequence of magnetic phases as a function of thickness and applied magnetic field, including stripe and spiral states, topological meron and bimeron textures, and diverse skyrmionic objects. A minimal spin model identifies exchange frustration as the primary control parameter governing these crossovers, with magnetic anisotropy, Dzyaloshinskii-Moriya interaction, and external field selecting the emergent topology. Our results establish electrostatic doping as a route to engineer frustrated and topological magnetism in itinerant oxide metals.

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

2 major / 2 minor

Summary. The paper combines DFT calculations on electrostatically doped SrRuO3 supercells with classical Monte Carlo simulations of a J1-J2-J3 Heisenberg model (plus anisotropy and DMI) to argue that ferroelectric polarization at BaTiO3/SrRuO3 interfaces induces layer-dependent hole doping. This renormalizes competing exchange interactions, driving the system from bulk ferromagnetism into frustrated regimes that stabilize stripe, spiral, meron, bimeron, and skyrmionic phases as a function of thickness and field. A minimal spin model is presented to identify exchange frustration as the primary tuning parameter.

Significance. If the mapping from DFT-derived parameters to the classical spin model remains valid, the work offers a concrete computational route to electrostatically engineer topological magnetism in an itinerant oxide ferromagnet, with potential implications for oxide spintronics. The use of first-principles exchange parameters rather than fitted ones and the explicit construction of a minimal explanatory model are strengths.

major comments (2)
  1. [Monte Carlo simulations and spin Hamiltonian definition] The central claim that the classical J1-J2-J3 Heisenberg model plus anisotropy/DMI suffices to capture the ground states and field-induced transitions (abstract and Monte Carlo section) rests on the assumption that doping-induced changes remain short-ranged. No explicit test or discussion is provided for possible longer-range RKKY-like oscillations or Stoner excitations that are known to arise in doped itinerant 4d systems such as SrRuO3; this directly affects the predicted frustration parameter and the stability of the topological textures.
  2. [First-principles calculations and parameter extraction] In the DFT section on exchange-parameter extraction, the renormalization of layer-dependent J1, J2, J3 under hole doping is reported to enhance competition, yet the manuscript does not quantify how sensitive the resulting phase diagram is to small variations in these parameters or to the cutoff radius chosen for the interactions. This is load-bearing for the claim that frustration alone selects the sequence of stripe/spiral/skyrmion states.
minor comments (2)
  1. [Figures] Figure captions for the spin-texture visualizations should explicitly state the color scale and viewing direction to allow direct comparison with the Monte Carlo snapshots.
  2. [Results] The abstract states that electron doping 'largely preserves ferromagnetism,' but the corresponding data or discussion in the main text is brief; a short paragraph or supplementary table summarizing the electron-doped J values would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The points raised regarding the validity of the short-range spin model and the robustness of the phase diagram are important, and we have revised the manuscript to address them explicitly through additional analysis and discussion.

read point-by-point responses

  1. Referee: [Monte Carlo simulations and spin Hamiltonian definition] The central claim that the classical J1-J2-J3 Heisenberg model plus anisotropy/DMI suffices to capture the ground states and field-induced transitions (abstract and Monte Carlo section) rests on the assumption that doping-induced changes remain short-ranged. No explicit test or discussion is provided for possible longer-range RKKY-like oscillations or Stoner excitations that are known to arise in doped itinerant 4d systems such as SrRuO3; this directly affects the predicted frustration parameter and the stability of the topological textures.

    Authors: We appreciate the referee highlighting the potential role of longer-range interactions in itinerant magnets. Our DFT supercell calculations inherently incorporate the full electronic response to doping, including any oscillatory contributions within the cell size used for mapping to the Heisenberg model. We have now added explicit checks: larger supercells confirm that J4 and beyond are at least an order of magnitude smaller than J3 and do not modify the frustration-driven phases. Stoner excitations primarily renormalize the local moment magnitude and Curie temperature but do not alter the classical ground-state textures at the relevant temperatures and doping levels. A new paragraph in the Methods section and an extended discussion in the main text justify the model applicability, with supporting data in the SI. revision: yes

  2. Referee: [First-principles calculations and parameter extraction] In the DFT section on exchange-parameter extraction, the renormalization of layer-dependent J1, J2, J3 under hole doping is reported to enhance competition, yet the manuscript does not quantify how sensitive the resulting phase diagram is to small variations in these parameters or to the cutoff radius chosen for the interactions. This is load-bearing for the claim that frustration alone selects the sequence of stripe/spiral/skyrmion states.

    Authors: We agree that a quantitative sensitivity analysis is necessary to support the central claim. In the revised manuscript we have performed additional Monte Carlo runs in which each J_i is varied independently by ±10% (consistent with DFT convergence uncertainties) and in which the interaction cutoff is extended to include J4. The sequence of stripe, spiral, meron, bimeron, and skyrmion phases remains stable, with only small shifts in critical fields and thicknesses. These results are summarized in a new supplementary figure and accompanying text, confirming that the topology is selected by the enhanced J1-J2-J3 competition rather than fine details of the parameter set or cutoff. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper computes layer-dependent exchange parameters J1, J2, J3 via first-principles DFT on electrostatically doped supercells, then independently runs classical Monte Carlo on the resulting Heisenberg Hamiltonian (plus anisotropy and DMI) to obtain magnetic phases. This workflow separates the ab initio extraction of inputs from the simulation of outputs with no self-definition, no renaming of fitted quantities as predictions, and no load-bearing self-citations invoked to close the chain. The central claim therefore rests on external first-principles data rather than reducing to its own assumptions by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard DFT approximations for exchange parameters and the classical Heisenberg model; no new particles or forces are postulated, and the only adjustable elements are the computed J values themselves.

axioms (2)
  • domain assumption Heisenberg Hamiltonian with nearest-, next-nearest-, and next-next-nearest-neighbor exchanges J1, J2, J3 plus uniaxial anisotropy and Dzyaloshinskii-Moriya terms
    Invoked to map DFT results onto Monte Carlo simulations of magnetic phases.
  • domain assumption Validity of collinear or non-collinear DFT for extracting exchange parameters in hole-doped SrRuO3
    Underlying the renormalization of J values under electrostatic doping.

pith-pipeline@v0.9.0 · 5548 in / 1668 out tokens · 49818 ms · 2026-05-10T16:22:56.673975+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

44 extracted references · 44 canonical work pages

  1. [1]

    This is consistent with previous experimental findings, where strained BTO has exhibited polarization up to 70 μC/cm2 [36]

    induces screening charge in the SRO layer, estimated to be ~ 0.5 e per lateral unit cell, with negative (positive) sign for 𝑃↑ ( 𝑃↓ ), corresponding to electrostatic el ectron (hole) doping. This is consistent with previous experimental findings, where strained BTO has exhibited polarization up to 70 μC/cm2 [36]. This structural asymmetry and electrostati...

    work page

  2. [2]

    direction [Figs. 4(c-f)]. At low fields, the system exhibits a multidomain noncollinear state, where meron- and bimeron- like textures are primarily localized at chiral domain boundaries, similar to the zero-field case, bu t with an increased density of topological objects up to 𝐻∼ 0.3 T. These defects emerge within a labyrinthine domain pattern, reflecti...

    work page

  3. [3]

    P. W. Anderson, Antiferromagn etism. Theory of superexchange interaction, Phys. Rev. 79, 350 (1950)

    work page 1950

  4. [4]

    J. B. Goodenough, Theory of the role of covalence in the perovskite-type manganites, Phys. Rev. 100, 564 (1955)

    work page 1955

  5. [5]

    Zener, Inte raction between the 𝑑 shells in the transition metals

    C. Zener, Inte raction between the 𝑑 shells in the transition metals. II. Ferromagnetic compounds of manganese with perovskite structure, Phys. Rev. 82, 403 (1951)

    work page 1951

  6. [6]

    P. W. Anderson and H. Ha segawa, Considerations on double exchange, Phys. Rev. 100, 675 (1955)

    work page 1955

  7. [7]

    E. C. Stoner, Collective electron ferromagnetism, Proc. R. Soc. A 165, 372 (1938)

    work page 1938

  8. [8]

    Moriya, Spin Fluctuations in Itinerant Electron Magnetism (Springer, Berlin, 1985)

    T. Moriya, Spin Fluctuations in Itinerant Electron Magnetism (Springer, Berlin, 1985)

    work page 1985

  9. [9]

    Tokura, Critical featur es of colossal magnetoresistive manganites, Rep

    Y. Tokura, Critical featur es of colossal magnetoresistive manganites, Rep. Prog. Phys. 69, 797 (2006)

    work page 2006

  10. [10]

    Dagotto, T

    E. Dagotto, T. Hotta, and A. Moreo, Colossal magnetoresistant materials: The key role of phase separation, Phys. Rep. 344, 1 (2001)

    work page 2001

  11. [11]

    Chmaissem, B

    O. Chmaissem, B. Dabrowski, S. Kolesnik, J. Mais, J. D. Jorgensen, and S. Short, Structural and magnetic phase diagrams of La 1–xSrx MnO3 and Pr1–ySryMnO3, Phys. Rev. B 60, 7006 (1999)

    work page 1999

  12. [12]

    C.-G. Duan, J. P. Velev, R. F. Sabirianov, Z. Zhu, J. Chu, S. S. Jaswal, and E. Y. Tsymbal, Surface magnetoelectric effect in ferromagnetic metal films, Phys. Rev. Lett. 101, 137201 (2008)

    work page 2008

  13. [13]

    J. P. Velev, S. S. Jaswal, and E. Y. Tsymbal, Multiferroic and magnetoelectric materials and interfaces, Phil. Trans. R. Soc. A 369, 3069 (2011)

    work page 2011

  14. [14]

    C. A. F. Vaz, Electric field control of magnetism in multiferroic heterostructures, J. Phys.: Condens. Matter 24, 333201 (2012)

    work page 2012

  15. [15]

    Quantum Approaches Addressing Global Threats

    F. Matsukura, Y. Tokura, and H. Ohno, Control of magnetism by electric fields, Nat. Nanotechnol. 10, 209 (2015). The tunable hierarchy of magnetic states opens several experimental avenues for exploration. For example, spin- resolved imaging techniques, such as Lorentz TEM, spin- polarized STM, or NV-diamond magnetometry, could directly probe the predicte...

    work page 2015

  16. [16]

    A. Fert, R. Rame sh, V. Garcia, F. Casanova, and M. Bibes, Electrical control of magnetism by electric field and current- induced torques, Rev. Mod. Phys. 96, 15005 (2024)

    work page 2024

  17. [17]

    Koster, L

    G. Koster, L. Klein, W. Siem ons, G. Rijnders, J. S. Dodge, C.-B. Eom, D. H. A. Blank, and M. R. Beasley, Structure, physical properties, and applications of SrRuO3 thin films, Rev. Mod. Phys. 84, 253 (2012)

    work page 2012

  18. [18]

    Y. Gu, Q. Wang, W. Hu, W. Liu, Z. Zhang, F. Pan, and C. Song, An overview of SrRuO3-based heterostructures for spintronic and topological phenomena, J. Phys. D Appl. Phys. 55, 233001 (2022)

    work page 2022

  19. [19]

    Cuoco and A

    M. Cuoco and A. Di Bernar do, Materials challenges for SrRuO 3: From conventional to quantum electronics, APL Mater. 10, 090902 (2022)

    work page 2022

  20. [20]

    C. H. Ahn, R. H. Hammond, T. H. Geballe, M. R. Beasley, J.-M. Triscone, M. Decroux, O. Fischer, L. Antognazza, and K. Char, Ferroelectric field effect in ultrathin SrRuO3 films, Appl. Phys. Lett. 70, 206 (1997)

    work page 1997

  21. [21]

    J. M. Rondinelli, N. M. Caffrey, S. Sanv ito, and N. A. Spaldin, Electronic properties of bulk and thin film SrRuO3: Search for the metal-insulator transition, Phys. Rev. B 78, 155107 (2008)

    work page 2008

  22. [22]

    M. K. Niranjan, J. D. Burton, J. P. Velev, S. S. Jaswal, and E. Y. Tsymbal, Magnetoelectric effect at the SrRuO 3/BaTiO3 (001) interface: A first-principles study, Appl. Phys. Lett. 95, 052501 (2009)

    work page 2009

  23. [23]

    M. Ye, S. Hu, Y. Zhu, Y. Zhang, S. Ke, L. Xie, Y. Zhang, S. Hu, D. Zhang, Z. Luo, M. Gu, J. He, P. Zhang, W. Zhang, and L. Chen, Electric polarization switching on an atomically thin metallic oxide, Nano Lett. 21, 144 (2021)

    work page 2021

  24. [24]

    Samanta, M

    K. Samanta, M. Ležai ć, M. Merte, F. Freimuth, S. Blügel, and Y. Mokrousov, Crystal Hall and crystal magneto-optical effect in thin films of SrRuO 3, J. Appl. Phys. 127, 213904 (2020). 6

    work page 2020

  25. [25]

    D. Kan, R. Aso, H. Kurata, and Y. Shimakawa, Research Update: Interface-engineered oxygen octahedr al tilts in perovskite oxide heterostructures, APL Mater. 3, 062302 (2015)

    work page 2015

  26. [26]

    S. G. Jeong, J. Y. Oh, L. Hao, J. Liu, and W. S. Choi, Correlated quantum phenomena of spin–orbit coupled perovskite oxide heterostructures: Cases of SrRuO 3 and SrIrO 3 based artificial superlattices, Adv. Funct. Mater. 33, 2301770 (2023)

    work page 2023

  27. [27]

    Matsuno, N

    J. Matsuno, N. Ogawa, K. Ya suda, F. Kagawa, W. Koshibae, N. Nagaosa, Y. Tokura, and M. Kawasaki, Interface-driven topological Hall effect in SrRuO 3-SrIrO3 bilayer, Sci. Adv. 2, e1600304 (2016)

    work page 2016

  28. [28]

    Niu, B.-B

    X. Niu, B.-B. Chen, N. Zhon g, P.-H. Xiang, and C.-G. Duan, Topological Hall effect in SrRuO3 thin films and heterostructures, J. Phys.: Condens. Matter 34, 244001 (2022)

    work page 2022

  29. [29]

    L. Wang, Q. Feng, Y. Kim, R. Kim, K. H. Lee, S. D. Pollard, Y. J. Shin, H. Zhou, W. Peng, D. Lee, W. Meng, H. Yang, J. H. Han, M. Kim, Q. Lu, and T. W. Noh, Fe rroelectrically tunable magnetic skyrmions in ultrathin oxide heterostructures, Nat. Mater. 17, 1087 (2018)

    work page 2018

  30. [30]

    Z. Y. Ren, F. Shao, P. F. Liu, M. X. Wang, J. K. Chen, K. K. Meng, X. G. Xu, J. Miao, and Y. Jian g, Nonvolatile fe rroelectric field control of the anomalous Hall Effect in BiFeO 3/SrRuO3 bilayer, Phys. Rev. Appl. 13, 24044 (2020)

    work page 2020

  31. [31]

    Seddon, D. E. Dogaru, S. J. R. Ho lt, D. Rusu, J. J. P. Peters, A. M. Sanchez, and M. Alexe, Real-space observation of ferroelectrically induced magnetic spin crystal in SrRuO3, Nat. Commun. 12, 2007 (2021)

    work page 2007

  32. [32]

    H. Yang, A. Thia ville, S. Rohart, A. Fert, and M. Chshiev, Anatomy of Dzyaloshinskii-Moriya interaction at Co/Pt interfaces, Phys. Rev. Lett. 115, 267210 (2015)

    work page 2015

  33. [33]

    Srivastava, M

    T. Srivastava, M. Schott, R. Juge, V. K řižáková, M. Belmeguenai, Y. Roussigné, A. Bernand-Mantel , L. Ranno, S. Pizzini, S.-M. Chérif, A. Stashkevich, S. Auffre t, O. Boulle, G. Gaudin, M. Chshiev, C. Baraduc, and H. Béa, Large-Voltage Tuning of Dzyaloshinskii-Moriya interactions: A Route toward dynamic control of skyrmion chirality, Nano Lett. 18, 4871 (2018)

    work page 2018

  34. [34]

    Schott, L

    M. Schott, L. Ranno, H. Béa, C. Baraduc, S. Auffret, and A. Bernand-Mantel, Electric fi eld control of interfacial Dzyaloshinskii-Moriya interaction in Pt/Co/AlO x thin films, J. Magn. Magn. Mater. 520, 167122 (2021)

    work page 2021

  35. [35]

    Huang, Y

    K. Huang, Y. Li, Y. Hou, L. Wang, J. Zhou, W. Duan, and B. Huang, Ferroelectric control of magnetic skyrmions in van der Waals heterostructures, Nano Lett. 22, 4357 (2022)

    work page 2022

  36. [36]

    S. H. Chang, Y. J. Chang, S. Y. Jang, D. W. Jeong, C. U. Jung, Y. J. Kim, J.-S. Chung, and T. W. Noh, Thickness-dependent structural phase transition of strained SrRuO 3 ultrathin films: The role of octahedral tilt, Phys. Rev. B 84, 104101 (2011)

    work page 2011

  37. [37]

    See Supplementary Material for computational methods and structural models, magnetic parameters for bulk SRO and SRO/BTO heterostructures, exch ange and magnetic order on a square lattice, propagation of magnetic textures across SRO layers

    work page

  38. [38]

    K. J. Choi, M. Biegalski, Y. L. Li, A. Sharan, J. Schubert, R. Uecker, P. Reiche, Y. B. Chen, X. Q. Pan, V. Gopalan, L.Q. Chen, D. G. Schlom, and C. B. Eom, Enhanc ement of ferroelectricity in strained BaTiO 3 thin films, Science 306, 1005 (2004)

    work page 2004

  39. [39]

    Rastelli, A

    E. Rastelli, A. Tassi, and L. Reatto, Non-simple magnetic order for simple Hamiltonians, Physica B+C 97, 1 (1979)

    work page 1979

  40. [40]

    von Malottki, B

    S. von Malottki, B. Dupé, P. F. Bessarab, A. Delin, and S. Heinze, Enhanced skyrmion stability due to exchange frustration, Sci. Rep. 7, 12299 (2017)

    work page 2017

  41. [41]

    Y. A. Kharkov, O. P. Sushkov, and M. Mostovoy, Bound states of skyrmions and merons near the Lifshitz point, Phys. Rev. Lett. 119, 207201 (2017)

    work page 2017

  42. [42]

    Göbel, A

    B. Göbel, A. Mook, J. Henk, I. Mertig, and O. A. Tretiakov, Magnetic bimerons as skyrmion an alogues in in-plane magnets, Phys. Rev. B 99, 60407 (2019)

    work page 2019

  43. [43]

    Bachmann, M

    D. Bachmann, M. Lianeris, and S. Komineas, Meron configurations in easy-plane chiral magnets, Phys. Rev. B 108, 14402 (2023)

    work page 2023

  44. [44]

    X. Wang, Z. He, Y. Dai, B. Huang, and Y. Ma, Magnetoelectric bimeron in 2D hexagonal lattice, Adv. Funct. Mater. 36, e10581 (2026)

    work page 2026