{"id":"be7dcc4c-1dd7-48f2-aee5-85514266e419","arxiv_id":"2501.09619","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"K2Co2(SeO3)3 shows multiple fractional magnetization plateaus and a finite-field BKT phase region in a phase diagram the authors reproduce with a weak-dimer triangular-lattice model.","lead":"Researchers measured the magnetic phases of the layered cobalt compound K2Co2(SeO3)3 and found magnetization plateaus at 1/3, 1/2, 2/3, and 5/6 of full saturation. Their experiments plus simulations point to a low-field phase diagram controlled by an emergent U(1) x S3 symmetry, including a Berezinskii-Kosterlitz-Thouless region at finite magnetic field.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"BKT labeling rests on an unvalidated classical-MC parameter set; experiment shows only broad thermodynamic anomalies, so the finite-field BKT claim is not directly evidenced.","rationale":"The reader's weakest assumption—that the classical Monte Carlo with hand-picked parameters faithfully represents the real quantum magnet—is exactly the hinge of the paper's central claim. I agree with that assessment and add that the lack of any direct experimental probe of BKT correlations (e.g., algebraic decay or universal jump) makes the label purely inferential. The paper itself admits the parameters are undetermined and that quantum effects are needed for the 5/6 plateau, weakening the classical-to-quantum mapping. These admissions are explicit in the Discussion and should be weighed in the verdict. Nevertheless, the experimental dataset (magnetization plateaus, specific heat, MCE) is solid and the U(1)⊗S3 symmetry framework is a plausible organizing principle; the concerns are addressable with neutron scattering and quantum many-body calculations. Therefore the conditional verdict is appropriate: the central claim is not yet established, but the paper is worth publishing with the caveat that the BKT identification is a theory-driven interpretation rather than a directly measured finding. No change to the reader's verdict is needed.","tokens_in":11143,"tokens_out":7489,"duration_ms":81732,"concrete_test":"Identify a single crystal of K2Co2(SeO3)3 and perform neutron diffraction / inelastic neutron scattering in the candidate BKT region (e.g., B≈4 T, 1 K<T<3 K) to measure the static spin structure factor S(q). A BKT phase would show a power-law broadening (no magnetic Bragg peak) with a diverging S(q_K) as T→Tc1 from below; the observation of magnetic Bragg peaks would prove long-range order and falsify the BKT label. In parallel, run quantum Monte Carlo simulations of the S=1/2 version of Eq. (1) with parameters fitted to the measured Curie-Weiss temperature and the 1/3- and 1/2-plateau fields; if no finite-field BKT phase survives or the fitted parameters are outside the range used in Fig. 4(a), the claimed identification loses its theoretical support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's headline claim—a finite-field BKT phase region in K2Co2(SeO3)3—is carried by the transfer of phase labels (BKT, tss, uud, Y) from the classical Monte Carlo phase diagram of Eq. (1) with the single parameter set J'_⊥/J=1, J⊥/J=0.2, Δ=3 (Fig. 4(a)) onto the experimental phase diagram (Fig. 3(d)). This transfer is load-bearing and insecure for three reasons. First, the authors state in the Discussion that the exact exchange parameters 'remain undetermined' and that the same classical model 'successfully reproduces the 1/3, 1/2, and 2/3 plateaus' but cannot produce the observed 5/6 plateau, which they attribute to quantum effects; thus the classical model is at best partially validated and is acknowledged to miss one of the four plateau phases. Second, the experimental evidence for phase boundaries consists only of broad specific-heat anomalies and magnetization features; no measurement directly probes the algebraic spin correlations, the order-parameter symmetry, or the universal BKT jump that would distinguish a BKT quasi-long-range-ordered phase from a conventional ordered phase or from a broad crossover. Third, the 'weak-dimer' classification itself is inferred from the close similarity of susceptibility and critical fields to the monolayer compound rather than from a measured dimer gap, and the chosen J'_⊥/J=1 is not obviously small; if the true intra-dimer coupling is stronger, the S3 symmetry framework and the BKT region derived from it would not apply. In short, the finite-field BKT phase is a plausible but unvalidated theoretical inference, not an experimentally established finding.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript combines high-field magnetization, specific heat, and magnetocaloric measurements on single-crystal K2Co2(SeO3)3 with classical Monte Carlo simulations of an easy-axis bilayer triangular Heisenberg model. The experiments show magnetization plateaus at 1/3, 1/2, 2/3, and 5/6 of saturation and a series of low-field thermodynamic anomalies from which the authors construct a B-T phase diagram. Guided by simulations, the low-field phases are labeled as a U(1)-broken 'Y' phase, an S3-broken tss phase, a Z3-broken uud 1/3-plateau phase, and an intermediate finite-field BKT phase. The authors argue that the low-field phase structure of the weak-dimer bilayer is governed by an emergent U(1)⊗S3 symmetry with S3 = Z3⊗Z2^d, and that the close agreement between experiment and simulation validates this symmetry description.","tokens_in":11436,"tokens_out":6100,"duration_ms":61675,"significance":"Should the finite-field BKT region and the S3 phase topology be confirmed, the paper would be a valuable experimental realization of nontrivial symmetry-enriched phase structure in a frustrated bilayer magnet, with clean plateau fractions and a promising platform for further study. The paper includes useful raw experimental data and a transparent symmetry-based framework. Its strengths include the observation of four fractional plateaus, the consistency between temperature- and field-dependent thermodynamic data, and the classical Monte Carlo analysis with Binder ratios and correlation-length crossings. However, the headline BKT claim is currently carried by a model-based labeling scheme whose parameters are not fixed by experiment and whose finite-field BKT phase is not directly probed; the significance therefore hinges on an inference that is not yet established.","major_comments":[{"comment":"The phase labels BKT, tss, uud, and Y in Fig. 3(d) are explicitly imported from the classical Monte Carlo phase diagram of Eq. (1) computed at one parameter set (J'_⊥/J=1, J⊥/J=0.2, Δ=3). The same simulation is then used in the Discussion to claim that the agreement 'validates the pivotal role' of S3 symmetry. This is circular for phase identification: the thermodynamic anomalies establish transition lines, but not the broken-symmetry character of the phases between them. Moreover, the parameters are stated to remain undetermined. For the central BKT claim to be load-bearing, the authors must either determine J'_⊥/J, J⊥/J, and Δ from experiment (for example by fitting the susceptibility, critical fields, or neutron data) or explicitly present the phase labels as a model-based scenario rather than an experimental identification. The current wording 'identify a BKT phase region' in the abstract exceeds the evidence.","section":"Theoretical phase diagram; Fig. 4(a); Fig. 3(d) caption"},{"comment":"The experimental case for a BKT phase rests on broad specific-heat anomalies and field-derivative features; no measurement directly establishes quasi-long-range order, a universal BKT jump, or algebraic correlation decay. Since the observed peaks are broad and the transitions are identified only by peak positions, the same data could be described by weak first-order transitions or crossovers. To support a BKT label, the authors need either a direct probe of in-plane correlations (e.g., neutron scattering), a quantitative finite-size/universality test of the transition, or a clear statement that 'BKT' is a theoretical interpretation inferred from the simulation rather than an experimental determination.","section":"Phase diagram below 14 T; Fig. 3(a)-3(d)"},{"comment":"The weak-dimer classification is supported mainly by the similarity of χ(T) and critical fields to the monolayer compound and the absence of an obvious spin gap at 2 K. These are suggestive but do not determine J'_⊥/J; the simulation uses J'_⊥/J=1, which is not a small dimer coupling. If the actual intra-dimer coupling is stronger, the S3=Z3⊗Z2^d symmetry reduction and the BKT phase derived from it would not apply. The manuscript should provide a quantitative estimate of J'_⊥ (e.g., from a fit of the full susceptibility or from the dimer gap scale) before using the weak-dimer regime as the basis for the phase diagram.","section":"Weak-dimer interaction and multiple magnetization plateaus; Eq. (1)"},{"comment":"The classical model used for the low-field labels is acknowledged not to reproduce the 5/6 plateau, which is attributed to quantum effects. While this concerns high fields, it weakens the 'remarkable agreement' argument used to validate the model and thereby the transferability of the low-field labels. Please specify the field and temperature range over which the parameter set is actually validated, and discuss how the missing plateau affects the confidence in the low-field phase identification.","section":"Discussions and Conclusions"}],"minor_comments":[{"comment":"In the sentence 'the exact experimental parameters of K2Co2(SeO3)2 remain undetermined', the formula appears to be a typo for K2Co2(SeO3)3; please correct it.","section":"Discussions and Conclusions"},{"comment":"The order parameters in Eq. (2) do not include a normalization by the number of dimer sites, and the quantity δq in the correlation-length expression is only loosely defined. Please make the definitions self-contained.","section":"Eq. (2)"},{"comment":"The phase-boundary lines are 'guides for the eye'; please indicate which symbols correspond to specific-heat, magnetization, and MCE data and add error bars or symbol-size information so the reader can judge the uncertainty in boundary positions.","section":"Fig. 3(d)"},{"comment":"The paper references fine-tuned simulations that reproduce the 1/3, 1/2, and 2/3 plateaus (Fig. S2) and the MCE data (Fig. S4), but the text does not give the corresponding parameter values or a link to the supplementary file. Please provide these details for reproducibility.","section":"Supplementary Material"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within scope for a condensed-matter journal, and the experimental data set is valuable. The editorial decision should hinge on whether the authors are willing to soften the BKT claim to a model-based prediction and to separate experimental facts from theoretical labels. I would not recommend acceptance with the current abstract; a major revision with clearly distinguished experimental results and theoretical interpretations would be appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one for the experimental phase diagram, not for the BKT claim. Fu et al. give us a clean high-field magnetization staircase for K2Co2(SeO3)3—1/3, 1/2, 2/3, 5/6 plateaus—plus a low-field thermodynamic phase diagram from specific heat, magnetization, and MCE. That is new, internally consistent, and worth having on the record. The comparison with the monolayer K2Co(SeO3)2 is also useful and supports the weak-dimer classification only in the broad sense that the two compounds behave similarly; it is not a direct measurement of the dimer coupling.\n\nThe theoretical side is where I part company with the abstract. The U(1)⊗S3 symmetry framework is a reasonable organizing principle for a weak-dimer bilayer triangular magnet, and it builds honestly on Chen's earlier work (ref. 26). The classical MC phase diagram of Eq. (1) is a nice illustration of what that symmetry can produce. But the central claim—a finite-field BKT region in this material—rests on transferring the phase labels BKT, tss, uud, Y from that simulation to the experimental phase diagram. The same simulation is then invoked as validation of the experiment. That is circular for phase identification, and the stress-test note has it right: the experiment shows only broad specific-heat anomalies and magnetization features, with no direct probe of algebraic decay or a BKT universal jump. The parameter set J'_⊥/J = 1, J⊥/J = 0.2, Δ = 3 is admitted to be undetermined, and the classical model does not produce the observed 5/6 plateau, which the authors attribute to quantum effects. So the BKT region should be labeled a plausible inference, not an established finding.\n\nThe paper also states in the Discussion that 'exact experimental parameters remain undetermined,' which is honest, and that acknowledgment is partly why I still take the work seriously. The authors are explicit about what the classical model can and cannot do. That does not rescue the headline claim, but it makes the paper a legitimate candidate for publication after revision.\n\nWho gets value: experimentalists working on frustrated triangular magnets, and theorists interested in bilayer S3 symmetry. The experimental dataset is citable; the BKT label should be softened. I would send this to a serious referee, with instructions to require a tempered abstract and either neutron scattering or a measured dimer gap before the BKT region is presented as fact. On balance, accept after major revision.","headline":"The new magnetization and heat-capacity data on K2Co2(SeO3)3 are a real experimental contribution; the finite-field BKT claim is a plausible but unproven interpretation resting on a circular classical-MC label transfer.","tokens_in":12079,"tokens_out":2637,"would_cite":true,"duration_ms":24768,"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":"K2Co2(SeO3)3 exhibits a finite-field Berezinskii-Kosterlitz-Thouless phase region whose transitions are governed by U(1) ⊗ S3 symmetry.","keywords":["Berezinskii-Kosterlitz-Thouless transition","triangular lattice antiferromagnet","magnetization plateaus","bilayer quantum magnet","spin-1/2 dimer","emergent symmetry","frustrated magnetism","classical Monte Carlo"],"falsifier":"A neutron-scattering measurement in the region labeled BKT, roughly $B \\simeq 1$–$3$ T and $T \\simeq 0.1$–$0.5$ K, that found Bragg peaks instead of algebraically decaying spin correlations with a temperature-dependent exponent would refute the BKT assignment.","tokens_in":10924,"feed_emoji":"🧲","tokens_out":7974,"duration_ms":74714,"temperature":0.7,"pith_summary":"The paper claims that the bilayer triangular antiferromagnet K2Co2(SeO3)3 is a weak-dimer, easy-axis magnet whose low-field phases are controlled by an emergent $U(1) \\otimes S_3$ symmetry, with $S_3 = \\mathbb{Z}_3 \\otimes \\mathbb{Z}_2^d$ combining threefold sublattice permutations and bilayer dimer exchange. On this basis, it assigns the low-field experimental phase diagram a finite-field Berezinskii-Kosterlitz-Thouless (BKT) region, together with tss, uud, and Y phases, and reports magnetization plateaus at 1/3, 1/2, 2/3, and 5/6 of saturation. A reader should care because a stable finite-field BKT phase in a real frustrated magnet has been elusive; if the identification holds, this compound becomes a concrete platform for studying two-dimensional topological phase transitions in a magnetic field.","feed_headline":"BKT phase region found at finite fields in a triangular magnet","feed_subtitle":"Magnetization plateaus and heat-capacity data tie the transitions to a U(1) ⊗ S3 symmetry.","key_machinery":"The model is a spin-1/2 (treated classically in simulation) easy-axis antiferromagnet on a bilayer triangular lattice, Eq. (1), with intralayer $J$, intra-dimer $J'_\\perp$, and crossed-layer $J_\\perp$ couplings and anisotropy $\\Delta$. The order parameters $\\psi^x_K$, $\\psi^z_K$, and $\\psi^z_\\Gamma$ defined in Eq. (2) measure respectively U(1) spin-rotation breaking, $\\mathbb{Z}_3$ sublattice-permutation breaking, and $\\mathbb{Z}_2^d$ bilayer-dimer breaking; Binder ratios and scaled correlation lengths from parallel-tempering Monte Carlo locate the transition temperatures. The emergent symmetry $U(1) \\otimes S_3$, continuous c-axis spin rotations combined with the discrete $\\mathbb{Z}_3$ and $\\mathbb{Z}_2^d$ factors, is the organizing principle that fixes the phase topology at low fields.","core_discovery":"Combining single-crystal magnetization, heat capacity, and magnetocaloric measurements with classical Monte Carlo simulations of a two-layer easy-axis triangular-lattice Hamiltonian, the paper establishes a low-field phase diagram in which the paramagnet gives way through an intermediate BKT phase to the 1/3-plateau uud phase, and at lower fields and temperatures to the S3-broken tss phase and a U(1)-broken Y phase. The central claim is that the sequence and universality of the transitions are dictated by the emergent $U(1) \\otimes S_3$ symmetry, where $S_3 = \\mathbb{Z}_3 \\otimes \\mathbb{Z}_2^d$ combines the threefold sublattice permutation of the triangular lattice with the swap symmetry of the weakly coupled Co2 dimers. The same symmetry argument explains why a BKT region survives at finite fields in this bilayer compound, whereas in the monolayer counterpart such a region is confined to zero field.","pith_inferences":["Beyond the paper, the same symmetry argument predicts that isostructural bilayer dimers such as K2Ni2(SeO3)3 should develop a finite-field BKT pocket, since the $S_3$ symmetry comes from bilayer geometry rather than from finely tuned couplings; heat-capacity scans at fields up to a few tesla would test this.","Beyond the paper, the 5/6 plateau that classical Monte Carlo cannot produce could be a signature of quantum order by disorder in the dimer-triangle lattice, and a tensor-network ground-state calculation at zero temperature could decide whether it is a fully quantum plateau.","Beyond the paper, a stable finite-field BKT region implies a universal jump in the spin stiffness at the BKT line, which is in principle measurable through low-frequency torque or specific-heat scaling; no such measurement is reported here."],"forward_implications":["A finite-field BKT region exists in K2Co2(SeO3)3, making it a concrete material where two-dimensional BKT physics survives in a magnetic field.","The low-field transitions are governed by the emergent $U(1) \\otimes S_3$ symmetry, so the phase sequence PM–BKT–uud, with tss and Y at lower temperature, is fixed by symmetry rather than by the precise exchange-coupling values.","The magnetization plateaus at 1/3, 1/2, and 2/3 of saturation are reproduced by the classical model with the chosen parameters, while the 5/6 plateau requires effects beyond that model.","The transition between the S3-symmetric and $\\mathbb{Z}_3$-broken regimes is not direct; it proceeds through an intermediate $\\mathbb{Z}_2^d$-symmetric uud phase or through the BKT phase.","The finite-field BKT region distinguishes this bilayer compound from its monolayer counterpart and establishes the material as a platform for studying symmetry-driven phase transitions in frustrated magnets."],"supporting_citations":[{"why":"Supplies the material K2Co2(SeO3)3, its crystal structure, and the easy-axis spin-1/2 dimer picture.","marker":"[20]"},{"why":"Provides the monolayer counterpart K2Co(SeO3)2 used for susceptibility and magnetization comparison that motivates the weak-dimer classification.","marker":"[13]"},{"why":"Gives the theoretical prediction of emergent BKT physics in the triangular bilayer cobaltate that this paper tests experimentally.","marker":"[26]"},{"why":"Establishes the zero-field intermediate BKT phase and the $\\mathbb{Z}_3 \\to \\mathbb{Z}_6$ symmetry enlargement in easy-axis triangular antiferromagnets, anchoring the finite-field BKT interpretation.","marker":"[31]"},{"why":"Supplies the renormalization-group treatment of BKT transitions with symmetry-breaking perturbations, the universality class used to label the phase.","marker":"[32]"},{"why":"Provides prior experimental evidence of a BKT phase in a frustrated magnet, the comparison benchmark for the finite-field BKT claim.","marker":"[38]"},{"why":"Introduces Binder-ratio finite-size scaling, the method used to locate the critical temperatures in the Monte Carlo simulations.","marker":"[47]"},{"why":"Gives the classical triangular-lattice Heisenberg phase diagram in a magnetic field, the basis for the simulated phase structure and correlation-length analysis.","marker":"[50]"}],"fun_headline_variants":["BKT phase appears at finite fields in a triangular magnet","Fractional plateaus map BKT phase in easy-axis triangular magnet","Emergent U(1)⊗S3 symmetry explains BKT phase in weak-dimer magnet","Triangular weak-dimer magnet shows BKT region and 1/3 plateaus","Magnetization steps tie BKT phase to U(1)⊗S3 in bilayer magnet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classical Monte Carlo simulation with parameters $J'_\\perp/J = 1$, $J_\\perp/J = 0.2$, and $\\Delta = 3$ represents the low-field phase structure of the real spin-1/2 material faithfully enough to name the experimental phases.","fun_headline_variants_meta":{"raw":{"variants":["BKT phase appears at finite fields in a triangular magnet","Fractional plateaus map BKT phase in easy-axis triangular magnet","Emergent U(1)⊗S3 symmetry explains BKT phase in weak-dimer magnet","Triangular weak-dimer magnet shows BKT region and 1/3 plateaus","Magnetization steps tie BKT phase to U(1)⊗S3 in bilayer magnet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00137,"raw_usage":{"total_tokens":5596,"prompt_tokens":1030,"completion_tokens":4566,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":4461}},"tokens_in":646,"tokens_out":4566,"duration_ms":32736,"temperature":1.0,"reasoning_tokens":4461,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:50:34.274782+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A neutron-scattering measurement in the region labeled BKT, roughly $B \\simeq 1$–$3$ T and $T \\simeq 0.1$–$0.5$ K, that found Bragg peaks instead of algebraically decaying spin correlations with a temperature-dependent exponent would refute the BKT assignment.","supporting_citations":[{"cited_title":"Kofu, J.-H","cited_arxiv_id":null,"evidence_quote":"Supplies the material K2Co2(SeO3)3, its crystal structure, and the easy-axis spin-1/2 dimer picture."},{"cited_title":"Sheng, J.-W","cited_arxiv_id":null,"evidence_quote":"Provides the monolayer counterpart K2Co(SeO3)2 used for susceptibility and magnetization comparison that motivates the weak-dimer classification."},{"cited_title":"Stre ˇcka, K","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical prediction of emergent BKT physics in the triangular bilayer cobaltate that this paper tests experimentally."},{"cited_title":"Miyashita, Phase transition in spin systems with various types of fluctuations, Proceedings of the Japan Academy, Se- ries B 86, 643 (2010)","cited_arxiv_id":null,"evidence_quote":"Establishes the zero-field intermediate BKT phase and the $\\mathbb{Z}_3 \\to \\mathbb{Z}_6$ symmetry enlargement in easy-axis triangular antiferromagnets, anchoring the finite-field BKT interpretation."},{"cited_title":"Melchy and M","cited_arxiv_id":null,"evidence_quote":"Supplies the renormalization-group treatment of BKT transitions with symmetry-breaking perturbations, the universality class used to label the phase."},{"cited_title":"Chatterjee and X.-G","cited_arxiv_id":null,"evidence_quote":"Provides prior experimental evidence of a BKT phase in a frustrated magnet, the comparison benchmark for the finite-field BKT claim."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces Binder-ratio finite-size scaling, the method used to locate the critical temperatures in the Monte Carlo simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the classical triangular-lattice Heisenberg phase diagram in a magnetic field, the basis for the simulated phase structure and correlation-length analysis."}],"review_version":1}