{"id":"94e8a2a4-838d-416e-a771-649b5f298280","arxiv_id":"2507.11391","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"NdZnAl11O19 shows a well-separated effective spin-1/2 doublet with moderate Ising anisotropy and persistent spin fluctuations down to 0.28 K, with no ordering to 50 mK, making it a candidate quantum spin liquid.","lead":"Researchers studied the magnetic behavior of NdZnAl11O19, a triangular-lattice magnet, and found signs of persistent spin motion down to very low temperatures, with no sign of magnetic freezing. The material is a new candidate for a quantum spin liquid, a state where electron spins remain correlated but never settle into an ordered pattern.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MuSR persistent-fluctuation plateau is the load-bearing evidence for QSL candidacy, but the analysis does not separate majority 2d-site Nd spins from the known ~20% 6h-site defect spins; if orphan 6h spins dominate, the QSL claim loses grounding.","rationale":"I read the manuscript as a candidate-nomination paper, not a proof of a QSL. Its strongest evidence is the combination of no ordering to 50 mK, a large CEF gap, and persistent muSR fluctuations. The CEF model has a real weakness (only two of four 2d transitions observed, third level set by a specific-heat decomposition), but that uncertainty affects the anisotropy parameters, not the existence of a Kramers doublet, and even a 3.5 meV gap would leave the system in the ground doublet at 0.28 K. The ac susceptibility frequency independence and ZF/FC agreement are independent support against spin glass. The decisive load-bearing assumption is therefore the muSR plateau: it is the only evidence that low-temperature dynamics are intrinsic and quantum, and it is exactly where known structural disorder (2d/6h Nd) can mimic the signal. The reader's verdict identified this same assumption. My concern is not that the authors are wrong, but that the paper does not yet provide a test separating 6h defect spins from 2d bulk spins; the planned single-crystal and neutron-continuum experiments are the natural route. Since the paper is explicitly a candidate nomination and the reader already conditions acceptance on this issue, the verdict should remain CONDITIONAL.","tokens_in":13502,"tokens_out":7095,"duration_ms":88399,"concrete_test":"Synthesize a series of polycrystalline Nd1-xLaxZnAl11O19 samples with the same 2d/6h site-occupancy distribution (verified by Rietveld refinement) and measure ZF- and LF-muSR down to 0.28 K. If the low-T lambda plateau persists or increases with the 6h-site occupancy and survives strong dilution of Nd on the 2d sites, the relaxation is defect-dominated; if it scales with the bulk 2d Nd concentration and disappears upon dilution, the intrinsic interpretation is confirmed. This directly tests whether the Eq. 4 single-component dynamic fit is an adequate description or a hidden static channel from 6h Nd is present.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To reach 'good QSL candidate,' the paper must show that the persistent low-temperature spin fluctuations are intrinsic to the triangular 2d Nd sublattice. The strongest evidence for this is the ZF-muSR relaxation plateau (lambda ~ 8.5 per microsecond below 15 K) plus the LF dependence at 0.28 K, fitted to the single-component Redfield form (Eq. 4) with nu = 85 MHz and Delta = 21.76 mT. The problem is that the same data would be produced if the relaxation is dominated by orphan or slowly fluctuating Nd moments at the 6h sites, whose occupancy is roughly 20% of the Nd content (Table S1: Nd2(6h), occ. 0.072; the paper assigns the broad 3.5/9 meV INS modes to these 6h Nd ions). The stretched exponent beta of about 0.6 already signals a distribution of relaxation rates, but Eq. 4 assumes one dynamic channel. A static or very slow 6h contribution would produce the observed plateau and LF decoupling, and the paper's counterarguments (resolution-limited CEF widths for 2d peaks, no ac spin-glass signature) do not exclude a separate magnetic population on 6h. Thus the central claim rests on an unseparated mixture: without confirming the 6h spins are inert, the plateau does not uniquely imply persistent quantum fluctuations on the triangular lattice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a combined ac-susceptibility, powder inelastic neutron scattering (INS), and zero/longitudinal-field muon spin relaxation (muSR) study of the triangular-lattice compound NdZnAl11O19. The authors extract a crystal-electric-field (CEF) scheme with a J_eff = 1/2 ground doublet dominated by |m_J = ±7/2>, with g_c = 4.54 and g_ab = 1.42, separated by about 9.5 meV from the first excited doublet. They observe no magnetic ordering or spin freezing down to 50 mK, a low-temperature Curie-Weiss temperature of -0.42 K, and a temperature-independent muSR relaxation rate lambda = 8.5 microsecond^-1 below about 15 K, which they interpret as persistent quantum spin fluctuations. The paper concludes that NdZnAl11O19 is a good candidate quantum spin liquid with moderate Ising anisotropy, analogous to recently proposed rare-earth triangular-lattice spin-liquid candidates.","tokens_in":13812,"tokens_out":6980,"duration_ms":81367,"significance":"If the persistent muSR relaxation is intrinsic to the triangular 2d Nd sublattice, this paper would add a useful new member to the small family of rare-earth triangular-lattice quantum spin liquid candidates, with a different anisotropy ratio from CeMgAl11O19 and NdTa7O19. The strengths of the paper are its combination of three complementary probes, the clear reporting of fit parameters with error bars, and the explicit acknowledgment of some ambiguities (e.g., only two of four CEF transitions observed). The central phenomenon, however, is not yet uniquely established because the muSR signal likely contains contributions from the minority Nd 6h sites, and the CEF model used to characterize the ground state is underdetermined. The paper is therefore of interest but requires additional analysis before the QSL candidacy claim is load-bearing.","major_comments":[{"comment":"The central claim that the persistent muSR relaxation is intrinsic to the triangular 2d Nd sublattice is not established, because the sample contains a substantial Nd population on the 6h site. With Nd1(2d) occupation 0.780 and Nd2(6h) occupation 0.072, the 6h sites carry roughly 22% of the total Nd content, and the manuscript itself assigns the broad 3.5 and 9 meV INS peaks to those 6h ions. The ZF plateau (lambda = 8.5 microsecond^-1) and the LF data at 0.28 K are fitted with a single-component Redfield form (Eq. 4), but a slow or static distribution of 6h spins would produce the same plateau, the same LF decoupling, and the stretched exponent beta = 0.6 observed below 15 K. The counterarguments given in Section III, namely resolution-limited CEF widths on the 2d site, absence of precession, and no ac spin-glass signature, are not site-selective and cannot exclude a separate magnetic population on 6h. The supplement itself concedes that a single-crystal study is required to better quantify the disorder. I ask for a two-component analysis of the ZF/LF spectra (for example, one dynamic Redfield component plus one static Kubo-Toyabe component) and an assessment of how the fast relaxation rate scales with the 2d-site fraction, ideally using an isostructural compound with a different 6h occupancy. Without this, the persistent-fluctuation plateau is ambiguous evidence for a quantum spin liquid on the triangular lattice.","section":"III, Fig. 5, Eq. (4); SM Table S1"},{"comment":"The CEF determination is underdetermined and partly circular. Only two of the four expected ground-state transitions are observed, and the supplementary specific-heat analysis (Eq. 6) assigns the 3.5 and 9 meV peaks to the 6h site and then fixes E1 = 9.5(1) meV for the 2d site; that value is then used as an input to the simultaneous INS/magnetization fit in the main text. Thus E1 is not an independent constraint, and the reported agreement of the calculated 100 K susceptibility is a consistency check of the same fitted model rather than an out-of-sample prediction. Because the quoted g factors (g_c = 4.54, g_ab = 1.42) are derived from the fitted CEF parameters, the 'moderate Ising anisotropy' characterization inherits the model ambiguity. Please report parameter uncertainties from the covariance matrix or bootstrap, and show that the conclusions are stable when the positions and intensities of the two unobserved doublets are varied over ranges consistent with the INS data.","section":"III, Fig. 3 and Table I; SM Eqs. (5)-(6)"}],"minor_comments":[{"comment":"The title contains 'tr iangular' and the abstract contains 'temprature'; the paper would benefit from a careful proofread.","section":"Title and Abstract"},{"comment":"Please specify the magnetic field and temperature used in the magnetization fit and the number of data points, and state clearly whether the 100 K data shown in Fig. 3(b),(c) are included in the fit or only used as a consistency check.","section":"III, Fig. 3(c)"},{"comment":"The claim that delta = 14.7(3) meV 'agrees reasonably well' with the first CEF excitation E1 = 9.5(1) meV is difficult to follow, since the values differ by a factor of about 1.5; please justify this statement or provide an alternative interpretation (e.g., a distribution of gaps or a different relaxation mechanism).","section":"III, Fig. 5 and text"},{"comment":"Reference [17] is given as 'see the supplementary materials'; it should be replaced with a proper citation or a DOI for the supplementary data.","section":"Reference [17]"},{"comment":"The text and figure caption should be reconciled: the weak features near 7.5 and 12.5 meV mentioned in Fig. 2(d) are not the same as the 3.5 and 9 meV excitations discussed in the text, and the reader should not have to infer which features are phonon remnants and which are assigned to the 6h site.","section":"Fig. 2(d) and Section III"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper is an incremental but reasonable extension of the CeMgAl11O19 template to a Nd analogue with different anisotropy. The main risk is that the 6h-site contribution to muSR is not separated, and the CEF fit is underdetermined; reviewers in the field will likely require the two-component muSR analysis and a more careful propagation of CEF uncertainties before accepting the QSL candidacy claim. No ethical concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. This is the first detailed experimental study of NdZnAl11O19, and it does the basics right: INS with phonon subtraction via a La analog, clear 23.8 and 70.7 meV CEF excitations with the right Q dependence, muSR down to 0.28 K showing no precession and a relaxation plateau, and ac susceptibility to 50 mK with no frequency dependence. The CEF analysis gives a J_eff = 1/2 doublet with moderate Ising anisotropy, and the authors are careful to call it a candidate rather than a verified QSL. That is the right level of claim.\n\nThe soft spots are real but not disqualifying. The CEF scheme is fit to only two observed transitions with two levels missing; the first excited state at 9.5 meV comes from a specific heat fit whose weighting already assumes the weak 3.5/9 meV modes belong to the 6h Nd. That's a bit circular, and the 100 K susceptibility is a consistency check of the same model, not an independent prediction. Still, the ground-state doublet and g factors are probably robust because the two observed peaks are well-resolved and the magnetization fits are decent.\n\nThe bigger issue is the muSR plateau. The paper interprets the temperature-independent lambda ~8.5 us^-1 below 15 K as persistent quantum fluctuations on the triangular 2d lattice. But the structure refinement puts ~20% of Nd on the 6h site, and the same weak broad modes assigned to those ions have first excitations at 3.5/9 meV. A population of orphan or slowly fluctuating 6h spins would produce just this kind of plateau and LF dependence. The stretched exponential beta ~0.6 already hints at a distribution of rates, yet Eq. 4 assumes a single dynamic component. The authors' counterarguments—2d peak widths are resolution-limited, no ac spin-glass signature—don't exclude a separate 6h magnetic population. So the central claim rests on an unseparated mixture.\n\nThat said, the paper is honest about its limitations and the candidate label is appropriately qualified. It deserves a serious referee, but the referee should push for single-crystal work, a direct search for the missing CEF levels, and a muSR analysis that explicitly models the 6h contribution (or experiments on a site-diluted/ordered variant). I would not cite it in my own work yet, but it's a useful addition to the family.","headline":"Solid first characterization of a new rare-earth QSL candidate, but the muSR plateau—the main evidence—does not cleanly separate the 2d triangular sublattice from ~20% 6h Nd disorder.","tokens_in":14421,"tokens_out":2367,"would_cite":false,"duration_ms":27135,"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":"NdZnAl11O19 is a candidate quantum spin liquid: its effective spin-1/2 moments keep fluctuating to 0.28 K, with no magnetic order or freezing down to 50 mK.","keywords":["quantum spin liquid","triangular lattice antiferromagnet","rare-earth hexaaluminate","muon spin relaxation","inelastic neutron scattering","crystal electric field","Ising anisotropy","NdZnAl11O19"],"falsifier":"A decisive test would be single-crystal neutron diffraction and inelastic scattering at millikelvin temperatures: finding long-range magnetic order, a spin-freezing transition, or a conventional gapped spin-wave spectrum would falsify the quantum spin liquid claim, as would a muon experiment showing that the low-temperature relaxation is static in origin rather than dynamically fluctuating.","tokens_in":2122,"feed_emoji":"🧲","tokens_out":2213,"duration_ms":96466,"temperature":0.7,"pith_summary":"The paper argues that the triangular-lattice antiferromagnet NdZnAl11O19 is a good candidate for realizing a quantum spin liquid state. It bases this on three combined observations: inelastic neutron scattering shows the Nd moments form a well-defined $J_{\\mathrm{eff}}=1/2$ ground doublet with moderate Ising anisotropy, muon spin relaxation shows persistent spin fluctuations down to at least 0.28 K, and ac susceptibility shows no magnetic ordering or spin freezing down to 50 mK despite an antiferromagnetic Curie-Weiss temperature of $-0.42$ K. If correct, this would add a member to the small family of rare-earth triangular-lattice quantum spin liquid candidates and sharpen the role of anisotropy in realizing such states.","feed_headline":"No order down to 50 mK: new quantum spin liquid candidate","feed_subtitle":"Muon and neutron data show persistent spin fluctuations to 0.28 K in NdZnAl11O19, a candidate quantum spin liquid.","key_machinery":"The central object is the crystal-electric-field Hamiltonian $H_{\\mathrm{CEF}} = \\sum B_l^m O_l^m$ acting on the $|J, m_J\\rangle$ basis of the Nd$^{3+}$ $^4I_{9/2}$ multiplet, with only $B_2^0$, $B_4^0$, $B_6^0$, and $B_6^6$ nonzero under the $\\bar{6}$ site symmetry. This Hamiltonian produces the effective $J_{\\mathrm{eff}}=1/2$ ground doublet and the measured $g$-factor anisotropy, and the same fit is constrained by the inelastic neutron scattering intensities, the powder magnetization, and the specific-heat excitation at 9.5 meV. The argument for persistent spin fluctuations is carried by the stretched-exponential muon relaxation analysis and by the modified Redfield formula $\\lambda(B_{\\mathrm{LF}}) = 2(\\gamma_\\mu \\Delta)^2 \\nu / [\\nu^2 + (\\gamma_\\mu B_{\\mathrm{LF}})^2] + \\lambda_0$, which separates the dynamic electron-spin contribution from static nuclear and disorder contributions.","core_discovery":"The paper establishes that NdZnAl11O19 hosts a magnetic ground state that resists order and never freezes, and interprets this as quantum spin liquid behavior. The inelastic neutron scattering data, fit with a Stevens-operator crystal-electric-field Hamiltonian, put the first CEF excitation at about 9.5 meV and yield a ground doublet dominated by $|m_J = \\pm 7/2\\rangle$ with a small $|m_J = \\mp 5/2\\rangle$ admixture, giving $g_c = 4.54$, $g_{\\mathrm{ab}} = 1.42$. Zero-field muon spin relaxation spectra are described by stretched exponentials whose rate $\\lambda$ saturates near $8.5\\,\\mu\\mathrm{s}^{-1}$ below about 15 K; longitudinal-field data at 0.28 K give a spin fluctuation rate $\\nu = 85.3$ MHz and an internal field width $\\Delta = 21.76$ mT, with $\\nu/(\\gamma_\\mu \\Delta) = 4.5$, consistent with fast fluctuations in the motional-narrowing regime. AC susceptibility shows no transition and no frequency dependence down to 50 mK, while the low-temperature Curie-Weiss temperature is $-0.42$ K, giving a frustration index above 8.4. The paper concludes that NdZnAl11O19 may host a quantum spin liquid state with dominant Ising anisotropy.","pith_inferences":["If the quantum spin liquid interpretation holds, the moderate Ising anisotropy of NdZnAl11O19 relative to CeMgAl11O19 may place it closer to an XY-like regime where quantum fluctuations are more effective; comparing the two systems could reveal whether the U(1) Dirac state survives over a range of anisotropy.","A testable extension is to dilute the Nd sublattice with nonmagnetic ions: if the low-temperature muon relaxation plateau is intrinsic to the spin liquid, it should be robust against moderate dilution, whereas a plateau driven by orphan spins or disorder should change markedly.","The mismatch between the muSR plateau onset near 15 K and the ac-susceptibility flattening near 0.1 K may reflect the difference between a local probe summing over all wavevectors and a bulk $Q=0$ probe; a wavevector-resolved neutron study below 1 K could directly test this.","Single-crystal growth would allow a direct search for the spin excitation continuum, the sharpest experimental fingerprint of quantum spin liquid behavior, and would also allow quantifying the role of the known 2d/6h site disorder that the powder study cannot fully resolve."],"forward_implications":["The material becomes a concrete testing ground for the $J_{\\mathrm{eff}}=1/2$ triangular-lattice quantum spin liquid scenario with Ising-type anisotropy, comparable to NdTa7O19 and CeMgAl11O19.","Because NdZnAl11O19 has moderate anisotropy ($g_c/g_{\\mathrm{ab}}\\approx 3.2$) while CeMgAl11O19 has strong Ising anisotropy ($g_c/g_{\\mathrm{ab}}\\approx 9.5$), the hexaaluminate family offers a natural comparison series for how anisotropy controls the stability of quantum spin liquid phases.","The large first CEF gap ($\\sim 110$ K) justifies modelling the low-temperature physics with an effective spin-1/2 Hamiltonian, so theoretical predictions for such a model become directly testable on this material.","Persistent muon fluctuations down to at least 0.28 K together with the absence of spin freezing down to 50 mK set an upper bound on any ordering or glassy transition energy scale, directing future searches for a spin excitation continuum.","The paper's suggestion that future single-crystal synthesis and neutron scattering could look for an excitation continuum identifies the next experimental step that would confirm or refute the quantum spin liquid interpretation."],"supporting_citations":[{"why":"Supplies the isostructural CeMgAl11O19 analogue and the proposed U(1) Dirac quantum spin liquid state that this paper extends and compares against.","marker":"[13]"},{"why":"Gives the NdTa7O19 quantum spin liquid candidate whose g-factor anisotropy, muon plateau onset, and Curie-Weiss temperature are the direct comparison points.","marker":"[9]"},{"why":"Provides the crystal-electric-field fitting software used to extract the CEF parameters, eigenvalues, and g factors from the INS and magnetization data.","marker":"[22]"},{"why":"Supplies the dipole-approximation neutron cross-section formula used to compute the powder-averaged INS intensities in the CEF fit.","marker":"[23]"},{"why":"Supplies the modified Redfield formula that the paper uses to extract the fluctuation rate and field width from the longitudinal-field muon data.","marker":"[34]"},{"why":"Provides the supplementary structural refinement, specific heat, 2d/6h site occupancies, and the 9.5 meV first CEF excitation used as constraints in the analysis.","marker":"[17]"},{"why":"Documents the analogous disorder and Ising anisotropy in PrMgAl11O19, supporting the treatment of site disorder in the hexaaluminate family.","marker":"[14]"},{"why":"Provides the Orbach spin-lattice relaxation process used to interpret the thermally activated rise of the muon relaxation rate toward the plateau.","marker":"[27]"}],"fun_headline_variants":["Quantum spin liquid candidate: no order to 50 mK","NdZnAl11O19: no magnetic order, spin fluctuations persist","Ising-anisotropic NdZnAl11O19 resists freezing to 50 mK","Persistent spin fluctuations in NdZnAl11O19: no ordering","Quantum spin liquid candidate: persistent fluctuations, zero freezing"],"cache_read_input_tokens":16384,"weakest_assumption_plain":"The load-bearing premise is that the persistent muon relaxation below 15 K comes from fast intrinsic electronic spin fluctuations rather than from the material's known Nd-site disorder, orphan spins, or nuclear moments; if disorder is the true source, the quantum spin liquid candidacy loses its footing.","fun_headline_variants_meta":{"raw":{"variants":["Quantum spin liquid candidate: no order to 50 mK","NdZnAl11O19: no magnetic order, spin fluctuations persist","Ising-anisotropic NdZnAl11O19 resists freezing to 50 mK","Persistent spin fluctuations in NdZnAl11O19: no ordering","Quantum spin liquid candidate: persistent fluctuations, zero freezing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000851,"raw_usage":{"total_tokens":3761,"prompt_tokens":1069,"completion_tokens":2692,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":2597}},"tokens_in":685,"tokens_out":2692,"duration_ms":22833,"temperature":1.0,"reasoning_tokens":2597,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:09:44.836576+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be single-crystal neutron diffraction and inelastic scattering at millikelvin temperatures: finding long-range magnetic order, a spin-freezing transition, or a conventional gapped spin-wave spectrum would falsify the quantum spin liquid claim, as would a muon experiment showing that the low-temperature relaxation is static in origin rather than dynamically fluctuating.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the isostructural CeMgAl11O19 analogue and the proposed U(1) Dirac quantum spin liquid state that this paper extends and compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the NdTa7O19 quantum spin liquid candidate whose g-factor anisotropy, muon plateau onset, and Curie-Weiss temperature are the direct comparison points."},{"cited_title":"Suter and B","cited_arxiv_id":null,"evidence_quote":"Provides the crystal-electric-field fitting software used to extract the CEF parameters, eigenvalues, and g factors from the INS and magnetization data."},{"cited_title":"Scheie, PyCrystalField: software for calculation, anal- ysis and fitting of crystal electric field Hamiltonians, J","cited_arxiv_id":null,"evidence_quote":"Supplies the dipole-approximation neutron cross-section formula used to compute the powder-averaged INS intensities in the CEF fit."},{"cited_title":"Yang, C.-Y","cited_arxiv_id":null,"evidence_quote":"Supplies the modified Redfield formula that the paper uses to extract the fluctuation rate and field width from the longitudinal-field muon data."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the supplementary structural refinement, specific heat, 2d/6h site occupancies, and the 9.5 meV first CEF excitation used as constraints in the analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the analogous disorder and Ising anisotropy in PrMgAl11O19, supporting the treatment of site disorder in the hexaaluminate family."},{"cited_title":"Scheie, M","cited_arxiv_id":null,"evidence_quote":"Provides the Orbach spin-lattice relaxation process used to interpret the thermally activated rise of the muon relaxation rate toward the plateau."}],"review_version":1}