Pulse-Driven Reconfiguration of Fractional Polar Topology in Zr-Substituted Barium Titanate
Pith reviewed 2026-05-21 03:09 UTC · model grok-4.3
The pith
Picosecond electric pulses on six vortex cores reconfigure fractional topological charges into 64 distinct metastable states inside a ferroelectric nanodomain.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Chemical doubling along the polar axis in 12.5 percent Zr-substituted BaTiO3 produces a coupled nanodomain consisting of alternating Q = -2 antiskyrmionic and Q = +4 skyrmionic slices whose local topological charge fragments into six -1/3 and six +2/3 contributions called topological quarks. Picosecond pulses applied to selected vortex-core columns, following a binary mask protocol on all six cores, drive the system to 64 distinct relaxed metastable configurations that are distinguishable by sector-resolved topological fingerprints and by their real-space polarization fields; the states remain stable over at least 1 ns of field-free evolution in the cryogenic simulation protocol.
What carries the argument
The binary pulse-mask protocol applied to the six vortex cores, which locally excites and relaxes the fractional topological texture inside the 2.6 nm nanodomain.
If this is right
- All 64 binary masks applied to the six vortex cores produce unique relaxed metastable states.
- The resulting states remain distinguishable through both sector-resolved topological charge distributions and real-space polarization vectors.
- Each reconfigured state persists without external fields for at least one nanosecond on the simulation timescale.
- Fractional polar topology thereby functions as a multistate configurational degree of freedom in the idealized low-temperature setting.
Where Pith is reading between the lines
- If the same pulse protocol works in real samples at higher temperature, the number of accessible states per nanodomain could enable dense multistate memory elements.
- The six-core addressing scheme might be extended to larger arrays of vortex cores to increase the number of reachable configurations.
- Similar reconfiguration could be tested in other ferroelectric compositions that support vortex or skyrmion textures.
- Direct imaging of polarization after targeted pulses would provide an experimental test of whether the fractional charges remain localized during switching.
Load-bearing premise
The effective-Hamiltonian molecular-dynamics model together with the specific cryogenic temperature protocol correctly captures the real-material stability and dynamics of the fractional topological textures over nanosecond timescales.
What would settle it
Apply the same binary pulse sequences to the six vortex cores in an actual Zr-substituted BaTiO3 sample at cryogenic temperature and determine whether the resulting polarization maps and topological charge distributions collapse to fewer than 64 distinct, long-lived states or match the full set predicted by the simulation.
read the original abstract
Polar topological textures in ferroelectrics can host internal structure beyond a single integer topological charge. Here, effective-Hamiltonian molecular-dynamics simulations are used to examine whether such internal fractional topology can be reconfigured by local electric excitation in ordered 12.5% Zr-substituted barium titanate. Chemical doubling along the polar axis stabilizes a coupled nanodomain texture consisting of alternating Q = -2 antiskyrmionic and Q = +4 skyrmionic slices, in which the local topological charge fragments into six -1/3 and six +2/3 localized contributions, denoted here as topological quarks, separated by Bloch-point-like singular conversion regions. Picosecond local electric-field pulses applied to selected vortex-core columns drive reconfiguration of the internal dipolar texture of a 2.6 nm nanodomain. Under a binary pulse-mask protocol addressing the six vortex cores, all 64 masks lead, within the chosen low-temperature simulation protocol, to distinct relaxed metastable configurations. The switching calculations are performed in a cryogenic regime, and the programmed states remain stable over at least 1 ns of field-free evolution on the simulation timescale. The resulting configurations are distinguishable both by sector-resolved topological fingerprints and by their real-space polarization fields. These results provide a computational proof of concept that fractional polar topology in a ferroelectric nanodomain can be locally reconfigured by ultrafast electric excitation and used as a multistate configurational degree of freedom in an idealized low-temperature setting.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses effective-Hamiltonian molecular-dynamics simulations to study 12.5% Zr-substituted BaTiO3, where chemical ordering stabilizes a nanodomain texture with alternating Q=-2 antiskyrmionic and Q=+4 skyrmionic slices. Local topological charge fragments into six -1/3 and six +2/3 contributions separated by Bloch-point-like singular regions. Picosecond electric-field pulses applied via binary masks to the six vortex cores are shown to produce 64 distinct relaxed metastable configurations, distinguishable by sector-resolved topological fingerprints and real-space polarization fields, that remain stable for at least 1 ns in field-free evolution under a cryogenic simulation protocol. The work is framed as a computational proof of concept for multistate configurational degrees of freedom based on fractional polar topology.
Significance. If the reported reconfiguration and distinguishability of the 64 states hold under more rigorous validation, the results would demonstrate a route to accessing a large number of metastable states in a ferroelectric nanodomain via local electric excitation, with potential implications for multistate logic or memory devices. The explicit mapping from 6 binary inputs to 64 outputs and the use of a concrete material system with Zr substitution are concrete strengths. The cryogenic regime and coarse-grained model, however, constrain immediate experimental translation.
major comments (2)
- Abstract, switching calculations paragraph: The central claim that all 64 binary pulse masks produce distinct relaxed metastable configurations rests on trajectories generated by the effective-Hamiltonian MD model. No quantitative metrics (e.g., distances between sector-resolved topological fingerprints), error bars, or convergence tests with respect to simulation length or ensemble averaging are reported, making it impossible to verify that the 64 states remain robustly distinguishable rather than partially overlapping or protocol-dependent.
- Abstract, paragraph on fractional topology: The fragmentation into six -1/3 and six +2/3 charges separated by Bloch-point-like singular conversion regions is a load-bearing feature for the reconfiguration claim. The effective-Hamiltonian lattice model with its Landau-type expansion and short-range interactions is not guaranteed to reproduce the correct energetics or stability of these atomic-scale singularities; lattice smoothing or artifacts could collapse distinctions among configurations. An atomistic benchmark or explicit test of singularity resolution is required.
minor comments (2)
- The term 'topological quarks' is introduced for the fractional charges without a precise definition or comparison to prior literature on fractional topology in ferroelectrics; this should be clarified in the main text.
- No information is provided on code or trajectory availability, which would aid reproducibility of the 64-mask protocol and stability runs.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable suggestions regarding our manuscript on pulse-driven reconfiguration of fractional polar topology in Zr-substituted BaTiO3. We address each major comment below, indicating where revisions will be made and providing substantive responses on the points raised.
read point-by-point responses
-
Referee: Abstract, switching calculations paragraph: The central claim that all 64 binary pulse masks produce distinct relaxed metastable configurations rests on trajectories generated by the effective-Hamiltonian MD model. No quantitative metrics (e.g., distances between sector-resolved topological fingerprints), error bars, or convergence tests with respect to simulation length or ensemble averaging are reported, making it impossible to verify that the 64 states remain robustly distinguishable rather than partially overlapping or protocol-dependent.
Authors: We agree that the absence of quantitative metrics limits the strength of the distinguishability claim. In the revised manuscript we will add pairwise distance calculations between the sector-resolved topological fingerprints of all 64 configurations, together with standard deviations obtained from an ensemble of independent trajectories. We will also report convergence tests with respect to simulation length to confirm that the observed distinctions are robust and independent of the specific protocol details. revision: yes
-
Referee: Abstract, paragraph on fractional topology: The fragmentation into six -1/3 and six +2/3 charges separated by Bloch-point-like singular conversion regions is a load-bearing feature for the reconfiguration claim. The effective-Hamiltonian lattice model with its Landau-type expansion and short-range interactions is not guaranteed to reproduce the correct energetics or stability of these atomic-scale singularities; lattice smoothing or artifacts could collapse distinctions among configurations. An atomistic benchmark or explicit test of singularity resolution is required.
Authors: The effective-Hamiltonian model has been validated in the literature for capturing the formation and stability of polar skyrmions and antiskyrmions in perovskites. The fractional fragmentation arises directly from the symmetry imposed by the Zr chemical ordering within the nanodomain. We will expand the discussion section to explicitly address the coarse-grained nature of the model and its limitations for atomic-scale singularities. However, a full atomistic benchmark lies outside the present scope and computational resources. revision: partial
- Requirement for an atomistic benchmark or explicit test of singularity resolution, which cannot be performed within the current effective-Hamiltonian framework and available resources.
Circularity Check
No circularity: results obtained directly from MD simulation trajectories
full rationale
The paper reports outcomes of effective-Hamiltonian molecular-dynamics runs in which binary pulse masks are applied to vortex cores and the system is allowed to relax under a fixed low-temperature protocol. The 64 distinct configurations are identified by inspecting the resulting polarization fields and sector-resolved topological charges after relaxation; these quantities are computed outputs of the trajectories rather than inputs that have been redefined or fitted to match themselves. No self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation chain appears in the derivation. The work is therefore a self-contained computational exploration whose central claim does not reduce to its own assumptions by construction.
Axiom & Free-Parameter Ledger
free parameters (2)
- Zr substitution concentration
- Simulation temperature
axioms (2)
- domain assumption Effective-Hamiltonian model parameters for BaTiO3 with Zr substitution accurately reproduce the dipolar energetics and topological textures.
- domain assumption The 1 ns field-free evolution timescale is sufficient to establish metastability of the reconfigured states.
invented entities (2)
-
topological quarks
no independent evidence
-
Bloch-point-like singular conversion regions
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
effective-Hamiltonian molecular-dynamics simulations... fractional topology... six -1/3 and six +2/3 localized contributions... Bloch-point-like singular conversion regions... binary pulse-mask protocol... all 64 masks lead... distinct relaxed metastable configurations
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Back et al., The 2020 Skyrmionics Roadmap, J
C. Back et al., The 2020 Skyrmionics Roadmap, J. Phys. D. Appl. Phys. 53, 363001 (2020)
work page 2020
-
[2]
A. Fert, N. Reyren, and V. Cros, Magnetic Skyrmions: Advances in Physics and Potential Applications, Nat. Rev. Mater. 2, 17031 (2017)
work page 2017
-
[3]
S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni, Skyrmion Lattice in a Chiral Magnet, Science 323, 915 (2009)
work page 2009
-
[4]
S. Seki and M. Mochizuki, Skyrmions in Magnetic Materials, 1st ed. (Springer International Publishing, Cham, Switzerland, 2015). 16
work page 2015
- [5]
-
[6]
J. Tang, Y. Wu, J. Jiang, L. Kong, W. Liu, S. Wang, M. Tian, and H. Du, Sewing Skyrmion and Antiskyrmion by Quadrupole of Bloch Points, Sci. Bull. 68, 2919 (2023)
work page 2023
-
[7]
J. Zang, M. Mostovoy, J. H. Han, and N. Nagaosa, Dynamics of Skyrmion Crystals in Metallic Thin Films, Phys. Rev. Lett. 107, 136804 (2011)
work page 2011
- [8]
-
[9]
A. K. Nayak, V. Kumar, T. Ma, P. Werner, E. Pippel, R. Sahoo, F. Damay, U. K. Rößler, C. Felser, and S. S. P. Parkin, Magnetic Antiskyrmions above Room Temperature in Tetragonal Heusler Materials, Nature 548, 561 (2017)
work page 2017
-
[10]
H. Aramberri and J. Íñiguez-González, Brownian Electric Bubble Quasiparticles, Phys. Rev. Lett. 132, 136801 (2024)
work page 2024
-
[11]
Junquera et al., Topological Phases in Polar Oxide Nanostructures, Rev
J. Junquera et al., Topological Phases in Polar Oxide Nanostructures, Rev. Mod. Phys. 95, 025001 (2023)
work page 2023
-
[12]
Das et al., Local Negative Permittivity and Topological Phase Transition in Polar Skyrmions, Nat
S. Das et al., Local Negative Permittivity and Topological Phase Transition in Polar Skyrmions, Nat. Mater. 20, 194 (2021)
work page 2021
-
[13]
Das et al., Observation of Room-Temperature Polar Skyrmions, Nature 568, 368 (2019)
S. Das et al., Observation of Room-Temperature Polar Skyrmions, Nature 568, 368 (2019)
work page 2019
-
[14]
A. K. Yadav et al., Observation of Polar Vortices in Oxide Superlattices, Nature 530, 198 (2016)
work page 2016
- [15]
-
[16]
L. Hu, Y. Wu, Y. Huang, H. Tian, and Z. Hong, Dynamic Motion of Polar Skyrmions in Oxide Heterostructures, Nano Lett. 23, 11353 (2023)
work page 2023
-
[17]
S. Prosandeev, S. Prokhorenko, Y. Nahas, J. Grollier, D. Talbayev, B. Dkhil, and L. Bellaiche, Ultrafast Activation and Tuning of Topological Textures in Ferroelectric Nanostructures, Adv. Electron. Mater. 8, 2200808 (2022)
work page 2022
-
[18]
Zhu et al., Dynamics of Polar Skyrmion Bubbles under Electric Fields, Phys
R. Zhu et al., Dynamics of Polar Skyrmion Bubbles under Electric Fields, Phys. Rev. Lett. 129, 107601 (2022)
work page 2022
-
[19]
Eremenko et al., Emergent Topological Polarization Textures in Relaxor Ferroelectrics, Nat
M. Eremenko et al., Emergent Topological Polarization Textures in Relaxor Ferroelectrics, Nat. Commun. 16, 7531 (2025)
work page 2025
-
[20]
P. Chen, Y. Nahas, S. Prokhorenko, and L. Bellaiche, Deterministic Generation of a Single- Byte Electric Skyrmion Bubble, Phys. Rev. Lett. 134, 256802 (2025)
work page 2025
-
[21]
M. A. P. Gonçalves, C. Escorihuela-Sayalero, P. Garca-Fernández, J. Junquera, and J. Íñiguez, Theoretical Guidelines to Create and Tune Electric Skyrmion Bubbles, Sci. Adv. 5, eaau7023 (2019)
work page 2019
-
[22]
Y. Qi, S. Liu, I. Grinberg, and A. M. Rappe, Atomistic Description for Temperature-Driven Phase Transitions in BaTiO3, Phys. Rev. B 94, 134308 (2016)
work page 2016
-
[23]
M. A. P. Gonçalves, M. Paściak, and J. Hlinka, Antiskyrmions in Ferroelectric Barium Titanate, Phys. Rev. Lett. 133, 066802 (2024)
work page 2024
-
[24]
F. Mayer and J. Hlinka, Thermal Stability and Topological Charge Fragmentation in Antiskyrmions of Rhombohedral Barium Titanate, Phys. Rev. B 111, 174106 (2025). 17
work page 2025
-
[25]
C. Halcrow and E. Babaev, Fractional Skyrme Lines in Ferroelectric Barium Titanate, Phys. Rev. Res. 6, L032011 (2024)
work page 2024
-
[26]
F. Gómez-Ortiz, L. Bastogne, S. Anand, M. Yu, X. He, and P. Ghosez, Switchable Skyrmion-- Antiskyrmion Tubes in Rhombohedral BaTiO3 and Related Materials, Phys. Rev. B 111, L180104 (2025)
work page 2025
-
[27]
V. Stepkova and J. Hlinka, Creation, Annihilation and Transport of Nonmagnetic Antiskyrmions within Ginzburg–Landau–Devonshire Model, J. Appl. Phys. 137, 034102 (2025)
work page 2025
-
[28]
Mayer, Topological Textures in Zr-Substituted Barium Titanate, Phys
F. Mayer, Topological Textures in Zr-Substituted Barium Titanate, Phys. Rev. B 113, 054111 (2026)
work page 2026
-
[29]
Y.-J. Wang, Y.-P. Feng, Y.-L. Tang, Y.-L. Zhu, Y. Cao, M.-J. Zou, W.-R. Geng, and X.-L. Ma, Polar Bloch Points in Strained Ferroelectric Films, Nat. Commun. 15, 3949 (2024)
work page 2024
-
[30]
D. Wolf, S. Schneider, U. K. Rößler, A. Kovács, M. Schmidt, R. E. Dunin-Borkowski, B. Büchner, B. Rellinghaus, and A. Lubk, Unveiling the Three-Dimensional Magnetic Texture of Skyrmion Tubes, Nat. Nanotechnol. 17, 250 (2022)
work page 2022
-
[31]
Y. Li, L. Pierobon, M. Charilaou, H.-B. Braun, N. R. Walet, J. F. Löffler, J. J. Miles, and C. Moutafis, Tunable Terahertz Oscillation Arising from Bloch-Point Dynamics in Chiral Magnets, Phys. Rev. Res. 2, 033006 (2020)
work page 2020
-
[32]
L. Bellaiche, A. García, and D. Vanderbilt, Low-Temperature Properties of Pb(Zr1-xTix)O3 Solid Solutions near the Morphotropic Phase Boundary, Ferroelectrics 266, 41 (2002)
work page 2002
-
[33]
L. Bellaiche, A. García, and D. Vanderbilt, Finite-Temperature Properties of Pb(Zr1-xTix)O3 Alloys from First Principles, Phys. Rev. Lett. 84, 5427 (2000)
work page 2000
-
[34]
U. V Waghmare and K. M. Rabe, Ab Initio Statistical Mechanics of the Ferroelectric Phase Transition in PbTiO3, Phys. Rev. B 55, 6161 (1997)
work page 1997
-
[35]
K. M. Rabe and U. V Waghmare, Localized Basis for Effective Lattice Hamiltonians: Lattice Wannier Functions, Phys. Rev. B 52, 13236 (1995)
work page 1995
- [36]
-
[37]
R. D. King-Smith and D. Vanderbilt, First-Principles Investigation of Ferroelectricity in Perovskite Compounds, Phys. Rev. B 49, 5828 (1994)
work page 1994
-
[38]
T. Nishimatsu, U. V. Waghmare, Y. Kawazoe, and D. Vanderbilt, Fast Molecular-Dynamics Simulation for Ferroelectric Thin-Film Capacitors Using a First-Principles Effective Hamiltonian, Phys. Rev. B - Condens. Matter Mater. Phys. 78, 104104 (2008)
work page 2008
-
[39]
T. Nishimatsu, M. Iwamoto, Y. Kawazoe, and U. V Waghmare, First-Principles Accurate Total Energy Surfaces for Polar Structural Distortions of BaTiO3, PbTiO3, and SrTiO3: Consequences for Structural Transition Temperatures, Phys. Rev. B 82, 134106 (2010)
work page 2010
- [40]
-
[41]
T. Nishimatsu, A. Grünebohm, U. V. Waghmare, and M. Kubo, Molecular Dynamics Simulations of Chemically Disordered Ferroelectric (Ba,Sr)TiO3 with a Semi-Empirical Effective Hamiltonian, J. Phys. Soc. Japan 85, 114714 (2016)
work page 2016
-
[42]
C. Mentzer, S. Lisenkov, Z. G. Fthenakis, and I. Ponomareva, Phase Evolution in the Ferroelectric Relaxor Ba(Ti1-x,Zrx)O3 from Atomistic Simulations, Phys. Rev. B 99, 064111 (2019). 18
work page 2019
-
[43]
A. Paul, J. Sun, J. P. Perdew, and U. V Waghmare, Accuracy of First-Principles Interatomic Interactions and Predictions of Ferroelectric Phase Transitions in Perovskite Oxides: Energy Functional and Effective Hamiltonian, Phys. Rev. B 95, 054111 (2017)
work page 2017
- [44]
- [45]
-
[46]
S. D. Bond, B. J. Leimkuhler, and B. B. Laird, The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics, J. Comput. Phys. 151, 114 (1999)
work page 1999
-
[47]
See Supplemental Material for details of the constrained energetic analysis, thermal-st ability calculations, and tests of alternative excitation protocols and pulse parameters
-
[48]
B. Berg and M. Lüscher, Definition and Statistical Distributions of a Topological Number in the Lattice O(3) σ-Model, Nucl. Phys. B 190, 412 (1981)
work page 1981
-
[49]
C. Heo, N. S. Kiselev, A. K. Nandy, S. Blügel, and T. Rasing, Switching of Chiral Magnetic Skyrmions by Picosecond Magnetic Field Pulses via Transient Topological States, Sci. Rep. 6, 27146 (2016)
work page 2016
-
[50]
S. Prosandeev, J. Grollier, D. Talbayev, B. Dkhil, and L. Bellaiche, Ultrafast Neuromorphic Dynamics Using Hidden Phases in the Prototype of Relaxor Ferroelectrics, Phys. Rev. Lett. 126, 027602 (2021). 1 Supplemental Material Pulse-Driven Reconfiguration of Fractional Polar Topology in Zr-Substituted Barium Titanate Florian Mayer1* 1Materials Center Leobe...
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.