{"id":"a52e92b0-e134-43cf-a800-c562d60dacab","arxiv_id":"2607.17109","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Static lattice disorder from a single thermal ionic configuration is enough to generate effective dephasing in real-time TDDFT; the apparent damping is dominated by population changes in an unfolded primitive-cell basis, not by off-diagonal coherence decay.","lead":"The paper shows that freezing the ions into a single thermally disordered supercell configuration makes electrons in real-time simulations lose phase coherence, producing Drude-like current damping in metals and clean high-harmonic spectra in diamond without an empirical dephasing time. A semiconductor-minded reader should read it because it offers a mechanism, grounded in standard DFT, for replacing the fitted T2 damping that most ultrafast models carry by hand.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-size convergence of the static-disorder dephasing is not established; the reported ~10 fs damping could be a finite-supercell artifact.","rationale":"The reader's weakest assumption—single-shot static-disorder equivalence and missing infinite-size extrapolation—is also the most load-bearing concern. The paper is internally consistent and the mechanism is physically plausible; the ten-configuration test (Fig. S1) and the prior single-shot validation in Ref. [15] provide real supporting evidence. However, the headline claim goes beyond what is shown: 'single configuration suffices' is a statement about statistical typicality, which can only be validated by convergence in system size and/or configuration count. The 50 fs propagation window is within the phonon period scale, so the static-disorder picture is nominally valid only in that window, and finite-size recurrences may be masked by the short window. A controlled supercell-size extrapolation would distinguish genuine dephasing from finite-size phase mixing. If the extrapolation fails, the title claim overreaches; if it converges, the conditional verdict can be upgraded. I also considered the basis-dependence of the unfolding decomposition, but the more fundamental issue is whether the damping itself is thermodynamically converged. Since the reader already conditioned on this assumption, no verdict change is needed.","tokens_in":8205,"tokens_out":4779,"duration_ms":50163,"concrete_test":"Run the same SALMON setup for Al at 300 K with 5×5×5 and 6×6×6 supercells (k-point grids 9^3 and 8^3 to keep total Bloch states comparable, same displacement prescription as SM). Fit σ(t) over 0–50 fs to A e^{-t/τ}+c and report τ and c for L=2,3,4,5,6. In addition, extend the 4×4×4 run to 200 fs and inspect for recurrence or revival of the current. The static-disorder claim is supported only if τ changes by <10% from 4×4×4 to 5×5×5/6×6×6 and no recurrence appears; if τ grows with L or revivals occur, the reported dephasing is a finite-size artifact rather than the proposed microscopic mechanism. The analogous 5×5×5 check for Si and diamond would further confirm whether the single-shot typicality claim extends beyond Al.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that static lattice disorder in a single frozen configuration is sufficient to generate effective dephasing—requires that the finite-supercell calculation be representative of the thermodynamic limit. This is not established. The Al evidence (Fig. 1) shows only 2×2×2, 3×3×3, and 4×4×4 supercells; the fitted τ is quoted as 'on the order of 10 fs' but no τ(L) values or L→∞ extrapolation are reported, and the long-time value of σ(t) visibly decreases with L. For Si (Fig. S2) and diamond (Fig. 3), no size convergence or configuration statistics are shown at all. Since the supercell Hamiltonian is time-independent and finite, the current after the step perturbation is a quasiperiodic sum over discrete eigenfrequencies; apparent exponential decay on a 50 fs window could be a finite-size/level-spacing artifact rather than true dephasing. The Conclusions state that effective dephasing 'does not originate from ensemble averaging' and emerges from a single configuration; this is exactly the load-bearing assumption that needs a controlled thermodynamic-limit test.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes that static lattice disorder, represented by frozen thermal and zero-point ionic displacements in a finite supercell, is sufficient to produce effective electronic dephasing in real-time TDDFT without an empirical T2. In the linear response of Al, supercell calculations show an apparent exponential decay of the time-domain conductivity σ(t) with τ ∼ 10 fs, whereas the primitive cell does not; similar damping is reported for Si. In the nonlinear regime, a 4×4×4 diamond supercell gives clean high-harmonic spectra, while the primitive cell shows irregular post-pulse oscillations. The authors introduce a time-dependent unfolding of supercell Bloch states onto primitive-cell Bloch states, and decompose the induced current into diagonal and off-diagonal contributions (Eqs. 6–10). They conclude that the dominant contribution to Drude-like damping is the diagonal/population part of the unfolded density matrix, not decay of off-diagonal coherence, and that dephasing emerges from phase mixing in a single disordered configuration rather than from ensemble averaging.","tokens_in":8523,"tokens_out":8265,"duration_ms":118816,"significance":"The paper addresses an important open problem—how to include dephasing in rt-TDDFT—and does so with a physically motivated, parameter-free prescription: ionic displacements are fixed by phonon amplitudes and temperature, and no relaxation time is fitted. The time-dependent unfolding diagnostic is a genuine conceptual contribution that could be useful beyond this work. If the finite-size and single-shot issues are resolved, the claim that static disorder alone gives Drude damping and HHG spectral cleaning would be a significant step. The current evidence is suggestive but not yet conclusive because the central mechanism is demonstrated on finite supercells without a thermodynamic-limit analysis.","major_comments":[{"comment":"The central result is the apparent Drude-like decay in Al, but the finite-size evidence is incomplete. Only 2×2×2, 3×3×3, and 4×4×4 supercells are shown; the fit giving 'τ ∼ 10 fs' is not described (fitting window, functional form, uncertainty), and no τ(L) values or L→∞ extrapolation are reported. The long-time value of σ(t) decreases with L, so the decay on the 50 fs window may still be a finite-size/level-spacing artifact. Please provide τ(L), fit details, and a larger-size or extrapolated result; this is essential to support the claim that a single static configuration represents the thermodynamic limit.","section":"Fig. 1 and Eq. (3)"},{"comment":"The Si result is offered as corroboration, but it shows only 10 fs of dynamics and no configuration statistics. The statement that σ(t) damps on a 'similar timescale' is not supported by a fit, and the 4×4×4 curve is not shown to be converged. Without these, the silicon case does not independently strengthen the central mechanism.","section":"Fig. S2 (Si)"},{"comment":"For the HHG claim, only a primitive cell and one 4×4×4 supercell are compared. Since the supercell Hamiltonian is finite and time-independent, the suppression of spurious high-frequency oscillations could be a consequence of the discrete spectrum rather than electron-phonon dephasing in the infinite system. Please show HHG spectra for at least two supercell sizes (and ideally a second independent configuration) to demonstrate that harmonic peak widths and relative intensities are converged with L.","section":"Fig. 3 (bottom) / §3"},{"comment":"The paper uses the special-displacement prescription of Ref. [15] with 'phonon amplitudes set equal to their standard deviations.' This is not a generic thermal ensemble; it fixes the amplitude and randomizes only phases. The authors validate it for Al linear response with ten configurations (Fig. S1), but no equivalent test is provided for Si or diamond, and no comparison with an explicit ensemble average over the amplitude distribution is made. Because the dephasing rate could be sensitive to the displacement statistics, this validation is load-bearing for the single-shot claim.","section":"after Eq. (1); SM 'Computational Details'"}],"minor_comments":[{"comment":"The legend uses '1x1x1' for the primitive cell; please relabel as 'primitive (1×1×1)' to avoid confusion with a genuine supercell.","section":"Fig. 1"},{"comment":"Define F_MK explicitly (initial occupations of supercell orbitals) and state how occupations are treated during time propagation, since in the velocity gauge the time-dependent KS orbitals carry the dynamics.","section":"Eq. (7)"},{"comment":"Give the explicit expression for J_diamag, or state that it is the A(t)-dependent term; currently it is only described verbally.","section":"Eq. (9)"},{"comment":"The 'similar timescale' claim should be supported by a fit or a quantitative criterion; otherwise the wording is ambiguous.","section":"Si discussion / Fig. S2"},{"comment":"Conductivity is reported in atomic units without a conversion factor; adding one (e.g., 1 a.u. of conductivity) would improve readability for a broader audience.","section":"Figures"}],"recommendation":"major_revision","confidential_remarks":"This is a potentially important computational paper, but the central claim is not yet fully established because the finite-size analysis is incomplete. I see no fundamental flaw in the conceptual framework; the main risk is that the reported dephasing is a finite-supercell artifact. I recommend requiring the additional convergence and configuration-statistics tests in revision. The time-dependent unfolding analysis is interesting and likely robust, but the paper should also state the phase convention used in the projection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a careful read. The genuinely new piece is the time-dependent unfolding analysis: mapping the supercell wavefunction onto a primitive-cell density matrix and decomposing the current into diagonal and off-diagonal contributions. The claim that Drude-like damping in Al comes mostly from diagonal population dynamics, not off-diagonal coherence, is a real diagnostic advance and it is not in the cited prior work. The broader package—single-shot static disorder, Williams–Lax, no empirical T2—is a believable and useful explanation for why rt-TDDFT supercells show dephasing, and the Al size trend in Fig. 1 is suggestive in the right direction.\n\nThe soft spots are real and concentrated where the stress-test says they are. The 10 fs timescale is read off from three supercell sizes with no fit range, no uncertainty, and no L→∞ extrapolation. Since the supercell Hamiltonian is time-independent and finite, the current is a quasiperiodic sum over discrete eigenfrequencies, and an apparent exponential decay on a 50 fs window can be a level-spacing or beating artifact. The monotonic decrease of the long-time asymptotic value with supercell size argues that the damping is physical, but it does not quantify convergence. Si and diamond get no configuration statistics and no size convergence at all; the diamond HHG is a single 4×4×4 run. The 10-configuration test in Fig. S1 only shows shot-to-shot reproducibility for Al, not thermodynamic-limit behavior.\n\nTwo further caveats, both minor. The diagonal/off-diagonal decomposition is basis-dependent: it uses the conventional cubic primitive cell as the unfolding reference, and a different unfolding choice could shift weight between the two terms. So \"populations dominate\" should be read as a useful diagnostic in a chosen basis, not a basis-invariant microscopic fact. And the title overreaches a bit: the paper isolates static disorder and explicitly says phonon dynamics may add more at long times, but the title promises the microscopic origin of dephasing in solids. The HHG spectral cleaning itself was already reported in Ref. [22]; what is new here is the Williams–Lax attribution and the connection to the linear-response dephasing.\n\nThe authors are honest about the unitary, single-configuration nature of the calculation, and the conclusions do not oversell beyond what the numerics show internally. The missing piece is a controlled thermodynamic-limit test, which would make the quantitative timescale solid. As it stands, I read the paper as a strong qualitative mechanism plus a promising diagnostic, with the exact tau still open.\n\nI would send this to peer review. The diagnostic and the framework deserve referee time, and the convergence question is exactly what a good referee should push on. I would not cite the 10 fs number as a material property, but I would cite the unfolding method for dephasing analysis.","headline":"A useful, clever rt-TDDFT dephasing mechanism with a genuinely new diagnostic, but the headline 10 fs is a finite-size number until proven otherwise.","tokens_in":8933,"tokens_out":1692,"would_cite":true,"duration_ms":19137,"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 claims that static lattice disorder from thermal and zero-point motion, captured through a single frozen ionic configuration in a supercell, is enough to produce effective dephasing in real-time TDDFT—yielding Drude-like current","keywords":["dephasing","electron-phonon interaction","real-time TDDFT","Williams-Lax formalism","Drude damping","high-harmonic generation","band unfolding","thermal disorder"],"falsifier":"A concrete falsifier would be to compute the diamond HHG spectrum (or the silicon conductivity) using an explicit ensemble average over many (e.g., 20 or more) independently sampled ionic configurations and compare it to the single-shot result. If the single-shot spectrum is not statistically representative—for instance, if the clean harmonic peaks wash out or the extracted damping time shifts by more than a few femtoseconds under averaging—the single-configuration typicality claim fails. Additionally, running the aluminum conductivity at a 5x5x5 supercell and checking whether the ~10 fs decay","tokens_in":8118,"feed_emoji":"⚛️","tokens_out":2514,"duration_ms":29199,"temperature":0.7,"pith_summary":"The authors aim to establish that electron-phonon interactions, treated as frozen static disorder in a supercell, provide a first-principles mechanism for dephasing in solids within real-time time-dependent density functional theory. They argue that a single thermally disordered ionic configuration, without ensemble averaging, reproduces Drude-like exponential current damping in aluminum and produces well-resolved high-harmonic spectra in diamond that otherwise require an ultrashort phenomenological dephasing time. By unfolding the supercell dynamics onto a primitive-cell density matrix, they show that the dominant contribution to this dephasing comes from diagonal population dynamics, not the decay of off-diagonal coherence assumed in conventional T2 pictures. If correct, this removes the need for empirical relaxation parameters in real-time simulations and provides a microscopic interpretation of dephasing as a coarse-graining effect of unitary evolution in a disordered environment.","feed_headline":"Frozen lattice disorder yields Drude damping and clean harmonics","feed_subtitle":"One thermal snapshot in supercell real-time TDDFT gives first-principles dephasing in aluminum and diamond, no empirical T2 needed.","key_machinery":"The key machinery is the time-domain Williams-Lax framework: ionic positions are frozen in a supercell with displacements sampled from harmonic phonon modes (including zero-point motion), and the time-dependent Kohn-Sham equation is solved for that single configuration. To interpret the dephasing, the authors introduce a time-dependent unfolding procedure that projects supercell wavefunctions onto primitive-cell Bloch states, yielding a primitive-cell density matrix ρ_nn'k(t). This unfolding acts as a spatial average over translational replicas and is what reveals the diagonal population dynamics as the dominant source of the apparent damping.","core_discovery":"The central discovery is that a single frozen ionic displacement pattern, generated from thermal and zero-point phonon amplitudes set to their standard deviations, is sufficient to induce realistic dephasing in real-time TDDFT. In aluminum, supercell calculations with this static disorder produce exponential decay of the time-domain conductivity with a characteristic time of about 10 fs, matching the expected Drude form. In diamond, the same approach suppresses spurious post-pulse current oscillations and yields clean high-harmonic peaks beyond the 20th order without introducing a T2 parameter or propagation effects. A time-dependent unfolding analysis maps the supercell density matrix onto","pith_inferences":["An explicit test of the single-shot equivalence would be to average the diamond HHG spectrum over many independently sampled ionic configurations and compare with the single-shot result; the paper only validates this robustness for aluminum linear response, not for nonlinear or dielectric cases.","The reported ~10 fs relaxation time for aluminum is likely supercell-size dependent because the authors do not provide an infinite-size extrapolation; a 5x5x5 or larger supercell calculation would indicate whether the damping time converges or keeps shifting with cell size.","The dominance of diagonal population dynamics suggests that a quantum master equation formulated in the primitive-cell basis might reproduce the same dephasing with a local-in-time dissipator, but this connection is not made in the paper and remains an inference.","The framework could be extended to finite-temperature ionic dynamics or quantum nuclear effects to test whether static disorder remains sufficient at longer time scales (beyond a few tens of femtoseconds) where phonon dynamics may contribute."],"forward_implications":["Real-time TDDFT can predict dephasing times and spectral broadening from first principles, without empirical T2 parameters or configuration averaging.","High-harmonic spectra in dielectrics can be computed reliably using supercell rt-TDDFT with thermal disorder, potentially eliminating the need for macroscopic propagation effects to explain spectral cleanliness.","The primitive-cell density-matrix picture suggests that phenomenological dephasing models should include population relaxation effects, not just off-diagonal coherence decay.","The single-shot result implies that statistical typicality may hold: one representative disordered configuration can stand in for a full thermal ensemble in ultrafast electron dynamics.","The same framework should apply to other response functions, such as phonon-assisted absorption and nonlinear photocurrents, where electron-phonon coupling is relevant."],"fun_headline_variants":["One frozen lattice snapshot mimics real dephasing in solids","Static disorder gives Drude damping and clean harmonics in TDDFT","Single thermal displacement yields first-principles dephasing","Dephasing without T2: frozen phonons in supercell TDDFT","Population decay, not coherence loss, drives solid-state dephasing"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that a single frozen snapshot of ionic positions, with phonon displacements set to their standard deviations, reproduces the dephasing that a full thermal ensemble or explicit phonon dynamics would produce—an assumption the authors validate only for aluminum's linear conductivity with ten configurations, not for silicon or diamond.","fun_headline_variants_meta":{"raw":{"variants":["One frozen lattice snapshot mimics real dephasing in solids","Static disorder gives Drude damping and clean harmonics in TDDFT","Single thermal displacement yields first-principles dephasing","Dephasing without T2: frozen phonons in supercell TDDFT","Population decay, not coherence loss, drives solid-state dephasing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":9.4e-05,"raw_usage":{"total_tokens":764,"prompt_tokens":607,"completion_tokens":157,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":351,"completion_tokens_details":{"reasoning_tokens":68}},"tokens_in":351,"tokens_out":157,"duration_ms":2486,"temperature":1.0,"reasoning_tokens":68,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T18:59:04.740203+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete falsifier would be to compute the diamond HHG spectrum (or the silicon conductivity) using an explicit ensemble average over many (e.g., 20 or more) independently sampled ionic configurations and compare it to the single-shot result. If the single-shot spectrum is not statistically representative—for instance, if the clean harmonic peaks wash out or the extracted damping time shifts by more than a few femtoseconds under averaging—the single-configuration typicality claim fails. Additionally, running the aluminum conductivity at a 5x5x5 supercell and checking whether the ~10 fs decay","supporting_citations":[],"review_version":1}