{"id":"2583a619-f578-43a5-ba43-d2fc203c7413","arxiv_id":"2412.04869","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Molecular dynamics shows that 4 nm antiskyrmion nanodomains in rhombohedral BaTiO3 are stable up to about 85 K, with larger domains fragmenting into -1/6 pre-quarks.","lead":"This paper simulates antiskyrmion nanodomains in barium titanate over a range of sizes and temperatures. It finds that 4 nm domains stay stable to about 85 K, while larger domains split into fractional charge fragments called pre-quarks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The -1/6 pre-quark fragmentation in §III.B rests on an undefined local topological-charge-density procedure with no discretization or convergence test, so the paper's most novel claim may be a lattice artifact.","rationale":"The reader correctly identified model sensitivity of T* as a weakness, and the authors honestly acknowledge it in Section IV. My independent concern is different and arguably closer to the paper's most novel claim: the pre-quark fragmentation is presented as a definitive finding, but the paper provides no definition of local topological charge density and no convergence test for the charge-resolved images. This is not an ad hominem or a disagreement with consensus; it is an internal missing-validation issue. The global charge -2 is robust because it follows from the Berg-Lüscher construction, but the fractional charges -1/3 and -1/6 are local labels whose physical meaning depends on a well-defined continuum limit or on an explicit lattice procedure with demonstrated insensitivity. The paper does give creditworthy evidence elsewhere: the effective Hamiltonian from Ref. [28] is independently parameterized, the agreement with the shell-model result for the 2.8 nm smallest stable antiskyrmion is a useful cross-check, and the near-threshold 1 ns simulations strengthen the T* determination. Those positives do not cure the missing test for the fragmentation claim, so the conditional verdict remains appropriate, with the condition extended to require a discretization/convergence analysis of the local topological charge.","tokens_in":10954,"tokens_out":4094,"duration_ms":48548,"concrete_test":"Recompute the local topological charge density for the 5.8 nm, 10.6 nm, and 13 nm cases using two independent estimators: (i) the Berg-Lüscher solid angle with the exact lattice triangles, and (ii) a finite-difference continuum estimator Q_density = (1/4π) m · (∂x m × ∂y m) applied to the same time-averaged dipole fields. Then vary the local integration radius around each candidate defect from 0.5 nm to 2.0 nm and count the number of spots and their integrated charges. If the 12 spots each integrate to -1/6 only for one estimator or for a narrow radius window, the fragmentation claim is a numerical artifact; if both estimators converge to the same 12-spot structure over the full radius range, the claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new observation of the paper, highlighted in the abstract and Section IV, is that each -1/3 quark splits into two -1/6 pre-quarks for nanodomains larger than about 4.5 nm. However, the local topological charge density that supports this claim is never defined. Only the global Berg-Lüscher formula (Eqs. 2 and 3) is given, and the local decomposition into six -1/3 or twelve -1/6 spots requires an arbitrary choice of integration region, lattice-triangle orientation, color/bin threshold, and plane-projection convention. On a discrete lattice, only the total topological charge is invariant; local fractional charges depend on the chosen discretization and on how the local solid angles are summed. No test is reported for convergence with respect to the triangulation, the 2 fs timestep, the 64x64x64 supercell, or the local charge integration radius. Figures 2d, 2f, and 2h are presented as visual evidence, but without a quantitative criterion for a 'pre-quark' this could reflect the stretching of a single extended -1/3 defect rather than a genuine split into two independent -1/6 objects. Because the pre-quark fragmentation is one of the two headline claims, this missing validation is load-bearing.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11231,"tokens_out":6096,"duration_ms":60377,"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":[{"comment":"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":"Section III.B and Figures 2b, 2d, 2f, 2h"},{"comment":"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":"Section III.D, Figure 4"},{"comment":"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.","section":"Section III.C"}],"minor_comments":[{"comment":"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":"Abstract and Section III.A"},{"comment":"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":"Section III.C and Section III.D"},{"comment":"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":"Section III.B"},{"comment":"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.","section":"Section II"},{"comment":"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.'","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a computational condensed-matter journal and uses a defensible simulation protocol: a published effective Hamiltonian, the open-source feram code, explicit checks of total simulation times, and reproduction of the 2.8 nm smallest stable antiskyrmion from a different model. The main obstacle to acceptance is the undefined local topological charge density behind the pre-quark claim, which is a central and novel result. This is fixable: the authors can define the local charge decomposition, state the integration and threshold criteria, and provide convergence tests. The time-averaging issue in Section III.D is also fixable and should be addressed. I do not see a fundamental circularity in using the authors' own effective Hamiltonian, as the comparison with Ref. [14] provides an independent check; however, the manuscript should more clearly acknowledge that the quantitative T* value depends on the chosen model and on the finite heating rate and grid spacing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new part is the finite-temperature map: a 4 nm antiskyrmion holds its shape and charge -2 up to roughly 85 K, and larger domains shrink to that diameter before collapsing. The authors also reproduce the 2.8 nm smallest stable domain from Gonçalves et al. using a completely different method (effective Hamiltonian vs. shell potential), which is a nice cross-check. The paper is clearly written, the heating protocol is sensible (800 ps thermalization, 200 ps averaging, direct near-threshold runs), and they acknowledge that T* is model-dependent. That is all worth having.\n\nThe soft spot is the pre-quark claim. The local topological charge density that shows six -1/3 quarks splitting into twelve -1/6 pre-quarks is never defined. The paper gives only the global Berg-Lüscher formula; the local decomposition depends on an arbitrary integration region, lattice triangle orientation, and color threshold. On a discrete lattice, only the total charge is invariant, and nothing is shown about convergence with respect to supercell size, timestep, or the local charge integration procedure. The figures look plausible, but without a quantitative definition and a convergence test, the split could be a lattice artifact or simply a stretching of a single extended defect. Since this is one of the two headline results, it is load-bearing.\n\nTwo smaller issues: the heating runs are single trajectories, so threshold temperatures have no error bars, and no code or input files are shipped. The empirical fits for Tm and Tthr are labeled empirical, so they are not the problem.\n\nOverall, this is a serious computational study that extends known ground-state results to finite temperature. It deserves a serious referee, but the referee should insist on a precise definition of the local topological charge density and a brute-force convergence check before the pre-quark language is accepted. I would cite the stability map and the 2.8 nm cross-method agreement even if the pre-quark result later proves fragile.","headline":"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.","tokens_in":11784,"tokens_out":1735,"would_cite":true,"duration_ms":20050,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A 4 nm antiskyrmion in BaTiO3 keeps its topological charge -2 up to about 85 K.","keywords":["antiskyrmion","barium titanate","topological charge","molecular dynamics","effective Hamiltonian","ferroelectric nanodomain","fractional charge","thermal stability"],"falsifier":"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.","tokens_in":10750,"feed_emoji":"🌀","tokens_out":5415,"duration_ms":52597,"temperature":0.7,"pith_summary":"This paper uses molecular dynamics simulations on an effective Hamiltonian to ask how thermally stable antiskyrmions are in rhombohedral barium titanate when no external field is applied. It reports that a 4 nm antiskyrmion with topological charge -2 is the most stable configuration, keeping its size, shape, and charge up to a characteristic temperature T* ≈ 85 K. Smaller domains, from 2.8 to about 4.5 nm, stay nearly unchanged until they collapse, while larger domains shrink as temperature rises and collapse once their inner diameter reaches about 4 nm. The paper also finds that the negative topological charge of large domains is carried by six -1/3 quarks below about 4.5 nm, and by twelve -1/6 pre-quarks above it. This matters because it sets a quantitative temperature benchmark for using barium titanate as a platform for topological nanodomain devices.","feed_headline":"Smallest antiskyrmions keep charge -2 up to 85 K","feed_subtitle":"MD simulations map when six -1/3 quarks split into twelve -1/6 pre-quarks in rhombohedral BaTiO3.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the antiskyrmion ground-state configuration in rhombohedral BaTiO3 with six -1/3 quarks that the 1 K simulations reproduce and benchmark against.","marker":"[14]"},{"why":"Supplies the extended effective-Hamiltonian parameterization with anharmonic couplings used in all molecular dynamics runs.","marker":"[28]"},{"why":"Establishes the first-principles effective-Hamiltonian approach for BaTiO3 on which the simulations are built.","marker":"[22]"},{"why":"Describes the fast molecular-dynamics simulation method used to integrate the effective Hamiltonian.","marker":"[24]"},{"why":"Gives the lattice definition of the topological charge used to compute Q from local dipole vectors.","marker":"[40]"},{"why":"Provides the spherical-triangle implementation of the topological charge formula applied in the simulations.","marker":"[41]"}],"fun_headline_variants":["4-nm antiskyrmion keeps charge -2 to 85 K","Six quarks split into twelve pre-quarks in BaTiO3","Thermal limit for antiskyrmions: 85 K","BaTiO3 antiskyrmions hold charge -2 to 85 K","Antiskyrmions stable to 85 K without bias"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["4-nm antiskyrmion keeps charge -2 to 85 K","Six quarks split into twelve pre-quarks in BaTiO3","Thermal limit for antiskyrmions: 85 K","BaTiO3 antiskyrmions hold charge -2 to 85 K","Antiskyrmions stable to 85 K without bias"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001435,"raw_usage":{"total_tokens":5837,"prompt_tokens":1046,"completion_tokens":4791,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":4698}},"tokens_in":662,"tokens_out":4791,"duration_ms":30791,"temperature":1.0,"reasoning_tokens":4698,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:10:23.192426+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Nishimatsu, U","cited_arxiv_id":null,"evidence_quote":"Describes the fast molecular-dynamics simulation method used to integrate the effective Hamiltonian."}],"review_version":1}