{"id":"497ad70f-255e-4581-9a92-a7d1e887b943","arxiv_id":"2502.00130","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Muon spin rotation and magnetization measurements on triangular-lattice YbZn2GaO5 find a dynamic non-magnetic ground state whose field-dependent spin fluctuations are consistent with a U1A01 Dirac quantum spin liquid.","lead":"This paper reports magnetic measurements on YbZn2GaO5, a material with triangular layers of ytterbium ions that may host a quantum spin liquid. The data show the spins remain dynamic down to 48 millikelvin, and the field dependence of the spin fluctuations matches a particular type of Dirac U(1) quantum spin liquid.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The U1A01 assignment rests on an unvalidated link: the fitted exponent n_D=0.08(3) in D2D is read as the Dirac-spinon signature, but no U1A01 calculation is shown to predict this value, and the J2D model is generic 2D diffusion.","rationale":"The reader's conditional verdict is appropriate, and my read does not change it. The empirical core of the paper, especially the absence of static order and persistent spin dynamics down to 48 mK, is well supported by the zero-field and longitudinal-field muon data, and the DFT-based muon-site analysis is a genuine independent check. The load-bearing concern is the specific classification step: the claim that the field-dependent relaxation is consistent with a U1A01 Dirac spin liquid depends on identifying the temperature exponent n_D = 0.08(3) of the fitted 2D spin-diffusion constant D2D with the spectral signature of a linear Dirac spinon dispersion. That identification is not derived anywhere in the manuscript. The model J2D(D2D, omega_e) is a generic 2D diffusive spectral density taken from the same group's earlier work on YbZnGaO4, and the paper does not show that this form is uniquely produced by the U1A01 state. A broader generic issue is that a nearly T-independent diffusion constant could arise in any gapless spin liquid or in a slowly fluctuating defect/spin-glass-like background. The fitted exponent is also fragile: the low-temperature 'quantum' region contains only four temperatures, and the total variation in D2D is comparable to the scatter of the global fit, so the 0.08(3) value should not be given the interpretive weight the paper places on it. The proposed direct calculation from the U1A01 parton theory would settle whether the observed muon relaxation is actually a prediction of the claimed state or instead a generic consistency argument. Since the paper carefully hedges the central claim as 'plausible description', a conditional acceptance remains the right outcome, with the requirement that this quantitative prediction be supplied before the title-level classification is taken as established.","tokens_in":47,"tokens_out":11106,"duration_ms":219424,"concrete_test":"Compute, from the U1A01 parton mean-field ansatz of Ref. [42] with J = 3.16 K and g_parallel = 3.436, the dynamic spin structure factor S(q, omega) at T = 0.05, 0.2, 0.4, and 1.6 K; then evaluate the longitudinal-field muon relaxation rate at the DFT site A, lambda(B) = gamma_mu^2 sum_q |A_d(q)|^2 S(q, omega_L), and compare its field dependence and temperature evolution with Fig. 3. If the predicted low-frequency diffusion constant follows D2D proportional to T^{0.08(3)}, the Dirac-QSL interpretation is quantitatively supported; if the predicted form or exponent differs, the J2D model is not specific to U1A01 and the classification remains unverified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive step in the Dirac-QSL classification is the interpretation of the low-temperature spin-diffusion parameter D2D (Fig. 4a): its nearly flat power law D2D proportional to T^{0.08(3)} is taken as the fingerprint of a linear spinon dispersion in the U1A01 state. This step is load-bearing but not demonstrated. The spectral density J2D in Eq. (2) is a generic 2D diffusive form imported from Ref. [26]; the paper does not derive J2D or the D2D(T) exponent from the U1A01 parton theory, nor does it show that competing states (U1A11, a spin glass, or a disorder-broadened paramagnet) predict a measurably different exponent. In addition, the 'quantum' regime T < J is covered by only four field scans (0.05, 0.2, 0.4, and 1.6 K), and the total change in D2D across this range is only about 30%; an exponent of 0.08(3) is therefore weakly constrained, and the systematic uncertainties of the 10-parameter global fit are not propagated into n_D. The empirical observation of persistent spin dynamics is solid, but the specific U1A01 assignment is a consistency argument, not a unique identification.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports zero-field and longitudinal-field muon spin rotation measurements on the triangular-lattice material YbZn2GaO5 down to 48 mK, together with pulsed-field magnetization at 0.5 K and DFT-based muon-site calculations. The zero-field data show exponential relaxation with no oscillations or missing asymmetry, establishing a dynamic non-magnetic ground state. Longitudinal-field relaxation data at eight temperatures are fit to a four-component model consisting of a 2D diffusive term, a 0D localized fluctuation term, a background, and a level-crossing resonance; the 2D diffusion rate D2D shows a crossover near J≈3.16 K and a weak power-law T^0.08(3) at low T, which the authors interpret as consistent with a gapless U(1) Dirac QSL of the U1A01 type. The paper also reports a 2.7 T level-crossing resonance, which it attributes to a muon-Yb-dimer singlet-triplet crossing, and derives entanglement length and Quantum Fisher Information from the fitted spectral density.","tokens_in":10515,"tokens_out":6108,"duration_ms":59244,"significance":"If the central identification were established, the paper would provide a showcase of a chemically ordered triangular-lattice material showing persistent spin dynamics down to milliKelvin temperatures, with a specific gapless U(1) Dirac QSL assignment and a constrained projective symmetry group. The paper's strongest asset is the raw muSR evidence: the base-temperature zero-field data show no static order, and the dynamic relaxation is present without invoking the global fit. The DFT muon-site assignment and the pulsed-field saturation moment are useful and appear carefully done. However, as the paper itself frames the QSL classification as 'consistent with', the specific U1A01 assignment is not uniquely determined: its decisive link, the exponent n_D of D2D(T) read from a generic 2D diffusive spectral density, is neither derived from the U1A01 parton theory nor compared with competing states. The stress-test concern that Eq. (2) is generic and that n_D is weakly constrained lands; the central claim is therefore plausible but currently under-supported.","major_comments":[{"comment":"The identification of YbZn2GaO5 as a U1A01 Dirac QSL rests on interpreting the nearly flat low-temperature power law D2D ∝ T^{0.08(3)} as the fingerprint of a linear spinon dispersion, but the manuscript never derives J2D or D2D(T) from the U1A01 parton theory. The spectral density in Eq. (2) is imported from Ref. [26] as a generic two-dimensional diffusive form, and no calculation is shown that a Dirac spinon spectrum in the U1A01 ansatz produces n_D = 0.08(3), nor that alternatives such as U1A11, a gapless Z2 state, or a disorder-broadened paramagnet lead to measurably different exponents. Without this, the consistency argument cannot distinguish the claimed state from other dynamic spin-liquid or glassy scenarios. Please provide the missing calculation or visibly soften the classification claim.","section":"Eq. (2), Fig. 4(a), Table II"},{"comment":"The low-temperature regime T<J is covered by only four field scans (0.05, 0.2, 0.4, and 1.6 K), and D2D changes by only about 30% over this range. With the global fit containing ten or more parameters (A, D, λBG, m, B0, Bwid and temperature-dependent D2D, ν, f), the reported n_D = 0.08(3) is weakly constrained, and the systematic uncertainties of the global fit are not propagated into n_D. Please add robustness tests, such as varying the 0D exponent m within its error, excluding individual temperatures, and using alternative forms for J2D, and report how n_D changes. If the exponent is not robust, the statement that the data are consistent with a Dirac dispersion should be withdrawn or replaced by an upper bound on the low-temperature variation of D2D.","section":"Fig. 4(a) and global fit (Eq. 1-4)"},{"comment":"The 2.7 T level-crossing resonance is explained post hoc by the singlet-triplet gap of a Yb pair using J=3.16 K and g∥=3.436 from the literature. This is an interesting observation, but it is not evidence for a QSL: it uses a pair picture whose relation to the proposed U(1) Dirac spin liquid is not explained. If the resonance is a property of the high-temperature paramagnetic regime, then its presence or absence below T∼J should not be used in the QSL classification without clarifying how it arises in the parton description. Please either provide a microscopic account of the level-crossing resonance within the proposed spin-liquid framework or explicitly remove it from the QSL evidence.","section":"End Matter, 'Level-crossing resonance', Fig. 7"},{"comment":"The entanglement length ξE and the Quantum Fisher Information FQ are computed from the fitted 2D spectral density J2D, not from an independent measurement. They therefore inherit all assumptions of Eq. (2) and of the global fit, and they cannot be used as independent confirmation of the QSL state. The text should state this model dependence explicitly and avoid presenting the large FQ as a direct experimental entanglement witness without the caveat that the witness is evaluated on the best-fit spectral function.","section":"Fig. 4(b,c) and End Matter, 'Quantum Fisher information'"}],"minor_comments":[{"comment":"The word 'low-freqency' should be 'low-frequency'.","section":"Final paragraph before Acknowledgments"},{"comment":"The shading that distinguishes the four spectral contributions is not legible in a monochrome version of the figure; please use distinct line styles or labels for the 2D, 0D, background, and LCR components.","section":"Fig. 3 caption"},{"comment":"The fitted exponent m=7(1) for the 0D spectral density is unusually large, and no physical motivation is given; a brief comment on what such a sharp cutoff implies for the local field distribution would help the reader.","section":"Eq. (3) and Table I"},{"comment":"The abstract says the relaxation is consistent with a single dominant muon site, whereas the DFT calculations in Table III find four sites; the text later argues that site A dominates, but a quantitative statement on the expected occupancy of sites B-D would make the assignment more transparent.","section":"Abstract and main text"}],"recommendation":"major_revision","confidential_remarks":"The paper's title and abstract make a stronger QSL claim than the evidence supports; the phrase 'consistent with' is appropriate, but the editors may wish to ensure that the revised version either supplies the missing theoretical link between the U1A01 parton theory and the measured D2D(T) exponent or visibly softens the classification. The work fits the journal's scope and would be of interest to the muSR and frustrated-magnetism communities."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my honest read of the YbZn2GaO5 paper. The empirical core is solid and new: this is the first μSR study of this clean triangular-lattice candidate, and the data convincingly show a dynamic, non-magnetic ground state down to 48 mK. The field-dependent relaxation is analyzed with a global multi-temperature fit that is careful, and the DFT+μ muon-site work is about as thorough as it gets. The contrast with YbZnGaO4 matters — this material avoids the site-disorder confound that has muddied earlier Yb triangle QSL claims.\n\nThe soft spot is the step from the fitted spectral density to a specific U1A01 Dirac QSL. The exponent D2D ∝ T^{0.08(3)} is presented as the fingerprint of a linear spinon dispersion, but that link is asserted, not derived. The J2D spectral-density model is imported from Pratt et al. and is a generic 2D diffusive form; the paper doesn't show that U1A01 predicts this exponent, or that plausible competitors (U1A11, a spin glass, a disorder-broadened paramagnet) would give a measurably different one. In the quantum regime below J there are only four temperatures, and D2D changes by only about 30%, so the exponent is weakly constrained. The entanglement length and QFI are derived from the fitted spectral density, so they inherit its assumptions; they're not independent witnesses.\n\nThat said, the authors deserve credit for hedging appropriately — they say \"consistent with\" and \"plausible description,\" not \"proves.\" The 2.7 T level-crossing feature is a nice addition: the crossing field predicted from literature J and g values matches, so it's a genuine consistency check, though its bearing on the ground state is indirect since it appears only at high T.\n\nOverall: the empirical finding and the analytical care merit publication. The specific PSG assignment is a consistency argument, and I wouldn't treat it as established. This paper deserves a serious referee, and the referee should push the authors to make the theoretical connection explicit — what does U1A01 actually predict for D2D(T), and what would the alternatives predict? I'd cite this for the μSR characterization of this material, and I'd take the dynamic-ground-state result as solid.","headline":"The muon data robustly establish a dynamic non-magnetic ground state in YbZn2GaO5; the U1A01 Dirac QSL assignment is a plausible but not uniquely forced interpretation.","tokens_in":11130,"tokens_out":4094,"would_cite":true,"duration_ms":40437,"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":"YbZn2GaO5's low-temperature spin dynamics match a gapless U(1) Dirac quantum spin liquid.","keywords":["quantum spin liquid","U(1) Dirac spin liquid","muon spin rotation","triangular lattice magnet","YbZn2GaO5","spin dynamics","entanglement witness","frustrated magnetism"],"falsifier":"Measure the zero-field muon asymmetry at temperatures below 48 mK on a high-quality crystal: resolved oscillations or a persistent 1/3 tail would show static magnetic order and falsify the dynamic spin-liquid picture. Independently, compute or measure whether the same $J_{2D}$ field dependence and weak power-law exponent are reproduced by a classical disordered triangular magnet; a match would remove the Dirac-spin-liquid specificity of the claim.","tokens_in":9909,"feed_emoji":"🧲","tokens_out":12217,"duration_ms":104929,"temperature":0.7,"pith_summary":"The paper sets out to establish that the triangular-lattice magnet YbZn$_2$GaO$_5$ has a quantum spin liquid ground state, and more specifically that its low-temperature spin dynamics match a gapless U(1) Dirac spin liquid of the U1A01 type. Zero-field $\\mu$SR measurements show no static magnetic order down to 48 mK, and longitudinal-field scans at eight temperatures are decomposed into a two-dimensional diffusive spectral density, a local fluctuation term, a background, and a level-crossing resonance. The extracted spin-diffusion rate follows $D_{2D}\\propto T^{0.08(3)}$ in the low-temperature regime, which the authors interpret as the dynamical fingerprint of linearly dispersing spinon excitations. A sympathetic reader would care because this compound, unlike its better-known siblings, is free of chemical site disorder, so the fractionalized-spin-liquid interpretation is not an artifact of random Mg/Ga occupancy.","feed_headline":"Muon data point to a Dirac spin liquid in YbZn2GaO5","feed_subtitle":"No magnetic order down to 48 mK; spin fluctuations match a gapless U(1) Dirac quantum spin liquid.","key_machinery":"The central object is the two-dimensional diffusive spectral density $J_{2D}(\\omega)$ entering Eq. (2), which relates the measured muon relaxation rate to a spin-diffusion constant $D_{2D}$. The longitudinal-field dependence is fit as the sum of this 2D term, a 0D localized fluctuation term, a constant background, and a Gaussian level-crossing resonance, with only $D_{2D}$, $\\nu$ and $f$ varying with temperature. The temperature dependence of $D_{2D}$ is the load-bearing quantity: its crossover at $T\\sim J$ separates classical from quantum regimes, and its weak power law $\\propto T^{0.08(3)}$ is read as the signature of the linear spinon dispersion of a U(1) Dirac spin liquid, classified as U1A01 in the projective symmetry group scheme, denoting a gapless state with zero flux per triangular plaquette and Dirac-cone spinon excitations.","core_discovery":"The central discovery is that the longitudinal-field muon relaxation rate $\\lambda(B_{\\rm LF})$ in YbZn$_2$GaO$_5$ is well described by a four-term model, and the dominant low-field term follows a two-dimensional diffusive spectral density $J_{2D}$ whose spin-diffusion constant falls steeply on cooling and then flattens into $D_{2D}\\propto T^{0.08(3)}$ below about $J$. This weak power law is the predicted dynamic signature of a U(1) Dirac spin liquid with a linear spinon dispersion, and the authors identify the state as U1A01 in the projective symmetry group classification. The same fits yield an entanglement length and a Quantum Fisher Information that grow strongly below $T\\sim J$, supporting a crossover from classical fluctuations to a quantum entangled regime. The magnetization saturating near 15 T with a corrected moment of $2.1(1)\\,\\mu_{\\mathrm{B}}$ and the absence of order down to 48 mK are consistent with this gapless, dynamic picture.","pith_inferences":["If the weak power-law scaling is genuine, nuclear magnetic resonance spin-lattice relaxation on the same crystals should show a correspondingly weak temperature dependence in the quantum regime, providing an independent check of the scenario.","Applying the same four-component decomposition to other clean triangular-lattice Yb materials could reveal whether $D_{2D}\\propto T^{0.08(3)}$ is a generic Dirac spin liquid signature or a compound-specific fit.","A decisive extension would be momentum-resolved inelastic neutron scattering under applied field, since the U1A01 theory predicts a gapless spinon continuum whose field evolution differs from that of a gapped or $\\pi$-flux state."],"forward_implications":["If the assignment is correct, YbZn$_2$GaO$_5$ becomes a clean, disorder-free platform for studying U(1) Dirac spinon physics on a triangular lattice, in contrast to the sister compounds where mixed-site disorder can mimic quantum spin liquid signals.","The scaling $D_{2D}\\propto T^{0.08(3)}$ gains status as a dynamical fingerprint for linear spinon dispersions, so future $\\mu$SR experiments on other candidate materials can test whether the same weak power law appears.","The growth of entanglement length and Quantum Fisher Information below $T\\sim J$ shows that longitudinal-field muon data can act as an entanglement witness in frustrated magnets.","The 2.7 T level-crossing resonance, present at high temperature and suppressed in the quantum regime, provides a pair-excitation probe whose disappearance marks the crossover into entangled spin-liquid behavior."],"supporting_citations":[{"why":"Supplies the two-dimensional diffusive spectral density $J_{2D}$ and the fitting strategy that converts muon relaxation into the spin-diffusion rate $D_{2D}$ used here.","marker":"[26]"},{"why":"Reports the synthesis, magnetic susceptibility, and the inelastic neutron scattering continuum of YbZn$_2$GaO$_5$ that identify it as a U(1) Dirac QSL candidate.","marker":"[27]"},{"why":"Provides the theoretical U(1) Dirac spin liquid description whose spinon spectrum and gapless excitations the present muon data are compared with.","marker":"[28]"},{"why":"Gives the projective symmetry group classification in which the U1A01 state is defined and distinguished from U1A11.","marker":"[42]"},{"why":"Supplies the muon spectroscopy formalism connecting longitudinal-field relaxation to the spectral density of spin fluctuations.","marker":"[30]"},{"why":"Introduces the Quantum Fisher Information expression used to derive an entanglement witness from the measured spectral density.","marker":"[32]"},{"why":"Defines the DFT+$\\mu$ method used to identify candidate muon stopping sites and support the assumption that only the lowest-energy site contributes.","marker":"[35]"},{"why":"Provides the plane-wave DFT code used to relax the muon-containing supercells and locate the stopping sites.","marker":"[38]"}],"fun_headline_variants":["Dirac spin liquid signature seen in muon spin relaxation","Muon study hints at gapless Dirac spin liquid in YbZn2GaO5","Weak power law reveals Dirac quantum spin liquid state","No order down to 48 mK: U(1) Dirac QSL dynamics in YbZn2GaO5","Spin fluctuations in YbZn2GaO5 match gapless Dirac QSL"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole conclusion rests on the assumption that the mathematical model used to fit the muon data, a two-dimensional diffusive spectral density carried over from the earlier YbZnGaO$_4$ study, is unique to a U(1) Dirac spin liquid; if the same model or the same weak power law $T^{0.08(3)}$ can arise from an ordinary or disordered magnet, the quantum spin liquid identification loses its force.","fun_headline_variants_meta":{"raw":{"variants":["Dirac spin liquid signature seen in muon spin relaxation","Muon study hints at gapless Dirac spin liquid in YbZn2GaO5","Weak power law reveals Dirac quantum spin liquid state","No order down to 48 mK: U(1) Dirac QSL dynamics in YbZn2GaO5","Spin fluctuations in YbZn2GaO5 match gapless Dirac QSL"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001193,"raw_usage":{"total_tokens":4909,"prompt_tokens":919,"completion_tokens":3990,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":3885}},"tokens_in":535,"tokens_out":3990,"duration_ms":28111,"temperature":1.0,"reasoning_tokens":3885,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T20:05:05.335407+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the zero-field muon asymmetry at temperatures below 48 mK on a high-quality crystal: resolved oscillations or a persistent 1/3 tail would show static magnetic order and falsify the dynamic spin-liquid picture. Independently, compute or measure whether the same $J_{2D}$ field dependence and weak power-law exponent are reproduced by a classical disordered triangular magnet; a match would remove the Dirac-spin-liquid specificity of the claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-dimensional diffusive spectral density $J_{2D}$ and the fitting strategy that converts muon relaxation into the spin-diffusion rate $D_{2D}$ used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the synthesis, magnetic susceptibility, and the inelastic neutron scattering continuum of YbZn$_2$GaO$_5$ that identify it as a U(1) Dirac QSL candidate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the theoretical U(1) Dirac spin liquid description whose spinon spectrum and gapless excitations the present muon data are compared with."},{"cited_title":"Li, Y.-M","cited_arxiv_id":null,"evidence_quote":"Gives the projective symmetry group classification in which the U1A01 state is defined and distinguished from U1A11."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the muon spectroscopy formalism connecting longitudinal-field relaxation to the spectral density of spin fluctuations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the DFT+$\\mu$ method used to identify candidate muon stopping sites and support the assumption that only the lowest-energy site contributes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the plane-wave DFT code used to relax the muon-containing supercells and locate the stopping sites."}],"review_version":1}