REVIEW 3 major objections 5 minor 41 references
Thermal Stability and Topological Charge Fragmentation in Antiskyrmions of Rhombohedral Barium Titanate
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A 4 nm antiskyrmion in BaTiO3 keeps its topological charge -2 up to about 85 K.
desk verdict A solid MD temperature study with one headline claim—the -1/6 pre-quark splitting—that is under-defined and needs validation before it carries the paper. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the polarization-field antiskyrmion, a cylindrical nanodomain whose dipole texture winds around a hexagon-like cross-section with six vortices and a net topological charge of -2. The topological charge is computed with the Berg-Lüscher lattice formula, which sums signed spherical-triangle areas formed by normalized local dipole vectors; this is what lets the authors identify -1/3 quarks and -1/6 pre-quarks as local fractional-charge concentrations. The dynamics come from an effective Hamiltonian for BaTiO3 whose anharmonic couplings were fitted to density functional theory and to experimental transition temperatures, with local-mode amplitudes directly representing dipoles. The load-bearing step is the heating protocol: starting from each 1 K relaxed nanodomain, the system is heated in 5 K steps with no bias field, so the temperature at which the charge and size change is an intrinsic property of the model rather than an imposed boundary condition.
What would settle it
An independent measurement or atomistic calculation that follows a 4 nm reversed domain in zero field and finds its topological charge leaving -2 below about 85 K would falsify the central stability claim. A practical check is to compare the effective-Hamiltonian threshold with a different BaTiO3 model or with an experiment that images the polarization texture of such domains as a function of temperature.
Extended reading notes
Core claim
At 1 K, induced cylindrical nanodomains in the rhombohedral phase relax into antiskyrmions with total topological charge Q = -2 for all diameters from 2.8 to 13 nm, matching the symmetry and quark structure of earlier shell-model results. The paper's central quantitative finding is that the 4 nm antiskyrmion is optimal: its diameter and shape remain constant up to T* ≈ 85 K, and it is the critical size toward which all larger domains shrink before thermally collapsing into the single-domain state. In domains larger than about 4.5 nm, each of the six -1/3 topological quarks splits into two -1/6 pre-quarks, so the charge -2 is shared among twelve fractional defects located near the hexagon vertices. Above the threshold, larger domains show increasing quark mobility, shape fluctuations, and occasional excursions of the total topological charge to integer values such as -3, -1, 0, or +1 before collapse.
Load-bearing premise
The simulation's effective Hamiltonian accurately represents real rhombohedral barium titanate at nanometer length scales, so the computed 85 K stability and the quark-to-pre-quark split reflect the actual material and not the model's approximations.
Editorial extensions
If this is right
- At T* ≈ 85 K, a 4 nm antiskyrmion is a stable, zero-field topological object, giving a concrete target for experiments and device designs.
- Domains larger than 4.5 nm do not fail by abrupt unwinding; they first shed area and fragment their fractional charges, so their decay path is observable.
- The appearance of integer total charges -3, -1, 0, and +1 at high temperature means thermal energy can create neighboring topological states close in energy to the -2 ground state.
- The empirical threshold and depinning laws reported for the diameter evolution give testable predictions for how stability scales with domain diameter.
Reading between the lines
- If T* is set by domain-wall energy rather than by the specific Hamiltonian, similar stability should appear in other rhombohedral ferroelectrics with the same 3m symmetry; this is my inference, not a claim in the paper.
- The pre-quark splitting suggests a hierarchy of fractional-charge condensation, with pre-quarks pairing into quarks as the domain shrinks; this could be tested by simulating intermediate diameters between 4.5 and 5.8 nm.
- Applying an electric bias field, which the paper deliberately omits, may raise T* well above 85 K; this is a natural next simulation that follows from the paper's zero-field benchmark.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports molecular dynamics simulations, based on the published feram effective Hamiltonian for BaTiO3, of the thermal stability of antiskyrmion nanodomains in rhombohedral BaTiO3. Nineteen initial diameters from 2.2 nm to 13 nm are heated from 1 K to 211 K in 5 K steps, with 800 ps thermalization and 200 ps averaging per temperature. The simulations reproduce the smallest stable antiskyrmion diameter of 2.8 nm from Ref. [14] and find that all stable nanodomains from 2.8 nm to 13 nm carry a total topological charge of -2 at 1 K. The central claims are: (i) the most robust antiskyrmion has a diameter of 4 nm and is stable up to a characteristic temperature T* ~ 85 K; (ii) for diameters between 2.8 and about 4.5 nm the topological charge is fragmented into six -1/3 quarks, while for larger diameters each quark splits into two -1/6 pre-quarks; and (iii) above T* all larger nanodomains shrink to about 4 nm before collapsing into a single-domain state. The paper also presents empirical fits for the depinning and threshold temperatures as functions of diameter.
Significance. If the results are robust, the paper provides a valuable computational benchmark for ferroelectric antiskyrmion stability and introduces a size-dependent fragmentation mechanism of fractional topological charges that is not present in the earlier Ref. [14] description. The study is strengthened by the use of a published effective Hamiltonian with improved transition temperatures, the open-source feram MD code, stated convergence checks on simulation times, and the quantitative reproduction of the 2.8 nm smallest stable domain from an independent shell-model simulation. The authors also honestly acknowledge at the end of Section IV that T* is likely sensitive to model details. However, the paper's most novel claim, the pre-quark fragmentation, is supported by a local topological charge density that is never defined, and the time-averaging procedure used for topological charges raises additional concerns. These issues prevent the central claims from being fully validated as presented.
major comments (3)
- [Section III.B and Figures 2b, 2d, 2f, 2h] The local topological charge density that underlies the pre-quark claim is never defined. Equations (2) and (3) define only the global Berg-Lüscher topological charge, but the paper states in Section III.B that 'integrating locally over these newly formed pre-quarks confirms a fractional charge of -1/6.' This requires a precise local assignment of the solid-angle contributions to lattice sites or plaquettes, a choice of integration regions around each spot, a binning or coloring threshold for the density plots, and a convention for projecting onto the (111) planes. None of these is specified, and no convergence test is reported with respect to the lattice triangulation, the 2 fs timestep, the 64x64x64 supercell, or the local integration radius. Without such a definition and convergence test, the visual blue spots in Figures 2d, 2f, and 2h could equally reflect a single stretched -1/3 quark rather than two independent -1/6 pre-quarks. Because the pre-quark fragmentation is highlighted in the abstract and in Section IV as the most intriguing observation, this missing validation is load-bearing for the paper's central claim.
- [Section III.D, Figure 4] The topological charge is stated to be computed from time-averaged local dipole moments. The topological charge of a time-averaged configuration is not generally equal to the time average of the instantaneous topological charges, especially near the onset T_Q where fluctuations are exactly the quantities of interest. The reported integer jumps of Q to -3, -1, 0, and +1 could therefore be artifacts of time-averaging applied to a nonlinear function of the polarization field. The authors should clarify whether Q is computed from instantaneous snapshots and then averaged, or from a single time-averaged texture, and should test the sensitivity of the reported Q(T) trajectories to the length of the averaging window.
- [Section III.C] The characteristic temperature T* is identified with the threshold temperature of the 4 nm domain, but the definition of T* as 'the configuration's maximum stability against temperature' is only as precise as the discrete sampling of diameters and temperatures. The heating runs use 5 K steps and diameters spaced by 0.6 nm, so the statement T* ~ 85 K should be accompanied by an uncertainty estimate or a discussion of how the 5 K grid and the diameter spacing affect this benchmark. The near-threshold 1 ns checks are useful, but they do not resolve the discretization issue.
minor comments (5)
- [Abstract and Section III.A] The abstract states that antiskyrmions with topological charge -2 are the most stable nanodomain state 'across all examined diameters,' but Section III.A reports that diameters below 2.2 nm relax to a single-domain state; please qualify this statement by specifying the examined diameter range (2.8-13 nm).
- [Section III.C and Section III.D] The empirical formulas for T_m, T_thr, and T_Q introduce several fitted parameters (d_m ~ 14 nm, d_0 = 2.8 nm, prefactors 100 K, 80 K, 12 K, and alpha ~ 36 K), but no fitting procedure, statistical uncertainties, or validation against the data are provided. Since these formulas are not used in any later derivation, please present them clearly as descriptive fits and show the quality of the fits, for example by adding the curves to Figure 3.
- [Section III.B] The term 'pre-quark' is introduced before it is defined; please give an explicit operational definition of a pre-quark when it first appears, and distinguish it from a merely stretched quark region.
- [Section II] The velocity scaling thermostat is said to be 'validated against Nose-Poincaré results,' but no reference to this validation is given; please cite the validation or show a brief comparison.
- [References] Reference [29] is a 2002 preprint from the particle-physics literature; consider replacing it with a standard reference on fractional topological charges or justify why this particular preprint is the appropriate source for the term 'pre-quark.'
Circularity Check
No significant circularity: the antiskyrmion stability and fragmentation results are MD outputs from an independently parametrized effective Hamiltonian; no prediction reduces to a fit or self-citation by construction.
full rationale
The paper's central claims (T* ~ 85 K, fragmentation of six -1/3 quarks into twelve -1/6 pre-quarks, size-dependent collapse) are molecular dynamics outputs produced by the feram code from the effective Hamiltonian of Ref. [28]. That Hamiltonian is not fitted to antiskyrmion data; it is a first-principles effective Hamiltonian parametrized by DFT calculations and validated against experimental BaTiO3 transition temperatures, so importing it is independent support rather than circular. The topological charge is computed with the standard Berg-Lüscher formula (Eqs. 2-3), and the total charge is independently checked against the known -2 value. The empirical expressions for T_m, T_thr, and T_Q are explicit post-hoc fits to the simulation curves in Figs. 3-4, presented as descriptive summaries rather than predicted laws, so they do not make the simulation outputs circular. The paper does note in Section IV that T* is likely sensitive to model details, and the local topological charge density supporting the pre-quark claim is not defined with a discretization or convergence test; however, these are validation/correctness concerns, not equation-level circular reductions. No Eq. X is defined in terms of Eq. Y, and no fitted parameter is renamed as a prediction. Therefore the analysis is self-contained in the circularity sense.
Assumptions & free parameters
free parameters (4)
- alpha (T_Q fit coefficient) =
36 K
- dm (depinning temperature scale) =
14 nm
- d0 (threshold temperature offset) =
2.8 nm
- Empirical prefactors (100 K, 80 K, 12 K) =
not unique
assumptions (4)
- domain assumption The effective Hamiltonian of Ref 28 accurately represents the potential energy surface of BaTiO3, including anharmonic couplings.
- domain assumption The 64x64x64 supercell with periodic boundary conditions is large enough to avoid finite-size effects on nanodomains up to 13 nm.
- standard math The Berg-Lüscher formula gives the topological charge of the discrete dipole configuration.
- domain assumption A cylindrical domain with inverted <111> polarization, relaxed at 1 K, samples the relevant metastable states.
invented entities (1)
-
Pre-quarks
Cite this review
Pith. "Pith review of Thermal Stability and Topological Charge Fragmentation in Antiskyrmions of Rhombohedral Barium Titanate." pith.science (2026). https://pith.science/paper/NG7GTVWT
@misc{pith2026241204869,
author = {Pith},
title = {Pith review of: Thermal Stability and Topological Charge Fragmentation in Antiskyrmions of Rhombohedral Barium Titanate},
year = {2026},
howpublished = {\url{https://pith.science/paper/NG7GTVWT}},
note = {Machine review of arXiv:2412.04869}
}
read the original abstract
Antiskyrmions, as topological quasi-particles, hold significant promise for spintronics and nanoscale data storage applications. Using molecular dynamics simulations based on effective Hamiltonians, we investigate the thermal stability of antiskyrmion nanodomains in rhombohedral barium titanate. At 1 K, antiskyrmions with a topological charge of -2 emerge as the most stable nanodomain state across all examined diameters. In our systematic study, the most robust antiskyrmion was found to have a diameter of 4 nm, maintaining its original size, shape, and topological charge up to the characteristic temperature T* ~ 85 K. Domains with diameters between 2.8 and approximately 4.5 nm exhibit fragmentation into six topological defects, termed quarks, each carrying a fractional skyrmion charge of -1/3. For domains larger than 4.5 nm, each topological quark splits into two pre-quarks, each with a charge of -1/6. These larger nanodomains demonstrate increased mobility and a growing tendency for shape and skyrmion charge fluctuations. Above the T* temperature, all larger nanodomains gradually shrink to a diameter of about 4 nm before collapsing into a single-domain state. These findings reveal the relatively high stability of antiskyrmions over a broad temperature range, even in the absence of a stabilizing bias field, and emphasize the pivotal role of topological quark dynamics. This establishes barium titanate as a key platform for exploring and applying topological phenomena.
Reference graph
Works this paper leans on
-
[14]
M. A. P. Gonçalves, M. Paściak, and J. Hlinka, Antiskyrmions in Ferroelectric Barium Titanate, Phys. Rev. Lett. 133, 066802 (2024)
2024
-
[1]
A. Fert, N. Reyren, and V. Cros, Magnetic Skyrmions: Advances in Physics and Potential Applications, Nat. Rev. Mater. 2, 17031 (2017)
2017
-
[2]
Back et al., The 2020 Skyrmionics Roadmap, J
C. Back et al., The 2020 Skyrmionics Roadmap, J. Phys. D. Appl. Phys. 53, 363001 (2020)
2020
-
[3]
S. Seki and M. Mochizuki, Skyrmions in Magnetic Materials, 1st ed. (Springer International Publishing, Cham, Switzerland, 2015)
work page 2015
-
[4]
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). 12
work page 2009
-
[5]
Heinze, K
S. Heinze, K. von Bergmann, M. Menzel, J. Brede, A. Kubetzka, R. Wiesendanger, G. Bihlmayer, and S. Blügel, Spontaneous Atomic-Scale Magnetic Skyrmion Lattice in Two Dimensions, Nat. Phys. 7, 713 (2011)
2011
-
[6]
U. K. Rößler, A. N. Bogdanov, and C. Pfleiderer, Spontaneous Skyrmion Ground States in Magnetic Metals, Nature 442, 797 (2006)
work page 2006
-
[7]
Schulz, R
T. Schulz, R. Ritz, A. Bauer, M. Halder, M. Wagner, C. Franz, C. Pfleiderer, K. Everschor, M. Garst, and A. Rosch, Emergent Electrodynamics of Skyrmions in a Chiral Magnet, Nat. Phys. 8, 301 (2012)
2012
Show all 41 references
-
[8]
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)
2017
-
[9]
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)
2019
-
[10]
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)
2019
-
[11]
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)
2021
-
[12]
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)
2022
-
[13]
Shafer et al., Emergent Chirality in the Electric Polarization Texture of Titanate Superlattices, Proc
P. Shafer et al., Emergent Chirality in the Electric Polarization Texture of Titanate Superlattices, Proc. Natl. Acad. Sci. 115, 915 (2018)
2018
-
[15]
Stepkova and J
V. Stepkova and J. Hlinka, Creation, Annihilation and Transport of Nonmagnetic Antiskyrmions within Ginzburg-Landau-Devonshire Model, ArXiv:2410.05906 [Cond- Mat.Mtrl-Sci] (2024)
2024 arXiv
-
[16]
Halcrow and E
C. Halcrow and E. Babaev, Fractional Skyrme Lines in Ferroelectric Barium Titanate, Phys. Rev. Res. 6, L032011 (2024)
2024
-
[17]
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)
2023
-
[18]
Aramberri and J
H. Aramberri and J. Íñiguez-González, Brownian Electric Bubble Quasiparticles, Phys. Rev. Lett. 132, 136801 (2024)
2024
-
[19]
Nahas, S
Y. Nahas, S. Prokhorenko, L. Louis, Z. Gui, I. Kornev, and L. Bellaiche, Discovery of Stable Skyrmionic State in Ferroelectric Nanocomposites, Nat. Commun. 6, 8542 (2015)
2015
-
[20]
Gómez-Ortiz, L
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, ArXiv:2411.16395 [Cond-Mat.Mtrl-Sci] (2024)
2024 arXiv
-
[21]
Hashimoto and H
T. Hashimoto and H. Moriwake, Structure of Amorphous BaTiO3 by Molecular Dynamics Simulations Using a Shell Model, Phys. B Condens. Matter 579, 411799 (2020)
2020
-
[22]
Zhong, D
W. Zhong, D. Vanderbilt, and K. M. Rabe, First-Principles Theory of Ferroelectric Phase Transitions for Perovskites: The Case of BaTiO3, Phys. Rev. B 52, 6301 (1995)
1995
-
[23]
V Waghmare and K
U. V Waghmare and K. M. Rabe, Ab Initio Statistical Mechanics of the Ferroelectric Phase Transition in PbTiO3, Phys. Rev. B 55, 6161 (1997)
1997
-
[24]
Nishimatsu, U
T. Nishimatsu, U. V. Waghmare, Y. Kawazoe, and D. Vanderbilt, Fast Molecular-Dynamics 13 Simulation for Ferroelectric Thin-Film Capacitors Using a First-Principles Effective Hamiltonian, Phys. Rev. B 78, 104104 (2008)
2008
-
[25]
Bellaiche, A
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)
2000
-
[26]
Mayer, M
F. Mayer, M. Deluca, and M. N. Popov, Hidden Phases in Homovalent and Heterovalent Substituted BaTiO3, Phys. Rev. B 107, 184307 (2023)
2023
-
[27]
Mentzer, S
C. Mentzer, S. Lisenkov, Z. G. Fthenakis, and I. Ponomareva, Phase Evolution in the Ferroelectric Relaxor Ba(Ti1-xZrx)O3 from Atomistic Simulations, Phys. Rev. B 99, 064111 (2019)
2019
-
[28]
Mayer, M
F. Mayer, M. N. Popov, D. M. Evans, S. Krohns, M. Deluca, and J. Spitaler, Improved Description of the Potential Energy Surface in BaTiO3 by Anharmonic Phonon Coupling, Phys. Rev. B 106, 064108 (2022)
2022
-
[29]
M. E. de Souza, The General Structure of Matter, ArXiv:Hep-Ph/0207301 (2002)
2002
-
[30]
R. D. King-Smith and D. Vanderbilt, First-Principles Investigation of Ferroelectricity in Perovskite Compounds, Phys. Rev. B 49, 5828 (1994)
1994
-
[31]
K. M. Rabe and U. V Waghmare, Localized Basis for Effective Lattice Hamiltonians: Lattice Wannier Functions, Phys. Rev. B 52, 13236 (1995)
1995
-
[32]
J. H. Lee, U. V Waghmare, and J. Yu, First-Principles Effective Hamiltonian for Ferroelectric Polarization in BaTiO3/SrTiO3 Superlattices, J. Appl. Phys. 103, 124106 (2008)
2008
-
[33]
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)
2017
-
[34]
Mayer, M
F. Mayer, M. N. Popov, P. Ondrejkovic, J. Hlinka, J. Spitaler, and M. Deluca, Finite- Temperature Investigation of Homovalent and Heterovalent Substituted BaTiO3 from First Principles, Phys. Rev. B 106, 224109 (2022)
2022
-
[35]
Nishimatsu, M
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)
2010
-
[36]
J. Paul, T. Nishimatsu, Y. Kawazoe, and U. V Waghmare, Ferroelectric Phase Transitions in Ultrathin Films of BaTiO3, Phys. Rev. Lett. 99, 077601 (2007)
2007
-
[37]
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)
1999
-
[38]
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)
2016
-
[39]
M. C. Oliveira, R. A. P. Ribeiro, E. Longo, M. R. D. Bomio, F. V Motta, and S. R. de Lazaro, Temperature Dependence on Phase Evolution in the BaTiO3 Polytypes Studied Using Ab Initio Calculations, Int. J. Quantum Chem. 120, e26054 (2020)
2020
-
[40]
Berg and M
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)
1981
-
[41]
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)
2016
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.