{"id":"f2da2f13-1981-4214-ac26-fa22435e5982","arxiv_id":"1908.07594","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In 2D disordered nonlinear lattices, the second moment of spreading wave packets grows as t^(1/5) in the weak and t^(1/3) in the strong chaos regime, while the maximum Lyapunov exponent decays as t^(-0.37) and t^(-0.46).","lead":"Wave packets in two-dimensional disordered nonlinear lattices spread slowly with a power law in time, and the chaos that drives the spreading weakens as a power law too. This paper measures those rates in two standard models and proposes a scaling rule that connects one- and two-dimensional behavior.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (10) is an assumed, post hoc validated scaling conjecture and the sole basis for the predicted α_Λ values; the predictive claim needs a direct test.","rationale":"The Reader's verdict is CONDITIONAL, and my analysis supports that disposition. The empirical results — spreading exponents near 1/5 and 1/3, MLE decay exponents near -0.37 and -0.46, and the observed wandering of chaotic seeds — are plausible and consistent with prior 1D studies; they are not the locus of my concern. The load-bearing soft spot is the new theoretical element, Eq. (10). The paper's own wording describes it as an assumption and a conjecture, and the only validation offered is a post hoc comparison of the exponent combination in Eq. (11) with fitted α_Λ values, at just two regimes and without stated uncertainties on α_Λ. This does not invalidate the paper: if Eq. (10) is false, the numerical findings remain as empirical facts, but the proposed universal scaling and its predictive claim do not survive. Since the paper itself acknowledges the conjecture and a direct test is feasible with existing data, the appropriate verdict is to keep the CONDITIONAL status rather than treat the scaling as established or reject the manuscript entirely. I therefore recommend UNCHANGED.","tokens_in":10734,"tokens_out":5395,"duration_ms":84312,"concrete_test":"Directly test Eq. (10) at the level of raw observables: using the 1D simulations of Refs. [17,19] and the 2D simulations of this paper, compute R(t) = Λ(t)/m2(t) for each realization over a common time window (e.g., t = 10^4 to 10^6), average in log space, and plot R_1D(t) and R_2D(t) on the same log-log axes for the weak and strong chaos cases of both DKG and DDNLS. If Eq. (10) holds, the curves should collapse onto a common power law; the degree of collapse (e.g., curves within a factor of 2 over at least two decades) would settle whether the scaling is real. Reporting local slopes with confidence intervals would also quantify the precision of the claimed 0.01 agreement from Eq. (11).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is the dimension-independent scaling law, Eq. (10): Λ1(t)/m2_1(t) = Λ(t)/m2(t). The paper introduces it in Sec. III as an assumption ('This can be quantified by assuming ...'), with no derivation from the equations of motion and no independent evidence. The headline claim that the proposed scaling 'allows the prediction of the obtained α_Λ values' rests entirely on this relation. Eq. (11) combines 1D measured exponents and theoretical a_m values to produce predicted exponents -0.38 (weak) and -0.47 (strong), which match the numerically fitted -0.37 and -0.46. But the 2D exponents were measured on the same data, no uncertainty is quoted for α_Λ, and the agreement at the 0.01 level is not quantified. If Eq. (10) is false, the measured am and α_Λ values remain valid empirical observations, but the proposed universal scaling and its predictive use collapse. The paper itself labels the relation a conjecture and defers a 3D test, so the scaling law is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports extensive numerical simulations of wave-packet spreading in two-dimensional disordered nonlinear lattices, specifically the 2D Klein-Gordon (DKG) model and the 2D discrete nonlinear Schrödinger equation (DDNLS), in both the weak and strong chaos regimes. The authors find that the second moment of the wave packet grows as m2 ∝ t^{a_m} with a_m ≈ 1/5 in the weak chaos case and a_m ≈ 1/3 in the strong chaos case, consistent with earlier theoretical predictions by Flach. They also compute the finite-time maximum Lyapunov exponent Λ and find Λ ∝ t^{α_Λ} with α_Λ ≈ −0.37 (weak chaos) and −0.46 (strong chaos). The deviation-vector distributions show localized chaotic seeds wandering within the excited part of the lattice. The central new proposal is a dimension-independent scaling relation, Eq. (10): Λ_1(t)/m2_1(t) = Λ(t)/m2(t), which the authors use to predict α_Λ values of −0.38 and −0.47. The paper concludes that the obtained α_Λ values are accurately predicted by this scaling.","tokens_in":11121,"tokens_out":3986,"duration_ms":476663,"significance":"If the empirical results hold, this paper provides a substantial advance: it extends the study of subdiffusive spreading and chaoticity in disordered nonlinear lattices from one to two dimensions, covering both the weak and strong chaos regimes, and it demonstrates universality across two different models with quartic nonlinearities. The numerical effort is a clear strength: integration times up to 10^8, lattice sizes up to 450×450, 50 disorder realizations, and several parameter sets per model. The measured spreading exponents a_m ≈ 1/5 and 1/3 match independent theoretical predictions with small quoted statistical errors, and the DVD analysis provides a plausible mechanism for thermalization. However, the paper's central new claim—a dimension-independent scaling that predicts α_Λ—is based on an assumed relation, Eq. (10), that is not derived and is only validated indirectly. The manuscript is therefore a valuable empirical study whose theoretical interpretation needs to be reframed or further supported.","major_comments":[{"comment":"The dimension-independent scaling relation Λ_1(t)/m2_1(t) = Λ(t)/m2(t) is introduced as an assumption ('This can be quantified by assuming...') and is the sole basis for the predicted α_Λ values in Eq. (11). Because the authors themselves label this a conjecture and defer a direct 3D test, the headline statement that the scaling 'allows the prediction of the obtained α_Λ values' overstates the evidence. The agreement between the predicted exponents (−0.38, −0.47) and the fitted ones (−0.37, −0.46) is not a convincing validation because the scaling relation was formulated after the 2D results were known and because no independent data are used to test it. I recommend either (a) providing a direct test of Eq. (10) by computing Λ(t)/m2(t) for the 1D and 2D time series over the same time window and showing that the ratio is independent of dimension, or (b) explicitly reframing Eq. (11) as a heuristic consistency check rather than a prediction. As written, the central claim of a universal scaling law is not established.","section":null},{"comment":"The exponents α_Λ are quoted as ≈ −0.37 and −0.46 without any uncertainty estimate, while the a_m values in footnotes [33] and [34] are given with standard errors. Moreover, the fits are performed over the last two decades of the evolution, a window that appears to be selected post hoc. This matters because the agreement between the predicted and observed α_Λ values is at the 0.01 level, and without error bars or a stated window-selection criterion that agreement cannot be quantitatively assessed. Please report uncertainties for α_Λ, for example by using a bootstrap regression over several fitting windows, and state the rule used to choose the fitting interval.","section":null},{"comment":"The claim that the measured exponents a_m ≈ 1/3 confirm the strong chaos regime relies on the asymptotic behavior visible in the last two decades of the data. Since the strong chaos cases have larger initial excitations and a crossover from a transient regime is possible, the manuscript should demonstrate that the fitted values are robust to the choice of the fitting window, for instance by showing that the local derivative a_m(t) is flat over at least two non-overlapping sub-windows within the fitted range. Without this, the distinction between weak and strong chaos rests on a fitting procedure that is not fully documented.","section":null}],"minor_comments":[{"comment":"The last entry in footnote [34] reads '0.3362±0.0002 (S2D)', but S2D was already listed earlier; this is presumably a typo for S3D.","section":null},{"comment":"In the list of strong chaos DKG cases, the case with W = 12.5, L = 15, h_{l,m} = 0.035 is labeled '(Case W 3K)', which is inconsistent with the S1K/S2K labels used for the other strong chaos cases; this should likely be S3K.","section":null},{"comment":"The notation r^{(D)}_{l,m} and \\bar{r}^{(D)}_{l,m} is confusing because the superscript D is used both for the deviation-vector distribution and for the dimensional label of the center; please use a clearer notation, e.g., \\bar{\\mathbf{r}}^{(D)} for the center of the distribution.","section":null},{"comment":"The notation 'max[0,t]{l^D(t)}' is not standard; please explain that it denotes the maximum over the time interval [0,t], or use a more explicit notation such as \\max_{0 \\le s \\le t} l^D(s).","section":null},{"comment":"The caption lists a_D^P = 0.0 for the strong chaos insets, while the text states 'a very slow increase P_D ∝ t^{0.045}'; please clarify whether the strong chaos participation number is strictly constant or grows with a very small exponent.","section":null}],"recommendation":"major_revision","confidential_remarks":"The empirical core of this paper—long-time 2D simulations of DKG and DDNLS in both chaos regimes—appears solid and is likely of genuine interest to the nonlinear disordered lattices community. My main reservation is the overinterpretation of Eq. (10) as a predictive scaling law when it is introduced as an assumption and validated only through the very data it is used to predict. This is fixable by reframing the claim and adding error bars and robustness checks for the α_Λ fits. I do not see grounds for rejection, but the revision should be substantial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is a careful numerical study that gives the first 2D measurements of MLE decay for these two standard models and confirms the expected spreading exponents (1/5 weak, 1/3 strong). The practical value is in the new data, not in the scaling conjecture.\n\nWhat's new: they simulate both DKG and DDNLS in 2D up to t=10^8, with 50 realizations and large lattices, and get clean power laws for m2 and Λ. The strong-chaos regime in 2D had not been done before. The DVD analysis shows localized chaotic seeds wandering, which is a nice visualization and matches the 1D picture. The measured spreading exponents match Flach's theory well, and the quoted errors are small.\n\nSoft spots: the α_Λ exponents are quoted without error bars, and the fitting window is the last two decades, chosen post hoc. The scaling law, Eq. (10), is explicitly introduced as an assumption, and the agreement of its predicted exponents (-0.38, -0.47) with the measured (-0.37, -0.46) is the kind of agreement that would be more convincing with a derivation or a direct test on independent data. The paper itself defers a 3D test, so it is a conjecture, not a result. That is fine as long as it is not oversold—and they do not oversell it.\n\nThe main empirical claims—the exponents and the chaoticity decay—are solid. I do not think the scaling conjecture undermines them. I would send this to a referee; it deserves a careful look. I would ask the authors to add error bars for α_Λ and to be clearer about the window selection. The lack of deposited code/data is a minor annoyance, not a blocker.\n\nRecommendation: accept for review. It is a useful paper for anyone working on disorder and nonlinearity.","headline":"Solid 2D numerics for spreading and MLE decay; the scaling conjecture is the soft spot but honestly labeled.","tokens_in":11582,"tokens_out":2671,"would_cite":false,"duration_ms":126962,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In 2D disordered nonlinear lattices, wave packets spread by t^{1/5} (weak chaos) or t^{1/3} (strong chaos), while the maximum Lyapunov exponent decays as t^{-0.37} or t^{-0.46}, and a proposed dimension-independent scaling law ties the…","keywords":["wave packet spreading","Anderson localization","disordered nonlinear lattices","Klein-Gordon lattice","discrete nonlinear Schrödinger equation","Lyapunov exponent","subdiffusion","chaotic seeds"],"falsifier":"Compute the ratio Λ(t)/m2(t) for both one- and two-dimensional DKG and DDNLS systems in the same dynamical regime and check whether the curves for different dimensions collapse onto each other; failure to collapse would disprove the dimension-independent scaling law. A complementary test would simulate a 3D disordered nonlinear lattice and compare the measured Lyapunov-exponent decay exponent with the value predicted by Eq. (11) under the same scaling assumption.","tokens_in":10523,"feed_emoji":"🌀","tokens_out":5616,"duration_ms":52777,"temperature":0.7,"pith_summary":"This paper asks what happens to Anderson localization in two-dimensional disordered lattices when nonlinearity is present, and tries to establish that the resulting wave-packet spreading and chaoticity obey universal power laws shared by two quite different models: the disordered Klein-Gordon (DKG) lattice and the disordered discrete nonlinear Schrödinger equation (DDNLS). Through long-time numerical simulations, it reports that the wave packet's second moment grows as m2 ∝ $t^{{1/5}}$ in the weak chaos regime and as $t^{{1/3}}$ in strong chaos, matching previously predicted exponents for two dimensions. It also finds that chaos persists but weakens in time, with the finite-time maximum Lyapunov exponent decaying as Λ ∝ $t^{{-0.37}}$ in weak chaos and $t^{{-0.46}}$ in strong chaos. The paper proposes a dimension-independent scaling relation between spreading and chaoticity, from which these Lyapunov-decay exponents can be predicted rather than merely measured. If correct, this would unify the subdiffusive spreading and thermalization of disordered nonlinear systems across spatial dimensions.","feed_headline":"Wave packets spread by t^1/5 and t^1/3 in 2D disorder","feed_subtitle":"Two different nonlinear disordered models show the same subdiffusive spreading and chaos decay, linked by a scaling law.","key_machinery":"The central objects are the normalized energy (or norm) density distribution ξ_{l,m}, its second moment m2, and the finite-time maximum Lyapunov exponent Λ computed from a deviation vector's evolution. The load-bearing identity is the proposed scaling law of Eq. (10), Λ1(t)/m2_1(t) = Λ(t)/m2(t), which asserts that the ratio of chaoticity to wave-packet extent is independent of lattice dimensionality; combined with known exponents it predicts Λ ∝ $t^{{a_m - a_m^1 + α_Λ^1}}$. The deviation-vector distributions and their wandering area A(t) supply the mechanism: moving chaotic seeds thermalize the wave packet interior.","core_discovery":"The paper claims that in both 2D DKG and DDNLS lattices with quartic nonlinearities, an initially compact wave packet spreads subdiffusively and destroys Anderson localization, with the same power laws in both models: m2 ∝ $t^{{1/(1+2d)}}$ = $t^{{1/5}}$ in the weak chaos regime and m2 ∝ $t^{{1/(1+d)}}$ = $t^{{1/3}}$ in the strong chaos regime. It further claims that chaos persists throughout the observed evolution, with the finite-time maximum Lyapunov exponent decaying as Λ ∝ $t^{{αΛ}}$ with αΛ ≈ −0.37 (weak chaos) and −0.46 (strong chaos), slower than the regular-motion decay $t^{{-1}}$, indicating no crossover to regularity. The deviation-vector distributions show localized 'chaotic seeds' wandering randomly through the excited part of the lattice, homogenizing chaos and supporting thermalization. Finally, the paper proposes the dimension-independent scaling assumption Λ1(t)/m2_1(t) = Λ(t)/m2(t) between one- and two-dimensional systems, which, combined with earlier 1D results and theoretical spreading exponents, predicts αΛ ≈ −0.38 and −0.47, in good agreement with the numerically observed values.","pith_inferences":["If the dimension-independent scaling law holds in three dimensions, Eq. (11) predicts αΛ ≈ −0.42 for weak chaos (with a_m = 1/7) and ≈ −0.55 for strong chaos (with a_m = 1/4), a testable extension the paper leaves for future work.","The scaling relation Λ/m2 being dimension-independent suggests that the ratio of the chaoticity rate to the effective wave-packet area is a dynamical invariant across lattice geometries, which could be probed directly by computing Λ(t) m2(t) for one- and two-dimensional systems in the same regime.","The contrast between the unweighted wandering area A(t) growing as t^{0.5}–t^{0.55} and the weighted second moment growing as t^{0.2}–t^{0.33} implies that the 'thermalization front' extends beyond the energy-carrying core of the wave packet; measuring this front separately in experiments could provide a finer test of the thermalization picture.","Because both models yield nearly identical exponents for every observable studied, the results support a universality class for disordered nonlinear wave-packet spreading that may extend beyond quartic nonlinearities to other power-law nonlinearities."],"forward_implications":["The destruction of Anderson localization by nonlinearity is generic in two dimensions, since both DKG and DDNLS show identical spreading exponents, t^{1/5} and t^{1/3}, in their respective regimes.","Chaos persists for at least the simulated time scales, as Λ decays more slowly than t^{-1}, so no crossover to regular or quasiperiodic dynamics is expected on those time scales.","Thermalization precedes spreading: the Lyapunov time remains shorter than the spreading time, with T_D/T_L growing as t^{0.43} (weak chaos) and t^{0.21} (strong chaos).","The wandering area of the deviation-vector center grows as t^{0.5} (weak) and t^{0.55} (strong), showing that chaotic seeds visit a widening region as the wave packet expands.","The proposed scaling law connects chaoticity to spreading across dimensions, allowing Lyapunov-exponent decay to be predicted from spreading exponents rather than fitted independently."],"supporting_citations":[{"why":"Supplies the theoretical predictions a_m = 1/(1+2d) and 1/(1+d) for weak and strong chaos that the measured 2D exponents are compared against.","marker":"[12]"},{"why":"Provides prior 2D DKG weak-chaos results and the 'localization volume' context that the present long-time simulations extend.","marker":"[16]"},{"why":"Earlier 2D DDNLS weak-chaos study that this work extends by using larger systems, more disorder realizations, and MLE computations.","marker":"[5]"},{"why":"Establishes the 1D MLE and deviation-vector-distribution methodology and provides the 1D chaoticity data used in the proposed scaling law.","marker":"[17]"},{"why":"Supplies the 1D DKG and DDNLS Lyapunov-exponent decay exponents and the wandering-chaotic-seeds analysis that the 2D results are compared with.","marker":"[19]"},{"why":"Provides the ABA864 symplectic integrator with tangent-map method used for the DKG model's evolution and variational equations.","marker":"[23]"},{"why":"Provides the s11ABC6 integration scheme used for the DDNLS system's evolution.","marker":"[24]"},{"why":"Defines the finite-time maximum Lyapunov exponent method that the paper adopts for detecting chaos and its decay.","marker":"[25]"}],"fun_headline_variants":["Two 2D disordered models share t^1/5 and t^1/3 spreading","Chaos decay and spreading linked by scaling in 2D disorder","Wave packets spread subdiffusively in 2D disordered lattices","Universal spreading exponents in 2D nonlinear disorder","2D disordered systems: same spreading, same chaos decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the dimension-independent scaling law Λ1(t)/m2_1(t) = Λ(t)/m2(t), which is assumed without derivation; if that ratio equality fails, the predicted Lyapunov-decay exponents (−0.38 and −0.47) do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Two 2D disordered models share t^1/5 and t^1/3 spreading","Chaos decay and spreading linked by scaling in 2D disorder","Wave packets spread subdiffusively in 2D disordered lattices","Universal spreading exponents in 2D nonlinear disorder","2D disordered systems: same spreading, same chaos decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000353,"raw_usage":{"total_tokens":1973,"prompt_tokens":1045,"completion_tokens":928,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":661,"completion_tokens_details":{"reasoning_tokens":836}},"tokens_in":661,"tokens_out":928,"duration_ms":8060,"temperature":1.0,"reasoning_tokens":836,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:01:44.450587+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the ratio Λ(t)/m2(t) for both one- and two-dimensional DKG and DDNLS systems in the same dynamical regime and check whether the curves for different dimensions collapse onto each other; failure to collapse would disprove the dimension-independent scaling law. A complementary test would simulate a 3D disordered nonlinear lattice and compare the measured Lyapunov-exponent decay exponent with the value predicted by Eq. (11) under the same scaling assumption.","supporting_citations":[{"cited_title":"Flach, Chem","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical predictions a_m = 1/(1+2d) and 1/(1+d) for weak and strong chaos that the measured 2D exponents are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides prior 2D DKG weak-chaos results and the 'localization volume' context that the present long-time simulations extend."},{"cited_title":"Garc´ ıa-Mata and D","cited_arxiv_id":null,"evidence_quote":"Earlier 2D DDNLS weak-chaos study that this work extends by using larger systems, more disorder realizations, and MLE computations."},{"cited_title":"Skokos, I","cited_arxiv_id":null,"evidence_quote":"Establishes the 1D MLE and deviation-vector-distribution methodology and provides the 1D chaoticity data used in the proposed scaling law."},{"cited_title":"Senyange, B","cited_arxiv_id":null,"evidence_quote":"Supplies the 1D DKG and DDNLS Lyapunov-exponent decay exponents and the wandering-chaotic-seeds analysis that the 2D results are compared with."},{"cited_title":"Senyange and Ch","cited_arxiv_id":null,"evidence_quote":"Provides the ABA864 symplectic integrator with tangent-map method used for the DKG model's evolution and variational equations."},{"cited_title":"Danieli, B","cited_arxiv_id":null,"evidence_quote":"Provides the s11ABC6 integration scheme used for the DDNLS system's evolution."}],"review_version":1}