{"id":"98782f6d-8724-48bd-b648-397855f61257","arxiv_id":"1908.04314","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For strong power-law disorder and strong bath mismatch, the thermal conductance of a one-dimensional harmonic chain scales as L^{-1}, and the conductance follows a universal curve in terms of the ratio of the localization cutoff frequency to temperature.","lead":"This paper calculates how much heat flows along a one-dimensional chain of masses connected by springs whose stiffness varies randomly. It finds that for very uneven spring strengths and strong coupling to heat baths, the heat flow obeys the usual inverse-length law, and that all cases collapse onto one universal curve.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Fourier-law prediction for strong power-law disorder with impedance mismatch rests on an unverified assumption that each transmission resonance has a frequency-independent area; tests for the ε<1 regime are missing.","rationale":"The reader identified the frequency-independence of the resonance area as the weakest assumption, and this is indeed the most load-bearing point for the paper's central novel claim: Fourier's law for strong power-law disorder under strong impedance mismatch. The entire derivation of β=1 for ε≤1 with k≪1 or k≫1 depends on Eq. (4), which requires the per-resonance area to be ω-independent. For strong coupling, this property is derived from a first-order perturbative scaling of the boundary amplitudes that is verified numerically only for uniform disorder (W=0.5), not for power-law disorder with ε<1, where the density of states diverges and rare soft springs can dominate the low-frequency mode structure. For weak coupling, the area independence is imported from Ref. [21] without a derivation for the power-law case. Because this assumption is directly testable and the numerical evidence in the paper does not cover the crucial regime, the concern is substantive. However, it is not a demonstrated internal inconsistency; it is a gap that a targeted numerical check would resolve. The reader's CONDITIONAL verdict already accounts for this kind of unresolved assumption, so I recommend no change to the verdict: UNCHANGED. I agree with the reader that the k≈1 exponents rely on a crude step-function approximation and that the ε=1 DOS divergence is unresolved; these are secondary to the Lorentzian-area issue but reinforce the need for conditions rather than full acceptance.","tokens_in":18046,"tokens_out":8041,"duration_ms":84857,"concrete_test":"For a power-law disordered chain with ε=0.5 (strong disorder), N=2000, and k=100, compute the transmission coefficient τ(ω) via the transfer matrix (or by diagonalizing the fixed-boundary dynamical matrix) for many disorder realizations. For every resonance with eigenfrequency ω_i < ω_L (where ξ(ω_L)=L), fit τ(ω) to the Lorentzian form of Eq. (S.14) and extract the area γ_i. Test whether the fitted areas are consistent with a constant independent of ω_i, or whether they follow γ_i∼ω_i^δ with δ significantly different from zero. A second, equivalent check is to compute the low-frequency eigenmodes of the fixed-boundary chain and fit |v1| and |vN| to A ω^p/N^q; the claim requires p=1. If δ≠0 or p≠1, Eq. (4) is invalid and the β=1 prediction for strong power-law disorder must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The prediction β=1 for ε≤1 and strong impedance mismatch follows from Eq. (4), G≈Σ∫_0^{ω_L}ρ(ω)dω, together with the assumption that each transmission resonance has a frequency-independent area Σ. For k≫1, Supp. Sec. I derives γ≈2π√(km)/(MN) using the first-order perturbative boundary-amplitude scaling |v1|,|vN|∼ω/(k√N) (Eq. S.24). That scaling is obtained in Supp. Sec. II by treating Ad as a perturbation and using ordered-chain sinusoidal modes (Eq. S.23) as zeroth-order eigenmodes. The paper asserts without verification that low-frequency modes of a strongly disordered power-law chain follow approximately the same dispersion relation as the ordered chain. For ε<1, where ρ(ω)∼ω^{(ε-1)/(ε+1)} diverges and soft springs are abundant, the zeroth-order modes are not guaranteed to be sinusoidal; boundary amplitudes can be controlled by rare soft springs near the ends, changing the power p in |v1|∼ω^p/N^q. If p≠1, the integrated area scales as γ∼ω^{2p-2}, and Eq. (5) becomes G∝ω_L^{s+2p-1}, so the exact cancellation yielding β=1 is lost. The weak-coupling case (k≪1) is equally exposed: the area independence is imported from Ref. [21] without a derivation for power-law disorder. No numerical check of the frequency dependence of the resonance area for the ε<1 power-law case is provided; Fig. S3 reports only W=0.5 uniform disorder. Thus the cornerstone of the Fourier-law claim is an assumption that has not been tested in the regime where it matters most.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies steady-state heat transport through one-dimensional disordered harmonic chains connected to ordered harmonic reservoirs. The authors start from the Landauer formula and combine scaling laws for the localization length ξ(ω) and the density of states ρ(ω) to derive the scaling of the thermal conductance with the chain length L. For near impedance matching (k≈1) they approximate τ(ω)=1 below a disorder cut-off ω_L defined by ξ(ω_L)=L, obtaining G∝L^{−1/α}; for strong impedance mismatch (k≪1 or k≫1) they treat transmission as non-overlapping Lorentzians with frequency-independent area, obtaining G∝L^{−(s+1)/α}. The paper then uses published and partially re-derived results for uniform and power-law disorder: ρ(ω)∼Dω^s and ξ(ω)∼ω^{−α}. The combination (s+1)/α=1 for power-law disorder with ε≤1 gives a Fourier-law prediction G∝L^{−1}. The results are collected in Table I, and numerics for τ, ρ, ξ, and G are presented, along with a claimed universal collapse in terms of ω_L/ω_T (k≈1) and ω_L^{s+1}/ω_T (impedance mismatch).","tokens_in":18419,"tokens_out":8844,"duration_ms":87535,"significance":"Proving the existence of a concrete parameter regime where Fourier's law holds in a one-dimensional disordered harmonic chain is of considerable interest, given the long history of anomalous transport in low-dimensional systems. If the central claim is correct, the paper also offers a compact unifying scaling variable and predictions that can be tested in numerical simulations. The analytical structure is transparent, and the numerical work for the uniform-disorder and ε>1 power-law cases corroborates the predicted exponents. A particular strength of the paper is that the asymptotic exponents are not obtained by fitting; they are derived from independent scaling inputs and then compared with numerics. My reservations concern the untested assumptions in the ε≤1 strong-disorder, strong-impedance-mismatch regime, which is exactly the case used to claim Fourier's law.","major_comments":[{"comment":"The derivation of the strong-coupling Lorentzian area (Eq. (S.15)), with |v1|,|vN|∼ω/(k√N) from Eq. (S.24), assumes that the low-frequency modes of the disordered chain are close to ordered-chain sinusoidal modes. For power-law disorder with ε<1, where ρ(ω) diverges at zero frequency and soft springs are abundant, this assumption is not self-evident; rare soft springs near the boundaries could change the exponent in |v1|∼ω^p/N^q. If p≠1, the Lorentzian area scales as γ∼ω^{2p−2}, and Eq. (5) becomes G∝ω_L^{s+2p−1}, destroying the exact cancellation that yields β=1. The numerical evidence cited for Eq. (S.24), namely Figs. S.2 and S.3, is for uniform disorder only (W=0.5 and W=1.95), and no frequency-dependence check of the resonance area or of v1(ω) is provided for the ε<1 power-law case. I ask the authors to either prove Eq. (S.24) for power-law disorder or supply a dedicated numerical test in this regime.","section":"Supplemental Material, Sec. II and Eq. (S.24); main text Eqs. (4)-(5)"},{"comment":"For weak coupling k≪1, the frequency independence of the individual Lorentzian areas is imported from Ref. [21] without a derivation for the power-law disorder model. The soft-spring dominance in the ε<1 regime is precisely the case in which the delocalized-mode structure underlying Ref. [21] is most questionable. The authors should state the precise conditions under which the Ref. [21] result carries over, or verify numerically that Σ(ω) is flat for ε<1. Until this is supplied, the Fourier-law row of Table I rests on an untested premise.","section":"Main text, Eqs. (4)-(5) and Table I"},{"comment":"The paper reports an unexplained weak divergence of ρ(ω) at ε=1 in a regime where theory predicts s=0. Since ε=1 lies on the boundary of the strong-disorder row used for the β=1 prediction, this discrepancy should be addressed explicitly, for example by a finite-size scaling study or by sharpening the regime boundaries. If the divergence persists in the thermodynamic limit, the exponent s used in Eq. (5) at ε=1 is not the correct input to the central derivation.","section":"Fig. 2(b) and the paragraph near Fig. 2"}],"minor_comments":[{"comment":"There is a typo: 'condactuance' should be 'conductance'; similarly, footnote [38] reads 'infinte' for 'infinite'. The manuscript would benefit from a careful proofread.","section":"Supplemental Material, Sec. I"},{"comment":"The prefactor in Eq. (2) is easy to misread because g_q=π²T/(3h) already contains T and h; please state explicitly that the expression is written in natural units k_B=1 and ħ=1 for the thermal frequency, so that the prefactor 3/π² is unambiguous.","section":"Main text, Eq. (2)"},{"comment":"The caption introduces G̃qm=(s+1)Gqm/D, but the symbols G̃qm and D are not defined in the main text before the figure; please define them in the main text near Eq. (5) so that the collapse shown in panels (c) and (d) is self-contained.","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is interesting and the proposed Fourier-law mechanism is plausible, but the central ε≤1, strong-impedance-mismatch claim is currently supported by an assumption that has only been tested for uniform disorder. This gap is fixable with additional numerical checks, so major revision rather than rejection seems appropriate. The authors should also clarify the provenance of the scaling inputs they quote from refs. [21,29,30], since those inputs are load-bearing for the derived exponents."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: This paper has one genuinely new idea — a scaling table for thermal conductance in 1D disordered harmonic chains with different disorder and bath coupling, plus a dimensionless scaling variable that collapses data across models — and one central soft spot: the claim that strong power-law disorder plus strong impedance mismatch restores Fourier's law depends on an assumption about transmission resonance areas that is not tested in the regime where it matters most.\n\nWhat's good: The paper takes known scaling for localization length and DOS (mostly from Ziman and from the authors' earlier work) and turns it into a systematic prediction for conductance exponents. The k≈1 case is straightforward but useful. The strong-coupling limit analysis in the supplement is nontrivial: they derive the Lorentzian area and width for transmission and show the area is frequency-independent for delocalized modes using degenerate perturbation theory. The universal collapse in terms of ω_L/ω_T or ω_L^{s+1}/ω_T is a nice organizing principle, and the numerics in Fig. 4 are extensive across disorder types and couplings. The paper is honest: it flags the unexplained DOS divergence at ε=1.\n\nSoft spots: The stress-test is right. The strong-coupling derivation assumes the zeroth-order eigenmodes of a strongly disordered chain are sinusoidal, which is exactly what fails for ε<1 where soft springs produce rare weak links. They verify the amplitude scaling |v1|∼ω/(k√N) only for uniform disorder W=0.5 (Fig. S3), not for power-law ε<1. If the area became frequency-dependent, the cancellation giving β=1 would be lost. The weak-coupling side is also imported from ref [21] rather than re-derived for power-law disorder. That said, the numerical G(L) data in Fig. 4(b) do show β≈1 for ε≤1 and k=0.01, which is an indirect check; but no error bars are reported, so we can't judge how clean the exponent is. The k≈1 exponents rely on a crude step-function transmission approximation, acknowledged by the authors; that's probably fine for exponents but worth stating more carefully. The citation pattern is fine; the reliance on refs [21,33] is explicit and appropriate.\n\nBottom line: This is for the low-dimensional heat-transport community; the scaling table and collapse variable are worth engaging with. The Fourier-law claim is plausible but not fully pinned down; a referee should ask for a direct numerical check of the resonance-area frequency independence for ε<1 power-law disorder, and for error bars on exponents. This deserves peer review.","headline":"A useful scaling framework for 1D disordered heat transport, but the headline Fourier-law claim for strong power-law disorder rests on an untested resonance-area assumption that needs a direct numerical check.","tokens_in":18926,"tokens_out":5258,"would_cite":true,"duration_ms":47952,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["44.10.+i","05.60.-k","63.20.Pw"],"model":"deepseek-v4-flash","headline":"One-dimensional disordered harmonic chains transport heat as G∝L^{-β}; strong heavy-tailed disorder with impedance mismatch gives β=1 (Fourier's law).","keywords":["thermal conductance","disordered harmonic chains","localization length","Fourier's law","power-law disorder","Landauer formula","phonon transport","universal scaling"],"falsifier":"For a chain with power-law disorder ϵ=0.5 and coupling k=0.01, compute the transmission peaks below ω_L and measure the area γ of each peak as a function of its center frequency ω. If γ∝ω^δ with δ≠0, the scaling becomes G∝$L^{{-(s+1+δ)/α}}$, contradicting the paper's claimed β=1 for strong disorder.","tokens_in":1730,"feed_emoji":"🔥","tokens_out":3482,"duration_ms":76241,"temperature":0.7,"pith_summary":"This paper works out how the thermal conductance G of a one-dimensional disordered harmonic chain shrinks as the chain length L grows. The central claim is that the exponent β in G∝$L^{{-β}}$ is fixed by two low-frequency scaling laws: the localization length ξ(ω)∝$ω^{{-α}}$ and the phonon density of states ρ(ω)∝$ω^{{s}}$. With a matched bath coupling the conductance counts all delocalized modes and gives β=1/α; with a large impedance mismatch it counts resonant modes and gives β=(s+1)/α. The striking consequence is that for power-law disorder with ϵ≤1, where the mean compressibility diverges, s and α cancel to leave β=1, i.e. Fourier's law in a one-dimensional disordered system. For uniform disorder the exponents give β=1/2, independent of the bath coupling.","feed_headline":"Strong disorder restores Fourier's law in 1D heat chains","feed_subtitle":"G scales as 1/L only when heavy-tailed disorder meets impedance mismatch; otherwise transport remains anomalous.","key_machinery":"The cutoff frequency ω_L, defined by the equality ξ(ω_L)=L, is the central object: it separates delocalized transmitting phonons from localized insulating ones. The argument uses the Landauer formula for phonons and replaces the full transmission coefficient by a step function below ω_L for k≈1 (yielding G∝ω_L) or by Lorentzian peaks of equal area for k≪1 and k≫1 (yielding G∝∫$_0^{{ω_L}}$ρ(ω)dω). The equal-area property for weak coupling is imported from an earlier work of the same group; for strong coupling it is derived in the supplementary material from first-order perturbation theory about a clean chain with fixed boundaries, where the end amplitudes of low-frequency modes scale as |v_1|,|v_N|∼ω/(k√N).","core_discovery":"For frequencies ω below a cutoff ω_L defined by ξ(ω_L)=L, phonons are delocalized and carry heat; above it they are localized and do not. Taking the transmission to be roughly one for ω<ω_L in the impedance-matched case k≈1, the Landauer integral gives G∝ω_L/ω_T, hence β=1/α because ω_L∝$L^{{-1/α}}$. In the strong-impedance-mismatch case k≪1 or k≫1, transmission consists of well-resolved Lorentzians centered on the chain's eigenfrequencies, and the integrated area of each resonance is frequency-independent (derived for strong coupling from the boundary-amplitude scaling |v_1|,|v_N|∼ω/(k√N)). The conductance then reduces to counting modes below ω_L, G∝∫$_0^{{ω_L}}$ dω ρ(ω)∝$ω_L^{{s+1}}$∝$L^{{-(s+1)/α}}$. The paper demonstrates numerically that for uniform disorder s=0, α=2, giving β=1/2; for power-law disorder with 1<ϵ≤2, α=ϵ, giving β=1/ϵ; and for ϵ≤1, s=(ϵ-1)/(ϵ+1) and α=2ϵ/(1+ϵ), whose combination yields s+1=α and therefore β=1, normal heat conduction.","pith_inferences":["A direct corollary, not stated in the paper, is that any disorder model whose exponents satisfy s+1=α will show Fourier-like β=1 even without impedance mismatch; strong heavy-tailed disorder is the paper's concrete realization.","The universal variable ω_L^{s+1}/ω_T suggests a practical experimental way to extract both α and s from a single conductance-versus-temperature curve, avoiding separate localization-length measurements.","The derivation's Lorentzian-area assumption could be tested directly by computing the area of each transmission peak below ω_L; if the area acquires a frequency dependence δ, the exponent becomes β=(s+1+δ)/α rather than (s+1)/α."],"forward_implications":["Uniform mass or spring disorder gives G∝L^{-1/2} for all bath couplings, confirming and extending earlier free-boundary results.","Power-law disorder with 1<ϵ≤2 gives an anomalous exponent β=1/ϵ, tunable by the disorder strength.","Strong heavy-tailed disorder (ϵ≤1) together with strong impedance mismatch (k≪1 or k≫1) satisfies Fourier's law, G∝L^{-1}.","All disorder and temperature data collapse onto one universal curve when the conductance is plotted against ω_L/ω_T for k≈1 or against ω_L^{s+1}/ω_T for k≪1 and k≫1.","At fixed length L, the conductance saturates at the quantum of thermal conductance for T≪ℏω_L and then decays as T^{-1} for T≫ℏω_L."],"supporting_citations":[{"why":"Supplies the weak-coupling result that transmission resonances of delocalized modes have equal, frequency-independent area; the key input for the mode-counting formula.","marker":"[21]"},{"why":"Provides the localization-length and density-of-states scalings for power-law disorder summarized in Table I.","marker":"[33]"},{"why":"Reported β=1/2 for mass-disordered chains with free boundary conditions, a baseline the paper extends.","marker":"[13]"},{"why":"Established the dependence of the scaling on the heat bath and defined the free-boundary setup used here.","marker":"[10]"},{"why":"Reported the fixed-boundary exponent 3/2, which the paper contrasts with its free-boundary results.","marker":"[12]"},{"why":"Provides the Landauer formula for phonons and the quantum of thermal conductance used for normalization.","marker":"[34]"},{"why":"Provides the asymptotic ξ∝ω^{-2} scaling for uniform disorder used in both matched and mismatched cases.","marker":"[27]"}],"fun_headline_variants":["Disorder restores Fourier's law in 1D heat chains","Heavy-tailed disorder makes 1D heat conduction normal","Impedance mismatch plus disorder yields Fourier's law","Universal scaling emerges in disordered 1D chains"],"cache_read_input_tokens":20992,"weakest_assumption_plain":"The derivation assumes that every transmission resonance below the cutoff ω_L contributes a frequency-independent integrated area to the conductance; if the resonance area varies with frequency, the mode-counting step G∝∫ρ(ω)dω and all derived exponents, including Fourier's β=1, collapse.","fun_headline_variants_meta":{"raw":{"variants":["Disorder restores Fourier's law in 1D heat chains","Heavy-tailed disorder makes 1D heat conduction normal","Impedance mismatch plus disorder yields Fourier's law","Universal scaling emerges in disordered 1D chains"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1337,"prompt_tokens":941,"completion_tokens":396,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":330}},"tokens_in":557,"tokens_out":396,"duration_ms":4097,"temperature":1.0,"reasoning_tokens":330,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:46:35.677248+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a chain with power-law disorder ϵ=0.5 and coupling k=0.01, compute the transmission peaks below ω_L and measure the area γ of each peak as a function of its center frequency ω. If γ∝ω^δ with δ≠0, the scaling becomes G∝$L^{{-(s+1+δ)/α}}$, contradicting the paper's claimed β=1 for strong disorder.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the weak-coupling result that transmission resonances of delocalized modes have equal, frequency-independent area; the key input for the mode-counting formula."},{"cited_title":"Alexander, J","cited_arxiv_id":null,"evidence_quote":"Provides the localization-length and density-of-states scalings for power-law disorder summarized in Table I."},{"cited_title":"Casher and J","cited_arxiv_id":null,"evidence_quote":"Reported β=1/2 for mass-disordered chains with free boundary conditions, a baseline the paper extends."},{"cited_title":"Mai and O","cited_arxiv_id":null,"evidence_quote":"Established the dependence of the scaling on the heat bath and defined the free-boundary setup used here."},{"cited_title":"Moghaddasi Fereidani and D","cited_arxiv_id":null,"evidence_quote":"Reported the fixed-boundary exponent 3/2, which the paper contrasts with its free-boundary results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Landauer formula for phonons and the quantum of thermal conductance used for normalization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the asymptotic ξ∝ω^{-2} scaling for uniform disorder used in both matched and mismatched cases."}],"review_version":1}