{"id":"c4b7099e-8b6b-4758-a2ac-1c101940761d","arxiv_id":"2412.17913","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"BaCe2ZnS5 is a quantum paramagnet whose intra-dimer exchange is XY-type, giving an entangled ground state (|↑↑> - |↓↓>)/√2 instead of a Heisenberg singlet.","lead":"BaCe2ZnS5, a Shastry-Sutherland lattice magnet, shows no magnetic order down to 73 mK, and a dimer model fitted to neutron, magnetization, and heat capacity data yields an XY-type exchange with an entangled dimer ground state. The result adds a new anisotropic interaction class to this canonical frustrated lattice and predicts a field-induced quantum critical point near 11 tesla.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unquantified inter-dimer coupling J' in Eq. (1) could bias the fitted intra-dimer exchange and the inferred entangled ground state; the paper's own high-field heat capacity shows the dimer model fails well below the predicted QCP.","rationale":"The reader's weakest_assumption correctly identifies the neglect of inter-dimer coupling as the most load-bearing premise. The zero-field central claim is supported by multiple datasets, and the observed two-band structure with a degenerate lower doublet at 0.77 meV and an upper mode at 1.48 meV is a strong, direct signature of an XY-type anisotropic dimer: a Heisenberg dimer would give a triplet at a single energy, and an Ising dimer would not produce the observed intensity patterns. The quasi-flatness of the bands bounds the inter-dimer coupling to roughly the resolution scale, so the localized dimer model is likely a good starting point. However, the paper never converts the qualitative flatness observation into a quantitative bound on J', nor does it test whether the fitted Jxx, Jyy, Jzz and the resulting wavefunction are stable when a small J' is included. The high-field heat capacity data in Fig. 3(c) are particularly troubling because the low-temperature peak turns upward at 7.5–9 T, far below the dimer-model crossing at 11 T; if this is due to inter-dimer coupling, the deduced J' is not obviously negligible, and the zero-field parameters could be renormalized by an amount comparable to their quoted uncertainties. The manuscript is honest about the discrepancy, but it leaves the concern unresolved. Given that the reader's verdict is already CONDITIONAL on quantifying J' or demonstrating robustness, our read does not change the recommended verdict; the condition is appropriate and should be applied before the entangled-dimer ground state is accepted as definitive.","tokens_in":15886,"tokens_out":17241,"duration_ms":163799,"concrete_test":"Using the existing CNCS single-crystal data, extract the energy centroids of the E1+2 and E3 modes as a function of momentum along [h,0,0] and [0,k,0] and fit them to a Shastry-Sutherland dimer model with an inter-dimer coupling J' between nearest-neighbor dimers (two dimers per unit cell). Compare the best-fit J' to the quoted uncertainties and to the zero-field gap (0.77 meV). Then re-fit the full model including this J' and check whether the intra-dimer parameters Jxx, Jyy, Jzz shift by more than the reported 0.06 meV error bars. If the shift is within errors, the localized dimer result stands; if not, the entangled ground-state wavefunction is not uniquely determined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that BaCe2ZnS5 realizes XY-type dimers with ground state |ψ0>=(|↑↑>−|↓↓>)/√2 rests on fitting the data with a Hamiltonian (Eq. 1) that contains only intra-dimer exchange and Zeeman terms. The authors justify the omission of inter-dimer coupling J' by the quasi-flatness of the observed INS bands, but they never quantify J' or propagate its possible size into the fitted parameters. The flatness argument is qualitative: with an energy resolution of 0.1 meV FWHM, a bandwidth of up to ~0.05 meV would be unresolved, and such a J' is not negligible compared with the quoted uncertainty of 0.06 meV on Jxx, Jyy, Jzz. Moreover, the paper's own specific-heat data (Fig. 3(c)) show the low-temperature Schottky peak turning upward between 7.5 and 9 T for H ∥ [1,1,0], well below the dimer-model level crossing at ~11 T; the authors attribute this to inter-dimer couplings lifting the quasi-degeneracy. If J' is large enough to shift the effective level crossing by several tesla, it could also renormalize the zero-field fitted exchanges and admix the dimer eigenstates, changing the local ground state away from the pure |ψ0> product state. The manuscript provides no estimate of J' or a stability analysis of the fit under inclusion of J'.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a combined experimental and theoretical study of the Shastry-Sutherland lattice compound BaCe2ZnS5. Magnetization, specific heat, and inelastic neutron scattering (INS) data are analyzed with a localized spin-dimer Hamiltonian containing intra-dimer exchange and Zeeman terms only (Eq. 1). The authors fit three exchange parameters and three g-tensor components to five quasi-flat INS band energies, magnetization, and heat capacity, obtaining Jxx = 0.06(6) meV, Jyy = Jzz = -1.48(6) meV, with gxx = 1.2(1), gyy = 2.4(2), gzz = 2.1(1). From these parameters they deduce that the zero-field ground state of each dimer is the entangled state |ψ0⟩ = (|↑↑⟩ - |↓↓⟩)/√2, rather than the Heisenberg singlet, and they predict a field-induced level crossing near 11 T for fields along [1,1,0]. The SI provides an analytic diagonalization of the dimer Hamiltonian, together with details of the CEF analysis, polarized neutron diffraction, and the extraction of INS peak positions.","tokens_in":16145,"tokens_out":9021,"duration_ms":82302,"significance":"If the conclusions hold, the paper identifies a new type of Shastry-Sutherland magnet in which the intra-dimer exchange has strong XY-type anisotropy and the elementary magnetic units are entangled |↑↑⟩-|↓↓⟩ dimer states, in contrast to the singlet dimers of SrCu2(BO3)2. The central claim is well constrained: the fit is overdetermined (five energies plus magnetization and heat capacity for six parameters), the g-tensor is independently supported by CEF analysis and polarized neutron diffraction, and the model reproduces the INS spectral-weight modulations, not just the band energies. The prediction of a field-induced quantum critical point near 11 T is falsifiable. The main weakness is that the model neglects inter-dimer coupling J' without a quantitative bound, and the paper's own high-field heat capacity data deviate from the dimer model well below the predicted crossing, an effect attributed to inter-dimer couplings.","major_comments":[{"comment":"The localized dimer Hamiltonian omits inter-dimer coupling J' without a quantitative bound. The quasi-flatness of the INS bands is cited as evidence that J' is weak, but with an energy resolution of 0.1 meV FWHM, a bandwidth up to about 0.05 meV would be unresolved; this is comparable to the reported uncertainty of 0.06 meV on Jxx. Moreover, the authors' own specific-heat data (Fig. 3(c)) deviate from the dimer model between 7.5 and 9 T, an effect they attribute to inter-dimer couplings. To support the central claim that the intra-dimer exchange is XY-type and that the zero-field ground state is the product of local |ψ0⟩ states, the authors should estimate an upper bound on J' from the INS band widths and demonstrate that the fitted Jxx, Jyy, Jzz and the ground state are stable when a representative inter-dimer term is included.","section":"Eq. (1) and SI §IV"},{"comment":"The predicted level crossing at about 11 T is computed for the isolated dimer model, but the experimental specific-heat peak in Fig. 3(c) already begins to move upward between 7.5 and 9 T, indicating that inter-dimer couplings become relevant well below the bare crossing. The paper should clearly distinguish the bare dimer level crossing from the actual quantum critical field and either estimate J' from this discrepancy or explain why the deviation does not affect the zero-field determination of the exchange parameters.","section":"Fig. 3(a) and surrounding text"}],"minor_comments":[{"comment":"The third row of the table labels the 1.48 meV (0 T) and 1.75 meV (4 T) modes as 'E2'; these should be labeled E3 to be consistent with the text.","section":"SI Table S3"},{"comment":"In the conclusion, 'Fig. 3 (c)' should be 'Fig. 3 (d)' when referring to the schematic phase diagram.","section":"Conclusion"},{"comment":"The phrase 'in-plane Ising-spin nature with moments orthogonal to the dimer bond' is confusing because the effective spin model is described as XY-type; clarify that this refers to the anisotropic local susceptibility from polarized neutron diffraction.","section":"SI §III"},{"comment":"The CEF analysis reports the fitted parameters but does not give the corresponding g-tensor principal values; providing these would allow a direct comparison with the dimer-model g-factors in Eq. (3).","section":"SI §II"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong candidate for publication if the authors can provide a quantitative estimate of the inter-dimer coupling or show that the fitted exchange parameters and the inferred ground state are robust to a finite J'. This can be done with existing data, e.g., by fitting the upper bound on the band dispersion and re-fitting with a J' term in the SI. The request is for a stability analysis rather than new experiments, so it is within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Ma et al. report something genuinely new: the first resolved XY-type anisotropic exchange on a Shastry-Sutherland lattice. BaCe2ZnS5 is characterized with CEF analysis, polarized neutron diffraction, single-crystal INS, magnetization, and heat capacity. The localized dimer Hamiltonian (three exchange parameters, three g-factors) fits five quasi-flat INS bands and the thermodynamics. The fit is overdetermined, and the INS intensity modulations in the hk0 plane constrain transition matrix elements, so the resulting wavefunction (|↑↑>-|↓↓>)/√2 is not a pure restatement of the fitted energies. The circularity concern is real but partial.\n\nThe paper is honest about its own tension: the H∥[1,1,0] heat capacity deviates from the dimer model between 7.5 and 9 T, well below the predicted 11 T crossing, and the authors attribute this to inter-dimer couplings. That is the genuine soft spot. J' is never quantified, and the quasi-flatness argument is qualitative—a 0.05 meV bandwidth would be invisible under 0.1 meV resolution, and 0.05 meV is comparable to the 0.06 meV uncertainty on Jxx. A referee should demand a stability analysis: re-fit with a small inter-dimer coupling and show the intra-dimer parameters and the ground state are robust. Without that, the XY-dimer claim is plausible but conditional on J' being negligible.\n\nThat said, the stress-test overstates the risk to the zero-field ground state. The condition for |ψ0> to be the ground state is Jyy+Jzz<0, independent of Jxx. Both are -1.48(6) meV, so even a J' shift of 0.05 meV on Jxx doesn't change the ordering. The bigger risk—a sizable J' admixing the dimer product states—is indeed unaddressed, but the flat bands make it unlikely, not excluded.\n\nMinor: no raw data or code deposited, which is increasingly expected but not fatal.\n\nThis paper is for the Shastry-Sutherland and anisotropic-exchange community. It deserves a serious referee, who should ask for the J' bound and data deposition. The core result is new and probably correct; I'd send it out.","headline":"First resolved XY-type dimers in the Shastry-Sutherland family, with an honest but unquantified inter-dimer-coupling caveat that a referee should push on.","tokens_in":16809,"tokens_out":4801,"would_cite":true,"duration_ms":45778,"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":"BaCe2ZnS5 is a Shastry-Sutherland quantum paramagnet whose intra-dimer exchange is XY-type, so each dimer's zero-field ground state is the entangled state $(|\\uparrow\\uparrow\\rangle - |\\downarrow\\downarrow\\rangle)/\\sqrt{2}$, not a…","keywords":["Shastry-Sutherland lattice","XY exchange anisotropy","quantum paramagnet","spin dimer entanglement","inelastic neutron scattering","quantum critical point","BaCe2ZnS5","crystal electric field"],"falsifier":"Measure the supposedly flat 0.77 and 1.48 meV excitation bands across the full Brillouin zone with higher energy resolution: a resolved bandwidth comparable to or larger than the fitted exchange values would falsify the isolated-dimer assumption and require the exchange parameters and ground-state wavefunction to be revised.","tokens_in":15638,"feed_emoji":"🧲","tokens_out":12862,"duration_ms":106397,"temperature":0.7,"pith_summary":"The paper reports that BaCe2ZnS5 shows no magnetic order down to 73 mK and argues that this Shastry-Sutherland material is a quantum paramagnet built from entangled dimer pairs rather than the usual Heisenberg singlet dimers. Fitting a localized spin-dimer Hamiltonian to inelastic neutron scattering, magnetization, and heat capacity data, the authors extract an intra-dimer exchange tensor with strong XY anisotropy: $J_{xx}=0.06(6)$ meV and $J_{yy}=J_{zz}=-1.48(6)$ meV, together with an easy-plane $g$-tensor. In zero field, that exchange places each dimer's ground state at $|\\psi_0\\rangle = (|\\uparrow\\uparrow\\rangle - |\\downarrow\\downarrow\\rangle)/\\sqrt{2}$, an entangled state rather than the singlet $(|\\uparrow\\downarrow\\rangle-|\\downarrow\\uparrow\\rangle)/\\sqrt{2}$. If the dimer model is correct, this introduces a new anisotropic flavor into Shastry-Sutherland physics and predicts a field-induced quantum critical point near 11 T for fields along [1,1,0].","feed_headline":"BaCe2ZnS5 is a quantum paramagnet of entangled XY dimers","feed_subtitle":"Neutron, magnetization, and heat capacity all point to an 'up-up minus down-down' dimer ground state, not a singlet.","key_machinery":"The load-bearing object is the localized spin-dimer Hamiltonian $H = H_{\\text{intra}} + H_{\\text{Zeeman}}$ of Eq. (1), defined on effective spin-1/2 operators at the two Ce sites of each dimer, with dimer B obtained from dimer A by a 90-degree rotation. Symmetry reduces the exchange tensor to three diagonal components per dimer, and exact diagonalization in the four-dimensional Hilbert space provides analytic eigenstates whose lowest member is $|\\psi_0\\rangle$. The same wavefunctions enter the neutron cross-section through matrix elements and interference factors, so the model is checked against momentum-dependent intensity patterns of the flat bands, not only their energies. The quasi-flatness of the observed bands is the experimental evidence that inter-dimer coupling is weak enough to be omitted from this minimal model.","core_discovery":"The paper's central claim is that BaCe2ZnS5 realizes a Shastry-Sutherland lattice whose intra-dimer exchange is of XY type. The authors show that all low-energy magnetic degrees of freedom of the Ce3+ ground doublets can be described by a Hamiltonian containing only intra-dimer exchange and Zeeman terms, with symmetry-constrained diagonal exchange tensors for the two orthogonal dimer orientations A and B. Fitting the five observed quasi-flat neutron bands together with magnetization and heat capacity yields $J_{xx}=0.06(6)$ meV, $J_{yy}=J_{zz}=-1.48(6)$ meV and $g_{xx}=1.2(1)$, $g_{yy}=2.4(2)$, $g_{zz}=2.1(1)$. Because the $x$-component of the exchange is near zero while the $y$ and $z$ components are equal and negative, the zero-field wavefunction is the entangled state $|\\psi_0\\rangle = (|\\uparrow\\uparrow\\rangle - |\\downarrow\\downarrow\\rangle)/\\sqrt{2}$; the near degeneracy of the first two excitation modes at zero field is the direct signature of this XY anisotropy. With a field along [1,1,0], the model predicts a level crossing near 11 T, and with weak inter-dimer couplings that crossing becomes a quantum critical point with a dome of field-induced magnetic order and an Ising-like transition.","pith_inferences":["If the dimer model is correct, the heat-capacity deviation already visible between 7.5 and 9 T near the expected crossing offers a way to estimate the inter-dimer coupling $J'$; a quantitative fit of that region would likely place the true critical field below the isolated-dimer value of 11 T.","The XY anisotropy originates in the cerium crystal-field doublet, so substituting the rare-earth or ligand ions is a natural way to tune the intra-dimer exchange across Heisenberg, XY, and Ising limits on the same Shastry-Sutherland lattice.","The entangled $|\\psi_0\\rangle$ state has zero total $z$-magnetization but carries quadrupolar character, so measuring two-spin correlation functions or spin-nematic susceptibilities could reveal whether this dimer entanglement leaves observable signatures beyond single-dimer thermodynamics."],"forward_implications":["BaCe2ZnS5 becomes the first member of the BaR2ZnX5 Shastry-Sutherland family with a fully resolved anisotropic exchange Hamiltonian, and its zero-field state is a quantum paramagnet of entangled XY dimers rather than singlet dimers.","For a field along [1,1,0], the model predicts a level crossing near 11 T; weak inter-dimer couplings should convert that degeneracy into a quantum critical point with a dome of field-induced order in the field-temperature plane.","The local order parameter of that field-induced phase is the $y$-$z$ plane magnetization of the A-dimers, and the transition is expected to belong to the Ising universality class in $D=d+1$ dimensions.","Because the fit reproduces the momentum-dependent neutron intensities, the dimer wavefunctions are determined by spectral-weight distributions, giving a direct experimental handle on the entangled ground state."],"supporting_citations":[{"why":"It supplies the exact-product ground state of the Shastry-Sutherland model that defines the lattice being studied.","marker":"[12]"},{"why":"It provides the sister-compound study and the polarized-neutron method used to fix the orientation of the anisotropic g-tensor.","marker":"[16]"},{"why":"It gives the canonical Heisenberg Shastry-Sutherland realization with singlet dimer ground state, the benchmark the XY result is contrasted with.","marker":"[22]"},{"why":"It provides the crystal-field fitting routine used to extract the CEF parameters and g-tensor from neutron and susceptibility data.","marker":"[33]"},{"why":"It supplies the polarized-neutron diffraction approach for distinguishing Ising from XY anisotropy, applied here to orient the g-tensor.","marker":"[34]"},{"why":"It establishes the singlet-dimer ground state and low-lying excitations of the Heisenberg Shastry-Sutherland model that the entangled XY ground state is compared with.","marker":"[37, 38]"}],"fun_headline_variants":["BaCe2ZnS5: entangled XY dimers, not singlet state","Quantum paramagnet BaCe2ZnS5: entangled XY dimers","XY-type dimers entangle in Shastry-Sutherland BaCe2ZnS5","Not a singlet: BaCe2ZnS5's ground state is XY-entangled"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that Ce dimers are independent, with no inter-dimer coupling strong enough to matter, and if that coupling is not actually small, the fitted exchange values and the claimed entangled ground state could change.","fun_headline_variants_meta":{"raw":{"variants":["BaCe2ZnS5: entangled XY dimers, not singlet state","Quantum paramagnet BaCe2ZnS5: entangled XY dimers","XY-type dimers entangle in Shastry-Sutherland BaCe2ZnS5","Not a singlet: BaCe2ZnS5's ground state is XY-entangled"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00034,"raw_usage":{"total_tokens":1908,"prompt_tokens":1012,"completion_tokens":896,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":628,"completion_tokens_details":{"reasoning_tokens":807}},"tokens_in":628,"tokens_out":896,"duration_ms":7840,"temperature":1.0,"reasoning_tokens":807,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:09:04.264426+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the supposedly flat 0.77 and 1.48 meV excitation bands across the full Brillouin zone with higher energy resolution: a resolved bandwidth comparable to or larger than the fitted exchange values would falsify the isolated-dimer assumption and require the exchange parameters and ground-state wavefunction to be revised.","supporting_citations":[{"cited_title":"Structural fluctuations in the spin-liquid state of Tb2Ti2O7,","cited_arxiv_id":null,"evidence_quote":"It supplies the exact-product ground state of the Shastry-Sutherland model that defines the lattice being studied."},{"cited_title":"Topological triplon modes and bound states in a shastry– sutherland magnet,","cited_arxiv_id":null,"evidence_quote":"It provides the sister-compound study and the polarized-neutron method used to fix the orientation of the anisotropic g-tensor."},{"cited_title":"Distinct magnetic ground states in Shastry- Sutherland lattice materials: Pr2Be2GeO7 versus Nd2Be2GeO7,","cited_arxiv_id":null,"evidence_quote":"It gives the canonical Heisenberg Shastry-Sutherland realization with singlet dimer ground state, the benchmark the XY result is contrasted with."},{"cited_title":"Magnetic frustrations in the shas- try–sutherland system ErB4,","cited_arxiv_id":null,"evidence_quote":"It provides the crystal-field fitting routine used to extract the CEF parameters and g-tensor from neutron and susceptibility data."},{"cited_title":"Multi- step magnetization plateaus in the Shastry-Sutherland system TbB4,","cited_arxiv_id":null,"evidence_quote":"It supplies the polarized-neutron diffraction approach for distinguishing Ising from XY anisotropy, applied here to orient the g-tensor."}],"review_version":1}