{"id":"e95385bb-1d3e-4d3c-a9a5-6d62106fdad2","arxiv_id":"2412.16887","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A single local impurity in an integrable lattice produces a singular Lyapunov spectrum with exponentially localized leading Lyapunov vectors, while leaving KPZ spin transport unchanged.","lead":"In a discrete space-time spin lattice with a single-site magnetic impurity, the authors find that most Lyapunov exponents vanish as system size grows, and the few surviving ones correspond to exponentially localized chaotic modes. The same impurity leaves the system's spin transport at the Kardar-Parisi-Zhang universality class unchanged in the thermodynamic limit, a result that may extend to other integrable systems with impurities.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivation of Eq. (7) assumes a sharp light-cone delay and unmodified exponential growth at the global Lyapunov rate after the front; this is unproven and gives only an upper bound, so the central ξ=1/λ relation could fail if the post-front growth rate differs.","rationale":"The reader's weakest_assumption identifies exactly the same step: the heuristic factorization of the perturbation into a pure time delay and a single global Lyapunov growth rate, with no derivation and no justification for the extension to subleading vectors. I agree that this is the most load-bearing concern because Eq. (7) is the paper's central novelty: the singular Lyapunov spectrum itself is supported by Fig. 1 and Fig. 2(c), but the quantitative localization-exponent relation is the content that makes the result physically meaningful and universal. The transport claim is important but less novel (it extends known KPZ behavior of the unperturbed model) and is supported by fairly direct correlation-function scaling plots. The heuristic for Eq. (7) is the only theoretical link between the local perturbation and the localization length; if it is wrong, the numerics might still hold for the specific parameters (τ=1, φ=π/2) but the claimed universality and the exact ξ=1/λ scaling would be unsupported. I credit the paper for honestly labeling the argument as heuristic and the universality as a conjecture, and for providing numerical evidence over a range of parameters (Fig. S1) and a second model (Figs. S3–S4). Nevertheless, the precise proportionality constant (1 vs. v/λ) is not derived, and the sharp-front assumption is testable. The proposed experiment directly measures the time-delayed growth law and would settle whether the factor d(0)e^{λ_1(t−|x*|/v)} is accurate. If it reproduces the predicted delay and slope for k=1 and k=2, the concern is resolved and the paper's central claim is strengthened; if not, Eq. (7) must be revised. Since the reader's verdict is already CONDITIONAL and this concern is a core reason, the verdict should remain unchanged: the paper needs either a derivation of the growth law or an explicit conjecture status for Eq. (7), plus error bars and reproducibility artifacts.","tokens_in":11239,"tokens_out":13295,"duration_ms":127713,"concrete_test":"Perform a direct growth experiment on a single localized tangent perturbation. Choose a reference trajectory from the canonical ensemble at µ=0, L≥2^12, τ=1, φ=π/2. Initialize δX^0 nonzero only at site x* (e.g. x*=50, 100, 200), tangent to the sphere. Integrate the linearized map up to t≈5×10^3 and record d(t)=||δX^t||. Fit the late-time slope after the apparent crossover and compare to the independently computed λ_1; check that the crossover time equals |x*|/v with v=1. Repeat for k=2 by evolving the same initial vector but first projecting out the ψ_1 component via Gram-Schmidt, and compare the late-time slope to λ_2. If the post-front slope deviates from λ_k by more than a few percent, Eq. (7) fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is Eq. (7): ρ_k(x) ≍ e^{-λ_k|x|}, i.e. ξ_k = 1/λ_k. The only theoretical support is the heuristic in the paragraph after Eq. (8). It assumes that once the perturbation front from x=0 reaches a site x*, the tangent-space deviation seeded there grows as d(t) ∼ d(0) e^{λ_1(t−|x*|/v)} for t > |x*|/v. This factorizes the perturbation influence into a pure time delay |x*|/v and an unchanged global growth rate λ_1. The argument yields only the bound ξ_1 ≤ v/λ_1, not the observed equality; the saturation to v/λ_1 (with v=1) is not derived. More importantly, the post-front growth rate is assumed to equal λ_1 without justification. The local perturbation at x=0 modifies the Jacobian at that site, so a tangent vector passing through the impurity can be scattered into other directions; the asymptotic growth rate of that component need not be λ_1. If the effective rate after the front is λ' ≠ λ_1, the comparison with Eq. (8) gives |ψ_1,x*| ∼ e^{-λ'|x*|/v}, and the observed ξ=1/λ relation would be an accident of parameters. The extension to subleading k>1 is even less secure: the argument is constructed for the leading tangent direction, and the orthogonalized subspace dynamics (which defines ψ_k) is not described by the same delay-growth picture. The manuscript itself hedges, saying the argument 'suggests' the scaling and later 'conjectures' universality, but Eq. (7) is presented as a main finding. Since this relation is the paper's most novel and load-bearing claim, the heuristic's unverified factorization is the key point of vulnerability.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the discrete space-time Landau-Lifshitz magnet with a single-site local rotation that breaks integrability. It reports three main findings: (i) the Lyapunov spectrum becomes singular in the thermodynamic limit, with a vanishing fraction of nonzero exponents and a density converging to a delta function; (ii) the Lyapunov vectors corresponding to the nonzero exponents are exponentially localized around the impurity, with localization lengths satisfying ξ_k = 1/λ_k; and (iii) spin-spin and current-current correlations retain KPZ scaling with dynamical exponent z=3/2 at zero magnetization and ballistic scaling at nonzero magnetization, so the local impurity does not change bulk transport universality. The theoretical support for the localization-length relation is a heuristic delay-growth argument, while the transport and spectral claims are supported by extensive numerical simulations, including a second classical unitary circuit model in the supplemental material.","tokens_in":11622,"tokens_out":4734,"duration_ms":47513,"significance":"If the results hold, they identify a new and potentially universal mechanism for local integrability breaking: a single-site impurity produces a singular Lyapunov spectrum and exponentially localized Lyapunov vectors, while leaving bulk transport in the integrable universality class unchanged. The numerical work is substantial and carefully executed: long integration times (up to 10^6–10^7), large ensembles (up to 5×10^7 initial states), and scaling collapses that are visually convincing. The paper also contains a useful consistency check: the localization decay rate is taken from the independently measured Lyapunov spectrum rather than fitted from the spatial profile. The universality conjecture is supported by a second model in the supplemental material. The main weakness is that the central quantitative relation ξ_k = 1/λ_k is not derived from the tangent-space dynamics and is presented with a degree of confidence that the heuristic argument does not fully justify.","major_comments":[{"comment":"The heuristic argument does not establish Eq. (7) as a derived result. For k=1, the argument yields only the upper bound ξ1 ≤ v/λ1, not the observed equality; the saturation to v/λ1 is assumed rather than explained. More importantly, the post-front growth of the tangent-space deviation is assumed to occur at the unperturbed global rate λ1, but the local rotation at x=0 modifies the Jacobian at that site, so a deviation passing through the impurity can be scattered into other directions and may grow at a different rate. If the effective post-front rate were λ' ≠ λ1, the same comparison would give ξ1 = v/λ', and the numerical agreement in Fig. 2(a) would be parametric. The extension to subleading vectors k>1 is even less secure because the orthogonalized subspace dynamics defining ψ_k is not described by the same delay-growth picture. I recommend either deriving the equality from the tangent-space dynamics or explicitly presenting Eq. (7) as a numerically supported conjecture, with the upper bound as the only proven statement.","section":"Lyapunov spectrum and Lyapunov vectors, Eq. (7) and paragraph after Eq. (8)"},{"comment":"The claim that lim_{L→∞} P(λ) = δ(λ) is an extrapolation from visual narrowing of the central spectral peak for L = 2^7, ..., 2^10. The paper does not provide error bars or a quantitative finite-size scaling analysis of the peak width, peak height, or the fraction of exponents inside a fixed window. The convergence of the five largest exponents in Fig. 2(c) does not by itself imply that the density becomes singular. Please add a scaling collapse for the central peak (e.g., width w(L) ∝ L^{-α} with α>0, or an equivalent measure) or soften the statement to a finite-size-supported conjecture.","section":"Fig. 1 and text after 'Lyapunov spectrum and Lyapunov vectors'"}],"minor_comments":[{"comment":"'Lyapunov exponense' should be 'Lyapunov exponents'.","section":"Paragraph after Eq. (6)"},{"comment":"The text says 'Fig. 2 (b) shows Ik for the first 100 Lyapunov vectors', but the IPR panel is Fig. 2 (d); please correct the reference.","section":"Fig. 2 caption and text"},{"comment":"Please clarify the normalization convention in the QR-based estimate of λ_k: the expression uses Rjj(2n) only, and it should be stated whether this corresponds to a full map step or a half step and how the time average is normalized.","section":"Eq. (S7) in the Supplemental Material"},{"comment":"The notation d(t) ≡ |δX^t| and the subsequent comparison with Eq. (8) would benefit from a statement that the norm is taken after projecting onto the tangent sphere; otherwise the argument is ambiguous about whether δX^t is the full tangent vector or its magnitude.","section":"Eq. (8) and following paragraph"}],"recommendation":"major_revision","confidential_remarks":"The main physics is interesting and the numerical evidence is strong, but the central relation ξ_k = 1/λ_k is presented as a finding while the derivation is only an upper-bound heuristic. This is a load-bearing point for the paper's novelty. I do not see circularity in the numerical consistency check, and the transport results are credible. The revision should either strengthen the derivation or clearly demote Eq. (7) to a numerically supported conjecture. The finite-size support for the δ-limit should also be quantified. These issues are fixable within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"My take in one sentence: this is a solid numerical paper with a genuinely new result, but the one theoretical explanation on offer is a heuristic that should be presented as a conjecture, and the absence of error bars matters more than usual because the central claim is a quantitative scaling law.\n\nWhat is new: prior Lyapunov-vector work on integrability breaking used global perturbations; this paper studies a single-site impurity. The singular spectrum — vanishing density P(λ)→δ(λ) in the thermodynamic limit with a finite sequence of non-zero exponents — is new, and the exponential localization of the corresponding vectors with ξ_k=1/λ_k is, as far as I know, not in the literature. The transport result, KPZ scaling z=3/2 unchanged by the impurity at zero magnetization, extends earlier integrable-model results and is worth having.\n\nWhat the paper does well: the numerics are beefy. Integration times up to 10^6, ensembles up to 5×10^7, many system sizes, parameter sweeps in φ and τ, and a second model (unitary circuit) with the same localization behavior. The scaling collapses in Figs. 2 and 3 look convincing. Importantly, the localization length is not fitted from the spatial profile; the exponent λ_k is measured independently from the spectrum and then compared with the profile. That is a genuine consistency check.\n\nWhere the soft spots are: the derivation of ξ=1/λ after Eq. (8) assumes a sharp light-cone delay and that the perturbation's influence grows at the unmodified global rate λ_1 once the front passes. That gives only an upper bound ξ_1 ≲ v/λ_1, not the observed equality, and the extension to subleading k>1 is asserted without an argument. The stress-test note worries that a different post-front growth rate would break the relation; I don't think that is demonstrated, but it is a real gap in the theory. The paper itself is honest about this ('suggests', 'conjectures'), which I credit. What worries me more is the absence of error bars on λ_k, IPR, and the scaling exponents. For a paper whose headline is a quantitative law, I want to see error bars or at least convergence data. Also, the vanishing-density result is an extrapolation from L up to 2^14; that is fine, but it should be flagged as an extrapolation. No code or data is provided, which makes independent checking harder.\n\nBottom line: the central claims are likely correct and the paper deserves serious refereeing. I would ask the authors to either tighten the theoretical argument or explicitly label Eq. (7) as a conjecture supported by numerics, and to add error bars and reproducibility artifacts. This is a conditional accept, not a reject.","headline":"Genuinely new numerical result on local integrability breaking, with extensive numerics and a plausible but unproven central relation; deserves serious refereeing, conditional on error bars and honest labeling of the conjecture.","tokens_in":12166,"tokens_out":2609,"would_cite":true,"duration_ms":23816,"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":"The paper shows that a single-site integrability-breaking perturbation in a discrete space-time spin lattice makes the Lyapunov spectrum singular in the thermodynamic limit: the density of exponents collapses to a delta function while a…","keywords":["Lyapunov spectrum","Lyapunov vectors","local integrability breaking","exponential localization","Kardar-Parisi-Zhang scaling","space-time lattice","Landau-Lifshitz magnet","dynamical exponent"],"falsifier":"Compute, for system sizes up to $L=2^{14}$ and long times, the logarithm of the spatial profile of the second and third backward Lyapunov vectors; if the slope $-\\log \\rho_k(x)/|x|$ is not numerically equal to $\\lambda_k$ for large $|x|$ (or deviates systematically as the perturbation strength $\\phi\\to 0$), the exponent of localization is not the Lyapunov exponent and Eq. (7) fails. Alternatively, measure $C(x,-x,t)$ at fixed small $x$ and increasing $t$: a crossover from $t^{-1/3}$ (KPZ) to $t^{-1/2}$ (diffusive) scaling would falsify the claim that the local impurity leaves the bulk transport universality class unchanged.","tokens_in":11049,"feed_emoji":"🧲","tokens_out":9129,"duration_ms":72252,"temperature":0.7,"pith_summary":"The paper studies what happens to deterministic chaos and transport when a discrete space-time lattice that is integrable is perturbed at a single site. It claims that in the thermodynamic limit the density of Lyapunov exponents becomes a delta function $\\delta(\\lambda)$ at zero, so that almost all exponents vanish, while a vanishing fraction of non-zero exponents survives and converges. The Lyapunov vectors belonging to those surviving exponents are exponentially localized around the impurity, with localization length equal to the inverse of the corresponding exponent. The authors also find that this local impurity does not change the bulk transport universality class: spin-spin and current-current correlations keep the Kardar-Parisi-Zhang scaling with dynamical exponent $z=3/2$ at zero magnetization, and ballistic scaling at non-zero magnetization. A sympathetic reader would care because it suggests a simple, universal structure for how local integrability breaking affects many-body chaos.","feed_headline":"Single-site impurity collapses most Lyapunov exponents","feed_subtitle":"Surviving chaos modes localize with length 1/λ; bulk transport keeps KPZ scaling.","key_machinery":"The machinery is the local perturbation combined with a light-cone argument: an infinitesimal deviation seeded at a site $x^*$ first grows with the maximal Lyapunov exponent $\\lambda_1$ until the influence of the perturbation at $x=0$ reaches it after a time $|x^*|/v$ (with sound speed $v=1$ here); the resulting time-delayed growth, compared with the free-growth formula, forces the overlap of the first Lyapunov vector with site $x^*$ to decay as $e^{-\\lambda_1 |x^*|}$, hence $\\xi_1 = 1/\\lambda_1$, and the argument is extended to subleading vectors $k>1$. The numerical workhorse is a standard QR-decomposition scheme for extracting Lyapunov exponents and backward Lyapunov vectors from long trajectories. For transport, the load-bearing objects are the spin-spin correlation $C(x_1,x_2,t)$ and the current-current correlation $C_{JJ}(t)$, whose scaling collapses under $t^{1/z}$ with $z=3/2$ ($\\mu=0$) or $z=1$ ($\\mu\\neq 0$).","core_discovery":"Consider the SO(3)-symmetric discrete space-time Landau-Lifshitz magnet with a local rotation applied to a single site at $x=0$. The paper's central discovery is that this local integrability breaking produces a Lyapunov spectrum whose thermodynamic limit $L\\to\\infty$ is singular: the spectral density $P(\\lambda)$ tends to $\\delta(\\lambda)$, so the overwhelming majority of Lyapunov exponents vanish, while a sub-extensive sequence of finite exponents $\\lambda_k$ survives. The corresponding backward Lyapunov vectors are exponentially localized around the impurity, with spatial profile $\\rho_k(x) \\asymp e^{-\\lambda_k |x|}$, that is, the localization length is $\\xi_k = 1/\\lambda_k$. The same perturbation leaves the transport properties of the bulk unchanged in the scaling limit: the two-point spin and current correlators obey KPZ scaling with $z=3/2$ at zero magnetization (and ballistic scaling at non-zero magnetization), with the spread of the perturbation's influence governed by the same dynamical exponent. The authors conjecture this structure is universal for local integrability breaking in locally interacting many-body systems.","pith_inferences":["A testable extension the paper leaves open: if the number of impurities grows as a finite density, the singular spectrum should break down, with the number of non-zero exponents becoming extensive and localization lengths shortening; this crossover could be probed numerically.","A natural quantum analogue: in spin chains with a single impurity, the response of out-of-time-order correlators or operator spreading might show a similar exponential localization with length set by the local Lyapunov rate, if such a rate can be defined.","The sharp light-cone assumption could be relaxed: in models with sub-ballistic front propagation, the relation $\\xi_k=1/\\lambda_k$ would acquire corrections set by the front's spreading exponent, a modification not covered by the paper's heuristic derivation.","The two-parameter scaling function is presented as a limit; whether its shape is universal or changes with the perturbation strength $\\phi$ and the coupling $\\tau$ is not tested in the paper and could be checked by data collapse at different parameters."],"forward_implications":["In the thermodynamic limit, the Lyapunov spectral density collapses to $\\delta(\\lambda)$; all but a sub-extensive set of exponents vanish.","Each surviving non-zero exponent $\\lambda_k$ corresponds to a Lyapunov vector exponentially localized around the impurity, with localization length $\\xi_k=1/\\lambda_k$; vectors of zero exponents become delocalized.","The local impurity does not change the bulk transport universality class: the KPZ dynamical exponent $z=3/2$ at zero magnetization and ballistic transport at non-zero magnetization persist in the thermodynamic limit.","A two-parameter scaling function for the spin correlator exists, whose single-parameter slice is the standard KPZ scaling function.","The same structure is observed in a locally perturbed classical unitary circuit map, suggesting universality across different locally perturbed integrable systems."],"supporting_citations":[{"why":"Introduces the unperturbed discrete space-time Landau-Lifshitz lattice and establishes its KPZ scaling and conserved charges, providing the baseline for the perturbed model.","marker":"[23]"},{"why":"Provides the QR-decomposition method used to compute Lyapunov exponents and backward Lyapunov vectors.","marker":"[31]"},{"why":"Supplies the one-dimensional scaling function of the KPZ universality class that the paper expects as a slice of its two-dimensional scaling function.","marker":"[32]"},{"why":"A classical unitary circuit map with global perturbation serves as the comparison model; the paper applies the same local-perturbation analysis and finds compatible localization.","marker":"[10]"},{"why":"Provides evidence of KPZ scaling in integrable spin chains, supporting the connection between integrability and $z=3/2$ transport.","marker":"[24]"},{"why":"Discusses KPZ scaling in the space-time lattice setting, reinforcing the unperturbed transport picture that the paper extends.","marker":"[25]"}],"fun_headline_variants":["Local impurity collapses most Lyapunov exponents","Impurity kills Lyapunov spectrum, survivors localize","Single-site break: Lyapunov spectrum vanishes, KPZ stays","Lyapunov exponents vanish locally, KPZ transport persists","Exponential localization of survivors of local impurity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument for $\\xi_k=1/\\lambda_k$ assumes that a small deviation near site $x$ grows at the unperturbed maximal rate $\\lambda_1$ both before and after the local impurity's influence reaches it, the influence acting only as a time delay $|x|/v$, and that this same single-rate picture carries over to every subleading vector $k$.","fun_headline_variants_meta":{"raw":{"variants":["Local impurity collapses most Lyapunov exponents","Impurity kills Lyapunov spectrum, survivors localize","Single-site break: Lyapunov spectrum vanishes, KPZ stays","Lyapunov exponents vanish locally, KPZ transport persists","Exponential localization of survivors of local impurity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1590,"prompt_tokens":909,"completion_tokens":681,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":604}},"tokens_in":525,"tokens_out":681,"duration_ms":7055,"temperature":1.0,"reasoning_tokens":604,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T06:00:53.033415+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, for system sizes up to $L=2^{14}$ and long times, the logarithm of the spatial profile of the second and third backward Lyapunov vectors; if the slope $-\\log \\rho_k(x)/|x|$ is not numerically equal to $\\lambda_k$ for large $|x|$ (or deviates systematically as the perturbation strength $\\phi\\to 0$), the exponent of localization is not the Lyapunov exponent and Eq. (7) fails. Alternatively, measure $C(x,-x,t)$ at fixed small $x$ and increasing $t$: a crossover from $t^{-1/3}$ (KPZ) to $t^{-1/2}$ (diffusive) scaling would falsify the claim that the local impurity leaves the bulk transport universality class unchanged.","supporting_citations":[{"cited_title":"Local integrability breaking and exponential localization of leading Lyapunov vectors","cited_arxiv_id":"2412.16887","evidence_quote":"Introduces the unperturbed discrete space-time Landau-Lifshitz lattice and establishes its KPZ scaling and conserved charges, providing the baseline for the perturbed model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the QR-decomposition method used to compute Lyapunov exponents and backward Lyapunov vectors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the one-dimensional scaling function of the KPZ universality class that the paper expects as a slice of its two-dimensional scaling function."},{"cited_title":"Danieli, T","cited_arxiv_id":null,"evidence_quote":"A classical unitary circuit map with global perturbation serves as the comparison model; the paper applies the same local-perturbation analysis and finds compatible localization."},{"cited_title":"Taniguchi and G","cited_arxiv_id":null,"evidence_quote":"Provides evidence of KPZ scaling in integrable spin chains, supporting the connection between integrability and $z=3/2$ transport."},{"cited_title":"Krajnik and T","cited_arxiv_id":null,"evidence_quote":"Discusses KPZ scaling in the space-time lattice setting, reinforcing the unperturbed transport picture that the paper extends."}],"review_version":1}