{"id":"5e4563b6-683f-498c-9959-b4875f9728ad","arxiv_id":"2501.13326","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Raman data indicate that phonons in Nd2Zr2O7 are scattered by crystal-field excitations near 100 K, while lower-energy phonons in semimetallic Nd2Ir2O7 are dominated by electron-phonon scattering.","lead":"This paper measures how lattice vibrations (phonons) in two pyrochlore magnets change with temperature using Raman scattering. It finds that in the insulator Nd2Zr2O7 phonons are scattered by crystal-field excitations of neodymium near 100 K, while in the semimetal Nd2Ir2O7 lower-energy phonons scatter off electrons.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Nd2Zr2O7 phonon–CEF claim rests on a small excess over a fitted Klemens baseline with no error bars or quantitative CEF-scattering model; overlapping CEF peaks at 34.4/35.2 meV may bias the key T2g(1) width.","rationale":"The reader's weakest assumption identifies the Klemens baseline and missing error bars as the key risk; I agree that this is the central fragile element. I add a concrete mechanism by which the deviation could be spurious: the CEF peaks at 34.4/35.2 meV overlap the T2g(1) phonon, and if not included in the spectral fits they can distort the extracted width. I also note that the paper provides no quantitative CEF-scattering model, so the claim is not tested against an expected temperature dependence. These strengthen the case for the conditional verdict but do not prove the claim false; the proposed checks would settle it. Therefore the reader's CONDITIONAL verdict remains appropriate, and I do not recommend changing it.","tokens_in":11320,"tokens_out":6048,"duration_ms":55415,"concrete_test":"Re-fit the raw Nd2Zr2O7 (x,y) spectra between 30 and 45 meV with a sum of Voigt profiles for T2g(1), Eg, and the two CEF peaks at 34.4 and 35.2 meV, with Gaussian widths fixed and CEF amplitudes free at each temperature. If the fitted T2g(1) width no longer shows the 100 K bump, the anomaly is a fitting artifact. As an independent control, measure the same Γ-point phonons in non-magnetic La2Zr2O7 (or Y2Zr2O7) under identical conditions and compare the normalized width versus temperature; if a 100 K deviation persists in the non-magnetic analog, it is not specific to Nd CEF depopulation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the Nd2Zr2O7 lattice reacts to CEF depopulation is inferred from two observations: a small excess of phonon linewidth around 100 K relative to Klemens fits (Fig. 2, right panel) and a correlation with a broad Schottky-like feature in heat capacity. Both legs are weaker than the conclusion requires. No uncertainties are reported for the fitted widths, so the excess could be within scatter. The Klemens form γ(T) = γ0 + A(2n_B(ω/2)+1) is a single-phonon-decay approximation; higher-order anharmonicity and thermal expansion can produce similar non-monotonic deviations, and the paper does not model or exclude these alternatives. Crucially, the paper never computes the expected temperature dependence of phonon linewidth from a CEF-phonon coupling term, so the mechanism is not demonstrated by the data. There is also a spectral-modeling hazard: the CEF doublet at 34.4 and 35.2 meV sits immediately below the T2g(1) phonon at 37.9 meV and appears as a growing 'wing' on cooling. If these CEF peaks are not explicitly included in the Voigt fits, the fitted T2g(1) width can change as the CEF peaks sharpen, producing a spurious bump near 100 K. Since the same normalization is applied across phonons, a common baseline error could affect all widths. The electron-phonon interpretation for Nd2Ir2O7 is less central but also has an internal inconsistency: the model in the text uses n_F(ω_a) − n_F(ω_a+ω), while Table V lists n_F(ω_a) + n_F(ω_a+ω).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports temperature-dependent Raman scattering of zone-center phonons and crystal-field (CEF) excitations in the pyrochlores Nd2Zr2O7 and Nd2Ir2O7, measured between ~10–290 K for the zirconate and ~50–270 K for the iridate. For Nd2Zr2O7, the authors observe phonon frequencies that harden on cooling but flatten or soften below about 100 K, and phonon linewidths that deviate slightly from a Klemens anharmonic baseline around 100 K. They interpret this as an additional phonon scattering channel caused by depopulation of the first excited CEF level, supported by a correlation with a Schottky-like anomaly in heat capacity. For Nd2Ir2O7, two low-frequency phonons show linewidths that decrease faster than the Klemens prediction on cooling; the authors model this with an electron-phonon interband scattering term and conclude that these phonons are broadened by scattering from electronic interband transitions. The paper also resolves a CEF doublet near 34–35 meV in Nd2Zr2O7, which was previously unresolved in neutron scattering.","tokens_in":11882,"tokens_out":3823,"duration_ms":35860,"significance":"If the central claims hold, Raman phonon linewidths and frequencies in rare-earth pyrochlores become a direct probe of crystal-field depopulation and electronic interband excitations, which would be a useful addition to the spectroscopy of frustrated magnets and semimetals. The paper has clear strengths: it presents high-quality temperature-dependent Raman data over a wide range, explicitly compares an insulator and a semimetal with ostensibly identical phonon symmetries, resolves a predicted CEF doublet, and uses standard fitting procedures for phonon parameters. The comparison between the two compounds is informative, and the data on Nd2Ir2O7 in the paramagnetic state complement existing studies of the ordered state. However, the zirconate claim rests on a small deviation from a Klemens baseline fitted to the same data, with no reported uncertainties and no quantitative model for the proposed CEF-phonon scattering contribution; the iridate claim uses a per-phonon fitted interband energy and contains an internal sign inconsistency between the text and Table V. These issues currently limit the strength of the conclusions.","major_comments":[{"comment":"The central evidence for the Nd2Zr2O7 claim is a deviation of the phonon linewidth from a Klemens baseline around 100 K, yet no uncertainties are reported for the fitted linewidths or for the Klemens parameters A and gamma0 in Table IV. Because the baseline is fitted to the same data, the excess could be within scatter. Please provide error bars or confidence intervals for the fitted widths, show fit residuals, and test whether the deviation is statistically significant relative to alternative baselines such as higher-order anharmonicity or thermal-expansion corrections.","section":"III.A.1, Fig. 2 (right panel), Table IV"},{"comment":"The CEF doublet at 34.4 and 35.2 meV sits immediately below the T2g(1) phonon at 37.9 meV and appears as a growing \"wing\" on cooling in the (x,y) channel. If these CEF peaks are not explicitly included in the Voigt fits used to extract phonon parameters, the fitted T2g(1) width can change as the CEF peaks sharpen, producing a spurious bump near 100 K. The manuscript does not state whether the CEF peaks were part of the fit model. Please report the full fitting model for each phonon and demonstrate that the 100-K anomaly in the linewidth is stable with respect to inclusion or exclusion of the CEF lines.","section":"III.A.1, Figs. 1 and 2, and III.B"},{"comment":"The text defines the electron-phonon linewidth as Gamma_ph-el(T) = Gamma0 + F (n_F(hbar omega_a, T) - n_F(hbar omega_a + hbar omega_ph, T)), but Table V and the surrounding discussion use n_F(omega_a) + n_F(omega_a + omega). This is an internal inconsistency: the relative sign determines whether the scattering decreases or increases on cooling. Please correct this and clarify which form was actually fitted. In addition, omega_a is fitted independently for each phonon and is not compared with band-structure calculations or other measurements, so the agreement in Fig. 2 is partly a fit rather than a parameter-free prediction. Please provide uncertainties for omega_a and discuss how many free parameters are used for each temperature curve.","section":"III.A.2, text before Table V, and Table V"},{"comment":"The interpretation that the lattice reacts to CEF depopulation is based on a correlation with a broad Schottky-like heat capacity feature, but the paper never computes the expected temperature dependence of the phonon linewidth from a CEF-phonon coupling term. Without such a model, the mechanism is not demonstrated: the deviation from Klemens behavior could equally arise from thermal expansion, higher-order anharmonicity, or spectral contamination. Please add a quantitative estimate of the expected linewidth anomaly, or at least a falsifiable prediction of its size and temperature dependence, so that the proposed scattering channel can be distinguished from these alternatives.","section":"IV, Discussion of Nd2Zr2O7"}],"minor_comments":[{"comment":"There are several typos and grammatical issues: \"possess\" should be \"possesses\", \"optimum\" and \"optimal\" are used inconsistently, \"parmagnetic\" should be \"paramagnetic\", and \"two orderes of magnitude\" should be \"two orders of magnitude\".","section":"Abstract and Introduction"},{"comment":"The caption states the subtraction is chi''(T=14 K) - chi''(300 K) for Nd2Zr2O7, but Fig. 1 shows spectra down to 4 K; please clarify the exact temperatures used for the subtraction and whether the 14-K value is a typo.","section":"Fig. 4 caption"},{"comment":"The background subtraction procedure for Nd2Ir2O7 is described in words but no example of the raw spectrum, the modelled background, or the subtracted spectrum is shown; a figure or supplementary material would help the reader judge the reliability of the weak phonon features.","section":"II. Experimental, background subtraction"},{"comment":"The unidentified feature P1 at 472 cm^-1 appears in both scattering channels and is described as not assignable to phonons or CEF excitations, but its temperature dependence and possible influence on the phonon fits are not discussed; please address this explicitly.","section":"III.A.1, P1 feature"},{"comment":"The statement that phonon hardening on cooling follows from lattice contraction is cited to Refs. [31,32], but Ref. [32] is an Applied Physics Letters paper and does not appear to report thermal expansion for Nd2Zr2O7; please verify the citation for the thermal-contraction data.","section":"III.A.1 and Discussion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the scope of the journal and reports useful data, but the load-bearing Nd2Zr2O7 claim needs statistical support (error bars on widths and fit parameters) and a demonstration that the CEF peaks are not contaminating the T2g(1) fit. The sign inconsistency in the electron-phonon model for Nd2Ir2O7 is straightforward to fix but must be resolved before the fit results can be interpreted. If the authors can provide these analyses, the paper would be substantially strengthened; in the current form the conclusions outrun the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is worth taking seriously but not worth taking at face value. The new Raman data on Nd2Zr2O7 and Nd2Ir2O7 are careful, the temperature range is well chosen, and the comparison between an insulator and a semimetal is a good experimental idea. The clearest new result is the resolution of the CEF doublet at 34.4 and 35.2 meV in Nd2Zr2O7, which neutron scattering had seen only as one broad band. That part is solid and useful.\n\nThe soft spot is the central claim that the lattice reacts to CEF depopulation. The evidence is a small excess of the T2g(1) linewidth over a Klemens fit near 100 K, plus a correlation with a Schottky-like anomaly in heat capacity. No error bars are given on the fitted widths, and the excess looks within plausible scatter. The paper also never computes a quantitative phonon–CEF scattering rate, so the mechanism is asserted rather than demonstrated. There is a real spectral-modeling hazard: the CEF doublet sits just below the T2g(1) phonon and grows as a wing on cooling. If those CEF peaks are not included in the Voigt fits, the phonon width can change artificially and produce a bump near 100 K. The heat-capacity correlation is suggestive, but it is not an independent confirmation.\n\nThe Nd2Ir2O7 electron-phonon interpretation is less central but also shaky. The model uses a per-phonon fitted interband energy omega_a, so the good agreement is partly a fit. On top of that, the text writes the Fermi function combination as n_F(omega_a) − n_F(omega_a+omega), while Table V lists n_F(omega_a) + n_F(omega_a+omega). Those are different functions with different temperature dependence, and the paper does not say which was actually used. That needs to be fixed and justified.\n\nWhat the paper does well: it is a careful experimental study with clean phonon assignments, a sensible polarization analysis, and a useful comparison to prior work on Pr-based pyrochlores. The citation pattern is fine. The absence of data/code makes independent checking harder, but the raw spectra appear to support the qualitative trends.\n\nThis paper deserves a serious referee, but it needs revision before publication: error bars on linewidths, a check of whether the CEF peaks bias the T2g(1) fit, a quantitative model for the phonon–CEF scattering channel, and a resolution of the sign inconsistency. I would send it to review, with the clear instruction that the phonon–CEF claim is not established yet.","headline":"Careful Raman study with a solid CEF doublet result, but the headline phonon–CEF scattering claim is under-supported by the data and the electron-phonon analysis has a fitting-consistency problem.","tokens_in":12518,"tokens_out":1316,"would_cite":true,"duration_ms":14596,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.30.-j","63.20.-e","71.38.-k"],"model":"deepseek-v4-flash","headline":"Raman phonon widths in Nd2Zr2O7 and Nd2Ir2O7 reveal scattering channels beyond phonon-phonon: crystal fields in the zirconate, electrons in the iridate.","keywords":["pyrochlore","Raman scattering","phonon linewidth","crystal-field excitations","electron-phonon scattering","Klemens model","Nd2Zr2O7","Nd2Ir2O7"],"falsifier":"Measure the temperature-dependent Raman phonon linewidths of a nonmagnetic structural analog such as Y2Zr2O7 or Lu2Zr2O7 under identical conditions; if the same 100 K deviation appears, the anomaly is not caused by Nd3+ crystal-field depopulation. Alternatively, apply a magnetic field strong enough to shift or split the Nd3+ crystal-field levels and check whether the phonon anomaly follows the shift.","tokens_in":11097,"feed_emoji":"🔬","tokens_out":7052,"duration_ms":55529,"temperature":0.7,"pith_summary":"This paper uses temperature-dependent Raman scattering to track zone-center phonons and Nd3+ crystal-field excitations in two pyrochlores, insulating Nd2Zr2O7 and semimetallic Nd2Ir2O7, in their paramagnetic temperature ranges. The authors argue that phonon linewidths and frequencies carry signatures of interactions beyond ordinary phonon-phonon scattering. In Nd2Zr2O7, a deviation near 100 K coincides with the depopulation of the first excited crystal-field level and a Schottky anomaly in heat capacity, so the lattice appears to scatter off crystal-field excitations. In Nd2Ir2O7, the two lowest phonons narrow far faster than the Klemens model predicts, and this is attributed to electron-phonon scattering as the thermal population of electronic levels decreases. If correct, Raman phonon measurements become a way to see how the lattice couples to magnetic and electronic degrees of freedom in frustrated magnets.","feed_headline":"Raman phonons reveal crystal-field and electron scattering in pyrochlores","feed_subtitle":"In Nd2Zr2O7 a 100 K anomaly links phonon damping to crystal-field depopulation; in Nd2Ir2O7, electrons dominate.","key_machinery":"The quantitative backbone is the comparison of measured phonon linewidths with the Klemens model, gamma(T) = gamma0 + A(2 n_B(omega/2) + 1), which encodes phonon-phonon scattering through a Bose factor; deviations from it signal other scattering channels. For the iridate, the alternative model gamma(T) = gamma0 + F(n_F(omega_a) - n_F(omega_a + omega_ph)) describes phonon scattering on interband electronic excitations via Fermi occupation factors. The discrimination between the two materials rests on the insulator/semimetal contrast between Nd2Zr2O7 and Nd2Ir2O7, on Raman polarization selection rules that assign phonon symmetries, and on the D3d crystal-field level scheme of Nd3+ that locates the CEF excitations.","core_discovery":"The central claim is that in both compounds the temperature dependence of Raman-active phonons cannot be explained by phonon-phonon scattering alone. For Nd2Zr2O7, the phonon width follows the Klemens model overall but deviates upward around 100 K, and phonon frequencies stop hardening or soften below that temperature; the same temperature marks a Schottky-like feature in heat capacity assigned to depopulation of the first excited crystal-field level at 23.5 meV. The paper concludes that the lattice responds to that depopulation, providing an additional scattering channel for phonons. For Nd2Ir2O7, the T2g(1) and Eg phonons lose width much more rapidly on cooling than Klemens behavior allows, and the data are described by scattering on interband electronic transitions with Fermi factors; this is interpreted as electron-phonon scattering controlled by depopulation of electronic levels. The measurements also resolve in Nd2Zr2O7 a crystal-field doublet at 34.4 and 35.2 meV that neutron scattering saw as one broad band.","pith_inferences":["If phonon-CEF scattering is the cause, the anomaly temperature should track the CEF gap: replacing Nd with a rare earth of different first-excited-level energy would shift the deviation accordingly.","A direct control experiment on a nonmagnetic structural analog, such as Y2Zr2O7 or Lu2Zr2O7, should show no 100 K phonon anomaly; its presence would refute the CEF-scattering assignment.","The electron-phonon interpretation implies that the anomalous phonon widths in Nd2Ir2O7 should respond to doping or pressure that moves the Fermi level relative to the interband threshold.","Extending the same analysis below the magnetic ordering temperature could connect the paramagnetic scattering channels identified here to the phonon anomalies already reported in the ordered state."],"forward_implications":["Phonon linewidth measurements in rare-earth pyrochlores can serve as a local probe of crystal-field level depopulation, complementing heat capacity and neutron scattering.","In pyrochlore iridates, the broadening of low-energy phonons at high temperature carries information about the population of electronic bands, not just anharmonicity.","The absence in Nd-based pyrochlores of the Eg phonon splitting seen in Pr-based ones suggests that the lattice-magnetic coupling mechanism depends on the rare-earth ground-state symmetry, consistent with theory for non-Kramers systems.","Resolving the Nd2Zr2O7 crystal-field doublet by Raman shows that optical spectroscopy can sharpen the CEF level scheme in these quantum spin-ice candidate materials."],"supporting_citations":[{"why":"Supplies the Klemens model form used as the phonon-phonon baseline for the linewidth fits.","marker":"[32]"},{"why":"Provides the Nd2Zr2O7 crystal-field level energies from neutron scattering and low-temperature lattice data used to identify the 100 K anomaly.","marker":"[31]"},{"why":"Provides the heat-capacity Schottky anomaly near 100 K that anchors the assignment to crystal-field depopulation.","marker":"[34]"},{"why":"Provides the Pr2Ir2O7 Raman study and electron-phonon model that the Nd2Ir2O7 analysis follows.","marker":"[29]"},{"why":"Shows the electronic excitations in the semimetallic bands with which the lower Nd2Ir2O7 phonons overlap in energy.","marker":"[16]"},{"why":"Supplies the theoretical prediction that phonon Raman scattering can indirectly probe spin-ice-type magnetic excitations, motivating the search for magnetic scattering channels.","marker":"[25]"},{"why":"Documents the temperature behavior of crystal-field Raman scattering used to distinguish CEF peaks from phonons.","marker":"[28]"},{"why":"Provides phonon assignments and the Pr2Zr2O7 Raman comparison, including the Eg phonon splitting.","marker":"[20]"}],"fun_headline_variants":["Phonon damping shows crystal-field and electron scattering in pyrochlores","Raman phonons trace crystal-field depopulation and electron scattering","Pyrochlore phonons respond to crystal-field and electronic excitations","Phonon width anomalies reveal crystal-field and electron scattering","Crystal-field depopulation and electrons scatter Raman phonons in pyrochlores"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The Nd2Zr2O7 conclusion assumes that the Klemens model, fitted to the same data, is the correct baseline for ordinary phonon-phonon scattering; the reported deviation near 100 K is small and no error bars are given, so anharmonicity, thermal expansion, or background artifacts could in principle absorb it.","fun_headline_variants_meta":{"raw":{"variants":["Phonon damping shows crystal-field and electron scattering in pyrochlores","Raman phonons trace crystal-field depopulation and electron scattering","Pyrochlore phonons respond to crystal-field and electronic excitations","Phonon width anomalies reveal crystal-field and electron scattering","Crystal-field depopulation and electrons scatter Raman phonons in pyrochlores"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000759,"raw_usage":{"total_tokens":3397,"prompt_tokens":996,"completion_tokens":2401,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":2307}},"tokens_in":612,"tokens_out":2401,"duration_ms":17144,"temperature":1.0,"reasoning_tokens":2307,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:15:14.515507+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the temperature-dependent Raman phonon linewidths of a nonmagnetic structural analog such as Y2Zr2O7 or Lu2Zr2O7 under identical conditions; if the same 100 K deviation appears, the anomaly is not caused by Nd3+ crystal-field depopulation. Alternatively, apply a magnetic field strong enough to shift or split the Nd3+ crystal-field levels and check whether the phonon anomaly follows the shift.","supporting_citations":[{"cited_title":"Watahiki, K","cited_arxiv_id":null,"evidence_quote":"Supplies the Klemens model form used as the phonon-phonon baseline for the linewidth fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Nd2Zr2O7 crystal-field level energies from neutron scattering and low-temperature lattice data used to identify the 100 K anomaly."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the heat-capacity Schottky anomaly near 100 K that anchors the assignment to crystal-field depopulation."},{"cited_title":"Zoghlin, J","cited_arxiv_id":null,"evidence_quote":"Provides the Pr2Ir2O7 Raman study and electron-phonon model that the Nd2Ir2O7 analysis follows."},{"cited_title":"Witczak-Krempa, G","cited_arxiv_id":null,"evidence_quote":"Shows the electronic excitations in the semimetallic bands with which the lower Nd2Ir2O7 phonons overlap in energy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical prediction that phonon Raman scattering can indirectly probe spin-ice-type magnetic excitations, motivating the search for magnetic scattering channels."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the temperature behavior of crystal-field Raman scattering used to distinguish CEF peaks from phonons."},{"cited_title":"Ramification of complex magnetism in Nd$_2$Ir$_2$O$_7$ observed by Raman scattering spectroscopy","cited_arxiv_id":"2302.00579","evidence_quote":"Provides phonon assignments and the Pr2Zr2O7 Raman comparison, including the Eg phonon splitting."}],"review_version":1}