{"id":"51525c3d-2555-4c22-b406-60fd74b8360f","arxiv_id":"2507.11213","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Magneto-THz spectroscopy on twin-free Na3Co2SbO6 reveals a field-induced spin continuum and pronounced in-plane anisotropy between the a and b axes.","lead":"Using magneto-terahertz spectroscopy on a twin-free crystal of Na3Co2SbO6, researchers observed a continuous band of magnetic excitations when an in-plane magnetic field suppresses the material's ordered spin state. The result adds a new data point in the search for Kitaev quantum spin liquids, materials whose spins never freeze and could one day be useful for quantum computing.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The field-induced 'continuum' could be a superposition of magnon modes from coexisting zigzag domains; the paper explicitly leaves this possibility open, so the central claim is not yet established.","rationale":"The reader's weakest assumption focuses on the Delta=0, diagonal-g-tensor, fully-polarized LSWT model used for high-field fits. That is a valid concern, but it is secondary to the headline: the quantitative parameters support the spin model, but the central observation of a field-induced continuum would stand even if those parameters shifted. The more load-bearing weakness is that the 'continuum' itself may not be an intrinsic continuum. The paper explicitly leaves the two-domain zigzag scenario open, and a superposition of two magnon modes from coexisting magnetic domains would mimic a continuum at q=0. Since THz spectroscopy probes only zone-center excitations, it cannot distinguish an intrinsic continuum from overlapping discrete modes without additional control such as field training or polarization-resolved measurements. The high-field model's fully-polarized assumption being contradicted by the paper's own two-magnon observation is a related internal tension: it means the fit used to extract g-factors and Curie-Weiss differences rests on a state that the data suggest is not fully polarized. Together these issues do not disprove the observations, but they leave the central interpretation conditional. The reader's CONDITIONAL verdict is therefore appropriate; no verdict change is needed.","tokens_in":19844,"tokens_out":6913,"duration_ms":94046,"concrete_test":"Field-train a fresh twin-free crystal: cool through TN in B=0.2 T along a, then measure THz transmission at 2 K for B from 1.0 to 2.0 T with hTHz at several in-plane angles (e.g., 0, 45, 90 degrees from a). If the trained sample shows a single sharp magnon that softens and then reappears, rather than a structureless continuum, the intermediate-field continuum is a domain-superposition effect; if a broad continuum persists with unchanged shape and no sharp peaks at any polarization, the intrinsic-continuum claim survives. This directly tests the one scenario the paper admits it cannot exclude.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the broad THz response between Bc1 and Bc2 is an intrinsic magnetic continuum. The paper's own discussion undermines this: after noting that the data can be reconciled by 'two equivalent zigzag domains tilted from b or a vector sum of them, leading to a double-q', it concedes 'our results have no obvious contradiction with a double-q structure, they don't definitively exclude the two-domain zigzag scenario.' A structurally twin-free crystal can still host two magnetic domains, and if those domains have different field-dependent magnon energies, their overlapping q=0 modes would produce a continuum-like envelope without any intrinsic continuum. The no-domain-repopulation observation is not a test of domain coexistence. Separately, the high-field LSWT fit (SM Eq. S15) assumes a fully polarized state and Delta=0, yet the paper states two-magnon features 'suggest the high-field phase has not been fully polarized'; if the state is partially polarized, Eq. S9 is not the correct dispersion and the fitted ga, gb, and Theta differences in Table I are unsupported. The raw anisotropy remains, but both the continuum interpretation and the quantitative model are conditional.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports magneto-THz spectroscopy on a twin-free single crystal of Na3Co2SbO6. It identifies a 0.47 THz magnon below the N\\'eel temperature, a continuum-like response above TN, and a broad low-energy field-induced response between the two critical fields for both B\\|\\|a (1.3\\u20131.7 T) and B\\|\\|b (0.5\\u20130.8 T). In the high-field regime, well-defined one-magnon and two-magnon modes are observed, and the field dependence of the one-magnon energies is fitted with a phenomenological linear spin-wave model that yields anisotropic g-factors and Curie-Weiss temperature differences. The authors interpret the intermediate-field response as a possible magnetic continuum and emphasize the strong a\\u2013b anisotropy as evidence for a strongly anisotropic spin model.","tokens_in":20031,"tokens_out":8750,"duration_ms":106608,"significance":"If the intermediate-field response is an intrinsic magnetic continuum, this would be a significant result: NCSO is a relatively disorder-free 3d7 honeycomb Kitaev candidate, and the continuum spans a broader field range than in BaCo2(AsO4)2, providing a potentially clean platform for studying field-tunable quantum magnetism. The paper contains solid experimental assets, including the use of a twin-free crystal, clear polarization- and field-dependent data, and a spectral-weight sum-rule check that matches the dc susceptibility. However, the fractionalized-excitation interpretation is not established: the two-domain magnon scenario is explicitly left open, and the high-field quantitative model relies on assumptions that the paper itself calls into question.","major_comments":[{"comment":"The central claim that the intermediate-field response is an intrinsic magnetic continuum is not established because the manuscript explicitly concedes that the two-domain zigzag scenario is not excluded ('our results have no obvious contradiction with a double-q structure, they don\\u2019t definitively exclude the two-domain zigzag scenario'). A structurally twin-free crystal can still host two magnetic domains, and if those domains have different field-dependent q=0 magnon energies, their overlapping modes can produce a broad, continuum-like envelope without any intrinsic continuum. The absence of domain repopulation in the B \\u22a5 hTHz configuration is not a test of domain coexistence. The authors need either a discriminating measurement (for example, the field-training polarimetry they mention) or a quantitative comparison of the predicted two-domain magnon envelope with the observed line shape; absent that, the abstract\\u2019s statement that spin waves 'transform to a magnetic continuum' overstates what the data demonstrate.","section":"Main text, 'To further investigate the in-plane anisotropy' paragraph; SM, 'Magnetic susceptibility with B\\|\\|b'"},{"comment":"The high-field parameter extraction assumes a fully polarized state, a g-tensor diagonal in the crystallographic basis, and the Delta = 0 condition that diagonalizes the q=0 interaction. The main text states that the observation of two-magnon excitations 'suggest[s] the high-field phase has not been fully polarized'; if that is true, Eq. (S9) is not the correct dispersion and the fitted values ga = 6.31, gb = 7.14, and the Curie-Weiss temperature differences in Table I do not have the claimed Hamiltonian interpretation. The Delta = 0 assumption is introduced ad hoc in the SM ('To simplify the discussion, we assume \\u2026 = 0') and is load-bearing: without it JAB is not diagonal and Eq. (S9) does not follow. The authors should either justify these assumptions from independent data or explicitly reframe the fit as a purely phenomenological description.","section":"Supplemental Material, Eqs. (S9), (S15); main text Table I"},{"comment":"The agreement between the fit curves and the measured high-field mode frequencies, and the statement that epsilon_{a,+} vanishes at about 1.72 T 'close to the measured polarizing field', are presented as support for the model. Since the parameters in Table I were obtained by minimizing the squared deviations of these same data in Eq. (S15), the visual agreement is a measure of fit quality rather than an independent verification, and the 1.72 T value is an extrapolation that inherits all of the fit\\u2019s assumptions. This should be stated explicitly, and the 1.72 T value should not be presented as an independent prediction.","section":"Main text, 'Field-dependent mode frequencies' paragraph; SM, text after Eq. (S16)"}],"minor_comments":[{"comment":"The caption says panels (b) and (c) show hTHz \\|\\| a and hTHz \\|\\| b, but the panel labels in the figure place hTHz \\|\\| b on panel (b) and hTHz \\|\\| a on panel (c); please correct the mismatch.","section":"Figure 1 caption and panel labels"},{"comment":"The caption reads 'for magnetic fields along the a axis and a axis'; the second axis should presumably be b.","section":"Supplemental Material, Fig. S1 caption"},{"comment":"Both subcaptions state '(a) hTHz \\u22a5 b. (b) hTHz \\u22a5 b.'; one of them should be hTHz \\|\\| b, or the geometry labels should be corrected to match the figure.","section":"Supplemental Material, Fig. S2 caption"},{"comment":"The phrase 'Temperatures and fields evolution' should be 'Temperature and field evolution' for grammatical clarity.","section":"Figure 2 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid experimental study with valuable data, but the headline interpretation is conditional. The unresolved two-domain scenario and the unsupported assumptions in the high-field fit are addressable through additional analysis or more cautious wording, so I do not see grounds for rejection; however, the abstract currently overstates the certainty of the continuum interpretation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good paper to know about, worth sending out. The genuinely new thing here is clean experimental data: twin-free Na3Co2SbO6, magneto-THz with field along a and b, showing a clear a-b anisotropy and a broad field-induced response in the intermediate field range (1.3–1.7 T for B||a, 0.5–0.8 T for B||b). The spectral weight of the zero-field magnon integrates to the dc susceptibility, which is a real consistency check. The data look careful and the presentation is honest—the authors flag that the continuum's microscopic origin is unresolved and that the two-domain zigzag scenario is not excluded.\n\nThe soft spots are real but not disqualifying. The central claim that the field-induced response is an intrinsic magnetic continuum rests on excluding overlapping magnon modes from two coexisting zigzag domains. The paper observes no domain repopulation, but that does not test coexistence, and the authors concede the point themselves. So the headline interpretation is conditional, even though the raw anisotropy is solid.\n\nThe high-field fit is shakier. The LSWT dispersion assumes a fully polarized state and a diagonal g-tensor, and the SM sets Delta = 0 to make the q = 0 interaction diagonal. That is an ad hoc simplification, and the paper also reports two-magnon features that suggest the state is not fully polarized at the fields used. Fitted ga, gb, and the Theta differences come without uncertainties, and the agreement is not an independent test of the model because the same model produces the fitted curves. Still, the model is presented as phenomenological, and the extracted values are roughly consistent with earlier susceptibility fits, so this is a minor-to-moderate weakness rather than a fatal one.\n\nWho is this for? Anyone working on cobalt honeycomb Kitaev candidates or field-tunable quantum magnets. The experimental results will be a useful reference even if the continuum interpretation later turns out to be domain-related. I would send this to a serious referee. The main revision pressure should be on the interpretation language and on adding uncertainties to the fits; the data themselves deserve publication.","headline":"Careful new THz data on twin-free Na3Co2SbO6 show robust a-b anisotropy and a field-induced broad response, but the continuum interpretation is not yet established—and the paper says so itself.","tokens_in":20639,"tokens_out":2377,"would_cite":true,"duration_ms":28669,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In twin-free Na3Co2SbO6, an in-plane magnetic field transforms the zero-field 0.5 THz magnon into a magnetic continuum over an intermediate field range, with strong a-b anisotropy.","keywords":["Na3Co2SbO6","magneto-THz spectroscopy","spin continuum","Kitaev quantum spin liquid","magnetic anisotropy","spin-wave theory","fractionalized excitations","twin-free crystal"],"falsifier":"Fit the high-field one-magnon energies without imposing the simplifying cancellation among exchange couplings and without assuming the g-tensor is diagonal: if a model with a nonzero value of that combination or a rotated g-tensor matches the data equally well, then the extracted g-factors and Curie-Weiss temperature differences do not uniquely support the proposed anisotropic Hamiltonian, and the description of the high-field state as a fully polarized spin-wave state would be falsified.","tokens_in":19575,"feed_emoji":"🧲","tokens_out":18166,"duration_ms":185894,"temperature":0.7,"pith_summary":"Using magneto-terahertz spectroscopy on a twin-free crystal of the honeycomb magnet Na3Co2SbO6, the paper traces what happens to the 0.5 THz (2 meV) spin-wave mode as an in-plane magnetic field is applied. In an intermediate field range, 1.3 to 1.7 T for field along the a axis and 0.5 to 0.8 T for field along the b axis, the sharp magnon disappears and is replaced by a broad magnetic continuum before the system enters a spin-polarized state at higher fields. The same measurements show a strong difference between the a and b crystallographic axes in both magnon intensity and field threshold, which the authors take as evidence that any successful spin model for this material must be strongly anisotropic. Because the crystals are twin-free, the continuum cannot be blamed on twinning domains, strengthening the case that it reflects intrinsic, possibly fractionalized, excitations in a Kitaev-candidate material. The microscopic origin of the continuum is left open.","feed_headline":"Field melts a honeycomb magnet's magnons into a broad continuum","feed_subtitle":"In twin-free Na3Co2SbO6 a 0.5 THz spin wave becomes a continuum between 1.3 and 1.7 T.","key_machinery":"The central experimental object is the imaginary part of the dynamic magnetic susceptibility, $\\chi_2(\\omega)$, obtained from time-domain THz transmission through a twin-free crystal; it gives direct access to $q = 0$ spin excitations with polarization selection rules controlled by the THz magnetic field. The theoretical engine is a $q = 0$ linear spin-wave treatment of a fully polarized honeycomb magnet with an anisotropic g-tensor, which produces the magnon energies $\\epsilon_{a,+}$ and $\\epsilon_{b,+}$ in terms of the applied field and the Curie-Weiss temperatures $\\Theta_a$, $\\Theta_b$, and $\\Theta_{c^*}$. To make the $q = 0$ interaction diagonal in the crystallographic basis, the paper assumes that the particular parameter combination $\\Delta = (K_X - K_Z) - (\\Gamma_X + \\Gamma'_X + \\Gamma'_Z) + (2\\Gamma''_X + \\Gamma_Z)$ vanishes; under that assumption the measured field dependence fixes $g_a$, $g_b$, and the Curie-Weiss temperature differences, which in turn constrain a strongly anisotropic Kitaev-Heisenberg-Gamma Hamiltonian.","core_discovery":"At zero field and below the Neel temperature of 6.6 K, the zone-center spin excitation spectrum of Na3Co2SbO6 consists of a single magnon at 0.5 THz (2 meV) whose spectral weight accounts for the full dc magnetic susceptibility, while above the ordering temperature a low-energy continuum appears. The central experimental discovery is that an in-plane field does not simply soften and close this magnon: over a wide intermediate range, 1.3 to 1.7 T for B along a and 0.5 to 0.8 T for B along b, the magnon is replaced by a magnetic continuum, and only above the upper critical field do sharp spin waves reappear in a spin-polarized state. The field evolution differs strongly for the two in-plane axes, and the magnon intensities for THz polarization along a and b differ by roughly a factor of two, which the paper argues rules out a simple zigzag state with spins along b and favors a double-q structure or two tilted zigzag domains. In the high-field phase the one-magnon modes fit a linear spin-wave model with anisotropic g-factors ($g_a = 6.31$, $g_b = 7.14$) and Curie-Weiss temperature differences ($\\Theta_b - \\Theta_a = 3.66$ K, $\\Theta_b - \\Theta_{c^*} = 9.78$ K), while the coexistence of two-magnon features shows the high-field state is not fully polarized.","pith_inferences":["Extension: a direct test of the fractionalized-excitation scenario would be to extend measurements below 0.2 THz; the paper's sum-rule analysis leaves room for a low-energy asymmetric peak near $\\omega/J_1 \\sim 0.1$ predicted by parton mean-field theory, and finding or failing to find it would discriminate between that scenario and a multi-magnon origin.","Extension: the high-field fit could be redone without setting the exchange combination $\\Delta$ to zero and without assuming a diagonal g-tensor; if a model with nonzero $\\Delta$ or a rotated g-tensor fits equally well, the extracted g-factors and Curie-Weiss differences would no longer pin down the proposed anisotropic Hamiltonian.","Extension: measuring the same spectra on deliberately twinned crystals of NCSO would quantify how much of the continuum is intrinsic rather than disorder-related, since the twin-free choice already removes twinning as a source of spectral continua."],"forward_implications":["If the observations are correct, Na3Co2SbO6 becomes a model system for field-tunable quantum magnetism in a clean twin-free honeycomb lattice, with a continuum window much wider than the one seen in BaCo2(AsO4)2.","The strong a-b anisotropy rules out the simple zigzag spin structure with moments along b, so any acceptable spin Hamiltonian must be strongly anisotropic and must reproduce a double-q structure or two tilted zigzag domains.","The continuum above the ordering temperature and inside the intermediate field range is consistent with fractionalized excitations of a Kitaev-type spin liquid, although an XXZ-$J_1$-$J_3$ spinon scenario remains open and the paper does not claim to distinguish them.","The high-field phase is not truly fully polarized because two-magnon excitations persist, so the linear spin-wave fits describe only the dominant one-magnon response.","The spectral-weight analysis shows the zero-field magnon exhausts the dc magnetic susceptibility, implying the THz continuum appears through field- and temperature-driven changes in the spectrum rather than coexisting with a sharp magnon in the ordered state."],"supporting_citations":[{"why":"Supplies the twin-free NCSO crystals, the in-plane C2 symmetry, and the magnetization evidence that an in-plane field suppresses magnetic order with the critical fields used throughout.","marker":"[48]"},{"why":"Provides the BaCo2(AsO4)2 continuum spectrum that the NCSO continuum is compared with in profile, energy scale, and temperature dependence.","marker":"[54]"},{"why":"Establishes the double-q magnetic structure or two equivalent zigzag domains and the in-plane anisotropy that the THz polarization results are interpreted against.","marker":"[45]"},{"why":"Proposes the zigzag spin structure with spins along b whose predicted selection rules the measured magnon intensities are used to exclude.","marker":"[41]"},{"why":"Gives the XXZ-$J_1$-$J_3$ model with low-energy spinon excitations and the parton mean-field prediction of a low-energy THz peak, the main alternative interpretation of the continuum.","marker":"[47]"},{"why":"Supplies the Kitaev spin-liquid theory in which a continuum of fractionalized excitations is the diagnostic used to motivate the quantum-spin-liquid interpretation.","marker":"[55]"},{"why":"Provides the time-domain THz magneto-spectroscopy method and the extraction of the imaginary magnetic susceptibility used in all measurements.","marker":"[59]"},{"why":"Supplies the spectral-weight sum rule used to show that the zero-field magnon accounts for the full dc magnetic susceptibility.","marker":"[56]"},{"why":"Gives the relation between Curie-Weiss temperatures and zone-center exchange parameters that converts fitted magnon energies into interaction constraints.","marker":"[69]"},{"why":"Defines the anisotropic Kitaev-Heisenberg-Gamma model whose exchange parameters the fitted g-factors and Curie-Weiss differences are meant to constrain.","marker":"[70]"}],"fun_headline_variants":["Field melts honeycomb magnet's magnon into spin continuum","Twin-free Na3Co2SbO6 shows anisotropic field-induced continuum","Magnon turns to continuum under in-plane field in honeycomb magnet","Anisotropic spin continuum emerges in field-tuned Na3Co2SbO6"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that in the high-field phase the spins are almost fully polarized along the field and the interaction tensor at zero momentum is diagonal in the crystal axes, so that a particular combination of exchange couplings can be set to zero; the observed two-magnon excitations show this polarized limit is not fully reached, and if the assumption is relaxed the fitted g-factors and Curie-Weiss differences may change.","fun_headline_variants_meta":{"raw":{"variants":["Field melts honeycomb magnet's magnon into spin continuum","Twin-free Na3Co2SbO6 shows anisotropic field-induced continuum","Magnon turns to continuum under in-plane field in honeycomb magnet","Anisotropic spin continuum emerges in field-tuned Na3Co2SbO6"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000843,"raw_usage":{"total_tokens":3758,"prompt_tokens":1115,"completion_tokens":2643,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":731,"completion_tokens_details":{"reasoning_tokens":2564}},"tokens_in":731,"tokens_out":2643,"duration_ms":25297,"temperature":1.0,"reasoning_tokens":2564,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:13:11.888517+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the high-field one-magnon energies without imposing the simplifying cancellation among exchange couplings and without assuming the g-tensor is diagonal: if a model with a nonzero value of that combination or a rotated g-tensor matches the data equally well, then the extracted g-factors and Curie-Weiss temperature differences do not uniquely support the proposed anisotropic Hamiltonian, and the description of the high-field state as a fully polarized spin-wave state would be falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the twin-free NCSO crystals, the in-plane C2 symmetry, and the magnetization evidence that an in-plane field suppresses magnetic order with the critical fields used throughout."},{"cited_title":"Zhang, Y","cited_arxiv_id":null,"evidence_quote":"Provides the BaCo2(AsO4)2 continuum spectrum that the NCSO continuum is compared with in profile, energy scale, and temperature dependence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the double-q magnetic structure or two equivalent zigzag domains and the in-plane anisotropy that the THz polarization results are interpreted against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proposes the zigzag spin structure with spins along b whose predicted selection rules the measured magnon intensities are used to exclude."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the XXZ-$J_1$-$J_3$ model with low-energy spinon excitations and the parton mean-field prediction of a low-energy THz peak, the main alternative interpretation of the continuum."},{"cited_title":"Knolle, D","cited_arxiv_id":null,"evidence_quote":"Supplies the Kitaev spin-liquid theory in which a continuum of fractionalized excitations is the diagnostic used to motivate the quantum-spin-liquid interpretation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the time-domain THz magneto-spectroscopy method and the extraction of the imaginary magnetic susceptibility used in all measurements."},{"cited_title":"Little, L","cited_arxiv_id":null,"evidence_quote":"Supplies the spectral-weight sum rule used to show that the zero-field magnon accounts for the full dc magnetic susceptibility."},{"cited_title":"Lampen-Kelley, S","cited_arxiv_id":null,"evidence_quote":"Gives the relation between Curie-Weiss temperatures and zone-center exchange parameters that converts fitted magnon energies into interaction constraints."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the anisotropic Kitaev-Heisenberg-Gamma model whose exchange parameters the fitted g-factors and Curie-Weiss differences are meant to constrain."}],"review_version":1}