{"id":"8746f2f5-407c-4e33-aa62-23363be077e5","arxiv_id":"1909.01224","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The strongly anisotropic g factors of Yb3+ in NaYbS2 and NaYbO2 are controlled mainly by the charge contrast between in-plane Yb3+ neighbors and inter-layer Na+ ions, with a more homogeneous cation environment yielding nearly isotropic g factors.","lead":"This paper uses quantum chemical calculations to show that the anisotropic magnetic response of Yb3+ ions in the layered triangular magnets NaYbS2 and NaYbO2 is governed mainly by the charge imbalance between cations inside the magnetic layer and cations between layers. It offers a design rule for tuning magnetic anisotropy in rare-earth layered compounds, candidates for exotic spin-liquid states.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MgYbX2 control changes both inter-layer and in-plane cation charges, so the claimed inter-layer charge-imbalance mechanism is not isolated; a single-axis control calculation is required.","rationale":"The paper's computational protocol is standard and the NaYbS2 g factors agree reasonably with ESR for gab, with a larger deviation for gc; this gives some independent support. However, the central mechanism is established through a single control calculation, MgYbX2, that changes both the inter-layer cation charge and the in-plane Yb nearest-neighbor charge simultaneously. That confound, plus the unexplained large error in NaYbO2 gc, means the claim that inter-layer charge imbalance specifically governs the anisotropy is not yet demonstrated. I do not think this merits rejection: the authors report a plausible, potentially general design rule and provide transparent data. It does merit a conditional verdict with the requested control calculation. The reader identified exactly this confound, so I agree with the weakest-assumption analysis and the conditional recommendation.","tokens_in":9657,"tokens_out":5246,"duration_ms":48159,"concrete_test":"Run embedded-cluster CASSCF/MRCI calculations for NaYbS2 (same embedding and basis as the paper) in two control geometries: (i) replace only the inter-layer Na+ cations with Mg2+, leaving the six in-plane Yb nearest neighbors as Yb3+, with charge neutrality maintained by a distant compensating background charge; (ii) replace only the six in-plane Yb3+ neighbors with Yb2+, leaving Na+ on the A site, again with distant compensation. Compute the sans-SOC CF splittings and the SO-MRCI g factors for both controls and compare with the NaYbS2 and MgYbS2 values in Tables I–V. If (i) alone shifts (gab, gc) from (3.19, 0.93) toward (2.73, 2.49) and (ii) leaves them nearly unchanged, the inter-layer attribution survives. If the isotropic result requires both substitutions or is mainly reproduced by (ii), the headline claim must be weakened to a combined cation-charge-asymmetry effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing comparison for the central claim is the hypothetical MgYbX2 model introduced in the 'Cation charge imbalance effects' section and reported in Tables III–V. In that model, the inter-layer Na+ ions are replaced by Mg2+ while the six in-plane Yb3+ nearest neighbors are simultaneously replaced by Yb2+, and an additional negative charge is placed nearby to preserve neutrality. The paper's conclusion that 'the crucial role' is played by inter-layer cation charge imbalance, and its abstract's rule 'less inter-layer positive charge is associated with stronger in-plane magnetic response', require that the effect be specifically attributable to the A-site charge. The reported calculation does not isolate that variable: the quasi-cubic f-level pattern and isotropic g factors in MgYbS2 could come from removing the in-plane 3+ versus 2+ contrast, from the inter-layer substitution, or from the combination. The NaYbO2 comparison also shows a large discrepancy in the out-of-plane component (MRCI gc = 0.87 versus ESR gc = 1.75), which is the component most sensitive to the inter-layer environment; this does not by itself falsify the mechanism, but it adds uncertainty to the quantitative attribution. Because the central claim is about which structural feature governs the anisotropy, the confounded control is the most load-bearing weakness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports embedded-cluster CASSCF and MRCI calculations of the 4f13 crystal-field levels, spin-orbit levels, and g tensors of Yb3+ in the layered triangular-lattice delafossites NaYbS2 and NaYbO2. The authors find that the trigonal distortion of the YbX6 octahedra alone does not explain the strongly anisotropic g factors and argue, through comparison with hypothetical MgYbX2 lattices in which all surrounding cation charges are made 2+, that the dominant mechanism is inter-layer cation charge imbalance. They conclude with a proposed general rule: less inter-layer positive charge leads to a stronger in-plane magnetic response. The calculations reproduce the low-energy inelastic-neutron-scattering and ESR excitation energies and the NaYbS2 in-plane g factor reasonably well.","tokens_in":9955,"tokens_out":3700,"duration_ms":40353,"significance":"If the central attribution holds, the paper offers a concrete structural tuning principle for 4f13 triangular-lattice magnets, a family currently of high interest as quantum spin-liquid candidates. The calculations are first-principles embedded-cluster quantum chemistry, with no fitted parameters for the g factors, and the agreement with experiment for NaYbS2 is a genuine strength. The main limitation is that the decisive control calculation changes the inter-layer and in-plane cation charges simultaneously, so the reported data do not uniquely isolate the inter-layer charge-imbalance mechanism. The proposed general rule is falsifiable and should motivate further calculations, which raises the significance if the confound can be resolved.","major_comments":[{"comment":"The central claim is that inter-layer cation charge imbalance governs the g-factor anisotropy, but the MgYbX2 model changes both the inter-layer Na+ to Mg2+ and the six in-plane Yb3+ nearest neighbors to Yb2+, and then adds one negative charge to preserve neutrality. The resulting near-cubic splittings and isotropic g factors in MgYbX2 could arise from removing the in-plane 3+/2+ contrast, from the A-site substitution, or from the combination. A single-axis control calculation is required: replace only the inter-layer Na+ by Mg2+ while keeping the in-plane Yb3+ neighbors, and separately replace only the in-plane Yb3+ neighbors by Yb2+ while keeping Na+ on the A site. Without such calculations, the abstract's statement that 'less inter-layer positive charge is associated with stronger in-plane magnetic response' is not directly supported by the reported data.","section":"Cation charge imbalance effects; Tables III-V"},{"comment":"The charge-compensation scheme introduces an unspecified free parameter. The authors write that 'we added one negative charge within the nearby crystalline surroundings' to maintain overall neutrality in the MgYbX2 lattices. The position and spatial extent of this added charge are not documented. Since the reference Yb3+ is itself a charge defect in a 2+/2-/2+ lattice, the compensating charge can produce an additional low-symmetry potential that affects the f-level splittings and g factors. The robustness of Tables III-V should be demonstrated by placing the compensating charge at several distinct, physically reasonable sites and showing that the quasi-cubic level pattern and g factors are unchanged.","section":"Cation charge imbalance effects; MgYbX2 model construction"},{"comment":"For NaYbO2, the MRCI out-of-plane g factor is gc = 0.87 whereas the ESR value is gc = 1.75, a factor-of-two discrepancy in the component most sensitive to the inter-layer environment, while gab is well reproduced (3.31 versus 3.28). The authors attribute the discrepancy to uncertainties in the excited-state energies, but the size of the error limits confidence in the quantitative attribution of the anisotropy mechanism and in the predicted general trend, since the MgYbX2 comparison uses the same computational scheme. The manuscript should discuss whether this discrepancy could affect the conclusions drawn from the hypothetical MgYbX2 calculations.","section":"Table V; NaYbO2 comparison"}],"minor_comments":[{"comment":"The text states that inelastic neutron scattering finds 'three intense peaks at 35, 58, and 83 eV' for NaYbO2; the unit should evidently be meV.","section":"Basic electronic structure; NaYbO2 INS paragraph"},{"comment":"The g-tensor formula is attributed to Gerloch and McMeeking but the citation [45] is to Bolvin; a direct citation to the original derivation would be more informative, though this is a stylistic point.","section":"Acknowledgments/Bibliography"},{"comment":"The text says the ESR low-energy excitation is about 27 meV while the calculated first excited doublet is at 39 meV; making this comparison explicit in the text (alongside the INS values) would help the reader assess the discrepancy.","section":"Basic electronic structure; NaYbO2 ESR comparison"},{"comment":"Figure 1 labels the six nearest-neighbor Yb ions as 'Yb2' and the caption explains they are treated as large-core pseudopotentials; a sentence in the main text clarifying that these sites do not carry active electrons would avoid possible confusion.","section":"Material model, computational scheme; Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The paper's interpretative claim is attractive and likely qualitatively correct, but the reported control calculation cannot discriminate the inter-layer mechanism because it changes multiple variables at once. The missing single-axis control and the sensitivity analysis for the compensating charge are straightforward within the same embedded-cluster setup, so I see no reason to reject. I do not see a novelty or scope problem for this journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful embedded-cluster MRCI study of the low-energy 4f states and g tensors in NaYbS2 and NaYbO2, with a proposed mechanism for the strong single-ion anisotropy. The calculations are standard and described thoroughly, and the benchmark against INS/ESR is reasonable: the f-level splittings and the in-plane g of NaYbS2 come out close to experiment. The new piece is the suggestion that the anisotropy is governed by the charge contrast between the 3+ in-plane Yb neighbors and the 1+ inter-layer Na ions, rather than by the local trigonal distortion alone. That is worth taking seriously.\n\nWhere the paper is soft: the key comparison is the hypothetical MgYbX2 lattice, where the inter-layer Na+ is replaced by Mg2+ and at the same time the six in-plane Yb3+ neighbors are replaced by Yb2+, plus an extra negative charge somewhere nearby for neutrality. The calculation therefore removes the inter-layer versus in-plane charge imbalance, but it also changes the in-plane cation charge and the overall charge scale. As a result, you cannot attribute the change in g factors specifically to the inter-layer species. The abstract's rule \"less inter-layer positive charge is associated with stronger in-plane magnetic response\" is not actually isolated by the data as reported. This is a fixable problem — a control that varies only the inter-layer charge (e.g., replacing Na+ with a 2+ or 3+ cation while keeping Yb3+ in-plane) would settle it. The NaYbO2 out-of-plane gc discrepancy (0.87 vs 1.75 ESR) is also left unexplained; the paper mentions it candidly, but it does add uncertainty to the quantitative side.\n\nOverall: the computational work is solid and the interpretation is plausible but not yet demonstrated to the level the abstract claims. The paper advances the subfield and deserves a serious referee. I would send it back for a control calculation before accepting the mechanism, or accept it conditional on that. For a reading group, it's a useful example of how ab initio crystal-field calculations can interrogate spin-liquid candidates, but the central claim should be discussed with the caveat.","headline":"Solid quantum-chemistry study of Yb delafossites, but the central claim about inter-layer charge imbalance rests on a control that changes both inter-layer and in-plane cation charges simultaneously.","tokens_in":10457,"tokens_out":3028,"would_cite":true,"duration_ms":28550,"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":"In NaYbS$_2$ and NaYbO$_2$, the anisotropy of the Yb$^{3+}$ g tensor is set by inter-layer cation charge imbalance rather than by the local ligand-cage distortion.","keywords":["g-tensor anisotropy","Yb 4f13","triangular-lattice magnet","delafossite","cation charge imbalance","crystal-field splitting","spin-orbit coupling","quantum chemical cluster calculations"],"falsifier":"A single-variable calculation that replaces only the inter-layer Na$^+$ by Mg$^{2+}$ (keeping the in-plane Yb$^{3+}$ charges and the YbX$_6$ geometry unchanged) and finds that the g tensor remains strongly anisotropic would disprove the paper's attribution; the same test could be done experimentally by measuring g factors in a delafossite with a different inter-layer cation.","tokens_in":9464,"feed_emoji":"🧲","tokens_out":10608,"duration_ms":90254,"temperature":0.7,"pith_summary":"The paper argues that in the layered triangular-lattice compounds NaYbS$_2$ and NaYbO$_2$, the strong single-ion magnetic anisotropy of Yb$^{3+}$ ($g_{ab}\\approx 3.2$ versus $g_c\\approx 0.9$) is controlled primarily by the inter-layer cation charge imbalance rather than by the sizable trigonal distortion of the YbX$_6$ octahedra. The central evidence comes from embedded-cluster quantum chemical calculations (CASSCF/MRCI with spin-orbit) on the real compounds and on hypothetical MgYbX$_2$ lattices in which all surrounding cation charges are made +2. In the latter, the Yb 4f crystal-field splitting approaches a quasi-cubic pattern and the g factors become nearly isotropic, even though the local geometry is unchanged. The authors conclude that, as a general trend in 4f$^{13}$ layered compounds, lower inter-layer positive charge should yield a stronger in-plane magnetic response. This matters because g-tensor anisotropy directly shapes the effective spin Hamiltonian relevant to spin-liquid candidates.","feed_headline":"Charge imbalance, not distortion, sets Yb magnet anisotropy","feed_subtitle":"Quantum chemistry shows inter-layer cation charge, not trigonal distortion, sets the Yb g-factor anisotropy.","key_machinery":"The central object is the single-ion 4f crystal-field level structure of Yb$^{3+}$, computed with embedded-cluster multireference quantum chemistry: CASSCF for the seven f orbitals, MRCI for dynamical correlation, and then spin-orbit coupling, with the g tensor obtained through the Gerloch-McMeeking formula. The decisive comparison is between the real NaYbX$_2$ systems and hypothetical MgYbX$_2$ lattices, in which both the inter-layer cations (Na$^+$ replaced by Mg$^{2+}$) and the in-plane Yb neighbors (Yb$^{3+}$ replaced by Yb$^{2+}$) are changed, with an added negative charge for neutrality. This construction makes the cation charge environment homogeneous and thereby isolates, in the authors' reading, the role of cation charge imbalance relative to the trigonal ligand-cage distortion.","core_discovery":"The central claim is that the highly anisotropic, noncubic g factors of Yb$^{3+}$ in the 4f$^{13}$ delafossites NaYbS$_2$ and NaYbO$_2$ are predominantly caused by the charge asymmetry experienced by each magnetic center: the six in-plane nearest-neighbor Yb ions carry 3+ charges while the inter-layer Na$^+$ cations carry 1+. Replacing that asymmetric environment by a homogeneous 2+ cation distribution (Mg$^{2+}$ between layers and Yb$^{2+}$ in plane, with an added negative charge) transforms the f-level spectrum into a quasi-cubic pattern and makes the g tensor nearly isotropic—for example, from $(g_{ab},g_c)=(3.19,0.93)$ in NaYbS$_2$ to $(2.73,2.49)$ in the hypothetical MgYbS$_2$—despite the unchanged trigonal compression of the ligand octahedra. On this basis the paper asserts a general design rule for 4f$^{13}$ layered magnets: less inter-layer positive charge correlates with a stronger in-plane magnetic response.","pith_inferences":["If the principle is general, chemical substitution of the inter-layer A cation (for example K$^+$, Ag$^+$, or a divalent ion) in AYbX$_2$ delafossites should be a direct experimental lever on the g-tensor anisotropy, testable by electron-spin resonance.","The same charge-imbalance mechanism may operate in other f-electron layered systems such as 5f materials or 4f compounds with different A-site charges, where it could be probed by single-variable calculations that change only the inter-layer charge.","Because g-tensor anisotropy feeds directly into the anisotropic exchange parameters extracted from fits to magnetization and ESR data, the proposed mechanism could affect the quantitative spin models deduced for Yb triangular-lattice spin-liquid candidates.","A clean single-variable calculation—replacing only the inter-layer Na$^+$ by Mg$^{2+}$ while leaving in-plane Yb$^{3+}$ charges fixed—would separate inter-layer and intra-layer contributions; the paper does not report such a decomposition."],"forward_implications":["Less inter-layer positive charge should generally enhance the in-plane g factor of 4f$^{13}$ layered magnets, a trend the paper claims holds across the family.","Combining ligand-cage distortion, longer-range structural anisotropy, and cation charge imbalance offers a practical range of tunability for the single-site $g_{ab}$ and $g_c$ values.","In compounds where trigonal compression and charge-imbalance effects nearly cancel, the f-level spectrum can become quasi-cubic despite a strongly distorted local cage; NaYbS$_2$ is presented as such a case.","The paper calls for further computations to determine how these same knobs affect intersite magnetic couplings, beyond the single-site g tensor."],"supporting_citations":[{"why":"Supplies experimental ESR g factors and INS excitations for NaYbS$_2$ that the calculations are benchmarked against, and identifies the compound as a spin-liquid candidate.","marker":"[11]"},{"why":"Provides experimental ESR g factors for NaYbO$_2$ used as the key comparison for the computed g-tensor anisotropy.","marker":"[12]"},{"why":"Gives INS excitation energies for NaYbO$_2$ against which the computed low-lying Kramers doublets are checked.","marker":"[13]"},{"why":"Establishes the computational scheme for g tensors and the earlier demonstration that cation charge imbalance affects g factors in layered transition-metal oxides, the precedent extended here.","marker":"[7]"},{"why":"Supplies the large-core pseudopotentials that include the f electrons for neighboring Yb ions, making the single-site embedded-cluster model possible.","marker":"[22]"},{"why":"Provides the crystallographic data for NaYbS$_2$ used to construct the embedded cluster.","marker":"[31]"},{"why":"Provides the crystallographic data for NaYbO$_2$ used to construct the embedded cluster.","marker":"[32]"},{"why":"Gives the internally contracted spin-orbit MRCI method used to obtain the spin-orbit-coupled f states.","marker":"[34]"},{"why":"Provides the Gerloch-McMeeking-type formula used to compute the g tensor from the spin-orbit wavefunctions.","marker":"[45]"}],"fun_headline_variants":["Charge mismatch, not distortion, sets Yb anisotropy","Inter-layer cation charge controls Yb g-factor shape","Less inter-layer positive charge means stronger in-plane Yb response","Charge imbalance, not trigonal distortion, drives Yb g anisotropy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central conclusion depends on the assumption that the hypothetical MgYbX$_2$ model—which simultaneously changes the inter-layer cations, the in-plane Yb charges, and adds a compensating negative charge—represents the removal of inter-layer charge imbalance specifically, so the observed g-tensor isotropization can be attributed to that imbalance rather than to the other simultaneous changes.","fun_headline_variants_meta":{"raw":{"variants":["Charge mismatch, not distortion, sets Yb anisotropy","Inter-layer cation charge controls Yb g-factor shape","Less inter-layer positive charge means stronger in-plane Yb response","Charge imbalance, not trigonal distortion, drives Yb g anisotropy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000288,"raw_usage":{"total_tokens":1716,"prompt_tokens":1000,"completion_tokens":716,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":648}},"tokens_in":616,"tokens_out":716,"duration_ms":7802,"temperature":1.0,"reasoning_tokens":648,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:24:01.441965+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single-variable calculation that replaces only the inter-layer Na$^+$ by Mg$^{2+}$ (keeping the in-plane Yb$^{3+}$ charges and the YbX$_6$ geometry unchanged) and finds that the g tensor remains strongly anisotropic would disprove the paper's attribution; the same test could be done experimentally by measuring g factors in a delafossite with a different inter-layer cation.","supporting_citations":[{"cited_title":"NaYbS 2: A planar spin-1/2 triangular-lattice magnet and putative spin liquid,","cited_arxiv_id":null,"evidence_quote":"Supplies experimental ESR g factors and INS excitations for NaYbS$_2$ that the calculations are benchmarked against, and identifies the compound as a spin-liquid candidate."},{"cited_title":"Field-induced in- stability of the quantum-spin-liquid ground state in the Jeff = 1/2 triangular-lattice compound NaYbO 2,","cited_arxiv_id":null,"evidence_quote":"Provides experimental ESR g factors for NaYbO$_2$ used as the key comparison for the computed g-tensor anisotropy."},{"cited_title":"Gapless spin-liquid state in the structurally disorder-free triangular antiferromagnet NaYbO$_2$","cited_arxiv_id":"1901.07810","evidence_quote":"Gives INS excitation energies for NaYbO$_2$ against which the computed low-lying Kramers doublets are checked."},{"cited_title":"Orbital reconstruction in nonpolar tetravalent transition-metal oxide layers,","cited_arxiv_id":null,"evidence_quote":"Establishes the computational scheme for g tensors and the earlier demonstration that cation charge imbalance affects g factors in layered transition-metal oxides, the precedent extended here."},{"cited_title":"Energy- adjusted pseudopotentials for the rare earth elements,","cited_arxiv_id":null,"evidence_quote":"Supplies the large-core pseudopotentials that include the f electrons for neighboring Yb ions, making the single-site embedded-cluster model possible."},{"cited_title":"Schleid and F","cited_arxiv_id":null,"evidence_quote":"Provides the crystallographic data for NaYbS$_2$ used to construct the embedded cluster."},{"cited_title":"Mag- netic properties of ternary sodium oxides NaLnO 2 (Ln = rare earths),","cited_arxiv_id":null,"evidence_quote":"Provides the crystallographic data for NaYbO$_2$ used to construct the embedded cluster."},{"cited_title":"Spin-orbit matrix elements for inter- nally contracted multireference configuration interaction wavefunctions,","cited_arxiv_id":null,"evidence_quote":"Gives the internally contracted spin-orbit MRCI method used to obtain the spin-orbit-coupled f states."},{"cited_title":"An alternative approach to theg-matrix: the- ory and applications","cited_arxiv_id":null,"evidence_quote":"Provides the Gerloch-McMeeking-type formula used to compute the g tensor from the spin-orbit wavefunctions."}],"review_version":1}