{"id":"d58e09f3-cf1d-442b-a3f0-af06a592c3c6","arxiv_id":"1908.01497","paper_version":1,"verdict":"UNVERDICTED","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A pedagogical review of the DNLS equation with Mathematica code for wave, soliton, and stability analysis.","lead":"This paper is a tutorial review of the discrete nonlinear Schrödinger equation, covering plane waves, modulational instability, and discrete solitons. It includes Mathematica code and exercises that could help lecturers teach nonlinear lattice dynamics to advanced students.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (2) states a conservation law for H that is false for Eq. (1); the intersite-coupling term has the wrong sign, so Problem 3 asks students to prove a false statement.","rationale":"The reader's weakest assumption was pedagogical transferability, which is secondary. The more load-bearing issue is internal correctness: Eq. (2) presents a non-conserved quantity as a conserved one, and Problem 3 asks students to prove it analytically. Because the paper is explicitly a teaching resource, a false invariant in a core equation is a substantive defect even though no new research claim is at stake. The dispersion, modulational-instability, and numerical-code components appear standard and usable, so the paper is salvageable with corrections. I therefore recommend CONDITIONAL rather than UNVERDICTED: the text should be revised to correct Eq. (2) and Problem 3, and ideally to reconcile the apparent factor-of-two inconsistency between the continuum equation (8) and the soliton width in Eq. (9).","tokens_in":10034,"tokens_out":38791,"duration_ms":322879,"concrete_test":"At t=0 take the periodic N=4 lattice with β=1, γ=1, (u1,u2,u3,u4)=(1,2i,0,0). From Eq. (1) the derivatives are u1'=-2+i, u2'=-8+i, u3'=-2, u4'=i. Evaluate d/dt of Eq. (2) H = Σ|u_k-u_{k-1}|^2 + ½Σ|u_k|^4: the result is 24, not 0. Repeating with the coupling sign reversed, H' = -Σ|u_k-u_{k-1}|^2 + ½Σ|u_k|^4, gives dH'/dt=0. This settles whether Eq. (2) is a true invariant for Eq. (1).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing weak point is the claimed conserved quantity H in Eq. (2). For Eq. (1) with β0=0, differentiating the proposed H = Σ[β|u_k-u_{k-1}|^2 + (γ/2)|u_k|^4] along the flow gives dH/dt = γ d/dt Σ|u_k|^4, which is nonzero for generic DNLS dynamics; the sign of the coupling term is incorrect. The correct invariant is, up to the conserved term 2βΣ|u_k|^2, H_c = Σ[-β|u_k-u_{k-1}|^2 + (γ/2)|u_k|^4]. This is not cosmetic: on a periodic N=4 lattice with β=1, γ=1 and initial data (1,2i,0,0), direct differentiation using Eq. (1) gives dH/dt=24 for the H of Eq. (2), while the sign-corrected expression has dH/dt=0. Since Problem 3 explicitly instructs students to prove that P and H in Eq. (2) are invariants, the paper asks readers to prove a false statement. The rest of the analytical core (dispersion relation Eq. (4), MI condition Eq. (7), stability analysis) appears standard, but this error undercuts the paper's claim of being a reliable self-contained introduction.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a pedagogical introduction to nonlinear discrete systems, using the discrete nonlinear Schrödinger (DNLS) equation as a universal model. It covers plane-wave dispersion, modulational instability, discrete bright solitons, the anti-continuum limit, numerical construction of stationary solitons, and linear stability analysis, and it supplies a Mathematica script together with twelve problems. The paper is aimed at university courses and claims to provide a self-contained analytical and numerical toolkit that can be extended to other discrete systems.","tokens_in":10286,"tokens_out":22502,"duration_ms":207294,"significance":"If corrected, this would be a useful teaching resource for Eur. J. Phys.: the analytical derivations are standard, the Mathematica code is provided as a working supplement, and the problems extend the methods to impurity modes, surface modes, cubic-quintic, and Ablowitz-Ladik/Salerno models. The modulational-instability condition and the linear stability equations are consistent with the literature, and the code appears to reproduce known results. The paper does not claim new research results; its value lies in accessibility and reproducibility, which are genuine strengths.","major_comments":[{"comment":"The Hamiltonian H in Eq. (2) is not an invariant of Eq. (1). With β0=0, differentiating H = Σ[β|u_k-u_{k-1}|² + (γ/2)|u_k|⁴] along the flow gives a generally nonzero result. For a periodic N=4 chain with β=γ=1 and initial condition (1,2i,0,0), direct differentiation using Eq. (1) gives dH/dt=24. The correct conserved quantity is H' = Σ[β0|u_k|² + β|u_k-u_{k-1}|² - (γ/2)|u_k|⁴]. Because Problem 3 explicitly asks students to prove that P and H in Eq. (2) are invariants, this is a load-bearing error that must be corrected.","section":"Eq. (2), §3 Problem 3"},{"comment":"The continuum limit is misstated by a factor of two. Using the intended phase transformation u_k = w(x_k,t) e^{2iβt} and expanding u_{k+1}+u_{k-1} ≈ 2u_k + h²∂²u/∂x² in Eq. (1) gives i w_t + βh² w_xx + γ|w|²w = 0, not 2βh². Equation (9) is consistent with the coefficient β (for h=1), so Eq. (8) should be corrected, and the missing imaginary unit in the transformation w ≈ u_k exp(−2βt) should also be fixed.","section":"Sec. 2.2, Eq. (8)"}],"minor_comments":[{"comment":"The phrase 'nonlinear distributed systems' should read 'nonlinear discrete systems'; the paper treats discrete lattices, not distributed (continuum) systems.","section":"Abstract and §1"},{"comment":"The statement 'There is no exact DS solutions of the DNLS equation' is too absolute and should be qualified, for example as 'no exact closed-form bright-soliton solutions in the general nonintegrable case', since integrable variants such as the Ablowitz-Ladik equation in Problem 11 do possess exact soliton solutions.","section":"Sec. 2.2, first sentence"},{"comment":"Typographical error: 'staggered solions' should be 'staggered solitons'.","section":"Problem 5, §3"},{"comment":"The soliton width in Eq. (9) is written without the lattice spacing h; if h is retained, the width should be ν = h√(2β)/(A√γ).","section":"Eq. (9), Sec. 2.2"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript fits the journal's teaching-oriented scope. The two technical errors, the sign of the nonlinear term in the Hamiltonian and the factor of two in the continuum limit, are fixable but must be corrected before publication; the rest of the analytical and numerical content appears standard, and the code is a valuable supplement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a teaching paper, not a research contribution, and the main derivations are mostly standard and correct. The one serious problem is Eq. (2): the Hamiltonian as written, with the +γ/2 term, is not conserved by Eq. (1). I checked this by direct differentiation on a simple four-site example with β=γ=1 and initial data (1,2i,0,0); dH/dt comes out to 24, not zero. The sign of the nonlinear term is flipped. That makes Problem 3, which asks students to prove P and H are invariants, a false statement. For a pedagogical paper, this is not cosmetic.\n\nWhat the paper does well: it walks through plane waves, the modulational instability condition, approximate solitons, the anti-continuum limit, and the linear stability analysis, all with a working Mathematica script that reproduces the figures. The derivations are self-contained, the cited literature is appropriate, and the numerical methods (Newton–Raphson, linearized eigenvalues) are standard. For a lecturer who wants a ready-made module on nonlinear discrete systems, this is a solid starting point.\n\nThe other soft spots are minor. The statement that \"there is no exact DS solution\" is a bit loose; there are no simple closed-form solutions, but the paper itself computes exact stationary solutions numerically. The claim that the approach transfers to other discrete systems is plausible but not demonstrated, which is acceptable for a tutorial.\n\nThe dispersion relation, MI condition, and stability equations check out against the literature. The citation pattern is normal; there are no self-citation issues.\n\nBottom line: with the sign in Eq. (2) fixed (the γ/2 should be −γ/2, or equivalently the β term should be −β), this would be a decent teaching resource. As it stands, I would not want students to work through Problem 3 without rewriting it. Because the pedagogical framing is good and the flaw is easy to repair, it deserves a serious referee rather than a desk rejection.","headline":"A useful teaching tutorial on the DNLS equation, but Eq. (2) contains a real sign error that makes Problem 3 ask students to prove a false invariant.","tokens_in":10742,"tokens_out":5851,"would_cite":false,"duration_ms":53458,"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":"This paper proposes a self-contained, code-backed introduction to nonlinear discrete systems built on the DNLS equation, and argues that the same analytical and numerical sequence transfers to other lattice models.","keywords":["discrete nonlinear Schrödinger equation","discrete solitons","modulational instability","anti-continuum limit","linear stability analysis","waveguide arrays","Bose-Einstein condensates","Mathematica teaching code"],"falsifier":"Run the supplied Mathematica code with $\\beta=0.5$, $\\gamma=1$, and $\\eta=-2$: the on-site soliton should stay stationary with maximum imaginary part of the modulation eigenvalues around $10^{-7}$, the inter-site soliton should be unstable, and a noise-modulated plane wave with $\\epsilon=10^{-3}$ should break into soliton-like pulses; failure of any of these outputs would show the tutorial does not deliver its central claim.","tokens_in":9852,"feed_emoji":"⚛️","tokens_out":9458,"duration_ms":82173,"temperature":0.7,"pith_summary":"This paper tries to show that the discrete nonlinear Schrödinger (DNLS) equation is a good universal model for teaching nonlinear discrete systems. It walks through the full analysis pipeline: plane-wave solutions and their dispersion relation, modulational instability, approximate and numerically exact discrete solitons, and linear stability. Every step is implemented in a Mathematica program, so a student can reproduce the figures and then modify the code for other systems. The intended payoff is that the basic toolkit of nonlinear lattice physics becomes accessible without needing specialist monographs.","feed_headline":"The DNLS equation can teach the full nonlinear-lattice toolkit","feed_subtitle":"Paper derives dispersion, instabilities, and soliton stability, then gives Mathematica code for every step.","key_machinery":"The load-bearing object is the discrete nonlinear Schrödinger (DNLS) equation, a nearest-neighbor lattice model with cubic onsite nonlinearity. The argument is carried by a small set of ansätze: the plane wave $u_k=a\\,e^{i(qk-\\omega t)}$ for dispersion, the modulated plane wave $u_k=(a+b_k(t))e^{i(qk-\\omega t)}$ for modulational instability, the stationary mode $u_k=U_k e^{-i\\eta t}$ for solitons, and the perturbed stationary mode $u_k=(U_k+w_k(t))e^{-i\\eta t}$ for stability. Substituting each ansatz and linearizing in the small quantities reduces the problem to algebraic or eigenvalue equations, and the anti-continuum limit provides the localized seed that the Newton-Raphson solver (Mathematica's FindRoot) needs. The stability verdict is read from the complex frequencies of the final eigenvalue problem: any nonzero imaginary part means instability.","core_discovery":"The authors claim that the main ideas and methods of nonlinear discrete-system theory can be presented through the DNLS equation $i\\,du_k/dt + \\beta(u_{k+1}+u_{k-1}) + \\gamma|u_k|^2u_k=0$, and that the analysis sequence is generic. They derive the plane-wave dispersion relation $\\omega = -(2\\beta\\cos q + \\gamma a^2)$, the modulational-instability condition $(\\Omega - 2\\beta\\sin Q\\sin q)^2 = 8\\beta\\sin^2(Q/2)\\cos q\\,[2\\beta\\sin^2(Q/2)\\cos q - \\gamma a^2]$, the stationary soliton equation $-\\beta(U_{k+1}+U_{k-1}) - \\gamma|U_k|^2U_k = \\eta U_k$, and a linear-stability eigenvalue problem for small perturbations. The paper also uses the continuous limit and the anti-continuum limit $\\beta=0$ to construct soliton approximations, and supplies Mathematica code covering every one of these steps.","pith_inferences":["The Mathematica code could be refactored into a reusable template: changing the coupling term and nonlinearity in MC.1 and re-running MC.2–MC.10 would give a working stability analysis for other lattice models without re-deriving the linear algebra.","The eigenvalue formulation suggests a natural extension to two- and three-dimensional arrays by altering the neighbor indexing in the code, a step the paper mentions but does not carry out.","The paper's implicit pedagogical claim is testable: comparing students who work through this DNLS pipeline against students who only see linear lattice theory would show whether the nonlinear toolkit transfers as claimed.","Because the anti-continuum limit provides the numerical seed, the method is especially well suited to strongly localized modes, which is where discrete effects differ most from continuum solitons."],"forward_implications":["A student who follows the paper and runs the Mathematica code can reproduce the dispersion relation, the modulational-instability growth, and the soliton stability spectra for the stated parameters.","The same five-step sequence (plane wave, dispersion, modulational instability, stationary soliton, linear stability) can be applied to other discrete systems, as the paper claims in Section 4 and demonstrates via the suggested problems.","The code is written for nearest-neighbor equations and can be extended to more general couplings, different nonlinearities, and higher-dimensional arrays.","Modulational instability of a nonlinear plane wave leads to break-up into soliton-like pulses, which is how the paper motivates the relevance of discrete solitons."],"supporting_citations":[{"why":"supplies the optical-soliton background and the standard classification of on-site versus inter-site discrete solitons","marker":"[1]"},{"why":"provides the integrable continuous NLS model and the discrete-system context for the continuous-limit approximation","marker":"[2]"},{"why":"supports the waveguide-array realization and the discrete-soliton phenomenology used in the paper","marker":"[4]"},{"why":"is the main DNLS reference, backing the modulational-instability relation, the anti-continuum limit, and the soliton stability results","marker":"[5]"},{"why":"justifies the DNLS model for a Bose-Einstein condensate in an optical lattice by derivation from the Gross-Pitaevskii equation","marker":"[13]"},{"why":"underlies the Newton-Raphson method that the Mathematica FindRoot command uses to solve the stationary soliton equations","marker":"[15]"},{"why":"supports the transformation between unstaggered and staggered solitons and the corresponding sign change of frequency and nonlinearity","marker":"[17]"}],"fun_headline_variants":["DNLS: one equation for all nonlinear lattices","Nonlinear lattices via DNLS, with code and stability","Master nonlinear discrete systems with DNLS","DNLS: plane waves, solitons, and stability in one","A practical intro to nonlinear discrete systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the toolkit demonstrated on the DNLS equation is generic enough to transfer to other discrete systems, which the paper asserts in Section 4 but does not demonstrate in detail.","fun_headline_variants_meta":{"raw":{"variants":["DNLS: one equation for all nonlinear lattices","Nonlinear lattices via DNLS, with code and stability","Master nonlinear discrete systems with DNLS","DNLS: plane waves, solitons, and stability in one","A practical intro to nonlinear discrete systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000176,"raw_usage":{"total_tokens":1275,"prompt_tokens":914,"completion_tokens":361,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":285}},"tokens_in":530,"tokens_out":361,"duration_ms":4369,"temperature":1.0,"reasoning_tokens":285,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:10:31.143719+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the supplied Mathematica code with $\\beta=0.5$, $\\gamma=1$, and $\\eta=-2$: the on-site soliton should stay stationary with maximum imaginary part of the modulation eigenvalues around $10^{-7}$, the inter-site soliton should be unstable, and a noise-modulated plane wave with $\\epsilon=10^{-3}$ should break into soliton-like pulses; failure of any of these outputs would show the tutorial does not deliver its central claim.","supporting_citations":[{"cited_title":"Each site is characterized by some varia bles, so that at a given time these variables specify a state of the system","cited_arxiv_id":null,"evidence_quote":"supplies the optical-soliton background and the standard classification of on-site versus inter-site discrete solitons"},{"cited_title":"Plane waves For discrete systems, it is instructive to start the analysis from a p lane wave solution","cited_arxiv_id":null,"evidence_quote":"provides the integrable continuous NLS model and the discrete-system context for the continuous-limit approximation"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supports the waveguide-array realization and the discrete-soliton phenomenology used in the paper"},{"cited_title":"}; eqn= I*u[k]’[t] + β*(u[k+1][t]+u[k-1][t]) + γ*(u[k][t])∧ 2*Conjugate[u[k][t]]","cited_arxiv_id":null,"evidence_quote":"is the main DNLS reference, backing the modulational-instability relation, the anti-continuum limit, and the soliton stability results"},{"cited_title":"standard","cited_arxiv_id":null,"evidence_quote":"justifies the DNLS model for a Bose-Einstein condensate in an optical lattice by derivation from the Gross-Pitaevskii equation"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"underlies the Newton-Raphson method that the Mathematica FindRoot command uses to solve the stationary soliton equations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supports the transformation between unstaggered and staggered solitons and the corresponding sign change of frequency and nonlinearity"}],"review_version":1}