{"id":"38634ac3-9d26-4b9d-b0f1-af204b940c94","arxiv_id":"1908.09510","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"CdYb2Se4 has a zero-field Γ5 antiferromagnetic ground state that crosses over to a splayed ferromagnet at 3 T, with fitted exchanges placing it near a phase boundary.","lead":"This paper shows that the frustrated ytterbium magnet CdYb2Se4 orders into an antiferromagnetic state at zero field and switches to a splayed ferromagnet when a magnetic field of 3 T is applied. The authors use neutron scattering and an anisotropic exchange model to argue the compound sits near a phase boundary, much like the Yb pyrochlore magnets.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed proximity to the Γ5/splayed-ferromagnet phase boundary is underconstrained by the data: the fitted exchange parameters depend on a non-unique CEF ground doublet and on qualitative matching of two powder Bragg curves.","rationale":"The experimental core of the paper is credible: powder neutron diffraction clearly shows a zero-field Γ5 state and a field-induced crossover near 3 T, and the low-energy inelastic signal at 0.72 meV is shown to be magnetic. The reader's CONDITIONAL verdict and the critique of underdetermination are appropriate. My stress-test focuses more specifically on the parameter-fitting chain: the exchange ratios in Eq. (4) are derived from a non-unique CEF ground doublet and from qualitative matching of powder Bragg intensities, with the paper itself conceding that quantitative agreement should not be expected. This undermines the strength of the phase-boundary proximity claim, but it does not invalidate the experimental findings or the qualitative model scenario. The proposed test—refitting with the alternative CEF parameterization—would settle whether the proximity claim is robust. Since the reader already assigned CONDITIONAL and I do not see grounds to move more than that, the verdict remains UNCHANGED.","tokens_in":14478,"tokens_out":7141,"duration_ms":78157,"concrete_test":"Re-run the exchange parameter fits using the alternative CEF parameter set of Higo et al. for CdYb2S4, quoted in Table I (g‖ = 2.67, g⊥ = 1.33), while keeping the same powder-averaging formulas (Eqs. (7)-(8)) and the same fitting protocol for D/J, K/J, and Γ/J. If the resulting best-fit parameters move significantly away from the Γ5/SF phase boundary, or if the region of (K/J, Γ/J, D/J) consistent with the powder data straddles the boundary, then the claimed proximity is an artifact of the non-unique CEF fit rather than a robust conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that CdYb2Se4 lies close to the Γ5/SF phase boundary rests on the fitted exchange ratios K/J = Γ/J = -0.07 and D/J = -0.3 in Eq. (4). Section IV.B.2 describes the fitting procedure: D/J is fixed by matching the position of the broad, resolution-broadened, powder-averaged intensity maximum near 0.7 meV (Eqs. (6)-(7)), while K/J and Γ/J are tuned together to reproduce 'qualitatively' the field dependence of the (200) and (111) Bragg intensities. The text explicitly warns that 'some details, such as the range of the crossover field or evolution of relative intensity of the peaks, are sensitive to the precise values of the g-factors, and the small symmetric anisotropies K and Γ', and that 'a quantitative agreement between theory and experiment should not be expected.' At the same time, Section IV.A admits the CEF parameter space fitting the INS data is 'rather wide' and that the ground-doublet composition is not unambiguous. Since the Zeeman term in Eq. (2) and the field-driven crossover depend on the g-tensor derived from this CEF fit, the exchange parameters are not uniquely pinned by the data. In particular, the proximity to the Γ5/SF boundary is essentially one particular point in a degenerate parameter region, so the claim of phase-boundary proximity is a plausible interpretation rather than an established result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a combined experimental and modeling study of the spinel CdYb2Se4, a pyrochlore-like magnet of Yb3+ ions. High-energy inelastic neutron scattering is used to determine crystal-electric-field (CEF) parameters, yielding a well-isolated ground doublet with a weakly Ising character (g||=3.67, g⊥=2.11). Powder neutron diffraction at 0.45 K identifies zero-field k=0 magnetic order in the Γ5 irreducible representation (ψ2/ψ3), and shows that an applied field drives a crossover near Bc=3 T to a Γ9 splayed ferromagnet. Low-energy inelastic neutron scattering reveals a weakly dispersive magnetic band near 0.72 meV that weakens with field and vanishes above 3 T. The authors model the low-energy physics with a nearest-neighbour anisotropic exchange Hamiltonian for effective S=1/2 Kramers doublets, fitting J, D, K, and Γ. They conclude that the estimated exchanges place CdYb2Se4 close to the phase boundary between the Γ5 and splayed-ferromagnet states, analogous to the Yb pyrochlores.","tokens_in":14731,"tokens_out":5236,"duration_ms":52636,"significance":"The experimental characterization is a valuable contribution: the identification of the Γ5 zero-field state and the Γ9 splayed ferromagnet at high field is directly supported by Rietveld refinements, including a March preferred-orientation model, and the 0.72 meV magnetic band is clearly observed in the inelastic data. The work extends the study of frustrated XY pyrochlore magnets to the spinel structure and suggests a new material for exploring Γ5/SF competition. However, the central quantitative claim of proximity to the Γ5/splayed-ferromagnet phase boundary is not established by the data: the exchange parameters are fitted to the same spectra and intensities that are later shown as agreement, and the CEF parameter space—and hence the g-tensor entering the Zeeman term—is admitted to be non-unique. The paper is strongest as an experimental report; the theoretical interpretation is plausible but underconstrained.","major_comments":[{"comment":"The exchange parameters in Eq. (4) are fitted to the same data that are later presented as agreement: D/J=-0.3 is fixed by the position of the intensity maximum in the zero-field inelastic spectrum, and K/J=Γ/J=-0.07 is tuned to the field dependence of the (200) and (111) Bragg intensities. The subsequent matching of the calculated spectra and intensities is therefore a consistency check, not an independent test of the model. Since the claimed proximity to the Γ5/splayed-ferromagnet phase boundary rests entirely on these fitted values, the central claim is underconstrained by the data.","section":"§IV.B.2, Eq. (4)"},{"comment":"The authors state in §IV.A that the CEF parameter space fitting the INS data is 'rather wide' and that the ground-doublet composition is not unambiguous; they further state in §IV.B.2 that 'a quantitative agreement between theory and experiment should not be expected.' Because the Zeeman term in Eq. (2) and the field-driven crossover depend on the g-tensor derived from this CEF fit, the exchange ratios in Eq. (4) are not uniquely pinned. The proximity to the Γ5/SF phase boundary is thus one particular point in a degenerate parameter region rather than an established result; the authors should either demonstrate robustness of the phase-boundary claim to variations in the CEF parameters or re-frame the claim as a plausible interpretation.","section":"§IV.A and §IV.B.2"},{"comment":"The determination of K/J and Γ/J is qualitative: K and Γ are assumed equal, the fit targets 'qualitatively' the field dependence of the upturns in two Bragg peaks, and the low-field absolute intensities are explicitly inconsistent with the model. This does not yield a quantitative estimate of the exchange parameters or their uncertainties. The authors should show how the deduced phase-boundary proximity depends on the K=Γ assumption and on the fitted g-factors, or soften the conclusion accordingly.","section":"§IV.B.2"}],"minor_comments":[{"comment":"The abstract uses the undefined macro '\\CYS'; it should read 'CdYb2Se4' explicitly.","section":"Abstract"},{"comment":"The word 'neiboughring' in the caption should be 'neighboring'.","section":"Table II caption"},{"comment":"The unit-vector notation in the powder-averaging integrals is not defined; please define \\hat{B} and \\hat{Q} explicitly and mention that the average over \\hat{B} corresponds to averaging over crystallite orientations in a powder.","section":"§IV.B.2, Eqs. (7)-(8)"},{"comment":"The figure caption does not label the left, middle, and right panels; please add panel labels for clarity.","section":"Figure 7"},{"comment":"The sentence 'R=Yb, Er and B = Sn, Ti, Ge' appears to mix the rare-earth site and the B site; it should be clarified, since for Yb2B2O7 the rare earth is fixed to Yb and B is the non-magnetic cation.","section":"§V Outlook"}],"recommendation":"major_revision","confidential_remarks":"The experimental results are sound and well presented, and the paper would be a worthwhile addition to the literature after the theoretical claims are appropriately limited. The main risk is overclaiming the phase-boundary proximity based on a fit with admitted qualitative character and non-unique CEF input. A revision that honestly frames the exchange parameters as 'consistent with' rather than 'predicting' the proximity, and that adds a robustness analysis, would make the paper acceptable. The scope of the journal is appropriate for this combined experimental-modeling study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the thing to know: this is a solid experimental paper from a strong group, and the new data are the field-driven Γ5-to-splayed-ferromagnet crossover near 3 T and the 0.72 meV excitation band. The zero-field Γ5 order was already reported in Refs. 13 and 14, so the genuinely new content is the field evolution, the CEF scheme with a weakly Ising doublet, and the fitted exchange parameters. On the experimental side, the neutron diffraction refinements look careful, and the paper is appropriately cautious about the powder-averaging limitations and the need for single-crystal confirmation.\n\nThe soft spot is the theoretical interpretation. The exchange parameters are not determined independently: D/J is fixed by matching the position of the powder-averaged inelastic maximum, K/J and Γ/J are tuned to reproduce qualitatively the field dependence of Bragg intensities, and the CEF fit that gives the g-tensor is admitted to be non-unique. So the claim that CdYb2Se4 sits close to the Γ5/splayed-ferromagnet phase boundary is plausible but not pinned down. The paper itself says quantitative agreement should not be expected, which is honest but also a sign that the fit is not a strong test.\n\nThe powder-averaging treatment in Eqs. (7)-(8), including only Q perpendicular to B for each sampled field direction, is a significant approximation; it works for the qualitative features but could misrepresent the crossover field. This needs single-crystal data to validate.\n\nNone of this kills the paper. The experimental results are likely correct and will be useful to the community working on Yb pyrochlores and spinels. The model comparison is reasonable as a consistency check, not as a determination. I'd recommend sending it to a specialist journal like PRB; a good referee should ask for a more systematic treatment of the parameter uncertainties and a clearer separation between what is measured and what is inferred. I would cite the experimental findings, though I'd be careful about quoting the exchange parameters as established.\n\nBring to reading group? Yes, it's a good example of how powder data plus a model can be both useful and underdetermined.","headline":"Solid experimental study of a Yb spinel with a new field-induced crossover, but the exchange parameter fit is underdetermined and the phase-boundary claim should be read as plausible, not proven.","tokens_in":15465,"tokens_out":2593,"would_cite":true,"duration_ms":25750,"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":"The paper argues that a 3-tesla field switches CdYb2Se4 from an XY antiferromagnet to a splayed ferromagnet.","keywords":["frustrated magnetism","pyrochlore lattice","ytterbium spinel","CdYb2Se4","anisotropic exchange","Gamma-5 antiferromagnet","splayed ferromagnet","neutron scattering"],"falsifier":"A single-crystal neutron diffraction experiment in an applied field could settle the claim without powder averaging: the model predicts a crossover in the Bragg intensities near 2-3 T and a specific angle-dependent evolution of the $\\Gamma_5\\to\\Gamma_9$ conversion. If the observed crossover field, moment directions, or the disappearance of the 0.72 meV band above 3 T disagreed with the calculated powder-averaged spectra for $K/J=\\Gamma/J=-0.07$, $D/J=-0.3$, the parameter set and the claimed proximity to the phase boundary would be ruled out.","tokens_in":14188,"feed_emoji":"🧲","tokens_out":9656,"duration_ms":79835,"temperature":0.7,"pith_summary":"The paper aims to establish that the low-temperature magnetism of the spinel CdYb2Se4 is controlled by a competition between two nearly degenerate ordered states. At zero field the Yb$^{3+}$ moments form a $k=0$ XY-type antiferromagnet of $\\Gamma_5$ symmetry, and an applied magnetic field of about 3 T crosses the system over to a splayed ice-like ferromagnet of $\\Gamma_9$ symmetry. The authors argue that both states are captured by a nearest-neighbour anisotropic exchange model for effective $S=1/2$ Kramers doublets on the pyrochlore lattice, with parameters that place CdYb2Se4 close to the phase boundary between the two orders. The weakly dispersive magnetic band at 0.72 meV is proposed as the experimental signature of that proximity, since it vanishes when the field suppresses the $\\Gamma_5$ state. A sympathetic reader would care because the same phase competition is believed to underlie the unusual behaviour of the ytterbium pyrochlores, and CdYb2Se4 offers a chemically distinct setting to test it.","feed_headline":"3-tesla field flips CdYb2Se4 between two magnetic orders","feed_subtitle":"A spinel with near-degenerate ordered states sits on the boundary where antiferromagnet meets splayed ferromagnet.","key_machinery":"The load-bearing object is the nearest-neighbour anisotropic exchange Hamiltonian for effective $S=1/2$ Kramers doublets on the pyrochlore lattice, $H=\\sum_{\\langle ij\\rangle} S_i^\\top J_{ij} S_j - \\mu_B B \\cdot \\sum_i g_i S_i$, with four couplings: Heisenberg $J$, Kitaev $K$, symmetric off-diagonal $\\Gamma$, and Dzyaloshinskii-Moriya $D$. The paper's quantitative conclusions rest on a powder-averaging procedure (Eqs. (7) and (8)) in which, for each sampled field direction, only wave-vectors perpendicular to the field contribute to the computed Bragg and inelastic intensities. Within that scheme, $D/J=-0.3$ is fixed by the position of the flat band at 0.72 meV, while $K/J=\\Gamma/J=-0.07$ is tuned to the field dependence of the (002) and (111) Bragg peaks; the resulting spectrum is computed in linear spin-wave theory and broadened to match the experimental resolution.","core_discovery":"The paper reports that CdYb2Se4 orders below $T_N = 1.8$ K into a $\\Gamma_5$ antiferromagnetic state with an ordered moment of $0.634(5)\\,\\mu_B$/Yb, much smaller than the single-ion value of $1.37\\,\\mu_B$, and that a field near $B_c=3$ T converts this into a $\\Gamma_9$ splayed ferromagnet whose ferromagnetic component reaches $0.94(1)\\,\\mu_B$/Yb at 5 T. The zero-field inelastic spectrum shows a weakly dispersive band at $\\hbar\\omega\\approx0.72$ meV that disappears above 3 T. The authors attribute this behaviour to an effective nearest-neighbour anisotropic exchange model with $J\\sim7$ K, $K/J=\\Gamma/J=-0.07$, and $D/J=-0.3$, which places the material close to the $\\Gamma_5$/splayed-ferromagnet phase boundary. They further show that the field dependence of the magnetic Bragg intensities and the disappearance of the 0.72 meV band follow from the same model once powder averaging over field orientations is accounted for.","pith_inferences":["If the model is right, CdYb2Se4 is a tunable laboratory for order-by-disorder: single-crystal scattering should reveal which $\\Gamma_5$ state (likely $\\psi_3$) is selected and whether the selection changes under field.","The discrepancy between the zero-field ordered moment (0.63 $\\mu_B$) and the single-ion expectation (1.37 $\\mu_B$) may quantify quantum fluctuations near the phase boundary; a similar reduction is expected in the spin correlations, which could be searched for in diffuse scattering.","The same parameter-extraction scheme could be applied to other $T$Yb$_2X_4$ spinels to map how the $\\Gamma_5$/splayed-ferromagnet boundary shifts with the crystal-field $g$-tensor, providing a chemical series in which to test the universality of the Yb-pyrochlore physics."],"forward_implications":["At zero field, the ground state is a $\\Gamma_5$ XY antiferromagnet with a strongly reduced ordered moment, consistent with strong quantum fluctuations.","A field of about 3 T (lower at elevated temperature) drives a crossover to a $\\Gamma_9$ splayed ice-like ferromagnet, and the $\\Gamma_5$ Bragg contribution vanishes at $B_c$.","The 0.72 meV excitation band, present at zero field, weakens with field and vanishes above 3 T, providing a direct spectroscopic marker of the phase competition.","The estimated exchanges $K/J=\\Gamma/J=-0.07$, $D/J=-0.3$ with $J\\sim7$ K place CdYb2Se4 in the same region of the phase diagram as the ytterbium pyrochlores.","Quantum (or other) order-by-disorder is expected to select one of the degenerate $\\Gamma_5$ states, likely $\\psi_3$, although powder data cannot distinguish $\\psi_2$ from $\\psi_3$."],"supporting_citations":[{"why":"It supplies the super-exchange argument used to expect $J>0$, $D<0$, and small $K$ and $\\Gamma$, placing the material near the $\\Gamma_5$/splayed-ferromagnet boundary.","marker":"[8]"},{"why":"It established the magnetic behaviour of TYb$_2$X$_4$ spinels and provided the crystal-field starting parameters used in the fit.","marker":"[13]"},{"why":"It identified the $k=0$ magnetic order and the $\\Gamma_5$ $\\psi_2/\\psi_3$ basis vectors that the zero-field refinement is compared with.","marker":"[14]"},{"why":"It gives experimentally determined exchange parameters in a related breathing pyrochlore that motivate and calibrate the chosen $K/J=\\Gamma/J$ values.","marker":"[17]"},{"why":"It describes the splayed ice-like ferromagnet state in Tb$_2$Sn$_2$O$_7$ used as the reference for the $\\Gamma_9$ high-field order.","marker":"[24]"},{"why":"It defines the parametrization of the pyrochlore nearest-neighbour anisotropic exchange model used in Eq. (2).","marker":"[26]"}],"fun_headline_variants":["CdYb2Se4 flips from AFM to splayed ferromagnet at 3 T","3 T field tunes CdYb2Se4 through magnetic phase competition","CdYb2Se4 near critical field: antiferro to splayed ferro","Quantum spinel CdYb2Se4 order changes at 3 T","Field-driven transition in CdYb2Se4: AFM to splayed ferromagnet"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The modelling stands on the assumption that the powder-averaged neutron intensities are faithfully represented by wave-vectors perpendicular to the field for each sampled field direction, and that the crystal-field parameters are accurate enough to constrain the anisotropic $g$-tensor; the paper itself notes that the crystal-field parameter space is wide and the ground-state doublet composition is not unambiguous.","fun_headline_variants_meta":{"raw":{"variants":["CdYb2Se4 flips from AFM to splayed ferromagnet at 3 T","3 T field tunes CdYb2Se4 through magnetic phase competition","CdYb2Se4 near critical field: antiferro to splayed ferro","Quantum spinel CdYb2Se4 order changes at 3 T","Field-driven transition in CdYb2Se4: AFM to splayed ferromagnet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001095,"raw_usage":{"total_tokens":4620,"prompt_tokens":1044,"completion_tokens":3576,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":3464}},"tokens_in":660,"tokens_out":3576,"duration_ms":23863,"temperature":1.0,"reasoning_tokens":3464,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:09:51.767629+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single-crystal neutron diffraction experiment in an applied field could settle the claim without powder averaging: the model predicts a crossover in the Bragg intensities near 2-3 T and a specific angle-dependent evolution of the $\\Gamma_5\\to\\Gamma_9$ conversion. If the observed crossover field, moment directions, or the disappearance of the 0.72 meV band above 3 T disagreed with the calculated powder-averaged spectra for $K/J=\\Gamma/J=-0.07$, $D/J=-0.3$, the parameter set and the claimed proximity to the phase boundary would be ruled out.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the super-exchange argument used to expect $J>0$, $D<0$, and small $K$ and $\\Gamma$, placing the material near the $\\Gamma_5$/splayed-ferromagnet boundary."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It established the magnetic behaviour of TYb$_2$X$_4$ spinels and provided the crystal-field starting parameters used in the fit."},{"cited_title":"Dalmas de R\\' e otier, C","cited_arxiv_id":null,"evidence_quote":"It identified the $k=0$ magnetic order and the $\\Gamma_5$ $\\psi_2/\\psi_3$ basis vectors that the zero-field refinement is compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives experimentally determined exchange parameters in a related breathing pyrochlore that motivate and calibrate the chosen $K/J=\\Gamma/J$ values."},{"cited_title":"Mirebeau, A","cited_arxiv_id":null,"evidence_quote":"It describes the splayed ice-like ferromagnet state in Tb$_2$Sn$_2$O$_7$ used as the reference for the $\\Gamma_9$ high-field order."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the parametrization of the pyrochlore nearest-neighbour anisotropic exchange model used in Eq. (2)."}],"review_version":1}