{"id":"000ed68b-35cb-414c-b392-269940be0e7b","arxiv_id":"2501.07871","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Scattering by isotropic disorder slows the relaxation of acoustic phonon polarization because the polarization oscillates between collisions and the oscillations average out, giving relaxation rates proportional to the scattering time.","lead":"This paper predicts that random disorder can actually lengthen the lifetime of acoustic phonon polarization, including phonon angular momentum, through motional narrowing. The result is an analog of the Dyakonov-Perel spin relaxation mechanism applied to phonons, with linear polarizations relaxing more slowly than circular ones.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The scalar, branch-averaged collision integral (Eq. 3) omits order-one L↔T branch-flip scattering by mass defects, so the claimed τ_AM ∝ 1/τ* enhancement is not yet established for the solid-argon disorder model.","rationale":"The reader's weakest assumption is exactly the load-bearing point. I agree: the model is self-consistent and the derivation of Eqs. (8)–(10) from the scalar-K kinetic equation appears correct. However, the physical system used for validation (mass disorder in solid argon) does not by itself satisfy the scalar-K assumption; mass defects scatter between L and T branches with matrix elements of order unity. The paper explicitly flags this in footnote [75] and imposes condition (ii), but never checks it. Therefore the central counterintuitive claim—stronger disorder enhances AM lifetime—is not a demonstrated property of disordered crystals; it is a property of a restricted model. The multipole truncation is a lesser concern because the neglected l≥3 components are fast and enter only at higher order in Ωτ*. The right response is to keep the CONDITIONAL verdict and require the branch-resolved calculation or a clear argument that the Elliott-Yafet channel is suppressed. The paper is honest about limitations, and the analytical framework is valuable; no rejection is warranted.","tokens_in":21294,"tokens_out":7541,"duration_ms":85697,"concrete_test":"Using the solid-argon IFCs of Ref. [90] and the standard Tamura mass-defect golden rule 1/τ_{qs→q's'} ∝ Γ_md ω_q ω_{q'} |e_{qs}·e_{q's'}|^2 δ(ω_{qs}−ω_{q's'}), compute the full branch-resolved relaxation matrix for 36Ar(1−x)136Xe. Diagonalize the kinetic equation for ρ^a_00 including both the double-commutator dephasing term and this collision integral. Compare 1/τ_AM with Fig. 3(b) at q a/(2π)=0.4 for α·(2π/a)^4=28 ps^-1. If the full rate is not below 1/τ* and does not decrease with increasing α, the motional-narrowing mechanism is masked by the Elliott-Yafet channel.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result rests on Eq. (3), in which K_{q,q'} is a polarization-independent scalar and energy conservation uses the averaged dispersion ℏ v q. For isotopic mass disorder, the exact golden-rule matrix element is proportional to e_{q,s}·e_{q',s'}; after angular averaging, |e_L(q)·e_T(q')|^2 is of order 1/3, not suppressed by (c_L−c_T)/v. Hence mass disorder generates an Elliott-Yafet-like branch-flip channel with rate 1/τ_EY comparable to the momentum relaxation rate 1/τ*. The paper's condition (ii), that polarization-conserving scattering dominates polarization-flip scattering, is asserted but not verified for the 36Ar/136Xe system. If 1/τ_EY ~ 1/τ*, the total AM decay rate is roughly A (c_L−c_T)^2 q^2 τ* + B/τ*; the second term grows with disorder and can overwhelm the motional-narrowing term at the q values and α values in Fig. 3(b). Footnote [75] concedes the approximation is 'not universally valid.' Thus the quantitative rates in Fig. 3 and the variational prediction are conditional on a calculation that is not provided.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a motional-narrowing mechanism for the relaxation of phonon polarization in disordered crystals. The authors introduce eight Gell-Mann polarization parameters for acoustic phonons, show that band splitting causes coherent polarization oscillations, and derive from a quantum kinetic equation with a scalar elastic collision integral that the circular-polarization (AM) and linear-polarization relaxation rates are 1/τ_AM = (2/3)(c_L−c_T)^2 q^2 τ_2 and 1/τ_LP = (2/5)(c_L−c_T)^2 q^2 τ_2 for isotropic elastic bodies, together with a general formula (10) applicable to cubic crystals. They apply the formula to solid argon with mass-defect disorder, predict 1/τ_Pol ∝ q^{-2} to q^{-4}, and conclude that stronger disorder lengthens polarization lifetimes and that linear polarizations decay more slowly than circular polarizations. The Supplemental Material provides the Gell-Mann structure constants, spherical-harmonic expansion coefficients, explicit submatrices, and the mass-disorder strength estimate.","tokens_in":21553,"tokens_out":4544,"duration_ms":45874,"significance":"If the result holds beyond the scalar-collision approximation, the paper would establish a genuinely counterintuitive and falsifiable effect: polarization lifetimes increasing with disorder, with a clean separation between LP and CP rates. The analytic prefactors 2/3 and 2/5 are derived without fitted parameters, the multipole elimination is explicit, and the solid-argon estimate gives concrete numbers for experiment, including the prediction that the AM and LP diffusion lengths are independent of impurity concentration. These are real strengths. The main risk is not internal consistency of the derivation but the physical completeness of the collision model in Eq. (3), which the authors themselves flag in footnote [75] as not universally valid.","major_comments":[{"comment":"The central quantitative claim rests on a collision integral in which K_{q,q'} is a scalar and energy conservation uses the averaged dispersion ℏvq. For the isotopic mass disorder used in Fig. 3, the exact golden-rule matrix element is proportional to e_{q,s}·e_{q',s'}, so longitudinal-to-transverse branch-flip scattering is not suppressed by (c_L−c_T)/v; after angular averaging the Elliott-Yafet-like channel has a rate 1/τ_EY comparable to 1/τ*. Since footnote [75] concedes that Eq. (3) is 'not universally valid,' the paper needs an explicit estimate of 1/τ_EY for the 36Ar/136Xe model, or a separate calculation with the exact matrix element, before the rates in Fig. 3(b) and the claim τ_AM ∝ 1/τ* can be accepted for mass disorder. This is load-bearing because the total AM decay rate would be roughly A(c_L−c_T)^2 q^2 τ* + B/τ_EY, and the second term grows with disorder.","section":"Eq. (3) and footnote [75]"},{"comment":"The authors assert that the polarization-conserving scattering rate 1/τ* is higher than the polarization-flip scattering rate, but no calculation of the latter is provided. For the solid-argon example the disorder model itself generates branch-flip scattering, so this condition is not automatically satisfied. The authors should either verify the inequality quantitatively for the mass-substitution parameters used in Fig. 3, or state the prediction conditionally on it; without this, the variational prediction and the experimental estimate based on 36Ar_{1−x}136Xe_x are incomplete.","section":"Condition (ii) after Eq. (10)"},{"comment":"The replacement of τ_2 by τ* ≡ max{τ_l} in Eq. (10) is an uncontrolled approximation when the l=1 relaxation time differs strongly from τ_2. For Rayleigh scattering from mass defects, τ_1 and τ_2 are generally different, and the numerical rates in Fig. 3(b) would change if the actual τ_2 from the mass-disorder model were used instead of τ*. The authors should justify this replacement or use the appropriate τ_2 in the numerical calculation; this is secondary to the branch-flip issue but affects the quantitative predictions.","section":"Eq. (10) and the definition of τ*"}],"minor_comments":[{"comment":"The symbol O is used for zero submatrices in Eq. (7) but is not defined; please define it or use bold zero notation.","section":"Eq. (7)"},{"comment":"The sentence 'This is an analog of precession of an electron spin in the three-dimensional space' is informal; specify the analogy as precession on the eight-dimensional hypersphere described in the Supplement.","section":"After Eq. (2)"},{"comment":"Please check the numerical prefactor and units in the mass-difference scattering rate in Eq. (S46); the surrounding text should specify the density-of-states convention used, since the conventional Tamura/Klemens result can differ by factors of π/2 or by the unit-cell volume factor.","section":"Eq. (S46)"},{"comment":"The three disorder strengths are described in the caption but the curves are not individually labeled in the figure; adding labels or a legend would improve readability.","section":"Fig. 3(b)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the journal's readership if the branch-flip scattering issue can be resolved. The referee report focuses on the collision integral in Eq. (3); the authors should be asked to provide a quantitative estimate of the Elliott-Yafet-like channel for the solid-argon model. If they can show that 1/τ_EY is negligible compared to 1/τ* in the regime of Fig. 3(b), the central claim would be substantially strengthened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nHere's my take on Suzuki-Murakami. The paper's core is a real contribution: it sets up the eight-component polarization structure of acoustic phonons, derives the Dyakonov-Perel-type relaxation rates 1/τ_AM = (2/3)(cL−cT)^2 q^2 τ2 and 1/τ_LP = (2/5)(cL−cT)^2 q^2 τ2, and shows linear polarization relaxes slower than circular. The supplement is thorough, with explicit Gell-Mann and spherical-harmonic coefficients and a calculation for solid argon. Within the stated model, the math checks out.\n\nThe soft spot is exactly where the reader put it: the kinetic equation (3) uses a scalar, branch-averaged collision integral. That is a standard approximation in electron Dyakonov-Perel, but for phonons the 'spin' is the polarization vector itself, and isotope-mass disorder has a golden-rule matrix element proportional to e_{q,s}·e_{q',s'}. For L to T branch flips that overlap is order one, not suppressed by (cL−cT)/v. So there is a plausible Elliott-Yafet-like channel that relaxes polarization at a rate comparable to 1/τ*. The authors assert condition (ii), that polarization-conserving scattering dominates, but they don't verify it for the 36Ar/136Xe example. Footnote [75] explicitly concedes the approximation 'is not universally valid.' If a branch-flip rate B/τ* is present, the total AM decay rate is roughly A Ω^2 τ* + B/τ*, and for the q and α in Fig. 3(b) the second term can win. That would overturn the quantitative prediction for that system, though the motional-narrowing idea could still survive if a cleaner realization is found.\n\nThe multipole truncation at l=0,2 is a smaller issue; it's standard and the authors note the couplings.\n\nNet: this is a serious, readable theory paper that deserves referee time. My recommendation is to send it out, with a referee asked to estimate the branch-flip scattering rate for mass disorder or otherwise justify condition (ii) for the specific example. If that can't be done, the quantitative claims should be softened to a model statement. As is, I'd cite it as the first derivation of this mechanism, with a caveat.\n\nBest.","headline":"Clean derivation of a phonon-polarization motional-narrowing mechanism, but the mass-disorder example may be dominated by the very branch-flip scattering the model leaves out.","tokens_in":22094,"tokens_out":3932,"would_cite":true,"duration_ms":43112,"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":"This paper proposes that elastic scattering by isotropic disorder slows the relaxation of acoustic-phonon polarization through motional narrowing, so stronger disorder lengthens polarization lifetimes rather than shortening them.","keywords":["phonon angular momentum","phonon polarization","motional narrowing","Dyakonov-Perel mechanism","disorder scattering","quantum kinetic equation","circular polarization","linear polarization"],"falsifier":"Measure the circular- and linear-polarization relaxation times in a family of isotopically mass-substituted crystals (e.g., $^{36}$Ar$_{1-x}^{136}$Xe$_x$) at fixed temperature and wavevector $qa/2\\pi\\simeq0.4$, varying the substitution fraction $x$. The paper predicts $1/\\tau_{\\mathrm{AM}}$ and $1/\\tau_{\\mathrm{LP}}$ to decrease with $x$ and the diffusion lengths $\\sqrt{D\\tau}$ to stay constant; observing the opposite—rates increasing with impurity content—would falsify the mechanism.","tokens_in":21046,"feed_emoji":"⏳","tokens_out":13424,"duration_ms":117160,"temperature":0.7,"pith_summary":"The authors argue that in disordered crystals, isotropic elastic scattering protects—rather than destroys—the polarization of acoustic phonons. The mechanism is motional narrowing: longitudinal–transverse band splitting makes polarization precess, and frequent collisions randomize the precession so its net effect averages away, leaving a slow decay. Their kinetic calculation gives a circular-polarization (phonon angular momentum) relaxation time $\\tau_{\\mathrm{AM}} \\propto 1/\\tau^{*}$, so more impurity scattering means longer-lived phonon angular momentum, and linear polarizations relax even more slowly. If correct, this provides a concrete route to extending phonon angular-momentum lifetimes in mass-disordered crystals such as solid argon.","feed_headline":"More disorder slows phonon polarization decay, not speeds it","feed_subtitle":"Motional narrowing: impurity scattering lengthens phonon polarization lifetimes, with linear lasting longest","key_machinery":"The engine of the argument is the quantum kinetic equation for the $3\\times3$ phonon density matrix, expanded in the eight Gell-Mann matrices $\\lambda_a$ (the SU(3) analogue of Pauli matrices): three antisymmetric components measure circular polarization, i.e. phonon angular momentum, and five symmetric components measure linear polarization. The collision integral is taken as elastic, isotropic, and polarization-independent, and the equation is expanded in spherical harmonics; the slowly relaxing $\\ell=0$ components couple to fast-decaying $\\ell=2$ multipoles, and adiabatic elimination produces the Hamiltonian-independent double-commutator formula $\\partial\\langle\\rho_q\\rangle/\\partial t \\simeq -(\\tau^{*}/\\hbar^2)\\langle[H_q,[H_q,\\langle\\rho_q\\rangle]]\\rangle$. That formula, evaluated analytically for isotropic elastic bodies and numerically for a cubic crystal, yields the relaxation rates and the linear-slower-than-circular ordering.","core_discovery":"The central discovery is a relaxation mechanism for all acoustic-phonon polarization degrees of freedom in which disorder impedes relaxation. In the absence of scattering, the polarization vector oscillates because the longitudinal–transverse splitting acts as a wavevector-dependent effective field; when the scattering rate exceeds the oscillation frequency, adiabatic elimination of the fast multipole components yields slow decay with rates $1/\\tau_{\\mathrm{AM}} = \\frac{2}{3}(c_L-c_T)^2 q^2 \\tau_2$ for circular polarization and $1/\\tau_{\\mathrm{LP}} = \\frac{2}{5}(c_L-c_T)^2 q^2 \\tau_2$ for linear polarization in isotropic elastic bodies. Because these rates are proportional to the impurity relaxation time $\\tau_2 \\propto 1/\\alpha$, increasing disorder strength lowers the rates, and the linear-polarization rates fall below the circular one. The paper validates the mechanism numerically for solid argon with Rayleigh-type mass-disorder scattering, reproducing the inverse proportionality to impurity concentration and the slower relaxation of linear polarizations.","pith_inferences":["Editorial extension: the same double-commutator structure suggests that any three-component bosonic polarization degree of freedom whose band splitting acts as a $\\mathbf{q}$-dependent field—electron orbital angular momentum, possibly magnon or optical-phonon polarizations—should exhibit disorder-protected relaxation under isotropic elastic scattering.","Editorial extension: the predicted linear-slower-than-circular hierarchy implies an experimental asymmetry test: preparing pure linear polarization states should yield longer-lived signals than circular ones, which could distinguish this mechanism from Elliott-Yafet-like direct depolarization.","Editorial extension: because polarization lifetimes grow with wavevector while heat is carried by small-$q$ phonons, the mechanism implies a spectral separation between heat transport and polarization transport that could be exploited in devices converting heat gradients into chiral-phonon currents."],"forward_implications":["Larger mass-disorder strength $\\alpha$ increases both $\\tau_{\\mathrm{AM}}$ and $\\tau_{\\mathrm{LP}}$, so isotope engineering can extend phonon polarization lifetimes without changing the phonon band structure.","Linear polarizations outlive circular ones in the same crystal, so phonon-angular-momentum probes see the shortest of the three rates while linear-polarization probes see longer decay.","Short-wavelength phonons dominate polarization transport: the rate $1/\\tau_{\\mathrm{Pol}}$ scales as $q^{-2}$ at small $q$ and $q^{-4}$ at large $q$, so injecting high-$q$ phonons should reveal the longest-lived polarizations.","The polarization diffusion length $l=\\sqrt{D\\tau_{\\mathrm{Pol}}}$ is independent of impurity concentration because $D\\propto\\tau^{*}$ and $\\tau_{\\mathrm{Pol}}\\propto 1/\\tau^{*}$, an observable signature that does not require time-resolved detection.","In solid argon with 6% xenon mass substitution, the estimated polarization relaxation times exceed the impurity relaxation time near $qa/2\\pi \\simeq 0.4$ (about 40 K), indicating a realistic regime for observing long-lived phonon polarization."],"supporting_citations":[{"why":"Dyakonov–Perel quantum kinetic equation for spin relaxation; the formalism this paper transplants to phonon polarization.","marker":"[43]"},{"why":"Dyakonov–Perel double-commutator formula for spin relaxation, which becomes Eq. (10) and yields the phonon rates (9a)-(9b).","marker":"[44]"},{"why":"Previous use of the same kinetic equation for electron orbital angular momentum; justifies the polarization-independent collision-integral approximation in Eq. (3).","marker":"[74]"},{"why":"Review contrasting Dyakonov–Perel motional narrowing with Elliott–Yafet direct spin-flip relaxation; frames why the inverse disorder dependence is the counterintuitive claim.","marker":"[45]"},{"why":"Klemens mass-defect scattering model that supplies the impurity relaxation rate used for the solid-argon estimates.","marker":"[46]"},{"why":"Rayleigh-scattering treatment from which the $\\alpha q^4$ form of $1/\\tau^{*}$ is taken for the numerical validation.","marker":"[48]"},{"why":"Interatomic force constants (Table II) for solid argon provide the phonon dispersions in the numerical calculation.","marker":"[90]"},{"why":"Gell-Mann matrices define the eight polarization parameters that reduce the polarization dynamics to coupled multipole equations.","marker":"[58]"}],"fun_headline_variants":["Motional narrowing slows phonon polarization decay","Disorder slows phonon polarization decay via motional narrowing","Impurities prolong phonon polarization lifetimes","Linear phonon polarization outlasts circular with disorder","More scattering slows phonon polarization decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that disorder scattering cares only about the phonon's energy and direction, never about which polarization or which acoustic branch it scatters; if real defects can tell the branches apart or flip polarization directly, the protective averaging can be bypassed and disorder would shorten lifetimes instead of lengthening them.","fun_headline_variants_meta":{"raw":{"variants":["Motional narrowing slows phonon polarization decay","Disorder slows phonon polarization decay via motional narrowing","Impurities prolong phonon polarization lifetimes","Linear phonon polarization outlasts circular with disorder","More scattering slows phonon polarization decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000811,"raw_usage":{"total_tokens":3516,"prompt_tokens":864,"completion_tokens":2652,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":2583}},"tokens_in":480,"tokens_out":2652,"duration_ms":17843,"temperature":1.0,"reasoning_tokens":2583,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:30:23.948854+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the circular- and linear-polarization relaxation times in a family of isotopically mass-substituted crystals (e.g., $^{36}$Ar$_{1-x}^{136}$Xe$_x$) at fixed temperature and wavevector $qa/2\\pi\\simeq0.4$, varying the substitution fraction $x$. The paper predicts $1/\\tau_{\\mathrm{AM}}$ and $1/\\tau_{\\mathrm{LP}}$ to decrease with $x$ and the diffusion lengths $\\sqrt{D\\tau}$ to stay constant; observing the opposite—rates increasing with impurity content—would falsify the mechanism.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Dyakonov–Perel quantum kinetic equation for spin relaxation; the formalism this paper transplants to phonon polarization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Dyakonov–Perel double-commutator formula for spin relaxation, which becomes Eq. (10) and yields the phonon rates (9a)-(9b)."},{"cited_title":"Kailasvuori, J","cited_arxiv_id":null,"evidence_quote":"Previous use of the same kinetic equation for electron orbital angular momentum; justifies the polarization-independent collision-integral approximation in Eq. (3)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Review contrasting Dyakonov–Perel motional narrowing with Elliott–Yafet direct spin-flip relaxation; frames why the inverse disorder dependence is the counterintuitive claim."},{"cited_title":"Tamura, Phys","cited_arxiv_id":null,"evidence_quote":"Rayleigh-scattering treatment from which the $\\alpha q^4$ form of $1/\\tau^{*}$ is taken for the numerical validation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Interatomic force constants (Table II) for solid argon provide the phonon dispersions in the numerical calculation."},{"cited_title":"However, in isotropic elas- tic bodies, the AM relaxation time τ AM remains unaf- fected by this correction term [64]","cited_arxiv_id":null,"evidence_quote":"Gell-Mann matrices define the eight polarization parameters that reduce the polarization dynamics to coupled multipole equations."}],"review_version":1}