{"id":"bf185e6f-ef3e-4c8b-a44c-eab499bef6cc","arxiv_id":"2608.11196","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"A time-dependent numerical renormalization group calculation reveals universal power-law and scaling behavior in the work distribution of quenched quantum impurity models.","lead":"This paper computes the full probability distribution of work done on two benchmark quantum impurity models, an Anderson impurity and a spin-boson model, after a sudden parameter quench. It finds a universal power-law edge at minimum work controlled by the Anderson orthogonality catastrophe, and shows the work distribution can serve as a probe of strong-coupling renormalization-group flow.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The power-law threshold and scaling-collapse claims rest entirely on the TDNRG expression Eq.","rationale":"For the central claim to hold, the TDNRG WDF of Eq. (5) must be correct and converged in the threshold region. I examined the steps available in the Letter: the work definition Eq. (1) is exact within the two-projective-measurement protocol; Eq. (3) is exact for local quenches because the environment beyond shell m is common to H_i and H_f, giving δ_{ee'}, and the environment energy cancels in the work difference E_r^{m,f}-E_q^{m,i}. Thus the reader's stated weakest assumption is partly mitigated for the overlap factorization itself. The genuinely load-bearing approximation is the NRG truncation together with the shell-energy assignment entering the phases and weights of Eq. (5): with K=860, Λ=4, and η=1/32 fixed and no convergence analysis, and with the derivation deferred to an unavailable companion paper, one cannot yet exclude that the broadened shell spectrum produces an effective power law with the claimed exponent by construction rather than by physics. The Crooks and moment checks in Fig. 3 constrain integral and bulk properties, not the threshold singularity; the observed match with OC exponents in Figs. 2 and 4 is a consistency check against the same OC physics the paper claims to test. The decisive tests are (i) the exact U=0 free-fermion WDF, which directly benchmarks Eq. (5) and the OC exponent in a noninteracting limit, and (ii) a sweep over Λ, n_keep, and η, which would expose any truncation-induced exponent shift. Because the concern is about verification rather than a demonstrated error, the verdict CONDITIONAL remains appropriate and I keep the reader's verdict unchanged.","tokens_in":17959,"tokens_out":44195,"duration_ms":385902,"concrete_test":"Compute the exact U=0 AIM work distribution for the same level quench from the free-fermion characteristic function G(u)=det[1-n_i+n_i e^{iu h_f} e^{-iu h_i}] (single-particle Hamiltonians h_i, h_f; initial Fermi occupation n_i), broaden it identically, and overlay it on the TDNRG U=0 curves of Fig. 2(a) over the full W-W_min range; if exponent and shape match within NRG resolution, Eq. (5) is validated in the noninteracting benchmark. Separately, re-run the central U=60Γ (β=1/2) and one SBM (α=0.7, Δ_0=0.05) case with Λ=2 (z-averaged) and Λ=6, and with n_keep=500 and 1400 at fixed η=1/32; if the extracted β and the T_K/Δ_r data collapse are unchanged within ~5%, the truncation concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the TDNRG work distribution, Eqs. (4)-(5), be quantitatively correct in the threshold region. Three conditions must hold: (i) the complete-basis bookkeeping, including the 'exclude all three kept' rule in the primed sum, must be correct; (ii) the NRG energy-scale-separation/truncation must be converged so that the retained K states faithfully give the overlaps S^m_{rq} and reduced density matrices eR^m for the large-amplitude quenches considered; (iii) the extracted exponent must be stable under the NRG controls Λ=4, n_keep=860, and broadening η=1/32. The Letter defers the derivation of Eq. (5) to Ref. [53] (a companion paper, not yet available), reports a single parameter set, and gives no convergence analysis. The internal checks shown (Crooks relation, first/second moment sum rules, Figs. 3(c-f)) validate global and integral properties of P(W), not the threshold exponent; they cannot distinguish a physical OC power law from an effective power law induced by truncation or broadening. The agreement of the numerical asymptotes with the OC exponents in Figs. 2 and 4 is encouraging but is a consistency check against the same OC physics being claimed, not an independent calibration. Note that Eq. (3) is exact for local quenches (the environment beyond shell m is common to H_i and H_f, giving δ_{ee'}), so the reader's stated factorization risk is mitigated at the level of peak positions; the operative approximation is the truncation plus shell-energy assignment (with the neglected shell-environment coupling), which is precisely what a convergence study must control.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the time-dependent numerical renormalization group (TDNRG) method to compute the full quantum work distribution function (WDF) following a sudden quench in two strongly coupled many-body open quantum systems: the Anderson impurity model (AIM) and the Ohmic spin-boson model (SBM). The central results are, at zero temperature, a power-law threshold P(W) ~ |W-W_min|^{-beta} as W approaches W_min, with beta = 1 - (n_f - n_i)^2/2 for the AIM and beta = 1 - alpha/2 for the SBM, and universal scaling collapse below the emergent low-energy scales T_K and Delta_r. At finite temperature, the authors verify the Crooks relation and exact first- and second-moment sum rules. The paper motivates the work as a non-perturbative, thermodynamic-limit treatment of work statistics in strongly correlated impurity models.","tokens_in":18297,"tokens_out":8664,"duration_ms":79468,"significance":"If the threshold exponents and scaling collapse are confirmed, this is a significant advance: it provides the first non-perturbative calculation of the full work distribution for strongly coupled many-body open systems in the thermodynamic limit, and it connects work statistics directly to the impurity renormalization-group flow via the Anderson orthogonality catastrophe. The exponents are not fitted: they are derived from the Friedel sum rule (AIM) and from the bosonization mapping to the interacting resonant level model (SBM), which gives independent grounding for the central predictions. The paper also includes useful internal checks (Crooks relation and exact moment sum rules) and a concrete, falsifiable scaling-collapse prediction. The main caveat is that the numerical evidence for the threshold exponent itself is less complete than the evidence for the global properties of P(W), and the derivation of the central formula is deferred to a companion paper.","major_comments":[{"comment":"The central claim that P(W) has a power-law threshold with exponent beta is not supported by a convergence analysis in the threshold region. The manuscript uses a single set of NRG controls (Lambda=4, n_kept=860, eta=1/32, stated in the End Matter), and the validation in Fig. 3 (Crooks relation and first/second moment sum rules) tests global or integral properties of P(W), which are insensitive to the local threshold exponent. Because the smooth WDF is obtained by Gaussian broadening of discrete delta peaks, the apparent power law near W_min can in principle be modified by the broadening, and the agreement with the OC prediction in Figs. 2 and 4 is a consistency check with the same OC physics being claimed rather than an independent calibration. I request a convergence study (varying Lambda, n_kept, and eta, and stating the fitting window in W-W_min) and, if possible, a comparison with the exact noninteracting U=0 WDF for the AIM, for which the threshold exponent is known analytically.","section":"TDNRG method and End Matter"},{"comment":"Equation (5), which is the foundation of all numerical results, is not derived in the Letter; the text states that the full derivation is provided in the companion paper Ref. [53], which is not yet available. In particular, the 'exclude all three kept' rule in the primed sum and the way the full-density-matrix reduced sectors are combined with the overlap matrices are nontrivial bookkeeping steps that determine the weight of each Wilson shell. For the Letter to be self-contained and checkable, this derivation (or at least the key identities leading to the primed sum) should be included in a supplemental appendix.","section":"Eq. (5)"},{"comment":"The manuscript does not report how the exponent beta is extracted from the numerical curves or how sensitive the extracted value is to the chosen fit range. The data in Fig. 2(d) show normalized curves for many beta values, but there is no quantitative measure of the deviation between the numerical slope and the predicted beta over the claimed scaling window. I ask for a quantitative extraction (for example, the logarithmic derivative of P(W) times (W-W_min)^beta) or a table of fitted exponents with fit ranges, so that the central exponent claim is testable.","section":"WDF of the quenched AIM, Fig. 2"}],"minor_comments":[{"comment":"The label 'U = 0Gamma' in Fig. 2(a) is a typographical error; it should read U=0.","section":"Fig. 2(a)"},{"comment":"The text immediately after Eq. (E.1) begins with a stray period; the equation and the following sentence should be separated cleanly.","section":"End Matter, Eq. (E.1)"},{"comment":"The expression 'sqrt(15)/4 approx 1.9' is ambiguous; it should be written as sqrt(15/4) approx 1.94 to match the intended value of n_f.","section":"Footnote 81"},{"comment":"Table E.1 contains duplicated rows for the alpha=0.7 entries; these should be removed or clearly marked as repeated rows for different bare parameters.","section":"Table E.1"},{"comment":"The x-axis labels in the upper panels of Fig. 3 appear to be mis-set; the independent variable should be W-W_min (in units of Gamma or T_K), not 'T_K Gamma'.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is likely correct in its main physical claim, but the published version needs to be more self-contained and needs quantitative evidence that the threshold exponent is converged. The companion paper Ref. [53] may provide the missing derivation and benchmarks; if so, the requested revisions could be addressed without changing the central conclusions. The scope and citation pattern are appropriate for a quantum-thermodynamics Letter."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the first non-perturbative computation of the full work distribution function for the interacting Anderson and Ohmic spin-boson models in the thermodynamic limit. That alone makes it worth a look. The headline results—a power-law threshold at the minimum work with an exponent set by the orthogonality catastrophe, and universal scaling collapse below T_K or Delta_r—are physically well motivated and not fitted. The AIM exponent comes from the Friedel sum rule, the SBM exponent from the IRLM bosonization mapping. Those are independent grounding, not circular. The numerics are checked against the Crooks relation at finite temperature and against exact first- and second-moment sum rules; those are meaningful validations of the global and integral properties of the distribution. The scaling collapse shown in Figs. 2(c,d) and 4(c,d) is compelling. The underlying TDNRG machinery is established, so the methodological leap is plausible.\n\nThe real soft spots are the ones the reader flagged. The derivation of the central expression, Eq. (5), is deferred to a companion paper, and no code or data are provided. More concerning, the threshold exponent is extracted from a single parameter set (Lambda=4, n_keep=860, eta=1/32) with no convergence analysis. The Crooks and moment checks validate integral properties, not the shape of the threshold; they cannot distinguish a true OC power law from an effective one induced by truncation or broadening. The agreement of the numerically extracted asymptotes with the predicted exponents is encouraging, but it is a consistency check against the same OC physics being claimed. The stress-test note is right that Eq. (3) is exact for local quenches, so the factorization risk is smaller than one might initially fear. The operative approximation is the truncation plus shell-energy assignment, and that is precisely what a convergence study would control. Given the paper's claim of exponentially fine resolution, the absence of even a brief convergence statement is a real gap.\n\nIf the companion paper delivers the promised derivation and the exponent is stable under NRG controls, this will be a solid result. As it stands, the threshold exponent is not yet established to the level the abstract implies. I would send this to serious peer review rather than desk reject. The questions are answerable with additional NRG data, and the physics is important enough to warrant the referee time. I would not yet cite the threshold exponent as established, but I would cite the method and the prediction with a caveat.","headline":"First non-perturbative full work distribution for strongly coupled impurity models, with independently grounded OC exponents; the main soft spot is the deferred derivation and missing convergence study, but the physics looks right.","tokens_in":18832,"tokens_out":2072,"would_cite":true,"duration_ms":19213,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.30.-d","05.70.Ln","72.15.Qm"],"model":"deepseek-v4-flash","headline":"Universal power-law threshold governs work distributions after quantum quenches.","keywords":["quantum work distribution","Anderson impurity model","spin-boson model","Anderson orthogonality catastrophe","time-dependent numerical renormalization group","quantum thermodynamics","strong coupling","work statistics"],"falsifier":"A decisive check is to compute $P(W)$ for the interacting Anderson model with an independent numerically exact method, such as real-time Monte Carlo or a process-tensor calculation, for the same level quenches, and see whether the exponent $1-(n_f-n_i)^2/2$ and the scaling collapse below $T_K$ are reproduced as $U$, $\\varepsilon_d$, and the quench size are varied; failure to reproduce the exponent or the collapse would show that the TDNRG result is an artifact of the energy-scale-separation approximation.","tokens_in":17781,"feed_emoji":"⚛️","tokens_out":9958,"duration_ms":85859,"temperature":0.7,"pith_summary":"The paper develops a non-perturbative, numerically exact method for computing the full quantum work distribution function of a strongly coupled many-body open system after a sudden parameter quench, and applies it to two paradigmatic models: the Anderson impurity model and the Ohmic spin-boson model. It claims that at zero temperature the work distribution near its minimum value $W_{\\min}$ is a universal power law, $P(W) \\sim |W-W_{\\min}|^{-\\beta}$, whose exponent is fixed by the change in low-energy scattering phase shifts between the initial and final states (the Anderson orthogonality catastrophe). In the Anderson model $\\beta = 1 - (n_f - n_i)^2/2$, where $n_i,n_f$ are the total occupations of the impurity level before and after the quench; in the spin-boson model $\\beta = 1 - \\alpha/2$, where $\\alpha$ is the Ohmic dissipation strength. Below an emergent interaction-renormalized scale, the Kondo temperature $T_K$ for the Anderson model and the renormalized tunneling amplitude $\\Delta_r$ for the spin-boson model, the distribution collapses onto a universal curve. If correct, this makes the work distribution a direct probe of the impurity renormalization-group flow and opens strongly coupled many-body systems to exact work-statistics calculations.","feed_headline":"Work distributions obey universal power law after quantum quenches","feed_subtitle":"At T=0, the Anderson and spin-boson models show a power-law edge in work, with exponents fixed by scattering phase shifts.","key_machinery":"The central object is the complete-basis time-dependent numerical renormalization group (TDNRG) representation of the characteristic function $G(u) = \\mathrm{Tr}[e^{iuH_f} e^{-iuH_i} \\rho_i]$. The bath is logarithmically discretized into a Wilson chain whose couplings decay exponentially; separate NRG runs for the initial and final Hamiltonians produce shell eigenstates, and a complete many-body basis is formed from the discarded states of all iterations. Assuming energy-scale separation, the characteristic function becomes a single sum over Wilson shells, Eq. (4), in which the overlap matrix $S^m_{rq}$ between initial and final eigenbases encodes the quench and the reduced full density matrix encodes the initial thermal weights; Fourier transformation gives the work distribution, Eq. (5). The machinery resolves the power-law threshold over many orders of magnitude because each Wilson shell contributes transitions at its characteristic energy scale, while the phase-shift difference entering through the Friedel sum rule fixes the exponent.","core_discovery":"The central claim is that the work distribution of a quenched quantum impurity model is governed at low energy by the Anderson orthogonality catastrophe. In the thermodynamic limit at $T=0$, the probability density of work has an edge singularity at the minimum work $W_{\\min} = E_0^f - E_0^i$, behaving as $P(W) \\sim |W-W_{\\min}|^{-\\beta}$ for $W \\to W_{\\min}^+$, with exponent $\\beta = 1 - \\alpha_{\\mathrm{OC}}$ where $\\alpha_{\\mathrm{OC}} = (\\Delta\\delta/\\pi)^2$ is the orthogonality exponent and $\\Delta\\delta$ is the change in scattering phase shift. For the Anderson impurity model the Friedel sum rule gives $\\Delta\\delta = \\pi(n_f - n_i)$, hence $\\beta = 1 - (n_f - n_i)^2/2$; for the spin-boson model the mapping to the interacting resonant level model gives $\\alpha_{\\mathrm{OC}} = \\alpha/2$, hence $\\beta = 1 - \\alpha/2$. The power law sets in only below the emergent low-energy scale of the final state, $T_K$ or $\\Delta_r$, and the normalized distributions collapse onto a single universal curve in that regime. At finite temperature the distributions satisfy the Crooks relation exactly, and their first and second moments match analytic sum rules.","pith_inferences":["An implication the authors do not draw is that the same orthogonality-catastrophe mechanism should control x-ray edge singularities in quantum-dot absorption spectra, making the work distribution and the optical spectrum two observable faces of the same infrared-singular response.","Near a quantum critical point such as the spin-boson localization transition at $\\alpha=1$, the threshold exponent may acquire corrections or a different scaling form; the present study is restricted to the tunneling regime $\\alpha<1$, so this is a testable extension.","The TDNRG characteristic function can in principle be differentiated to obtain full counting statistics of other conserved quantities such as charge or spin transfer, which would generalize the work-statistics result beyond energy exchange.","In finite-size realizations, the orthogonality catastrophe is cut off at a scale set by the inverse system size, so the thermodynamic-limit power law implies observable finite-size rounding of the work edge in qubit-embedding or cold-atom experiments."],"forward_implications":["The full work distribution, not just its first moments, can be computed non-perturbatively for arbitrary quench amplitudes and at zero or finite temperature, directly in the thermodynamic limit.","For any quench in the Anderson model, the threshold exponent is set by the occupation change and the same power law appears across mixed-valence and Kondo regimes, with the Kondo temperature $T_K$ controlling where the universal regime begins.","For the spin-boson model, the dissipation strength $\\alpha$ alone fixes the threshold exponent, while the renormalized tunneling amplitude $\\Delta_r$ sets the scale below which the universal power law appears.","Below $T_K$ or $\\Delta_r$, work distributions for different microscopic parameters collapse onto one universal curve, so the low-energy many-body physics is encoded in a single scale plus the orthogonality exponent.","The method extends to other impurity models, to quenches of interaction or system-bath coupling parameters, and, via sequences of quenches, to finite-time driving."],"supporting_citations":[{"why":"Supplies the original complete-basis TDNRG formulation that the work-distribution calculation is built on.","marker":"[49]"},{"why":"Provides the finite-temperature and full-density-matrix generalization of TDNRG used for the initial thermal states.","marker":"[51]"},{"why":"Introduces the Anderson infrared catastrophe that the predicted power-law threshold is attributed to.","marker":"[65]"},{"why":"Gives the exact x-ray edge singularity solution that motivates the asymptotic form of the work distribution near threshold.","marker":"[66]"},{"why":"Defines the two-projective-measurement work distribution used in Eq. (1).","marker":"[6]"},{"why":"Provides the Crooks fluctuation relation used to validate the finite-temperature work distributions.","marker":"[8]"},{"why":"Supplies the Friedel sum rule and phase-shift description of the Anderson impurity model used to express the orthogonality exponent in terms of occupations.","marker":"[44]"},{"why":"Supplies the spin-boson model, its Ohmic spectral density, and the emergent renormalized tunneling scale $\\Delta_r$.","marker":"[45]"},{"why":"Supplies the mapping from the Ohmic spin-boson model to the interacting resonant level model that the SBM calculations use.","marker":"[76]"}],"fun_headline_variants":["Quench work distribution follows universal power law","Orthogonality catastrophe sets work power law","Exact work power law from quenched open systems","Quench-induced power law in work distribution","Universal work power law below orthogonality scale"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is a clear separation of energy scales: every NRG state kept at a given step is assigned a shell energy that does not depend on the untouched part of the environment, and if that separation fails for the quenches considered, the predicted power-law threshold would be a numerical artifact rather than the physical work distribution.","fun_headline_variants_meta":{"raw":{"variants":["Quench work distribution follows universal power law","Orthogonality catastrophe sets work power law","Exact work power law from quenched open systems","Quench-induced power law in work distribution","Universal work power law below orthogonality scale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0008,"raw_usage":{"total_tokens":3576,"prompt_tokens":1062,"completion_tokens":2514,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":2459}},"tokens_in":678,"tokens_out":2514,"duration_ms":16731,"temperature":1.0,"reasoning_tokens":2459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:16:47.199736+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to compute $P(W)$ for the interacting Anderson model with an independent numerically exact method, such as real-time Monte Carlo or a process-tensor calculation, for the same level quenches, and see whether the exponent $1-(n_f-n_i)^2/2$ and the scaling collapse below $T_K$ are reproduced as $U$, $\\varepsilon_d$, and the quench size are varied; failure to reproduce the exponent or the collapse would show that the TDNRG result is an artifact of the energy-scale-separation approximation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the finite-temperature and full-density-matrix generalization of TDNRG used for the initial thermal states."},{"cited_title":"Nozi `eres and C","cited_arxiv_id":null,"evidence_quote":"Gives the exact x-ray edge singularity solution that motivates the asymptotic form of the work distribution near threshold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Friedel sum rule and phase-shift description of the Anderson impurity model used to express the orthogonality exponent in terms of occupations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the mapping from the Ohmic spin-boson model to the interacting resonant level model that the SBM calculations use."}],"review_version":1}