{"id":"92637924-5baa-4d8a-a4e5-253973568621","arxiv_id":"1908.00803","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Thermal conductivity of the Kondo insulator CeRu4Sn6 is phonon-dominated and nearly isotropic from 80 mK to 80 K, with a low-temperature T^2 increase attributed to phonon-electron scattering.","lead":"Scientists measured heat flow in single crystals of the Kondo insulator CeRu4Sn6 from 80 millikelvin to 80 kelvin and found that heat moves almost identically along all crystal directions, carried mostly by lattice vibrations rather than electrons. The result is a useful constraint on this strongly correlated material, where electrical current is highly directional while heat flow is not.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stated Callaway form cannot produce the claimed T^2 law: with τ_pe^-1=Aω^2 the phonon-electron-dominated asymptote is κ∝T, and Table I makes boundary scattering negligible across the whole measured range.","rationale":"The reader's weakest_assumption already points at the Callaway phonon-electron scattering term, so I agree with that location. However, the sharper problem is that the model as written cannot actually produce the T^2 law it is used to explain. The derivation in the attack follows directly from inserting τ_pe^-1=Aω^2 into Eq. (1): the T^3 prefactor combines with the T^-2 from τ to give T, while a T^2 law requires τ_pe^-1∝ω. The tabulated A also makes B negligible at every experimental temperature, contradicting the claimed T^3 boundary regime and the extracted mean free path. This is not a disagreement with consensus about Pippard theory; it is an internal consistency check. If the exponent or A value is a typo, the correction is easy, but the fit parameters and inferred mean free paths would change, so the conclusion should remain conditional on that correction. The direct experimental observation of essentially isotropic, phonon-dominated κ is plausible and is not itself invalidated, but the load-bearing inference to isotropic carrier concentration depends on the corrected model. Since the reader's verdict was already CONDITIONAL, I mark the verdict unchanged while noting that the required condition is different and more specific than the one the reader identified.","tokens_in":5003,"tokens_out":12866,"duration_ms":133775,"concrete_test":"Evaluate Eq. (1) numerically at T=80 mK, 0.2, 0.5, 1, and 2 K using the Table I parameters for the c direction (A=2.2×10^7 s, B=2.8×10^6 s^-1, v_g=2316 m/s, θ_D=250 K), and compute the local log-log slope d ln κ_ph/d ln T over 0.2–2 K. If the slope is not 2, or is ~1 because Aω^2 fully dominates B, then the stated Callaway term cannot produce the claimed T^2 law. As a second part of the same check, refit the same data with the Pippard form τ_pe^-1=A'ω and verify whether a boundary-to-pe crossover with T^3 below ~0.2 K and T^2 up to 2 K emerges with a reasonable A'.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central mechanism identification rests on the Callaway fit with τ_pe^-1=Aω^2 (Eq. (2) and Table I). This is internally inconsistent with the claimed low-T power laws. In the phonon-electron-dominated regime, inserting τ≈1/(Aω^2) into Eq. (1) gives κ_ph ∝ T^3 ∫ x^4 / (A(k_BT/ħ)^2 x^2) e^x/(e^x-1)^2 dx ∝ T, not T^2. To obtain κ∝T^2 one needs τ_pe^-1∝ω, as in the usual Pippard treatment, not ω^2. Moreover, with the tabulated A≈2.2×10^7 s and B≈2.8×10^6 s^-1, Aω^2 exceeds B by more than 20 orders of magnitude at 80 mK, so the boundary-scattering term B can never dominate in the measured range; the claimed T^3 regime below ~0.2 K and the extracted boundary mean free paths of ~0.8 mm are not reproducible from the stated parameters. Thus the attribution of the T^2 range to phonon-electron scattering is not secured by the model as written. If the ω^2 form or the A value is a typographical error, the fit must be redone and the parameters re-reported before the conclusion about isotropic carrier concentration can be drawn.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports thermal conductivity measurements on single crystals of the Kondo semiconductor CeRu4Sn6 along the c axis and perpendicular to it, covering 80 mK to 80 K. The authors find that the thermal conductivity is essentially isotropic, that the electronic contribution estimated via the Wiedemann-Franz law is small, and that the phononic part dominates over the whole range. Below 2 K the data are reported to follow a T^2 law, which is attributed to phonon-electron scattering, with a T^3 boundary-scattering regime at the lowest temperatures. The data are fitted to a Callaway model including phonon-electron, boundary, defect, and Umklapp scattering. The conclusion is that the strongly anisotropic electrical resistivity is due to anisotropic charge-carrier mobility rather than anisotropic carrier concentration.","tokens_in":5384,"tokens_out":5106,"duration_ms":53127,"significance":"If the conclusions hold, the work provides a useful transport characterization of a non-centrosymmetric Kondo insulator and supports a scenario in which carrier concentration is isotropic while mobility is strongly anisotropic, with implications for the interpretation of resistivity, thermopower, and optical conductivity in CeRu4Sn6. The experimental strengths are the wide temperature range, the single-crystal measurements along two crystallographic directions, and the direct visibility of phonon dominance and approximate isotropy in the raw data. However, the specific mechanism identification—the T^2 law attributed to phonon-electron scattering—rests on a Callaway fit whose stated parameters are internally inconsistent, as detailed in the major comments. Because this attribution is the load-bearing link to the isotropic-carrier-concentration conclusion, the central interpretive claim is not yet secured.","major_comments":[{"comment":"The Callaway model as specified cannot yield the claimed low-temperature T^2 law. With τ_pe^-1 = A ω^2 substituted into Eq. (1), the phonon-electron-dominated conductivity is κ_ph ∝ T^3 ∫ x^4 e^x/(e^x-1)^2 dx / [A(k_B T/ħ)^2 x^2] ∝ T, not T^2. A T^2 asymptote requires τ_pe^-1 ∝ ω rather than ω^2. Therefore the attribution of the 0.2–2 K T^2 range to phonon-electron scattering is not supported by the model as written; the fit must be redone with the correct frequency dependence (or the table corrected) before this conclusion can be drawn.","section":"Eq. (2), Table I, and paragraph after Fig. 3"},{"comment":"The returned parameters are mutually inconsistent with the claimed boundary-scattering T^3 regime. Taking A ≈ 2.2×10^7 s and B ≈ 2.8×10^6 s^-1, at T = 80 mK, ω ≈ k_B T/ħ ≈ 10^10 s^-1, so Aω^2 ≈ 10^27 s^-1, which exceeds B by more than 20 orders of magnitude. Thus the boundary term never dominates in the measured range, and the T^3 behavior below ~0.2 K and the extracted mean free paths 0.77–0.83 mm cannot be reproduced from the stated parameters.","section":"Table I and paragraph after Fig. 3"},{"comment":"The identification of the T^2 range as phonon-electron scattering is not tested against alternatives with the same temperature dependence, such as scattering by structural disorder or two-level systems, or by magnetic excitations. Because the isotropic-carrier-concentration conclusion in the Conclusion depends directly on this identification, the authors should provide a quantitative comparison with at least one alternative scattering channel or otherwise justify its neglect.","section":"Physical properties, after Fig. 3"}],"minor_comments":[{"comment":"The figures do not show error bars; please state the measurement uncertainty and the number of samples measured.","section":"Figures 1–3"},{"comment":"The text refers to 'solid red and blue lines' in Fig. 3, but the caption and main text should be made consistent with the line styles used in Fig. 2; ensure that color and line labels match.","section":"Figure 3 caption and text"},{"comment":"Reference [19] is a PhD thesis; please cite the published Hall-effect data if available and give numerical values for the carrier concentration at high and low temperatures.","section":"Reference [19]"},{"comment":"Around 10 K the data show a hump in κ⊥c; the claim that κ is 'essentially isotropic' should be quantified with the maximum relative difference between κc and κ⊥c.","section":"Isotropy claim, around 10 K"},{"comment":"There are several typographical and formatting issues, such as 'cryst allographic' in the Experimental section and inconsistent spacing in Table I; please check against the journal style.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"No additional remarks beyond the report."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. It is the first thermal conductivity dataset for CeRu4Sn6, measured along both c and ⊥c from 80 mK to 80 K. The headline empirical result is that κ is essentially isotropic and phonon-dominated across the whole range, in sharp contrast to the anisotropic resistivity and thermopower. That part is directly visible in the data and is a genuinely useful contribution.\n\nThe soft spots start with presentation: no error bars, no raw data, and a small but real mismatch between Table I and the text about which boundary-scattering rate belongs to which direction. Those are fixable.\n\nThere is a more serious problem with the Callaway analysis. The paper attributes the 200 mK–2 K T^2 regime to phonon-electron scattering with τ_pe^-1 = Aω^2. That cannot produce κ ∝ T^2; inserting that τ into the Callaway integral gives κ ∝ T. The T^2 power law requires τ_pe^-1 ∝ ω, not ω^2. On top of that, the tabulated A and B values make phonon-electron scattering dominate boundary scattering by twenty orders of magnitude at 80 mK, so the claimed T^3 boundary regime and the ~0.8 mm mean free paths are not reproducible from the stated parameters. This is not a nitpick: the conclusion about isotropic carrier concentration rests on this mechanism attribution. If the ω^2 form or the A value is a typo, the fit needs to be redone and the parameters re-reported before that conclusion can be drawn.\n\nThe direct isotropy result probably survives: two independent measurements of κ_c and κ_⊥c are similar, and the electronic correction is small. But as written, the physical interpretation is not secured.\n\nBottom line: the data are worth refereeing. This is a plausible first measurement that a serious editor should send out, but the manuscript needs major revision—at minimum a corrected Callaway analysis, error estimates, and raw data release. If I were working on Kondo insulators, I would cite the isotropic phonon-dominated κ, but I would not trust the carrier-concentration inference until the fit issue is resolved.","headline":"New isotropic thermal conductivity data for CeRu4Sn6, but the Callaway analysis as reported cannot produce the claimed power laws—worth refereeing with major revisions.","tokens_in":5916,"tokens_out":3985,"would_cite":true,"duration_ms":35310,"reading_group":"maybe","serious_thinker":"unclear","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The thermal conductivity of CeRu4Sn6 is essentially isotropic and phonon-dominated, with a low-temperature T-squared law pointing to an isotropic carrier concentration.","keywords":["Kondo semiconductor","CeRu4Sn6","thermal conductivity","phonon-electron scattering","Callaway model","anisotropic resistivity","charge carrier mobility","Hall effect"],"falsifier":"Remeasure $\\kappa$ below 2 K on crystals with different smallest dimensions: if the $T^2$ regime is true phonon-electron scattering, its prefactor should stay fixed while the $T^3$ boundary-scattering tail moves to lower temperatures in larger crystals, whereas a geometry-dependent $T^2$ term would expose an additional scattering channel. Alternatively, tune the carrier concentration by doping or pressure and check whether the fitted phonon-electron prefactor $A$ tracks the Hall carrier density as that theory requires.","tokens_in":4823,"feed_emoji":"","tokens_out":15168,"duration_ms":128246,"temperature":0.7,"pith_summary":"This paper reports thermal conductivity measurements on single crystals of CeRu4Sn6, a Kondo semiconductor in which a narrow gap opens in the electronic density of states at low temperature and the electrical resistivity is strongly anisotropic. The central finding is that heat flow is essentially isotropic along and perpendicular to the c-axis between 80 mK and 80 K, and that phonons carry almost all of the heat in this entire range. Below 2 K the phonon part follows a $T^2$ law, which the authors attribute to phonon-electron scattering, with a $T^3$ boundary-scattering regime at the lowest temperatures. Because the phonon-electron scattering rate is expected to scale with the charge carrier concentration, this isotropic heat transport supports the earlier conjecture that the large anisotropy in electrical resistivity comes from anisotropic carrier mobility rather than an anisotropic carrier concentration.","feed_headline":"CeRu4Sn6 conducts heat almost identically in every direction","feed_subtitle":"Phonons, not electrons, carry the heat, so the crystal's strong electrical anisotropy hides an isotropic carrier density.","key_machinery":"The central object is the Callaway model for phonon thermal conductivity, in which the total phonon relaxation rate is the sum of independent scattering channels: boundary scattering with constant rate $B$, phonon-electron scattering with rate $A\\omega^2$, point-defect scattering with rate $C\\omega^4$, and Umklapp scattering with rate $D\\omega^2 T \\exp(-\\theta_D/(3T))$. The decisive channel is the phonon-electron term: a standard theory of phonon-electron scattering connects its prefactor $A$ to the charge carrier concentration $n$, so the near-equality of $A$ for heat flow along and across the c-axis is what turns an isotropic thermal conductivity into evidence for an isotropic carrier density. The low-temperature $T^2$ dependence of $\\kappa_{\\mathrm{ph}}$ is the observable signature of that $A\\omega^2$ channel.","core_discovery":"On its own terms, the paper establishes that the thermal conductivity $\\kappa$ of CeRu4Sn6 is phonon-dominated over the whole measured range and is essentially isotropic: $\\kappa$ measured along the c-axis and perpendicular to it agree closely, even though resistivity, thermopower, and optical conductivity are strongly direction-dependent. The phononic part $\\kappa_{\\mathrm{ph}} = \\kappa - \\kappa_{\\mathrm{WF}}$ behaves as $T^2$ from roughly 200 mK to 2 K, identified as scattering of phonons by electrons, and as $T^3$ below that, identified as boundary scattering with a mean free path near 0.8 mm, matching the sample dimensions. A Callaway fit that combines phonon-electron, boundary, point-defect, and Umklapp scattering reproduces the data over the full range. The near-equality of the fitted phonon-electron scattering prefactors in the two directions, combined with Hall-effect results showing an isotropic carrier concentration, is used to conclude that the strongly anisotropic resistivity of CeRu4Sn6 is mainly a mobility anisotropy.","pith_inferences":["A direct test of the paper's central attribution would be to measure $\\kappa$ at low temperatures on the same compound with the carrier concentration tuned by doping or pressure: if the fitted $T^2$ prefactor follows the Hall density, the phonon-electron assignment is confirmed, and if not, another $T^2$ channel is present.","The same logic suggests that in other Kondo semiconductors the low-temperature $T^2$ phonon conductivity could serve as a contact-free probe of mobile carrier density, useful where Hall contacts are difficult to make.","Because the bulk heat flow is phonon-dominated, future searches for topological surface transport in CeRu4Sn6 should expect thermal conductivity to remain insensitive to surface carriers; electrical transport on thin or surface-dominated samples would be the relevant probe."],"forward_implications":["If the paper is right, the strongly anisotropic resistivity of CeRu4Sn6 is a mobility anisotropy: the carrier concentration is essentially the same in the tetragonal plane and along the c-axis.","The small enhancement of $\\kappa_{\\perp c}$ near 10 K becomes a direct fingerprint of the anisotropic Kondo gap: electrons in the plane freeze out sooner, removing phonon scattering partners.","Below roughly 200 mK, thermal conductivity is set by the sample boundaries, with a mean free path of about 0.8 mm matching the smallest crystal dimension.","In the 0.2-2 K window, the $T^2$ law gives a quantitative measure of the phonon-electron scattering prefactor $A$, which can be compared across directions and with Hall carrier densities.","Thermal conductivity is not a useful probe of the electronic anisotropy in this material, because the phonon background dominates at all measured temperatures."],"supporting_citations":[{"why":"Supplies the electrical resistivity used to compute the electronic part $\\kappa_{\\mathrm{WF}}$ via the Wiedemann-Franz law and to contrast with the measured heat flow.","marker":"[13]"},{"why":"Optical conductivity and band-structure calculations that locate the gap opening in the tetragonal plane, used to explain the small anisotropy around 10 K.","marker":"[15]"},{"why":"Provides the Callaway model integral and the relaxation-time framework for the phonon thermal conductivity fits.","marker":"[17]"},{"why":"Supplies the relation between the phonon-electron scattering rate and the carrier concentration, the step that connects isotropic heat flow to isotropic carrier density.","marker":"[18]"},{"why":"Hall effect measurements on the same material finding an isotropic carrier concentration at high and low temperatures.","marker":"[19]"},{"why":"Source of the Debye temperature $\\theta_D = 250$ K used to compute the sound velocity and the boundary-scattering mean free paths.","marker":"[11]"},{"why":"Earlier observations of similar gap-anisotropy influences on thermal conductivity in CeNiSn and CeRhSb, supporting the same interpretation here.","marker":"[20–22]"}],"fun_headline_variants":["Heat flows same in all directions in CeRu4Sn6","CeRu4Sn6: phonons carry heat, electricity anisotropic","Isotropic heat transport in CeRu4Sn6, phonon-driven","CeRu4Sn6 heat: phonons dominate, no direction bias","Heat sees no anisotropy in CeRu4Sn6, phonons do it"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on assuming that the low-temperature $T^2$ rise in thermal conductivity is caused by phonons scattering off electrons with a rate proportional to frequency squared; if some other scattering channel, such as magnetic excitations, disorder, or two-level systems, also produces a $T^2$ law, then the inference of an isotropic carrier concentration loses its footing.","fun_headline_variants_meta":{"raw":{"variants":["Heat flows same in all directions in CeRu4Sn6","CeRu4Sn6: phonons carry heat, electricity anisotropic","Isotropic heat transport in CeRu4Sn6, phonon-driven","CeRu4Sn6 heat: phonons dominate, no direction bias","Heat sees no anisotropy in CeRu4Sn6, phonons do it"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001189,"raw_usage":{"total_tokens":4866,"prompt_tokens":860,"completion_tokens":4006,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":3912}},"tokens_in":476,"tokens_out":4006,"duration_ms":23496,"temperature":1.0,"reasoning_tokens":3912,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:32:05.282374+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Remeasure $\\kappa$ below 2 K on crystals with different smallest dimensions: if the $T^2$ regime is true phonon-electron scattering, its prefactor should stay fixed while the $T^3$ boundary-scattering tail moves to lower temperatures in larger crystals, whereas a geometry-dependent $T^2$ term would expose an additional scattering channel. Alternatively, tune the carrier concentration by doping or pressure and check whether the fitted phonon-electron prefactor $A$ tracks the Hall carrier density as that theory requires.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the electrical resistivity used to compute the electronic part $\\kappa_{\\mathrm{WF}}$ via the Wiedemann-Franz law and to contrast with the measured heat flow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Optical conductivity and band-structure calculations that locate the gap opening in the tetragonal plane, used to explain the small anisotropy around 10 K."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Callaway model integral and the relaxation-time framework for the phonon thermal conductivity fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the relation between the phonon-electron scattering rate and the carrier concentration, the step that connects isotropic heat flow to isotropic carrier density."},{"cited_title":"thesis Technische Universit¨ at Wien","cited_arxiv_id":null,"evidence_quote":"Hall effect measurements on the same material finding an isotropic carrier concentration at high and low temperatures."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the Debye temperature $\\theta_D = 250$ K used to compute the sound velocity and the boundary-scattering mean free paths."}],"review_version":1}