{"id":"b82f394c-50be-45c4-9ce7-1e1b733dc50a","arxiv_id":"2505.00971","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In ultraclean α-RuCl3, the low-temperature specific heat along the Ru-Ru bond direction gives a gapless coefficient matching the Kitaev Majorana prediction, supporting a quantum spin liquid.","lead":"Specific heat measurements on ultraclean α-RuCl3 crystals show that the field-induced disordered state hosts gapless excitations whose linear-temperature coefficient approaches the value predicted for Majorana Dirac cones in the Kitaev spin liquid. The result strengthens the case that this material realizes a Kitaev quantum spin liquid in the clean limit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative Majorana match is not yet established: the T→0 intercept of C/T^2 at H||b could still contain disorder-induced or small-gap contributions, and the clean-limit extrapolation in Fig. 4 is uncontrolled.","rationale":"The reader flagged the assumption that no other T-linear background contaminates α. My stress-test sharpens this in two ways: the clean-limit extrapolation itself is uncontrolled, and even a small residual gap or disorder-induced term can mimic the intercept. These concerns reinforce the CONDITIONAL verdict rather than overturn it; the paper's angle-resolved gap data and sample-quality trends are genuine supporting evidence, but the quantitative claim needs lower-temperature measurements and an explicit disorder-extrapolation analysis. I do not find a reason to reject the paper, only to require additional evidence before the 'quantitative' language is accepted.","tokens_in":10241,"tokens_out":13906,"duration_ms":166372,"concrete_test":"Measure C/T^2 on #S1 for H||b down to 0.2 K in a dilution refrigerator, and also in 0.5° field-angle steps around b. If the intercept is constant below 0.5 K and collapses sharply for angles off b, the Dirac-node interpretation is supported. If C/T^2 bends downward at low T or remains finite over a finite angular range, a small gap or background (e.g., nuclear Schottky) is present and α is not intrinsic. In parallel, refit the existing 0.7–2 K data with C/T^2 = α + βT + A/T^4 + B exp(−Δ/kBT) to test whether α remains robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that α = 0.3 mJ mol−1 K−3 in #S1 at 10 T (H||b) quantitatively matches the Kitaev Majorana Dirac prediction (0.15–0.61), showing the clean FIQD state is a KQSL. This is not yet load-bearing. First, α is extracted as an intercept of C/T^2 over 0.7–2 K; a residual Majorana gap of about 1 K, a nuclear Schottky tail, or impurity-induced in-gap states can produce a nearly constant intercept over that window. Second, the paper's own electron-irradiated #B2 sample shows disorder creates low-energy T-linear excitations for H||b, and the measured α decreases monotonically with increasing TN; part of α in #S1 may therefore still be disorder-induced. Third, the inset of Fig. 4 extrapolates α versus l_ph^-1 without reported uncertainties or a specified functional form, so the clean-limit value could plausibly be zero or far below 0.15. Fourth, the theoretical band is about a factor of 4 wide because J is taken from the literature range 5–10 meV, making the observed 0.3 consistent with a broad set of scenarios. Without lower-temperature data, a controlled background subtraction, or a quantified disorder extrapolation, the agreement is not quantitative evidence for Majorana dispersions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports high-resolution specific heat measurements on ultraclean α-RuCl3 crystals grown by a two-step sublimation method, focusing on the field-induced quantum disordered (FIQD) state. The authors observe sixfold in-plane anisotropy of C/T, a field-dependent Majorana-like gap for H || a with T_max^2 Δ_M ∝ H^3 scaling, and a residual gapless coefficient α = lim C/T^2 for H || b. In the cleanest sample #S1 they measure α ≈ 0.3 mJ mol−1 K−3 at 10 T and compare it with a theoretical estimate of 0.15–0.61 mJ mol−1 K−3 derived from the reported Kitaev coupling J ≈ 5–10 meV. They conclude that the finite gapless coefficient, its approach to the theoretical band as sample quality improves, and the preserved H || a gap constitute quantitative evidence that the FIQD phase of bulk α-RuCl3 is close to an ideal Kitaev quantum spin liquid with bulk Majorana Dirac dispersions.","tokens_in":10482,"tokens_out":5811,"duration_ms":62569,"significance":"If the central claim holds, this would be an important step: it would connect bulk thermodynamic signatures of Majorana fermions in ultraclean crystals to the half-integer quantized thermal Hall edge response reported by the same group, supporting the bulk-edge correspondence in the Kitaev spin liquid candidate. The strengths of the paper are the new ultraclean crystals (TN = 7.7 K), the systematic field-rotation study, the internally consistent H^3 scaling, and the explicit comparison across samples of different quality. However, the quantitative claim is currently weakened by the factor-of-four theoretical range, the lack of error bars on α, the uncontrolled clean-limit extrapolation, and the only partially documented background subtraction. These issues are load-bearing because the central conclusion is specifically the claim of quantitative agreement with the Majorana prediction.","major_comments":[{"comment":"The central quantity α is read as the zero-temperature limit of C/T^2 over the measured window 0.7–2 K. This identification assumes that no other contribution to the specific heat has the same T^2 form or a sufficiently flat C/T^2 in this window. A sub-Kelvin gap in the Majorana spectrum, a nuclear Schottky tail, or a disorder-broadened in-gap continuum could each mimic an apparent intercept, and the phonon background is not subtracted in the displayed data. Because the decomposition formula C(T,H)/T = β(H)T^2 + C_M/T + C_flux/T is only sketched in the main text, the cleanliness of the extracted intercept is not established. Data below 0.7 K or a quantitative background model would be needed to support the claim that the intercept is intrinsic.","section":"§3, Fig. 2(c)"},{"comment":"The authors show that α decreases monotonically with increasing TN, from about 2.5 mJ mol−1 K−3 in the electron-irradiated #B2 sample to 0.3 in #S1. This is exactly the trend expected for a disorder-induced low-energy contribution, so the claim that the #S1 value is intrinsic rests on the extrapolation to the clean limit. The inset plots α versus l_ph^{-1} with only four points, no reported uncertainties, and a dashed line described only as a guide to the eyes. With no specified functional form or confidence interval, the clean-limit value could plausibly be zero or fall below the theoretical band. This uncontrolled extrapolation is load-bearing for the central 'quantitative evidence' claim and needs to be put on a firmer footing.","section":"§4, Fig. 4"},{"comment":"The theoretical range α_theory = 0.15–0.61 mJ mol−1 K−3 is generated from the literature range of J ≈ 5–10 meV, a factor-of-four spread. A measured value of 0.3 lies inside this window, but the window is so wide that the comparison has low discriminative power: essentially any value between 0.15 and 0.61 would be called quantitative agreement. The 'quantitative' language is therefore overstated unless the authors either provide a single computed α(J) curve or justify a narrower range of J. As it stands, the comparison is consistent with the Majorana prediction but does not strongly exclude other gapless quasiparticle scenarios with a similar density of states.","section":"§4, gray band in Fig. 4"},{"comment":"The manuscript relies on the Supplemental Material for the definitions of βph, the phonon mean free path estimate, and the theoretical calculation of α, but this material is not available with the manuscript. Since these details are essential for reproducing Fig. 4 and the central comparison, the authors should include the relevant formulas in the main text or supply the SM during review. In addition, no uncertainties are reported for TN, α, or the specific heat data in Figs. 2–4, which makes it impossible to judge whether the difference between #S1 and #S2, or between #S1 and the lower edge of the theoretical band, is statistically significant.","section":"§2, Table I and reference [35]"}],"minor_comments":[{"comment":"Sample #S2 is described as slightly more disordered than #S1, yet it shows a larger C/T at TN (2.1 versus 1.5 J mol−1 K−2). The text should clarify whether the jump height or the presence of the shoulder-like anomaly is the decisive quality criterion, or the statement may appear internally inconsistent.","section":"Table I"},{"comment":"The inset in Fig. 4 would be more informative with a linear fit to the four α versus l_ph^{-1} points, including the extrapolated intercept and its confidence interval, rather than an unspecified dashed guide to the eyes.","section":"Fig. 4 inset"},{"comment":"The abstract states that the gapless excitations 'quantitatively match' theoretical predictions, while the main text says 'in quantitative agreement'. Given the factor-of-four theoretical range and the absence of error bars, the wording should be softened to 'consistent with' or the analysis should be strengthened to justify 'quantitative'.","section":"Abstract and summary"},{"comment":"The sixfold oscillation amplitude in Fig. 2(a) is small (of order 0.2–0.4 mJ mol−1 K−2); representative error bars on a few data points would help the reader assess the significance of the anisotropy and the residual intercept in Fig. 2(c).","section":"§2, Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript builds directly on the authors' previous work (refs. 27, 29, 33, 34), and while the ultraclean-crystal data and the sample-quality trend are new, the novelty relative to those papers is incremental. The main risk is overclaiming 'quantitative' evidence: the measured α is consistent with the Majorana prediction, but the wide theoretical band and the uncontrolled disorder extrapolation prevent a quantitative conclusion at this stage. The paper would be publishable if the analysis is strengthened with error propagation, a documented background decomposition, and a more careful statement of what the comparison can and cannot establish."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news is the sample-quality trend: α for H∥b decreases monotonically as TN rises and as the phonon mean free path grows, landing at 0.3 mJ/mol·K³ in the cleanest crystal, inside the theoretical Kitaev window of 0.15–0.61. That is a genuinely useful observation, and it is the paper's main contribution over earlier work, which had already reported the six-fold C/T anisotropy and the H³ scaling. The data on four samples, including the electron-irradiated one, make a coherent case that the gapless direction is intrinsic and not a disorder artifact. I also appreciate that the authors show the raw C/T² curves and let the reader see the intercept being extracted.\n\nWhere I agree with the stress-test: the 'quantitative agreement' claim is not yet load-bearing. The theoretical band is set by J ranging from 5 to 10 meV, a factor of four, so a single measured α sitting inside that band is at best consistency, not a precision test. The extraction of α as the T→0 intercept of C/T² over 0.7–2 K could in principle pick up a small-gap activated term or a nuclear Schottky tail, and the paper's own electron-irradiated data show that disorder can create low-energy T-linear contributions for H∥b. The inset of Fig. 4 is a guide to the eye; without error bars or a specified extrapolation form, the clean-limit value could plausibly be anywhere from near zero to 0.5. The supplement apparently contains the decomposition details and theoretical derivation, but it was not available for checking — and for a claim built on 'quantitative evidence,' that is a real omission.\n\nThat said, I do not think the central qualitative conclusion collapses. The anisotropy, the H³ scaling, and the clear decrease of α with sample quality all point toward an intrinsic gapless direction in the FIQD state. The disorder story is actually self-consistent: dirtier samples show more low-energy weight, and the cleanest sample still shows a finite α that matches the approximate Kitaev prediction. The citation pattern looks normal, and the self-citations are to the prior specific-heat and crystal-growth work where credit is due.\n\nThis paper deserves a serious referee. I would send it out, then ask the authors to (1) release the supplement with the fitting details and theory estimate, (2) report uncertainties on α and Δ_M, and (3) soften 'quantitative' to 'consistent with' given the wide J range. The trend and the clean-limit data will still stand after those changes.","headline":"A careful specific-heat study on ultraclean α-RuCl3 showing that the residual T-linear term for H∥b shrinks toward the Kitaev prediction as sample quality improves; the trend is real, but the 'quantitative agreement' is weaker than claimed because the theory band is factor-of-four wide and the clean-limit extrapolation is uncontrolled.","tokens_in":11103,"tokens_out":2375,"would_cite":true,"duration_ms":28663,"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 the cleanest α-RuCl3 crystals, the gapless T-linear specific heat along the Ru-Ru bond direction quantitatively matches the Kitaev model's prediction for Majorana Dirac cones.","keywords":["α-RuCl3","Kitaev quantum spin liquid","Majorana fermions","specific heat","field-induced quantum disordered state","Dirac cones","field-angle dependence","ultraclean crystals"],"falsifier":"Measure $C/T^2$ in the cleanest sample at $H \\parallel b$ down to $0.3\\,\\mathrm{K}$ in a dilution refrigerator. If $\\alpha$ changes as the temperature range is extended, or if the field-angle pattern of the residual term deviates from the $|\\cos 3\\phi|$-driven Dirac behavior, the assignment of the $T$-linear term to Majorana Dirac fermions would be falsified. Alternatively, a crystal with an even longer phonon mean free path that gives $\\alpha$ outside the $0.15$–$0.61\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ window would break the quantitative agreement.","tokens_in":10000,"feed_emoji":"🧲","tokens_out":15206,"duration_ms":119398,"temperature":0.7,"pith_summary":"This paper tries to establish that the low-energy bulk excitations of the Kitaev spin-liquid candidate $\\alpha$-RuCl$_3$ are genuinely the Majorana quasiparticles of the Kitaev honeycomb model, not artefacts of sample disorder. Using ultraclean single crystals and specific-heat measurements under in-plane magnetic-field rotation, the authors find a sixfold-oscillating excitation gap, a fully gapped response when the field points along the $a$ axis, and gapless, linearly dispersing excitations when it points along the Ru-Ru bond direction. The coefficient of the $T$-linear term in the bond direction is $0.3\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ in the cleanest crystal, inside the $0.15$–$0.61\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ range predicted for Majorana Dirac cones with the reported Kitaev coupling $J \\approx 5$–$10\\,\\mathrm{meV}$. If the interpretation is right, the field-induced quantum disordered phase of $\\alpha$-RuCl$_3$ is close to an ideal Kitaev quantum spin liquid in the clean limit.","feed_headline":"Ultraclean α-RuCl3 bulk excitations match Majorana prediction","feed_subtitle":"The gapless bond-direction term matches the Kitaev spin-liquid prediction.","key_machinery":"The object doing the work is the field-angle-dependent Majorana gap of the Kitaev honeycomb model. In a magnetic field, the itinerant Majorana fermions acquire a gap $\\Delta_M$ that oscillates as $|\\cos 3\\phi|$ with the in-plane field angle $\\phi$, vanishing when the field lies along the Ru-Ru bond direction ($H \\parallel b$) and leaving a two-dimensional linear Dirac dispersion $E = v|k|$. The paper measures this gap through the specific heat: a gapped Majorana contribution plus a $Z_2$-flux term plus a phonon background for $H \\parallel a$, and a gapless Dirac coefficient $\\alpha = \\lim_{T\\to 0} C/T^2$ for $H \\parallel b$. The quantitative claim is carried by comparing $\\alpha$ with the range computed from the Dirac-cone density of states using the reported Kitaev coupling $J$ in the $5$–$10\\,\\mathrm{meV}$ window, and by the sample series in which $\\alpha$ tracks the inverse phonon mean free path toward the clean limit.","core_discovery":"On the paper's own terms, the central discovery is that the bulk specific heat of ultraclean $\\alpha$-RuCl$_3$ in the field-induced quantum disordered state reproduces, quantitatively, the two characteristic features of the Kitaev model's Majorana sector. For $H \\parallel a$, where the Majorana gap is maximal, a clear gap $\\Delta_M \\approx 18.8\\,\\mathrm{K}$ at $10\\,\\mathrm{T}$ is observed, and it follows the predicted $H^3$ scaling when combined with the vison-peak temperature. For $H \\parallel b$, the specific heat takes the gapless form $C/T = \\alpha T$, and $\\alpha$ extrapolates to about $0.3\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ in the cleanest sample, inside the theoretical window $0.15$–$0.61\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ obtained from the reported Kitaev coupling $J = 5$–$10\\,\\mathrm{meV}$. Across samples with different disorder levels, the gap is essentially unchanged while $\\alpha$ decreases as crystal quality improves, which the paper reads as evidence that the anisotropic Majorana excitations are intrinsic and that the field-induced quantum disordered state approaches the ideal Kitaev quantum spin liquid as disorder is removed.","pith_inferences":["A testable extension would be to push specific-heat measurements below $0.5\\,\\mathrm{K}$ in the same crystals: if the Dirac-cone assignment is right, the coefficient $\\alpha$ should stay constant and the field-angle oscillations of $C/T$ should sharpen as $T$ drops, whereas a nuclear Schottky or other background term would bend $C/T^2$ upward.","Because $\\alpha$ is set by the Dirac velocity, a precise clean-limit value could become an independent thermodynamic estimate of the Kitaev coupling $J$, cross-checking neutron-scattering and theoretical ranges.","The paper leaves open how dilute vacancies affect the gapless $H \\parallel b$ Dirac cone; a theory for that field direction, which lies in a different symmetry class, would predict whether ultraclean growth is necessary or merely helpful.","The same field-rotation specific-heat analysis could be applied to other Kitaev candidates, where the sixfold anisotropy pattern would help discriminate against non-Kitaev ordered states."],"forward_implications":["If the interpretation is correct, the field-induced quantum disordered phase of $\\alpha$-RuCl$_3$ is a bulk realization of the Kitaev quantum spin liquid, not a disorder-stabilized state.","The observation that $\\alpha$ decreases toward the clean-limit value while the $H \\parallel a$ gap stays fixed separates intrinsic Majorana physics from defect-induced in-gap states.","Combined with the reported half-integer quantized thermal Hall plateau in similarly grown crystals, the bulk gapless Dirac excitations support the predicted bulk–edge correspondence of the topological Kitaev state.","Specific heat under in-plane field rotation becomes a bulk thermodynamic signature for identifying Kitaev spin liquids and for ranking candidate materials by disorder."],"supporting_citations":[{"why":"Supplies the exactly solvable Kitaev model whose Majorana fermion and vison excitations are the subject of the comparison.","marker":"[1]"},{"why":"One source of the reported ferromagnetic Kitaev coupling J in the 5–10 meV range used to predict α.","marker":"[5]"},{"why":"Calculates the effective Kitaev model parameters for α-RuCl3, anchoring the other end of the J window.","marker":"[6]"},{"why":"Reports the half-integer quantized thermal Hall plateau in similar ultraclean crystals, supporting the bulk–edge correspondence.","marker":"[21]"},{"why":"Introduced the field-rotation specific-heat method and the sixfold Majorana-gap anisotropy that the paper reproduces.","marker":"[27]"},{"why":"Previous observation of gapless C/T = αT for H ∥ b in a Bridgman crystal and its Majorana-fermion interpretation.","marker":"[29]"},{"why":"Describes the two-step growth method that produced the ultraclean samples #S1 and #S2.","marker":"[33]"},{"why":"Establishes the electron-irradiated disordered sample #B2, providing the high-disorder comparison point.","marker":"[34]"},{"why":"Contains the theoretical α estimate from the Majorana Dirac cones and the phonon mean-free-path analysis used in the clean-limit extrapolation.","marker":"[35]"}],"fun_headline_variants":["Majorana dispersions confirmed in ultraclean α-RuCl3","Quantitative match: α-RuCl3 excitations follow Kitaev prediction","Ultraclean crystals reveal Majorana gap and Dirac modes","α-RuCl3 bulk excitations: direct evidence for Majorana bands","Gapless Majorana modes in α-RuCl3 match Kitaev theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the residual $T$-linear specific heat seen for $H \\parallel b$ comes purely from the predicted two-dimensional Majorana Dirac cones, with no additional disorder-independent constant term from nuclear spins, surfaces, or other quasiparticles entering the $T \\to 0$ intercept.","fun_headline_variants_meta":{"raw":{"variants":["Majorana dispersions confirmed in ultraclean α-RuCl3","Quantitative match: α-RuCl3 excitations follow Kitaev prediction","Ultraclean crystals reveal Majorana gap and Dirac modes","α-RuCl3 bulk excitations: direct evidence for Majorana bands","Gapless Majorana modes in α-RuCl3 match Kitaev theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00092,"raw_usage":{"total_tokens":3992,"prompt_tokens":1037,"completion_tokens":2955,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":2859}},"tokens_in":653,"tokens_out":2955,"duration_ms":20978,"temperature":1.0,"reasoning_tokens":2859,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:30:37.463984+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $C/T^2$ in the cleanest sample at $H \\parallel b$ down to $0.3\\,\\mathrm{K}$ in a dilution refrigerator. If $\\alpha$ changes as the temperature range is extended, or if the field-angle pattern of the residual term deviates from the $|\\cos 3\\phi|$-driven Dirac behavior, the assignment of the $T$-linear term to Majorana Dirac fermions would be falsified. Alternatively, a crystal with an even longer phonon mean free path that gives $\\alpha$ outside the $0.15$–$0.61\\,\\mathrm{mJ\\,mol^{-1}\\,K^{-3}}$ window would break the quantitative agreement.","supporting_citations":[{"cited_title":"Suzuki, H","cited_arxiv_id":null,"evidence_quote":"One source of the reported ferromagnetic Kitaev coupling J in the 5–10 meV range used to predict α."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Calculates the effective Kitaev model parameters for α-RuCl3, anchoring the other end of the J window."},{"cited_title":"Tanaka, Y","cited_arxiv_id":null,"evidence_quote":"Introduced the field-rotation specific-heat method and the sixfold Majorana-gap anisotropy that the paper reproduces."},{"cited_title":"Imamura, S","cited_arxiv_id":null,"evidence_quote":"Previous observation of gapless C/T = αT for H ∥ b in a Bridgman crystal and its Majorana-fermion interpretation."},{"cited_title":"Namba, K","cited_arxiv_id":null,"evidence_quote":"Describes the two-step growth method that produced the ultraclean samples #S1 and #S2."},{"cited_title":"Imamura, Y","cited_arxiv_id":null,"evidence_quote":"Establishes the electron-irradiated disordered sample #B2, providing the high-disorder comparison point."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the theoretical α estimate from the Majorana Dirac cones and the phonon mean-free-path analysis used in the clean-limit extrapolation."}],"review_version":1}