{"id":"16d761e2-03aa-4b94-a547-409645b0b3ad","arxiv_id":"2411.19905","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In non-Hermitian disordered lattices, wave-packet spreading follows universal laws set by the tail of the imaginary-part density of states: t/(ln t)^{1/2}, t^{1/(d+1)}, or t^{2/(d+2)}.","lead":"This paper derives universal scaling laws for how waves spread in disordered non-Hermitian systems, where random gain and loss localize all eigenstates yet propagation persists. The scaling exponents depend only on the tail shape of the imaginary-part density of states and the spatial dimension, and are testable in photonic and dissipative quantum platforms.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gaussian ImDOS tail for weak bounded disorder is cut off at W; the t/(ln t)^{1/2} law is an intermediate regime, so the claimed universal long-time scaling is not the true t→∞ asymptotic.","rationale":"The paper's optimization mechanism is coherent and the strong-disorder numerical comparisons (Figs. 2a and 3) support the uniform and linear ImDOS scaling laws in the accessible time window. The reader's conditional verdict is appropriate. My concern sharpens the reader's weakest assumption: the issue is not merely that the CLT fails for band-edge states, but that for any bounded disorder the Gaussian ImDOS has compact support, so the extreme-value tail near the upper edge is non-Gaussian. Consequently, the weak-disorder universal scaling t/(ln t)^{1/2} cannot be the true thermodynamic-limit asymptotic; it is an intermediate regime controlled by the CLT until times exponentially large in the localization volume. This does not invalidate the paper's main dynamical-optimization framework, nor its strong-disorder results, but it does require the universality claim to be qualified. The proposed test isolates the tail of the ImDOS that enters Eq. (2), so it would settle whether the objection lands. No code or error bars are provided by the paper, but they are secondary to this analytical gap.","tokens_in":20529,"tokens_out":16186,"duration_ms":172558,"concrete_test":"Importance-sample the law of λ = (1/ξ)∑_{i=1}^ξ V_i for a box-like eigenfunction with V_i uniform on [-W,W], and measure log P(λ > W - Δ) versus log Δ for Δ ≪ W. If the slope is ξ (power law) rather than the Gaussian prediction −ξ(W-Δ)^2/(2σ^2) plus a constant, substitute this measured tail into Eqs. (2)-(3) and confirm that the asymptotic exponent becomes 1/(1+d/ξ) instead of 1/2 for d-dimensional spreading. This directly settles whether the Gaussian tail is the controlling extreme-value input.","verdict_should_be":"UNCHANGED","load_bearing_attack":"For Eq. (1) with bounded Vx in [-W,W], every eigenvalue satisfies Im(E_x) = ψ†Vψ/ψ†ψ ∈ [-W,W], so the ImDOS has compact support. The weak-disorder Gaussian form (SM Eq. (19)) is a central-limit statement about typical eigenstates; it cannot describe the extreme-value tail near λ = W that Eq. (2) and Table I invoke to set the long-time scaling. For an eigenstate that effectively averages ξ bounded variables, P(λ > W - Δ) behaves as Δ^ξ, a power law, not as the Gaussian exp(-ξ(W-Δ)^2/(2σ^2)). Inserting this exact tail into Eqs. (2)-(3) gives λmax ≈ W - C r^{-d/ξ} and r ∼ t^{1/(1+d/ξ)} in the strict t→∞ limit, an exponent depending on ξ, not the universal t/(ln t)^{1/2}. Thus for any bounded disorder, the Gaussian scaling is at best an intermediate asymptotic valid for ln t ≪ ξ; the statement that the tail of ImDOS dictates universal long-time behavior is not literally correct. The SM concedes the CLT argument fails near the band edge, but the finite-support truncation affects even band-center states and is the sharper obstruction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an extreme-value optimization theory for wave-packet spreading in non-Hermitian disordered lattices with purely imaginary disorder. It argues that although eigenstates are exponentially localized, the imaginary parts of their energies provide temporal amplification, and the spreading distance is controlled by the eigenstate with the largest growth factor within a spatial volume. The resulting scaling laws are determined by the tail of the imaginary-part density of states (ImDOS): Gaussian ImDOS gives |x_c| ~ t/(ln t)^{1/2}, uniform ImDOS gives |x_c| ~ t^{1/(d+1)}, and linear ImDOS gives |x_c| ~ t^{2/(d+2)}. The theory is supported by 1D and 2D numerical simulations, and a dynamical renormalization group analysis is added to predict |x| ~ t^{3/5} for the 3D delocalized phase. A mapping to Lindblad master equations shows that the same scaling applies to open systems with random losses.","tokens_in":20839,"tokens_out":4652,"duration_ms":46922,"significance":"If correct, this work would establish a novel and potentially universal class of dynamical scaling in non-Hermitian disordered systems, connecting the tails of the non-Hermitian density of states to transport exponents and explaining experimental observations of dynamical delocalization despite spectral localization. The optimization framework is physically transparent and the numerical support for the uniform and linear ImDOS cases (Figs. 2 and 3) is reasonably convincing. The paper also makes a falsifiable prediction for 3D weak-disorder dynamics. However, the central claim that the Gaussian ImDOS tail dictates the long-time limit is not the true asymptotic for bounded disorder, which is a load-bearing issue that requires revision.","major_comments":[{"comment":"The claim that the Gaussian ImDOS tail determines the long-time scaling is undermined by the compact support of bounded disorder. For V_x in [-W,W], every eigenvalue satisfies Im(E) in [-W,W], so the ImDOS is exactly zero beyond W. The Gaussian form in SM Eq. (19) is a central-limit statement about typical eigenstates, but Eq. (2) of the main text requires the extreme-value tail near the upper edge. For an eigenstate effectively averaging xi bounded variables, P(Im(E) > W - Delta) behaves as Delta^xi, which gives lambda_max ~ W - C r^{-d/xi} and hence |x_c| ~ t^{1/(1+d/xi)} in the strict t -> infinity limit. The t/(ln t)^{1/2} law can hold only in the intermediate regime ln t << xi, not as the true asymptotic. The paper should state this crossover explicitly and either reframe the Gaussian result as an intermediate asymptotic or justify why the bounded-support tail is irrelevant in the thermodynamic limit.","section":"Universal scaling under weak disorder; SM Sec. II B"},{"comment":"The one-loop dynamical RG flow of Eq. (46) is used to predict a 3D delocalization transition and the scaling |x| ~ t^{3/5}. This claim is not supported by numerical simulation or an independent non-perturbative argument, and the derivation follows the Hermitian KPZ/directed-polymer literature without addressing the complex-valued nature of the field Phi in the non-Hermitian problem. The fixed point g*_2 = 1/(2K_d) at d=3 is presented as an exact statement, but the one-loop beta function may not be reliable for the original lattice model with purely imaginary disorder. This result should be explicitly labeled as a speculative extension, or accompanied by numerical validation.","section":"Renormalization group analysis; Eq. (46)"}],"minor_comments":[{"comment":"The notation is difficult to parse because of OCR artifacts (e.g., \"ïxð\" for the delta function, \"À\" for the localization length, \"Ã\" for the standard deviation). A careful proofreading and a notation table would improve readability.","section":"Throughout"},{"comment":"The crossover from t^{1/2} to t^{2/3} is inferred from overlapping curves for different system sizes, but no quantitative comparison with the predicted crossover scale is provided. A finite-size scaling analysis of the crossover time would strengthen the claim.","section":"Fig. 2(a)"},{"comment":"The derivation of the linear tail in the strong-disorder limit is heuristic and relies on perturbative arguments; the statement following Eq. (13) that the coefficient in front of Delta is not explicitly given makes the derivation incomplete. The numerical agreement with the predicted linear-tail scaling, however, partially compensates for this.","section":"SM Sec. II A"},{"comment":"The mapping to the Lindblad master equation is clearly presented, but the claim that the same scaling applies to purely dissipative systems (SM Sec. IV B) is stated without a derivation of the effective Hamiltonian H and its ImDOS; adding a brief derivation would make the mapping more transparent.","section":"Liouvillian dynamics"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and interesting problem, and the optimization framework is attractive. The main concern is the treatment of the Gaussian ImDOS tail for bounded disorder: the claimed long-time universal behavior is not the true asymptotic, and the paper should either restrict the claim to an intermediate regime or provide a more careful extreme-value analysis. The RG section also appears under-supported. With appropriate revisions, the paper could be suitable for publication; I recommend reconsideration after major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague —\n\nRead this one. It is the first quantitative theory I know that turns the 'jumpy' non-Hermitian transport seen in photonic experiments into concrete scaling exponents. The extreme-value optimization is simple, and the resulting Table I — Gaussian ImDOS gives t/(ln t)^{1/2}, uniform gives t^{1/(d+1)}, linear gives t^{2/(d+2)} — is genuinely new. The strong-disorder 1D and 2D numerics support the uniform/linear predictions; the crossover from t^{1/2} to t^{2/3} in 1D is visible and the size-effect discussion is careful. The RG section is speculative but clearly marked, and the Lindblad mapping is a useful bridge to experiments. The ImDOS is measured in the same systems, but since the exponents are not fitted back to dynamics, I do not see a damaging circularity.\n\nThe soft spot is the weak-disorder Gaussian tail. For bounded V in [-W,W], every eigenvalue has Im E in [-W,W]; the ImDOS has compact support. The Gaussian form from the central limit theorem describes the bulk of the distribution, not the extreme tail that the optimization actually selects. An eigenstate averaging ξ sites has a power-law tail (W - λ)^{ξ-1} near W, and using that exact tail in Eqs. (2)-(3) gives |x| ~ t^{1/(1+d/ξ)}, not t/(ln t)^{1/2}. The Gaussian scaling survives only for ln|x| << ξ, i.e. as an intermediate regime. The SM admits the CLT argument fails near the band edge, but the compact-support cutoff is sharper and applies to band-center states too. This does not kill the strong-disorder results — the linear tail there is independently derived from rare-event perturbation theory — but it does falsify the paper's 'universal weak-disorder limit' claim as a strict t→∞ statement.\n\nMinor points: no code, no error bars, and the 2D late-time linear-tail prediction is not really checked.\n\nWho should read it: anyone working on non-Hermitian disordered lattices or disorder in dissipative photonics. It deserves a serious referee, but the referee should push for a clear distinction between intermediate and true asymptotic scaling, and for a numerical test at large ξ with bounded disorder.\n\nRecommendation: send to peer review with major revision.","headline":"New and largely convincing scaling theory for non-Hermitian disorder transport, but the weak-disorder Gaussian-tail universality is an intermediate regime, not the true long-time asymptote, for bounded disorder.","tokens_in":21285,"tokens_out":3868,"would_cite":true,"duration_ms":38083,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B44","81Q12"],"pacs":[],"model":"deepseek-v4-flash","headline":"Even when disorder localizes every eigenstate, a wave packet still spreads; the spreading distance follows universal scaling set by the imaginary spectrum's tail.","keywords":["non-Hermitian Anderson localization","imaginary disorder","dynamical delocalization","imaginary-part density of states","universal scaling","sub-ballistic transport","subdiffusion","Lindblad master equation"],"falsifier":"Measure or compute exactly the density of imaginary eigenvalues near the top of the spectrum for a 1D lattice with $V_x\\in[-W,W]$ and $W\\gg t_0$; if the edge density scales as $\\Delta^p$ with $p\\neq 1$, where $\\Delta=W-\\lambda$, then the predicted late-time exponent $t^{2/(d+2)}$ is wrong. A direct dynamical test is the predicted early-to-late crossover $|x_c|\\sim t^{1/2}\\to t^{2/3}$ in that same model, which finite-size-controlled simulations or an optical-lattice experiment could confirm or reject.","tokens_in":20356,"feed_emoji":"🌊","tokens_out":11153,"duration_ms":90419,"temperature":0.7,"pith_summary":"Disorder that makes every eigenstate of a lattice exponentially localized normally freezes wave propagation; the paper's point is that this logic fails when the disorder is non-Hermitian, i.e. when the on-site potential is a random imaginary gain or loss. In that setting a wave packet can still spread, because localized eigenstates with larger imaginary energies amplify in time and come to dominate the wavefront. The paper derives a universal law for the spreading distance: it depends only on the tail of the imaginary-part density of states and on spatial dimension, not on microscopic details. Concretely, Gaussian imaginary-part density of states gives $|x_c|\\sim t/(\\ln t)^{1/2}$, a uniform one gives $|x_c|\\sim t^{1/(d+1)}$, and a linear tail gives $|x_c|\\sim t^{2/(d+2)}$. These results explain earlier experiments and predict new slower-than-ballistic and subdiffusive behaviors that have no Hermitian counterpart.","feed_headline":"Localized eigenstates don't stop the wave: new universal scaling","feed_subtitle":"Spreading distance follows a universal curve set by the imaginary spectrum's tail.","key_machinery":"The load-bearing object is the imaginary-part density of states, $\\rho(\\lambda)=\\int d\\varepsilon\\,\\rho(\\varepsilon+i\\lambda)$, the density of the imaginary parts of the complex eigenvalues; its tail, not its bulk, sets the late-time exponents. The optimization principle is the weight-factor competition $W=-|x|/\\xi+\\lambda_x t$ combined with the extreme-value condition $\\int_{\\lambda_{\\max}}^\\infty \\rho\\,d\\lambda\\sim |x|^{-d}$. For the weak-disorder limit the central limit theorem supplies a Gaussian $\\rho$; for strong bounded disorder a short perturbative argument supplies a linear tail. A separate renormalization-group calculation, based on a Cole-Hopf transformation to a nonlinear equation, determines the $d=3$ delocalized-phase exponent $3/5$.","core_discovery":"The central claim is that in the eigenstate-localized regime of a non-Hermitian disordered lattice, the average spreading distance is controlled by a simple optimization. Writing the weight of a localized eigenstate centered at $x$ as $W(x,t)=-|x|/\\xi+\\lambda_x t$, where $\\xi$ is the localization length and $\\lambda_x=\\mathrm{Im}\\,E_x$, the wavefront is dominated by the eigenstate with the largest such weight. Extreme-value statistics over the $\\sim |x|^d$ available centers gives $\\int_{\\lambda_{\\max}(|x|)}^\\infty \\rho(\\lambda)d\\lambda \\sim |x|^{-d}$, which ties the maximal imaginary part to distance through the tail of the imaginary-part density of states $\\rho(\\lambda)$. Solving $\\partial W/\\partial |x|=0$ then yields the universal scalings: $t/(\\ln t)^{1/2}$ for Gaussian $\\rho$, $t^{1/(d+1)}$ for uniform $\\rho$, and $t^{2/(d+2)}$ for a linear tail, with late-time dynamics set by the tail. The paper validates the 1D and 2D scalings numerically, shows the weak-disorder limit is universally Gaussian by the central limit theorem, and maps Lindblad open systems with random losses onto the same problem.","pith_inferences":["Editorial extension: because the spreading exponent encodes the tail shape, time-resolved measurement of $|x_c|(t)$ could serve as a practical probe of the imaginary-part density of states near the spectral edge, without resolving individual eigenvalues.","Editorial extension: the predicted crossover from the uniform bulk scaling to the linear-tail scaling, for instance $t^{1/2}\\to t^{2/3}$ in 1D, should occur at a time that grows with system size; locating this crossover numerically could discriminate the linear-tail mechanism from finite-size artifacts.","Editorial extension: the RG fixed point at $d=3$ suggests the same setup may host a disorder-driven delocalization transition; tuning disorder strength across it and measuring the $t^{3/5}$ exponent would test a prediction that lies outside the eigenstate-localized regime covered by the main optimization argument."],"forward_implications":["Even with every eigenstate exponentially localized, a non-Hermitian disordered lattice does not stop transport: the wave packet spreads indefinitely, with a distance that grows as a power of time or as $t/(\\ln t)^{1/2}$.","The spreading exponent is fixed by the tail of the imaginary-part density of states and by the spatial dimension, so different disorder distributions give different universality classes of transport.","In the strong-disorder limit with a uniform bounded potential, 1D spreading should cross over from $t^{1/2}$ to $t^{2/3}$, and 2D from $t^{1/3}$ to $t^{1/2}$, as the linear edge tail takes over.","The same dynamics appears in Lindblad open systems with random local losses, because the two-point correlation matrix obeys the same non-Hermitian evolution.","In $d=3$, weak disorder flows to a renormalization-group fixed point giving $|x|\\sim t^{3/5}$, an exponent independent of the disorder statistics."],"supporting_citations":[{"why":"Experiment showing coexistence of dynamical delocalization and spectral localization; the paper's scaling theory is designed to explain and extend this observation.","marker":"[33]"},{"why":"Provides the variable-range-hopping optimization picture, with time playing the role of inverse temperature, on which the weight-factor optimization is modeled.","marker":"[41]"},{"why":"Establishes that a random imaginary refractive index exponentially localizes eigenstates, supplying the localization length used in the weight factor.","marker":"[42]"},{"why":"Supplies the dynamical renormalization-group procedure used to obtain the flow equation and the $d=3$ fixed point.","marker":"[55]"},{"why":"Provides the extreme-value statistics relation used to justify $\\int_{\\lambda_{\\max}}^\\infty \\rho\\,d\\lambda\\sim |x|^{-d}$.","marker":"[68]"},{"why":"Classical extreme-value statistics reference that underpins the same $\\lambda_{\\max}$ scaling argument.","marker":"[69]"},{"why":"Defines the Hermitian Anderson localization phenomenon that the paper contrasts with non-Hermitian spreading.","marker":"[1]"}],"fun_headline_variants":["Universal transport laws from disordered non-Hermitian systems","Imaginary spectrum tail drives universal wave spreading","Localized eigenstates don't block transport: universal scaling","Non-Hermitian disorder: three universal transport classes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the imaginary-part density of states has exactly the tail shapes it quotes: Gaussian all the way out for weak disorder and linear at the band edge for strong bounded disorder, so a different tail shape would change every quoted exponent.","fun_headline_variants_meta":{"raw":{"variants":["Universal transport laws from disordered non-Hermitian systems","Imaginary spectrum tail drives universal wave spreading","Localized eigenstates don't block transport: universal scaling","Non-Hermitian disorder: three universal transport classes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1506,"prompt_tokens":950,"completion_tokens":556,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":502}},"tokens_in":566,"tokens_out":556,"duration_ms":5795,"temperature":1.0,"reasoning_tokens":502,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:41:51.089654+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or compute exactly the density of imaginary eigenvalues near the top of the spectrum for a 1D lattice with $V_x\\in[-W,W]$ and $W\\gg t_0$; if the edge density scales as $\\Delta^p$ with $p\\neq 1$, where $\\Delta=W-\\lambda$, then the predicted late-time exponent $t^{2/(d+2)}$ is wrong. A direct dynamical test is the predicted early-to-late crossover $|x_c|\\sim t^{1/2}\\to t^{2/3}$ in that same model, which finite-size-controlled simulations or an optical-lattice experiment could confirm or reject.","supporting_citations":[],"review_version":1}