{"id":"8c1665f6-12b8-43ee-9fba-daa67c9e8ae3","arxiv_id":"2608.03365","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Nd3Li3W2O12 is a frustrated hyperkagome antiferromagnet with a Kramers doublet effective spin-1/2 ground state and a crystal electric field gap of about 9 meV, as determined by neutron scattering and bulk thermodynamics.","lead":"This paper characterizes a new neodymium-based magnet, a garnet with a frustrated hyperkagome lattice, and finds it does not order magnetically down to 0.1 K. It maps the crystal electric field levels with neutron scattering and shows the low-energy physics is an effective spin-1/2 doublet.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CEF parameter set selected using the same bulk data later claimed as validation; uniqueness and uncertainties are not established.","rationale":"The reader's weakest assumption identifies the selection of one CEF parameter set from nine equally good INS fits using the same thermodynamic data that are later presented as confirmation. This is exactly the concern I find most load-bearing. The paper's central claim has two parts: the experimentally observed CEF excitations (which are directly visible in the INS data and are robust) and the specific CEF Hamiltonian and derived single-ion properties (which depend on the selected parameter set). The selection process is described only qualitatively ('best overall agreement'), and the other eight solutions are not provided, making it impossible to judge whether the bulk data actually break the degeneracy. The absence of uncertainties on the reported CEF parameters compounds this issue. My proposed test would quantify whether the thermodynamic data uniquely select a single parameter set; if they do not, the confirmation is circular and the specific wavefunctions/g-tensor components cannot be considered definitive. The energy-level scheme and the J_eff=1/2 ground-state gap are still well supported, so the verdict remains CONDITIONAL, unchanged from the reader's assessment.","tokens_in":21273,"tokens_out":8542,"duration_ms":97711,"concrete_test":"Re-run the CEF fit to the INS spectra from many random starting points to map all local minima that reproduce the INS data to similar quality. For each distinct minimum, compute χ_CEF(T), M_CEF(H), and C_CEF(T) and compare them to the experimental data using a quantitative metric (e.g., reduced χ²). If more than one parameter set falls within the experimental uncertainty, then the bulk data do not uniquely select the reported set, and the validation in Sec. IIID is circular. As a corollary, report the spread in the ground-state g-tensor components and the first-excitation gap across all equally good sets.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the underdetermined selection of the CEF parameter set. Sec. IIID reports that nine distinct parameter sets fit the INS spectra equally well, but only one is presented (Table I). This set is chosen by comparing the calculated magnetic susceptibility, isotherms, and heat capacity to the experimental data. Because the same thermodynamic data are used for selection, their subsequent agreement (Figs. 5–7) is not an independent confirmation; it is a discrete fit. The paper does not report the other eight sets, so one cannot assess whether the bulk data actually discriminate among them. If several sets reproduce the thermodynamic data within error, the reported CEF parameters, wavefunctions, and g-tensor components (Tables I–II, Eq. 7) are not unique. This does not threaten the directly observed CEF energy levels (9.3, 21.9, 25.6, 91.8 meV) or the conclusion that the ground state is a J_eff = 1/2 Kramers doublet with a ~9 meV gap, but it weakens the stronger claim that the fitted Hamiltonian is the correct single-ion model. The absence of parameter uncertainties further prevents assessing the significance of the reported g-values.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports bulk magnetization, heat-capacity, and powder inelastic-neutron-scattering (INS) measurements on the previously unexplored hyperkagome garnet Nd3Li3W2O12. Four non-dispersive CEF excitations are observed at about 9.3, 21.9, 25.6, and 91.8 meV and are attributed to transitions from a Kramers ground doublet to excited doublets. A nine-parameter D2-symmetric Stevens CEF Hamiltonian is fitted simultaneously to INS spectra at 6, 100, and 250 K. The authors report that nine distinct parameter sets reproduce the INS data equally well, and they select one set by comparing calculated susceptibility, magnetization, and heat capacity with experiment. They then present the selected set as reproducing those same bulk data, provide CEF eigenfunctions and g-tensor components, and conclude that the low-temperature physics is governed by a J_eff = 1/2 Kramers doublet with a ~9 meV gap. No magnetic long-range order is found down to 0.1 K.","tokens_in":21639,"tokens_out":4293,"duration_ms":49700,"significance":"If the CEF model is considered established, the paper is a useful single-ion reference for a new Nd-based frustrated garnet and adds to the small family of hyperkagome magnets with a well-defined low-energy spin-1/2 description. The directly observed INS level energies and the temperature-dependent intensity evolution are valuable, and the bulk probes (low-T Curie-Weiss, saturation magnetization, Schottky gap, entropy) independently corroborate a J_eff = 1/2 ground state. However, the specific CEF parameter set, the associated wavefunctions, and the g-tensor anisotropy are not independently determined, because the same thermodynamic data used for selection are later presented as validation. The nine-fold ambiguity and the absence of parameter uncertainties therefore weaken the stronger single-ion-model claims while leaving the main physical conclusion—an effective spin-1/2 low-temperature description—largely intact.","major_comments":[{"comment":"The CEF parameter set is underdetermined by INS: the text states that '9 distinct sets of CEF parameters reproduce the INS data equally well', but only one set (Table I) is shown. The selection among the nine sets is made by comparing calculated χ(T), M(H), and Cmag(T) with the experimental bulk data in Sec. IIID. The later statement that the model 'reproduces' these data (Figs. 5–7, abstract) is therefore a discrete selection artifact, not an independent validation. The authors should report all nine parameter sets, or at least the spread of their predicted bulk properties, and show explicitly whether the bulk data discriminate among them. Without this, the uniqueness of the wavefunctions, the g-tensor components (Eq. 7, Table II), and the claim that the fitted Hamiltonian is the correct single-ion model are not established. This concern does not affect the directly observed level energ","section":"Sec. IIID"},{"comment":"No uncertainties are reported for the CEF parameters, and the derived g-values are quoted with error bars (gx = 2.66(1), gy = 1.99(1), gz = 1.32(1)) without any stated propagation. Given that nine parameter sets fit the INS data equally well, the true uncertainty in the g-tensor components is likely much larger than the quoted last-digit errors. The authors should provide parameter uncertainties from the fitting procedure or from the distribution over the nine acceptable solutions, and avoid over-precise digits unless justified. This is needed to assess the significance of the reported anisotropy.","section":"Table I and Sec. IIID"},{"comment":"The magnetic heat capacity Cmag is obtained by subtracting a phonon background modeled with four Debye temperatures, with no nonmagnetic analog to constrain the fit. Since Cmag(T) is one of the bulk quantities used to select among the nine INS-equivalent CEF parameter sets, an error or bias in the phonon subtraction could change the selected solution and hence the derived wavefunctions and g-tensor. The authors should include a sensitivity analysis (e.g., varying the Debye temperatures within their fit uncertainties) or otherwise demonstrate that the selection among the nine sets is robust to the phonon-model choice.","section":"Sec. IIIB, Eq. (4)"},{"comment":"The claim that the CEF model 'reproduces' the experimental χ(T), M(H), and Cmag(T) is overstated for the reason given in the first major comment. The paper should distinguish clearly between (i) the directly measured CEF level energies and their temperature evolution, which are robust, and (ii) the specific parameter-set validation, which is a selection procedure using the same data. The wording in the abstract and Section IV should be tempered accordingly.","section":"Secs. IIID and IV"}],"minor_comments":[{"comment":"The title 'Steven operators' should read 'Stevens operators' (also in the main text where the same misspelling appears).","section":"Appendix A"},{"comment":"'Gd chopper frequency of 400 Hz' is presumably a typo for the Fermi/bandwidth chopper frequency used on MARI; please check and correct.","section":"Sec. II"},{"comment":"The symbol g in MCEF = NA g μB Σ ... is not defined. It should be the Landé g-factor gJ, and the powder averaging over orientations should be stated explicitly, since a single orientation expression is shown.","section":"Appendix B, Eq. (9)"},{"comment":"The denominator '(ZkBT)^2' is typeset ambiguously as 'ZkBT);' it should be '(Z k_B T)^2'. Also define the indexing convention for Em and En.","section":"Appendix B, Eq. (10)"},{"comment":"Equation (2) would benefit from a parenthetical note that μ_eff,0 and μ_eff,1 are powder-averaged effective moments, since the text otherwise mixes single-ion g-tensor language with the powder expression.","section":"Sec. IIIA"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is suitable in scope for a condensed-matter physics journal. The main issue is the underdetermined CEF selection; the central low-energy conclusion is sound but the single-ion model claim needs reworking. I would require the authors to document the nine CEF solutions and provide uncertainty estimates before publication, possibly as supplementary material."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first magnetic and CEF characterization of Nd3Li3W2O12, and the directly observed INS excitations make a convincing case for a J_eff=1/2 Kramers doublet ground state with a ~9 meV gap and no magnetic order down to 0.1 K. The one caveat worth remembering is that the CEF parameter set in Table I is not unique; the paper picks it from nine equally good INS fits using the same bulk data that are later presented as confirmation.\n\nThe new data are real: powder XRD shows a clean garnet phase; magnetization, heat capacity, and INS on MARI cover the right temperature and energy ranges. The four CEF excitations at 9.3, 21.9, 25.6, and 91.8 meV are directly visible, and the thermodynamic probes independently support a low-energy doublet: low-T Curie-Weiss, saturation magnetization ~1.1 mu_B/Nd, Schottky gap, and entropy release to R ln 2. The two-level CW fit gives Δ/k_B ~93 K, consistent with the INS gap. That is a coherent story. I also appreciate that the authors state plainly that nine parameter sets fit the INS equally well and that they used the bulk data to break the tie. That is honest, and it is the right instinct; the problem is that the same bulk data then cannot be used as independent validation.\n\nWhere it gets soft: the selection procedure in Sec. IIID is a discrete fit, not a confirmation. The other eight parameter sets are not reported, so we cannot judge whether the thermodynamic data actually discriminate among them or whether a different choice would change the wavefunctions and g-tensor components. The reported g-values (gx=2.66, gy=1.99, gz=1.32) and the ground-state wavefunction in Eq. 7 are therefore tied to one arbitrary choice among nine equally good INS fits. The paper gives no uncertainties on the CEF parameters; given the degeneracy, that omission is load-bearing. This does not threaten the energy levels themselves, which are direct observables, nor the J_eff=1/2 conclusion, which is robust across bulk probes. But it means the detailed single-ion model, Tables I and II, should be treated as plausible rather than definitive.\n\nWho this is for: anyone working on rare-earth frustrated magnets, hyperkagome lattices, or CEF analysis of powder INS data. A serious referee can handle it; the paper deserves revision, not rejection. The authors should be asked to report the other eight fits, quantify how much the bulk data disfavor them, and add uncertainties or at least a sensitivity statement. If they can show that only one set survives the thermodynamic filter within error, the paper becomes much stronger. Even as it stands, the basic physical picture is likely right, and the new compound is a worthwhile addition to the small family of rare-earth hyperkagome magnets with a well-characterized low-energy description. I would send it to review and ask for the missing analysis rather than desk-reject.","headline":"A genuinely new data set for a previously uncharacterized hyperkagome magnet, with a solid J_eff=1/2 assignment but a CEF parameter set that is less unique than the paper suggests.","tokens_in":22104,"tokens_out":2272,"would_cite":true,"duration_ms":24679,"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":"Inelastic neutron scattering maps the five crystal-field doublets of Nd3Li3W2O12 and shows the low-temperature ground state is an effective spin-1/2 Kramers doublet.","keywords":["hyperkagome lattice","crystal electric field","inelastic neutron scattering","Kramers doublet","effective spin-1/2","frustrated magnetism","garnet","Nd3Li3W2O12"],"falsifier":"Grow single crystals and measure the low-temperature magnetization along three orthogonal axes: the paper predicts gx=2.66, gy=1.99, and gz=1.32, so a mismatch in the observed anisotropy would rule out the selected crystal-field set. A neutron diffraction search below 0.1 K would also directly test whether magnetic long-range order actually appears at lower temperatures.","tokens_in":21195,"feed_emoji":"🧲","tokens_out":8115,"duration_ms":85651,"temperature":0.7,"pith_summary":"This paper studies Nd3Li3W2O12, a garnet whose Nd3+ ions form a three-dimensional frustrated hyperkagome network. It claims that the crystal electric field splits the Nd3+ J=9/2 multiplet into five Kramers doublets, with the lowest doublet at least 9.3 meV below the first excited one, so at low temperatures the material behaves as an effective spin-1/2 magnet with weak antiferromagnetic interactions. The evidence is a simultaneous fit to inelastic neutron scattering data at 6, 100, and 250 K, producing a crystal-field Hamiltonian that also reproduces the measured susceptibility, magnetization isotherms, and field-dependent heat capacity. If correct, the compound is a clean single-ion spin-1/2 hyperkagome system with no magnetic order down to 0.1 K, making it a candidate for exploring frustration-driven quantum ground states.","feed_headline":"Nd hyperkagome magnet's ground state is effective spin 1/2","feed_subtitle":"A 9.3 meV gap separates the ground doublet, so below 100 K this frustrated lattice behaves as coupled spin-1/2 moments.","key_machinery":"The central object is the crystal-electric-field Hamiltonian written in Stevens operators, with nine parameters ($B2^{0}$, $B2^{2}$, $B4^{0}$, $B4^{2}$, $B4^{4}$, $B6^{0}$, $B6^{2}$, $B6^{4}$, $B6^{6}$) appropriate for D2 point symmetry. It is fitted simultaneously to the energy and intensity of neutron scattering excitations at three temperatures. Because nine parameter sets fit the neutron data equally well, the measured susceptibility, magnetization, and heat capacity are used to select one set. Diagonalizing the chosen Hamiltonian gives the five doublet energies and their wavefunctions, which in turn generate the crystal-field contributions to the thermodynamic quantities.","core_discovery":"The paper establishes that in Nd3Li3W2O12 the Nd3+ ion sits in a D2-symmetric crystal-field environment that splits the J=9/2 multiplet into five Kramers doublets (pairs of degenerate states protected by time-reversal symmetry) at energies 0, 9.3, 21.9, 25.6, and 91.8 meV. The lowest doublet is well separated and behaves as an effective spin 1/2 with anisotropic g-values gx=2.66, gy=1.99, and gz=1.32. A nine-parameter crystal-field Hamiltonian fitted to the neutron spectra reproduces the measured magnetic susceptibility, isotherms, heat capacity, and magnetic entropy, including the field-induced Schottky anomalies. No magnetic long-range order is observed down to 0.1 K, and the authors inter","pith_inferences":["Single-crystal neutron or optical measurements could resolve the nine-fold degeneracy in the crystal-field fit, since powder data alone cannot distinguish the parameter sets; the paper's thermodynamic tie-break is reasonable but not unique.","If the anisotropic g-tensor is confirmed, Nd3Li3W2O12 becomes a direct test bed for anisotropic exchange and quantum tunneling effects on a three-dimensional frustrated lattice, not just another effective Heisenberg spin-1/2 magnet.","A lanthanum-based nonmagnetic analog would provide a cleaner phonon subtraction than the multi-Debye fit used here and could sharpen the heat-capacity comparison in the 10-100 K range.","Measurements below 0.1 K, such as muon-spin rotation or neutron diffraction, could distinguish a true quantum spin liquid from a state with a very small ordered moment or a very low ordering temperature, a distinction the current data leave open."],"forward_implications":["Below temperatures of order 100 K, the 9.3 meV gap to the first excited doublet means the magnetic properties are governed by a single Kramers doublet, so the hyperkagome lattice can be modeled as coupled effective spin-1/2 moments.","The same crystal-field scheme reproduces both the high-temperature Schottky anomaly and the field-dependent low-temperature heat capacity, so it can predict thermodynamic behavior at fields and temperatures not directly measured.","The absence of magnetic long-range order down to 0.1 K, despite antiferromagnetic correlations developing near 2 K, indicates strong frustration on the hyperkagome network and makes the compound a candidate for a disordered or exotic ground state.","The computed g-tensor is strongly anisotropic, so any future low-energy spin model must include substantial anisotropy rather than simple Heisenberg exchange."],"supporting_citations":[{"why":"Provides the neutron spectrometer on which the inelastic neutron scattering data were collected.","marker":"[38]"},{"why":"Supplies the data-reduction and analysis framework used to fit the INS spectra.","marker":"[39]"},{"why":"Supplies the crystal structure and lattice parameters of Nd3Li3W2O12, the compound under study.","marker":"[36]"},{"why":"Supports the low-Q magnetic versus high-Q phonon separation used to isolate the crystal-field excitations.","marker":"[50]"},{"why":"Introduces the Stevens operator formalism in which the crystal-field Hamiltonian is written.","marker":"[51]"},{"why":"Provides the standard Stevens operator equivalents used to construct the nine-term Hamiltonian.","marker":"[52]"},{"why":"An analogous Nd-based square-lattice study that supplies the analysis route for crystal-field levels, g-values, and the effective spin-1/2 picture.","marker":"[9]"},{"why":"Supplies the method for computing crystal-field contributions to susceptibility, magnetization, and heat capacity used to validate the selected parameter set.","marker":"[40]"}],"fun_headline_variants":["Nd hyperkagome magnet: ground state is spin-1/2 doublet","Frustrated Nd garnet: crystal field yields effective spin-1/2","No long-range order in hyperkagome Nd magnet down to 0.1 K","Anisotropic g-values reveal effective spin-1/2 in Nd hyperkagome","Crystal field splits Nd3+ into five Kramers doublets in hyperkagome"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The central claim depends on choosing one of nine crystal-field parameter sets that fit the neutron data equally well; if the way the heat capacity and magnetization are corrected for phonon and non-spin magnetic backgrounds is biased, a different set would be selected and the wavefunctions, g-values, and effective-spin-1/2 conclusion would change.","fun_headline_variants_meta":{"raw":{"variants":["Nd hyperkagome magnet: ground state is spin-1/2 doublet","Frustrated Nd garnet: crystal field yields effective spin-1/2","No long-range order in hyperkagome Nd magnet down to 0.1 K","Anisotropic g-values reveal effective spin-1/2 in Nd hyperkagome","Crystal field splits Nd3+ into five Kramers doublets in hyperkagome"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000399,"raw_usage":{"total_tokens":2014,"prompt_tokens":925,"completion_tokens":1089,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":975}},"tokens_in":669,"tokens_out":1089,"duration_ms":11155,"temperature":1.0,"reasoning_tokens":975,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T20:13:09.935337+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Grow single crystals and measure the low-temperature magnetization along three orthogonal axes: the paper predicts gx=2.66, gy=1.99, and gz=1.32, so a mismatch in the observed anisotropy would rule out the selected crystal-field set. A neutron diffraction search below 0.1 K would also directly test whether magnetic long-range order actually appears at lower temperatures.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the data-reduction and analysis framework used to fit the INS spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the crystal structure and lattice parameters of Nd3Li3W2O12, the compound under study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the low-Q magnetic versus high-Q phonon separation used to isolate the crystal-field excitations."},{"cited_title":"Boothroyd, Principles of Neutron Scattering from Condensed Matter (OUP Oxford, 2020)","cited_arxiv_id":null,"evidence_quote":"Introduces the Stevens operator formalism in which the crystal-field Hamiltonian is written."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard Stevens operator equivalents used to construct the nine-term Hamiltonian."},{"cited_title":"Sibille, N","cited_arxiv_id":null,"evidence_quote":"An analogous Nd-based square-lattice study that supplies the analysis route for crystal-field levels, g-values, and the effective spin-1/2 picture."},{"cited_title":"Arnold, J","cited_arxiv_id":null,"evidence_quote":"Supplies the method for computing crystal-field contributions to susceptibility, magnetization, and heat capacity used to validate the selected parameter set."}],"review_version":1}