{"id":"911cedce-71d5-474a-8095-78e39d2567e4","arxiv_id":"2411.15052","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"In cage compounds the effective Grüneisen parameter is temperature-dependent and exhibits two approximate regimes, with measured crossovers at about 11 K in LaB6 and 10.5 K in LaPt4Ge12.","lead":"The authors extend a two-frequency phonon model to predict how thermal expansion behaves in cubic cage compounds, then test it on LaB6 and LaPt4Ge12. They find that the Grüneisen rule, normally a single constant, switches between two regimes at low temperature, which could help separate crystal vibrations from magnetic effects in rare-earth cage materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The LaPt4Ge12 confirmation is not decisive: the observed 10.5 K crossover is interpreted with a two-branch model that the paper's own specific-heat analysis shows is already failing near 15 K, so the inferred ħω0 ≈ 7.5 meV may reflect neglected phonon branches rather than a successful prediction.","rationale":"The paper is a transparent extension of the authors' earlier two-frequency model, and the LaB6 analysis is credible: the model fits the specific heat over a wide range and the predicted crossover temperature returns a frequency consistent with neutron data. The vulnerability is concentrated in LaPt4Ge12, the second example named in the title. There the measured Tco = 10.5 K is converted into ħω0 ≈ 7.5 meV using the model's Dph/Cph curve, but the model's own Cp fit gives ħω0 ≈ 5 meV and fails above about 15 K (Sec. 4.1.2). Since additional phonon branches are already perturbing the thermodynamics near the crossover temperature, the two-branch Dph/Cph curve is not a reliable transfer function for this material. The proposed check—extending the phonon DOS and repeating the Grüneisen analysis—would settle whether the crossover is genuinely produced by the two-branch mechanism. This is an addressable limitation, not a fundamental flaw, and the authors explicitly flag the model's inadequacy for LaPt4Ge12. I therefore agree with the reader's weakest assumption and would keep the CONDITIONAL verdict; no change to the reader's assessment is needed.","tokens_in":16805,"tokens_out":10053,"duration_ms":102626,"concrete_test":"Reanalyze LaPt4Ge12 with a phonon model that includes the measured higher branches: add Einstein oscillators at the neutron-derived peaks near 7.5, 10.5, 12.5 and 15.5 meV (or use the measured neutron-weighted DOS) on top of the two-branch model, refit the specific heat, and recompute α(T)/Cp(T) with the same Dph weighting. If the crossover near 10.5 K and the recovered ħω0 ≈ 7.5 meV survive with the extended DOS, the concern is resolved; if the crossover shifts or disappears, the two-branch interpretation for LaPt4Ge12 is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central validation step uses the measured crossover Tco together with the model's Dph/Cph curves (Fig. 5) and Eq. (25) to recover ħω0: Sec. 4.2.1 uses kBTco/ħω0 ≈ 0.07 for LaB6 and Sec. 4.2.2 uses kBTco/ħω0 ≈ 0.12 for LaPt4Ge12, giving 13.5 meV and 7.5 meV. For LaPt4Ge12 this inference is load-bearing and insecure. The paper's own specific-heat analysis (Sec. 4.1.2) states that the two-frequency fit fails above about 15 K, with fitted ħω0 ≈ 5 meV that conflicts with the neutron peak at 7.5 meV, and Sec. 5 attributes the failure to neglected vibrations within the Pt4Ge12 cage. Because the observed crossover at 10.5 K is only about 5 K below the temperature where unmodeled branches become influential, the kink in α/T versus Cp/T can originate from the same neglected branches that distort the Cp fit, rather than from the Dph/Cph step of the two-branch model. In that case the claimed consistency with neutron data for LaPt4Ge12 is not an independent confirmation: it selects the neutron value after the fact while the Cp fit points to a different value. The LaB6 example is considerably safer because the same model fits the specific heat up to roughly 50 K, but the paper's joint claim relies on both examples.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends a previously proposed two-frequency phonon model for rare-earth cage compounds (Ref. [7]) to thermal expansion within the quasi-harmonic approximation. The central result is an effective Grüneisen function γeff(T) = Γ0{1 + (1 - γ0/Γ0) Dph(T)/Cph(T)} (Eq. 25), which predicts that the Grüneisen rule is violated at low temperatures and approximately restored at intermediate temperatures where the flattened phonon branches dominate Dph/Cph. The authors compare this prediction with specific heat and thermal expansion measurements on LaB6 and LaPt4Ge12. For LaB6 the model fits the specific heat below about 50 K and the observed crossover at ~11 K yields ħω0 ≈ 13.5 meV, consistent with neutron data. For LaPt4Ge12 the specific heat fit is only acceptable below about 15 K, and the crossover at ~10.5 K is used to estimate ħω0 ≈ 7.5 meV, again claimed to be consistent with neutron spectroscopy. The paper concludes that the two-frequency model predicts distinct Grüneisen regimes in both compounds, which is intended to help separate phonon and magnetic contributions to thermal expansion in the broader RB6 and RPt4Ge12 series.","tokens_in":17252,"tokens_out":8316,"duration_ms":73209,"significance":"If the central claim holds, the paper offers a simple analytical framework for the low-temperature thermodynamics of cage compounds, showing in a transparent way why the standard Grüneisen rule fails and how an approximate rule is restored at intermediate temperatures. The derivation of Eq. (25) is internally consistent, the model is tested against independent thermal expansion data rather than fitted to those data, and the LaB6 example provides a genuine quantitative prediction (ħω0 ≈ 13.5 meV) that matches neutron spectroscopy. The authors are also explicit about the model's limitations, particularly for LaPt4Ge12. However, the LaPt4Ge12 confirmation is substantially weaker: the model is known to fail above ~15 K, and the crossover-based frequency estimate is entangled with the neutron data that were used to motivate the model's parameter choices. Because the paper's joint claim rests on both examples, this weakness is load-bearing, not merely cosmetic.","major_comments":[{"comment":"The inference of ħω0 ≈ 7.5 meV for LaPt4Ge12 from the crossover at Tco = 10.5 K is load-bearing and insecure. Section 4.1.2 states that the two-frequency model fails above about 15 K, that the specific-heat fit yields a peak slightly above 5 meV against the neutron peak at 7.5 meV, and that below 20 K the fit compensates for missing branches by lowering the energies of the included branches. Since the observed crossover lies only about 4.5 K below the model's known failure temperature, the kink in α/T versus Cp/T (Fig. 11) could originate from the same neglected branches that distort the Cp fit, rather than from the Dph/Cph step of the two-branch model. Moreover, case c) was selected in Sec. 2.2 partly because its density of states resembles the neutron spectrum, so the agreement between the crossover-derived ħω0 and the neutron value is not an independent confirmation. Please either extend the analysis to include the additional branches or provide a quantitative assessment of the systematic error in ħω0 before claiming consistency with neutron data.","section":"Sec. 4.2.2"},{"comment":"The Grüneisen regime boundaries (13–45 K for LaB6, 11–21 K for LaPt4Ge12) and the crossover temperatures (11 K and 10.5 K) are identified by visual inspection of Figs. 10 and 11, and the reported slopes χ0γeff/3 are given without uncertainties. Since the central claim is that the model predicts the location of the step in γeff(T) through kBTco/ħω0, the analysis should include a quantitative criterion (e.g., segmented linear regression with confidence intervals) to establish the linear regions, to assign errors to Tco and to the slopes, and to test whether the observed slopes are compatible with a constant γeff within each regime. Without such estimates, the claimed quantitative agreement with kBTco/ħω0 ≈ 0.07 and 0.12 cannot be assessed.","section":"Secs. 4.2.1 and 4.2.2"}],"minor_comments":[{"comment":"The synthesis reference for LaB6 is missing: the text cites 'reference [ ?]' in Sec. 3.1; please supply the correct citation.","section":"Sec. 3.1"},{"comment":"The sentence introducing Eq. (21) refers to 'the averaged top frequency Γ0', but the symbol Γ0 is defined as the Grüneisen parameter of Ω0; the text should say 'the averaged top frequency Ω0'.","section":"Eq. (21)"},{"comment":"There are several typographical errors: 'skuterrudites' should be 'skutterudites', 'againts' in the caption of Fig. 10 should be 'against', 'strait line' should be 'straight line', 'arrises' should be 'arises', and 'distincts' should be 'distinct'.","section":"Throughout"},{"comment":"The text uses the formats 'γef f' and 'γeff' interchangeably; please use a single consistent notation for the effective Grüneisen function.","section":"Secs. 4.2.1 and 4.2.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a condensed-matter physics journal, and the LaB6 analysis is a solid contribution. The LaPt4Ge12 confirmation, however, needs either a more complete phonon model or a carefully quantified caveat; as written, the claimed cross-check is too close to circularity because the neutron data informed the model selection. The authors should be asked to address the uncertainty-quantification issue as well, since the central claim depends on the visual identification of linear regimes."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper takes the two-frequency phonon model for cage compounds and extends it to thermal expansion in the quasi-harmonic approximation, then tests the resulting effective Grüneisen function against new specific heat and thermal expansion data for LaB6 and LaPt4Ge12. The extension is genuinely new: the authors get an expression for α_ph containing C_ph plus an extra D_ph term, which gives a temperature-dependent γ_eff that steps up from Γ0 to a plateau. That is a concrete, falsifiable prediction.\n\nThe LaB6 part works well. The model fits the specific heat below 50 K, the α/T versus C_p/T plot shows two linear regimes with a crossover near 11 K, and the crossover-derived ħω0 ≈ 13.5 meV lands on the neutron and specific heat values. That is a real success and makes the paper useful for the RB6 family.\n\nThe LaPt4Ge12 analysis is the soft spot. The model's own specific heat fit is inadequate above about 15 K, and its fitted ħω0 ≈ 5 meV disagrees with the neutron peak at 7.5 meV. The crossover in the Grüneisen plot sits at 10.5 K, only a few kelvin below where unmodeled branches start to matter. Using that crossover to recover ħω0 ≈ 7.5 meV is not an independent confirmation; it selects the neutron value while the Cp fit points elsewhere. The authors are candid about the model failure, but the inference from Tco is more fragile than they acknowledge.\n\nMinor issues: the Grüneisen regime boundaries are chosen by eye, the fitted slopes have no error bars, and the raw data are not deposited. These are addressable.\n\nThe paper deserves a serious referee. The central LaB6 claim seems solid, and the model extension is worth having in the literature for people working on caged compounds. My recommendation: send it to review, ask for error bars on the slopes and a more careful treatment of the LaPt4Ge12 crossover, including a discussion of what additional branches do to the D_ph/C_ph ratio.\n\nBest,","headline":"The quasi-harmonic extension of the two-frequency cage model gives a real prediction, and LaB6 confirms it; the LaPt4Ge12 test is much thinner than the paper claims.","tokens_in":17796,"tokens_out":2596,"would_cite":true,"duration_ms":24472,"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":"A two-frequency cage model shows that the Grüneisen rule—constant proportionality between thermal expansion and heat capacity—is violated at low temperatures but approximately restored at intermediate temperatures, with the two regimes…","keywords":["cage compounds","rattling modes","Grüneisen parameter","thermal expansion","quasi-harmonic approximation","specific heat","LaB6","LaPt4Ge12"],"falsifier":"Measure $\\alpha/T$ versus $C_p/T$ on a cage compound with an independently known guest frequency $\\omega_0$; the two-frequency model predicts the slope-change crossover at $k_BT_{\\mathrm{co}}/\\hbar\\omega_0 \\approx 0.07$ for a heavy guest and $\\approx 0.12$ for a light guest. If the observed $T_{\\mathrm{co}}$ deviates from these values by more than the experimental uncertainty, the central quantitative prediction is wrong.","tokens_in":16595,"feed_emoji":"🌡️","tokens_out":11673,"duration_ms":97686,"temperature":0.7,"pith_summary":"The paper sets out to show that in rare-earth cage compounds—crystals where a guest atom sits loosely inside a rigid cage—the usual Grüneisen rule, a constant proportionality between thermal expansion and heat capacity, is not a valid description even at low temperatures. Using a two-frequency quasi-harmonic model of the two lowest phonon branches, the authors derive an effective Grüneisen function that switches between two nearly constant regimes: one set by the cage lattice at very low temperatures and another set by the flattened guest-atom branches at intermediate temperatures. They then measure the thermal expansion and specific heat of LaB6 and LaPt4Ge12 and find the predicted change of regime, with crossovers near 11 K and 10.5 K that give guest vibration frequencies consistent with neutron data. The result matters because a reliable phonon background is the first step toward isolating crystal-field and magnetic contributions to thermal expansion in magnetic members of these families.","feed_headline":"Cage compounds break the Grüneisen rule—then restore it","feed_subtitle":"LaB6 and LaPt4Ge12 show the predicted two-regime crossover, fixing guest vibration frequencies without neutron data","key_machinery":"The central object is the two-frequency phonon model for a cubic lattice of cages, which keeps only the two lowest phonon branches. A guest of mass $m$ oscillates in a cage of mass $M$ with natural frequency $\\omega_0$, while the empty-cage lattice has an averaged zone-boundary frequency $\\Omega_0 = \\omega_0 Y_0$; the mass ratio $\\rho = m/M$ sets the gap between the acoustic branch, which flattens near $\\omega_0$, and the optical branch starting at $\\omega_0\\sqrt{1+\\rho}$. Volume dependence enters through two Grüneisen parameters, $\\gamma_0$ for the guest frequency and $\\Gamma_0$ for the cage frequency, linked mode-by-mode through the weight function $f_\\rho(y)=\\rho/((1-y^2)^2+\\rho)$. The function $D_{\\mathrm{ph}}(T)$ collects the specific-heat-like weight of modes close to $\\omega_0$, and its ratio to $C_{\\mathrm{ph}}(T)$ is what makes $\\gamma_{\\mathrm{eff}}(T)$ temperature dependent.","core_discovery":"Within the quasi-harmonic approximation, the phonon contribution to thermal expansion in a cubic cage compound is not proportional to the phonon specific heat. The two-frequency cage model gives $\\alpha_{\\mathrm{ph}}(T) = (\\chi_0 \\Gamma_0/3)[C_{\\mathrm{ph}}(T) + (1-\\gamma_0/\\Gamma_0)D_{\\mathrm{ph}}(T)]$ and hence $\\gamma_{\\mathrm{eff}}(T)=\\Gamma_0[1+(1-\\gamma_0/\\Gamma_0)D_{\\mathrm{ph}}(T)/C_{\\mathrm{ph}}(T)]$. The ratio $D_{\\mathrm{ph}}/C_{\\mathrm{ph}}$ vanishes at low temperature and then plateaus at intermediate temperatures, so the effective Grüneisen function takes two quasi-constant values separated by a crossover near $k_BT_{\\mathrm{co}}/\\hbar\\omega_0 \\approx 0.07$ for a heavy guest (LaB$_6$) and $\\approx 0.12$ for a light guest (LaPt$_4$Ge$_{12}$). The measured $\\alpha/T$ versus $C_p/T$ graphs for LaB$_6$ and LaPt$_4$Ge$_{12}$ show the predicted slope changes, at about 11 K and 10.5 K, from which the inferred $\\omega_0$ values agree with inelastic-neutron determinations. The model also describes the specific heat of LaB$_6$ below 50 K, while for LaPt$_4$Ge$_{12}$ the specific-heat fit fails above about 15 K because additional low-energy branches of the heavy cage are not included.","pith_inferences":["A practical extension: in magnetic cage compounds where neutron data are unavailable, the $\\alpha/T$ versus $C_p/T$ crossover could be used to locate the guest frequency, since only thermal expansion and specific heat are needed.","The sign of the $\\gamma_{\\mathrm{eff}}$ step—higher $\\gamma_{\\mathrm{eff}}$ below versus above the crossover in LaB6, and lower in LaPt4Ge12—points to a systematic dependence on the guest-to-cage mass ratio; testing this across other cage families would give a rule for how $\\gamma_0$ and $\\Gamma_0$ compare.","If the two-regime Grüneisen signature is generic, then clathrates and skutterudites with well-separated flat phonon branches should show the same slope change in $\\alpha/T$ against $C_p/T$, providing a thermodynamic fingerprint of rattling.","For magnetic rare-earth cage compounds, subtracting a non-magnetic reference by rescaling a Debye temperature will not remove the phonon contribution cleanly, because the two-regime Grüneisen background cannot be represented by a single effective parameter."],"forward_implications":["A constant Grüneisen parameter cannot be assumed when analysing thermal expansion of cage compounds; the phonon background itself has a two-regime structure that must be subtracted before magnetic contributions are interpreted.","The crossover temperature between the two regimes provides a direct estimate of the guest vibration frequency $\\omega_0$, giving values consistent with neutron scattering for both LaB6 and LaPt4Ge12.","For LaB6, the two-frequency model accurately reproduces the specific heat up to about 50 K (roughly $\\hbar\\omega_0/3k_B$), so thermodynamic analyses of magnetic hexaborides below that scale can rely on it.","For filled skutterudites like LaPt4Ge12, the same model is only reliable at the lowest temperatures; by about 15 K, neglected cage vibrations dominate, so quantitative fits to specific heat become misleading."],"supporting_citations":[{"why":"Supplies the two-frequency phonon model of rare-earth hexaborides that the paper extends to thermal expansion.","marker":"[7]"},{"why":"Provides the measured phonon dispersion of LaB6 from which the guest frequency near 13.5 meV is taken and the flattened-branch picture is drawn.","marker":"[16]"},{"why":"Gives the inelastic neutron spectra of LaPt4Ge12 with peaks at 7.5 and 12.5 meV, the reference for comparing fitted and true guest frequencies.","marker":"[19]"},{"why":"Supplies the quasi-harmonic phonon-pressure formalism connecting volume derivatives of phonon frequencies to thermal expansion.","marker":"[20]"},{"why":"Defines the Grüneisen rule whose validity in cage compounds the paper tests.","marker":"[8]"},{"why":"Provides the low-temperature elastic constants of LaB6 needed to convert the measured $\\alpha/T$–$C_p/T$ slope into an effective Grüneisen parameter.","marker":"[24]"},{"why":"Documents the superconducting transition and electronic specific heat coefficient of LaPt4Ge12, used to set $\\gamma_e$ and exclude the superconducting anomaly.","marker":"[9]"}],"fun_headline_variants":["Grüneisen rule broken by cage phonons","Two-regime Grüneisen from cage phonons","Cage compounds defy the Grüneisen rule","Cage phonons cause a Grüneisen crossover","Thermal expansion reveals cage Grüneisen steps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on assuming that two phonon branches, described by one averaged zone-boundary frequency, dominate the low-temperature thermodynamics; the paper itself shows this fails for LaPt$_4$Ge$_{12}$ above about 15 K, where missing branches pull the fitted guest frequency below the neutron value.","fun_headline_variants_meta":{"raw":{"variants":["Grüneisen rule broken by cage phonons","Two-regime Grüneisen from cage phonons","Cage compounds defy the Grüneisen rule","Cage phonons cause a Grüneisen crossover","Thermal expansion reveals cage Grüneisen steps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000728,"raw_usage":{"total_tokens":3332,"prompt_tokens":1088,"completion_tokens":2244,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":2184}},"tokens_in":704,"tokens_out":2244,"duration_ms":16887,"temperature":1.0,"reasoning_tokens":2184,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:33:51.868034+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\alpha/T$ versus $C_p/T$ on a cage compound with an independently known guest frequency $\\omega_0$; the two-frequency model predicts the slope-change crossover at $k_BT_{\\mathrm{co}}/\\hbar\\omega_0 \\approx 0.07$ for a heavy guest and $\\approx 0.12$ for a light guest. If the observed $T_{\\mathrm{co}}$ deviates from these values by more than the experimental uncertainty, the central quantitative prediction is wrong.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-frequency phonon model of rare-earth hexaborides that the paper extends to thermal expansion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the inelastic neutron spectra of LaPt4Ge12 with peaks at 7.5 and 12.5 meV, the reference for comparing fitted and true guest frequencies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quasi-harmonic phonon-pressure formalism connecting volume derivatives of phonon frequencies to thermal expansion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the low-temperature elastic constants of LaB6 needed to convert the measured $\\alpha/T$–$C_p/T$ slope into an effective Grüneisen parameter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the superconducting transition and electronic specific heat coefficient of LaPt4Ge12, used to set $\\gamma_e$ and exclude the superconducting anomaly."}],"review_version":1}