{"id":"c3a8c348-c977-45d0-9176-50117920ce8a","arxiv_id":"2504.18754","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Nonthermal fixed points support hydrodynamic excitations described by time-dependent shear viscosity and relaxation time that both grow linearly with time for the direct energy cascade.","lead":"This paper proposes that spatial ripples around out-of-equilibrium 'nonthermal fixed points' behave like a fluid, with a time-dependent shear viscosity. It derives this from kinetic theory and finds qualitative support in QCD simulations, suggesting new ways to test far-from-equilibrium fluidity in cold atom gases.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 14-moment closure is untested for the inhomogeneous modes central to the paper's claim; absent a spectral-gap separation, the predicted shear and sound dispersions rest on an assumption the paper itself defers to future work.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing point: the 14-moment truncation around an isotropic nonthermal fixed point is not justified by a demonstrated separation of scales, and the paper's numerical tests concern only homogeneous dynamics. My stress-test pass finds no independent reason to overturn that assessment. The derivation is internally coherent as far as the stated approximations go, and the homogeneous isotropization comparison provides meaningful support for the time-dependence structure. However, the central claim is about hydrodynamic excitations of spatial inhomogeneities, and the paper itself acknowledges that testing those requires going beyond the 14-moment approximation and performing ab initio simulations with broken homogeneity. A conditional verdict is therefore the right outcome: the paper is a well-motivated proposal with partial support, not a fully validated derivation of far-from-equilibrium hydrodynamics.","tokens_in":19355,"tokens_out":30241,"duration_ms":313104,"concrete_test":"In the same QCD kinetic theory framework (or in the O(N) scalar kinetic theory), initialize a small shear perturbation with finite spatial momentum, e.g. δf(k,x3) = A f0(t0,k) (k1 k3/|k|²) e^{i q x3}, at a time t0 in the scaling regime, choosing q so that q(t0−t∗)/t_ref ≪ 1. Extract the decay of δπ13(t,q) from the full Boltzmann evolution. Compare the q-dependence and time dependence with the solution of Eq. (32) using η(t) and τπ(t) from Eq. (19), in particular the predicted late-time behavior e^{−(3η̄/4E) q² ((t−t∗)/t_ref)²}. Agreement confirms the 14-moment closure for inhomogeneous perturbations; the appearance of additional comparably slow modes or a different effective diffusion coefficient would invalidate the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that spatial inhomogeneities around a nonthermal fixed point obey the hydrodynamics of Eq. (18), with transport coefficients from Eq. (19) and dispersion relations from Eqs. (32)-(34). The load-bearing step is the 14-moment truncation of Eq. (B8), which closes the exact but unclosed moment equation (B5) by assuming a single tensor form for δfk. In the near-equilibrium setting this truncation is controlled by a separation between slow hydrodynamic modes and fast collisional relaxation. No such separation is demonstrated for a nonthermal fixed point: the background f0 is itself time-dependent, and Ref. [31] finds a tower of power-law decaying modes at q = 0 rather than an exponential gap. The 14-moment ansatz selects one relaxation channel, and nothing in the paper shows that this channel dominates inhomogeneous perturbations. The numerical support in Sec. IV is homogeneous: Fig. 2 tests Eq. (29), which concerns isotropization, not the q-dependent shear and sound modes that constitute the proposed hydrodynamic description. The paper explicitly identifies this gap in Sec. VI, listing inhomogeneous ab initio simulations as future work. Until such a test is performed, the quantitative predictions of Eq. (18) — most importantly the time-dependent diffusion law e^{-(3η̄/4E) q̄² ((t−t∗)/t_ref)²} — are plausible but unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that long-wavelength spatial perturbations around isotropic nonthermal fixed points are governed by relativistic hydrodynamic equations with intrinsically time-dependent transport coefficients. Starting from the Boltzmann equation and taking the background f_0k to be an isotropic nonthermal fixed-point distribution rather than local equilibrium, the authors use the 14-moment approximation to close the moment hierarchy and derive a relaxation equation for the shear stress, Eq. (18), with shear viscosity eta(t) and relaxation time tau_pi(t) given in Eq. (19). For the direct energy cascade these coefficients grow linearly with time, Eqs. (26), while the ratio eta/(tau_pi(E+P)) equals 1/5, Eq. (20). The paper derives scaling forms for the transport coefficients, analyses shear and sound dispersion relations in the rescaled variables of Eq. (31), and compares the predicted homogeneous isotropization and the scaling of eta and tau_pi against QCD kinetic theory simulations. It also extends the scaling analysis to the nonrelativistic direct energy cascade relevant to cold atoms.","tokens_in":19629,"tokens_out":7151,"duration_ms":75410,"significance":"If the central claim holds, the paper establishes a genuinely new bridge between nonthermal fixed points and fluid dynamics, with falsifiable consequences for cold-atom experiments and possible implications for early-universe and heavy-ion physics. The derivation is self-contained in kinetic theory: the transport coefficients are evaluated from the linearized collision kernel, the ratio in Eq. (20) follows from the definitions without fitting, and the scaling analysis in Appendix C is explicit. The paper is also appropriately conservative in several places, noting that the QCD comparison is qualitative and that fully inhomogeneous ab initio tests remain future work. The main weakness is that the load-bearing 14-moment closure is assumed rather than derived for the nonthermal background, and the numerical support in Sec. IV is limited to homogeneous observables. The reported work is a plausible and interesting proposal whose central quantitative predictions require additional validation.","major_comments":[{"comment":"The central derivation closes the exact moment equation (B5) by assuming the 14-moment form delta f_k = f_0k (2 p^mu p^nu / (15 N_4(t))) pi_mu nu for perturbations around a nonthermal fixed point. Near thermal equilibrium this truncation is controlled by a separation between slow hydrodynamic modes and rapidly relaxing nonhydrodynamic modes; the manuscript's only stated justification for its nonthermal counterpart is the isotropy of f_0k, which is necessary but not sufficient. The background is itself time dependent, and Ref. [31] is cited as finding a tower of power-law decaying modes at q=0 rather than an exponential gap, so the dominance of the single pi^mu nu channel is not established. Because Eq. (18), the transport coefficients (19), and the dispersion relations (32)-(34) all rest on this closure, the paper should either supply a test of the truncation for inhomogeneous perturbations, for example a linearized full kinetic-theory computation along the lines suggested in Sec. VI, or explicitly frame the main claim as a conjecture pending such a test.","section":"§II.B, Eq. (15)/(B8), Appendix B"},{"comment":"The ab initio support in Sec. IV does not directly test the central claim about spatially inhomogeneous hydrodynamic modes. Fig. 2 compares homogeneous isotropization with Eq. (29), which is the q=0 sector of the approximated equations, and Fig. 1 checks only the scaling of the ingredients entering Eqs. (19); neither exercise tests the q-dependent shear and sound modes that carry the paper's main prediction. Moreover, the QCD simulation includes inelastic 1-to-2 processes and a different transition matrix element than the elastic 2-to-2 scalar collision kernel used in the derivation, a point the authors acknowledge, so the comparison is qualitative. A quantitative validation of Eqs. (32)-(34) with a broken-homogeneity kinetic-theory simulation would substantially strengthen the paper.","section":"§IV, Figs. 1-2"}],"minor_comments":[{"comment":"Please verify the exponent of the transient shear mode: from the homogeneous solution (29), the transient decay should scale as ((t-t*)/t_ref)^(-t_ref/tau_bar_pi), whereas the text appears to print the inverted exponent -tau_bar_pi/t_ref.","section":"§III.B, after Eq. (32)"},{"comment":"The gain-loss terms in the collision integrals have unbalanced parentheses; for example Eq. (10) reads \"(f_p f_p' (f_k + f_k') - f_k f_k' (f_p + f_p')\" with an opening parenthesis that is never closed.","section":"Eqs. (10) and (37)"},{"comment":"The sentence \"delta C is evaluated using 2x10^9 MC samples at every 100th time step (Q dt = 0.1, or every t Q = 10 steps)\" is ambiguous; please specify the time-step size and the sampling interval precisely.","section":"Fig. 1 caption"},{"comment":"The exponent t_ref/tau_bar_pi in Eq. (29) is dimensionless, but the notation invites confusion with the dimensionful bar_tau_pi; a brief definition of the dimensions of bar_tau_pi and bar_eta would improve readability.","section":"Eq. (29) and Sec. III.A"},{"comment":"Reference [22] contains an \"https://\" URL inside the title field and would benefit from standard journal citation formatting.","section":"Reference [22]"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the hep-th and heavy-ion theory communities, and the conceptual proposal is attractive. The key issue is that the 14-moment closure is used to promote an unclosed kinetic-theory equation into a closed hydrodynamic system, but the expected separation of scales that justifies this truncation near equilibrium has not been demonstrated for the nonthermal fixed point; the paper itself defers the decisive inhomogeneous test to future work. A revised version that either provides such a test or explicitly downgrades the status of the main claim would be much easier to evaluate. The minor typographical issues are straightforward to fix."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on Berges-Denicol-Heller-Preis. The new thing is real: they formulate linearized hydrodynamics around nonthermal fixed points, with shear viscosity and relaxation time that grow linearly in time for the direct energy cascade, and they derive a time-independent ratio η/(τπ(E+P)) = 1/5. The derivation is clean and self-contained: the transport coefficients come from the linearized collision kernel, not from fitting, and the scaling analysis in Appendix C matches the known exponents. The QCD kinetic theory comparison, while only for homogeneous isotropization, does test the time-dependent τπ in a qualitatively honest way, and the finite-time corrections via t* are handled carefully. The paper is also fair about its limitations: inhomogeneous ab initio simulations are explicitly listed as future work.\n\nThe soft spot is the load-bearing one. The 14-moment truncation (Eq. B8) closes the moment hierarchy by assuming δf has a single tensor form. Near equilibrium that is justified by separation between slow hydrodynamic modes and fast collisional relaxation. Around a nonthermal fixed point no such separation is demonstrated; in fact, Ref. [31] finds a tower of power-law decaying modes at q=0, not an exponential gap. So the central predictions for shear diffusion (the e^{-... q^2 t^2} law) and sound propagation rest on an assumption that is controlled near equilibrium but not here. The numerical test in Sec. IV is homogeneous only: it validates Eq. (29) for ΔP(t), not the q-dependent modes that are the heart of the proposal. This is a genuine gap, and not a small one; the authors acknowledge it, but it means the paper should be read as a well-motivated proposal rather than a demonstrated result for inhomogeneous dynamics.\n\nI think this deserves serious peer review. The question is sharp, the derivation is transparent, and even if the truncation fails for inhomogeneous modes, the framework will likely be useful and the paper clearly advances the discussion. Would I cite it? Yes, as the derivation of time-dependent transport coefficients at nonthermal fixed points. Would I bring it to reading group? Yes, precisely because the closure problem is a good discussion topic. My recommendation: send it to referees, and make sure they focus on whether the 14-moment approximation can be justified or replaced for inhomogeneous perturbations.","headline":"A serious, transparent proposal for hydrodynamics around nonthermal fixed points; the untested 14-moment closure for inhomogeneous modes is the key soft spot.","tokens_in":20168,"tokens_out":3490,"would_cite":true,"duration_ms":32128,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Nonthermal fixed points support hydrodynamic excitations, with a shear viscosity that grows linearly in time.","keywords":["nonthermal fixed points","far-from-equilibrium hydrodynamics","shear viscosity","relaxation time","14-moment approximation","relativistic kinetic theory","direct energy cascade","QCD kinetic theory"],"falsifier":"A spatially inhomogeneous simulation of the full Boltzmann equation around the relativistic direct-energy-cascade fixed point, measuring the shear-channel response at fixed momentum $q$, would settle the claim: if the rescaled dispersion relation $\\bar\\Omega(\\bar q)$ from Eqs. (31)–(34) is not reproduced, or if the long-time shear amplitude decays exponentially rather than as $\\exp[-c q^2 (t-t_*)^2]$, the hydrodynamic description fails.","tokens_in":19106,"feed_emoji":"🌊","tokens_out":16501,"duration_ms":134203,"temperature":0.7,"pith_summary":"This paper proposes that nonthermal fixed points—the self-similar, far-from-equilibrium cascades seen in quark-gluon plasma, early-universe cosmology, and cold atomic gases—support hydrodynamic excitations. Working in relativistic kinetic theory for massless particles with binary collisions, the authors expand about an isotropic nonthermal fixed point rather than about local equilibrium and derive a closed, relaxation-type equation for shear-stress perturbations, with a shear viscosity $\\eta(t)$ and relaxation time $\\tau_\\pi(t)$ that are intrinsically time dependent. For the direct energy cascade both coefficients grow linearly in time, $\\eta(t),\\tau_\\pi(t)\\propto (t-t_*)/t_{\\rm ref}$, while the dimensionless combination $\\eta/[\\tau_\\pi(E+P)]$ stays fixed at $1/5$. If correct, this gives nonthermal fixed points transport coefficients, predicts power-law rather than exponential relaxation of pressure anisotropy, and opens the way to measuring far-from-equilibrium viscosity in cold-atom experiments.","feed_headline":"Shear viscosity grows linearly in time at nonthermal fixed points","feed_subtitle":"Spatial ripples around far-from-equilibrium attractors obey fluid equations, giving a time-dependent shear viscosity.","key_machinery":"The load-bearing object is the generalized 14-moment approximation taken around an isotropic nonthermal fixed point instead of local equilibrium: writing the one-particle distribution as $f_k=f_{0k}+\\delta f_k$ with $\\delta f_k = f_{0k}\\,(2p^\\mu p^\\nu)/(15 N_4(t))\\,\\pi_{\\mu\\nu}(t)$ and $N_4(t)=\\int dK\\,|\\mathbf{k}|^4 f_{0k}$ closes the exact but unclosed shear-stress equation of motion through a linearized collision kernel $\\delta C(t)$. That kernel determines the time-dependent transport coefficients $\\eta(t)=\\frac{8}{15^2} E N_4(t)/\\delta C(t)$ and $\\tau_\\pi(t)=\\frac{2}{15} N_4(t)/\\delta C(t)$. A scaling analysis of $\\delta C(t)$ on the self-similar fixed-point distribution then yields the linear time growth $\\eta(t),\\tau_\\pi(t)\\propto(t-t_*)/t_{\\rm ref}$ for the direct energy cascade.","core_discovery":"The paper claims that spatial inhomogeneities near a nonthermal fixed point obey the same kind of hydrodynamic equations that describe small departures from thermal equilibrium, provided the background distribution used for the expansion is the self-similar nonthermal fixed point rather than local equilibrium. Concretely, the shear-stress perturbation satisfies $\\tau_\\pi(t)\\dot{\\pi}^{\\langle\\mu\\nu\\rangle}=-\\pi^{\\mu\\nu}+2\\eta(t)\\sigma^{\\mu\\nu}+\\tau_\\pi(t)(-\\tfrac{4}{3}\\pi^{\\mu\\nu}\\theta-\\tfrac{10}{7}\\pi^{\\langle\\mu}_{\\ \\ \\lambda}\\sigma^{\\nu\\rangle\\lambda}-2\\pi^{\\langle\\mu}_{\\ \\ \\lambda}\\omega^{\\nu\\rangle\\lambda})$, with $\\eta(t)$ and $\\tau_\\pi(t)$ intrinsically time dependent. On the relativistic direct energy cascade (scaling exponents $\\alpha=-4/7$, $\\beta=-1/7$), both coefficients grow linearly with time, $\\eta(t),\\tau_\\pi(t)=\\bar\\eta,\\bar\\tau_\\pi\\,(t-t_*)/t_{\\rm ref}$, so the ratio $\\eta(t)/[\\tau_\\pi(t)(E+P)] = 1/5$ is time independent. The authors corroborate this with QCD kinetic-theory simulations showing the same scaling of the transport coefficients and a power-law, rather than exponential, relaxation of the pressure anisotropy after an anisotropic perturbation, and they show the same linear growth for the non-relativistic direct energy cascade.","pith_inferences":["If the 14-moment truncation is relaxed in the full Boltzmann equation, the exact dispersion relations may shift, but the noncommutativity of the large-time and small-momentum limits should persist, implying a finite time window of hydrodynamics at any fixed nonzero $q$; a large-volume cold-atom experiment could look for this breakdown.","The same construction should apply to the anisotropic nonthermal fixed point relevant to early-time heavy-ion collisions, potentially linking this far-from-equilibrium hydrodynamics to anisotropic-hydrodynamics descriptions used there.","A direct experimental test would be to imprint an anisotropic perturbation on a cold-atom system realizing the non-relativistic direct energy cascade and measure the pressure-anisotropy relaxation; a power-law rather than exponential decay with the predicted exponent would confirm the picture.","The time-independence of $\\eta/[\\tau_\\pi(E+P)]$, in contrast to the time dependence of $\\eta/s$, suggests that ratios of transport coefficients, rather than entropy-normalized viscosity, may be the more robust universal observables for comparing far-from-equilibrium systems across couplings."],"forward_implications":["Pressure anisotropy near the fixed point decays as a power law, $\\Delta P(t)\\sim [(t_0-t_*)/(t-t_*)]^{t_{\\rm ref}/\\bar\\tau_\\pi}$, instead of the exponential decay seen near equilibrium.","Shear perturbations at small spatial momentum $q$ become arbitrarily long-lived as $q\\to 0$, diffusing as $\\exp[- (3\\bar\\eta/4E)\\, q^2 ((t-t_*)/t_{\\rm ref})^2]$ rather than with the usual $e^{-D q^2 t}$.","Sound waves still propagate at speed $1/\\sqrt3$ around the fixed point, but their attenuation is governed by the same time-dependent transport coefficients.","The time-independent ratio $\\eta/[\\tau_\\pi(E+P)] = 1/5$ provides a dimensionless transport benchmark that connects far-from-equilibrium behavior to the near-equilibrium and strong-coupling values for the same underlying theory.","In QCD kinetic theory, the extracted $\\eta(t)\\sim t$ and $\\tau_\\pi(t)\\sim t$ scaling, together with power-law isotropization, are consistent with the fixed-point predictions even though inelastic $1\\leftrightarrow 2$ processes are present."],"supporting_citations":[{"why":"Provides the generalized 14-moment construction for far-from-equilibrium backgrounds that the derivation follows.","marker":"[46]"},{"why":"Supplies the method-of-moments equation of motion for the shear stress that is closed by the 14-moment approximation.","marker":"[56]"},{"why":"Gives the original 14-moment closure of relativistic hydrodynamics, which this work adapts to a nonthermal fixed-point background.","marker":"[50]"},{"why":"Provides the scaling exponents and QCD kinetic-theory analysis used to predict and test the linear time growth.","marker":"[47]"},{"why":"Supplies the QCD kinetic-theory simulations used as the ab initio comparison for the transport coefficients.","marker":"[15]"},{"why":"Reports power-law relaxation in the linearized spectrum around the same fixed point, consistent with the derived isotropization.","marker":"[31]"},{"why":"Establishes stability and the linearized perturbation spectrum of the nonthermal fixed point.","marker":"[30]"},{"why":"Gives the causality bound on the ratio of shear viscosity to relaxation time against which the value 1/5 is compared.","marker":"[53]"},{"why":"Provides the holographic relaxation time used to estimate the strong-coupling value of the same ratio.","marker":"[32]"}],"fun_headline_variants":["Viscosity grows linearly at nonthermal fixed points","Nonthermal fixed points get time-dependent viscosity","Time-dependent viscosity from nonthermal fixed points","Hydrodynamics with growing viscosity at nonthermal fixed points"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation hinges on assuming that the 14-moment truncation—writing the deviation from the isotropic nonthermal distribution solely in terms of the shear-stress tensor—remains an accurate closure around a nonthermal fixed point, even though no small-parameter separation of scales has been demonstrated and only the homogeneous isotropization channel, not the inhomogeneous modes, has been tested directly.","fun_headline_variants_meta":{"raw":{"variants":["Viscosity grows linearly at nonthermal fixed points","Nonthermal fixed points get time-dependent viscosity","Time-dependent viscosity from nonthermal fixed points","Hydrodynamics with growing viscosity at nonthermal fixed points"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000585,"raw_usage":{"total_tokens":2748,"prompt_tokens":943,"completion_tokens":1805,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":559,"completion_tokens_details":{"reasoning_tokens":1746}},"tokens_in":559,"tokens_out":1805,"duration_ms":15547,"temperature":1.0,"reasoning_tokens":1746,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:10:15.298490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spatially inhomogeneous simulation of the full Boltzmann equation around the relativistic direct-energy-cascade fixed point, measuring the shear-channel response at fixed momentum $q$, would settle the claim: if the rescaled dispersion relation $\\bar\\Omega(\\bar q)$ from Eqs. (31)–(34) is not reproduced, or if the long-time shear amplitude decays exponentially rather than as $\\exp[-c q^2 (t-t_*)^2]$, the hydrodynamic description fails.","supporting_citations":[],"review_version":1}