{"id":"f5d262e5-d3f3-4f35-bdf4-b707c6c8a9dd","arxiv_id":"2501.06967","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the helicoidal Peyrard-Bishop DNA model, kink solitons are stable only when viscosity is present and the wave travels below the sound speed.","lead":"This paper analyzes the stability of kink-shaped solitary waves in a DNA dynamics model. It reports that stable waves require both viscosity and subsonic speed, but the broader claim that no stable wave exists without viscosity is not proven.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The stability conclusion is based on eigenvalues of the travelling-wave ODE in ξ, not on the spectrum of the linearized time evolution; this does not establish temporal stability even in the viscous case.","rationale":"I read the paper as claiming physical, temporal stability of kink solitons in the helicoidal Peyrard–Bishop model. The derivation of the travelling-wave solutions from Eq. (6) is a routine application of the METHF method and appears internally consistent. The decisive weakness is earlier than the reader located it: the stability test in §3 substitutes perturbations into the ODE (6), whose independent variable is ξ, and then treats the sign of Re λ in the ξ-dynamics as the stability criterion. This does not settle whether perturbations of the PDE (5) grow or decay in time. Because Eq. (5) contains y_tt, a legitimate proof must produce the spectrum of the linearized time-evolution operator in a moving frame; no such spectrum is computed. The variable coefficient α+2ψ_i(ξ) makes the constant eigenvalues (22)–(23) especially problematic, since they are at best asymptotic/ local rather than global. The reader's identified overgeneralization for ρ=0 is real and remains a problem, but it is secondary to the conceptual mismatch between ODE-in-ξ stability and PDE-in-time stability. The REJECT verdict is therefore appropriate, and I would strengthen the grounds by pointing to §3's stability notion rather than only to the no-viscosity extrapolation.","tokens_in":8121,"tokens_out":8934,"duration_ms":95217,"concrete_test":"Perform a temporal linear stability analysis: set y=ψ_i(ξ)+ε e^{λτ}g(ξ) with ξ=x−Vt in Eq. (5), derive the eigenvalue problem for λ, and compute its spectrum numerically (e.g. by Chebyshev collocation or an Evans-function method) for the parameter values used in Figs. 1–7, including ρ=0. If the subsonic kink ψ_2 has any eigenvalue with Re λ>0, or if the ρ=0 spectrum contains unstable modes, the central claim fails; if all Re λ<0 for ρ>0, the viscous claim survives and the remaining objection is the unsupported no-viscosity generalization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing flaw is that the stability analysis in §3 does not analyze temporal stability of the solitary waves. The authors substitute ψ_i+f into the travelling-wave ODE (6), not into the time-dependent PDE (5), and study the resulting linear system (19) in the independent variable ξ. The eigenvalues (22)–(23) therefore describe decay or growth of profile perturbations along ξ in the ODE, not decay in time under the PDE. For a second-order-in-time PDE, temporal stability requires linearizing (5) in a moving frame, e.g. y=ψ_i(ξ)+ε e^{λτ}g(ξ), and computing the spectrum in λ; the ODE eigenvalues are not that spectrum. Moreover, the linearized coefficient α+2ψ_i(ξ) depends on ξ, so the constant eigenvalues (22)–(23) are only asymptotic or local and are not justified as global stability criteria. Separately, the no-viscosity conclusion in §3 is extrapolated from numerical integration of (27) for only four values of α and a few initial conditions to the universal statement that no wave is stable. Both deficiencies bear directly on the abstract's headline claims.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the helicoidal Peyrard-Bishop model of DNA and claims to demonstrate, from a continuum approximation, that traveling kink solitons are stable only when they are subsonic and viscosity is present, and that no solitary wave is stable in the absence of viscosity. The solutions themselves are taken from the authors' prior work [10] and are presented as explicit tanh-type kinks. The stability analysis in Section 3 proceeds by substituting a perturbed solution into the traveling-wave ODE, deriving a linear system, and computing eigenvalues; for the no-viscosity case, the paper integrates a nonlinear perturbation system numerically for a few parameter values and initial conditions. The conclusions state that viscosity is essential for stable solitary waves.","tokens_in":8417,"tokens_out":5942,"duration_ms":56318,"significance":"If the claims were established, they would be relevant to the nonlinear dynamics of DNA and to the broader literature on solitary waves in damped lattice models. The paper has the virtue of working with explicit analytic traveling-wave solutions and of making a falsifiable prediction about the role of viscosity. However, the stability argument is not connected to the temporal evolution of the original PDE, and the no-viscosity conclusion is supported only by a small number of numerical phase portraits. The paper does not provide machine-checked proofs or reproducible code, and the central stability claims are therefore not currently supported.","major_comments":[{"comment":"The stability analysis is carried out on the traveling-wave ODE (6), not on the time-dependent PDE (5). Substituting ψ_i + f into Eq. (6) and computing eigenvalues λ of the linearized system (19) yields rates of growth or decay in the moving coordinate ξ, not in time. For the second-order-in-time PDE (5), temporal stability requires linearizing in a moving frame, for example y = ψ_i(ξ) + e^{λτ} g(ξ), and examining the spectrum of the resulting operator in λ. The eigenvalues (22)–(23) are not that spectrum. Consequently, the abstract's central claim that subsonic kinks are stable and supersonic ones are not is not established by the arguments in the paper.","section":"§3, Eqs. (19)–(23)"},{"comment":"The linearized coefficient depends on ξ through ψ_i(ξ), yet the calculation leading to (22)–(23) replaces ψ_1 by the constant value 1/4, which is its value at ξ = 0. This is a frozen-coefficient approximation of a non-autonomous linear system. Even within the ODE framework, the constant eigenvalues are only a local approximation and do not determine the stability of the full non-autonomous linearization. The paper does not justify this freezing step, nor does it provide Floquet, asymptotic, or spectral arguments for the actual variable-coefficient system.","section":"§3, Eq. (20)"},{"comment":"The claim that no solution is stable in the absence of viscosity is inferred from numerical phase portraits of system (27) at only four values of α (0.8, 1.5, 0.5, and −0.8) and a few initial conditions. No error-controlled integration, systematic parameter sweep, Lyapunov function, or analytical stability criterion is provided. A phase portrait whose trajectory does not encircle the origin, or that blows up for one negative initial condition at ξ ≈ 12, is not a proof of instability of the original solitary wave; bounded non-periodic orbits or slow drift could produce similar portraits. The universal statement in the abstract therefore goes beyond the evidence presented.","section":"§3, Eqs. (24)–(27), Figs. 1–7"},{"comment":"The paper never states which notion of stability is meant: Lyapunov stability, asymptotic stability, or orbital stability. For solitary waves, orbital stability with respect to the PDE dynamics is the standard notion, and one must account for neutral modes associated with translation invariance. The eigenvalue calculation in ξ arguably describes asymptotic behavior of a fixed point of the spatial ODE, which is a different object. The manuscript should define the stability notion and connect the calculation to that definition.","section":"§3"}],"minor_comments":[{"comment":"Equation (6) is garbled in the submitted text; the displayed formula appears to be missing explicit terms, which makes it difficult to verify the subsequent linearization.","section":"Eq. (6)"},{"comment":"The word \"soluitons\" should be \"solitons.\"","section":"§3 heading"},{"comment":"The conclusions refer to \"MTs\" (microtubules), although the paper concerns DNA; this appears to be a leftover from another manuscript and should be corrected.","section":"§4"},{"comment":"The figures are only represented by captions in the text provided; ensure that the actual phase portraits are included and legible, with axes labeled and initial conditions indicated.","section":"Figs. 1–7"}],"recommendation":"reject","confidential_remarks":"The manuscript has a clear model and explicit analytic traveling-wave solutions, but the stability analysis does not address temporal stability of the PDE solutions. The eigenvalue calculation in §3 concerns the spatial ODE, and the no-viscosity claim rests on a small number of numerical phase portraits. These are load-bearing flaws for the stated conclusions. If the authors reframed the paper as a study of the linearized spatial dynamics, or performed a genuine spectral or Lyapunov stability analysis of the time-dependent problem, a resubmission could be considered, but as it stands the central claims are not supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper's central stability conclusions do not hold. The analysis in §3 linearizes the traveling-wave ODE, not the time-dependent PDE, so the eigenvalues in ξ say nothing about temporal stability. The no-viscosity claim is also overgeneralized from a few numerical runs.\n\nWhat's good: the paper is concise and the derivation of the two kink solutions from Eq. (6) is clean. The identification of α with V²-c² is a useful way to organize the subsonic/supersonic cases. If the stability claim were properly established, the subsonic-plus-viscosity condition would be a practical guide for biophysical modeling.\n\nThe problem: substituting ψ_i+f into the ODE (6) and studying the linear system in ξ gives information about how a profile perturbation grows as you move along the traveling wave coordinate, not about how the solution evolves in time under Eq. (5). For a second-order-in-time PDE, you need to linearize around the traveling wave in a moving frame and look for eigenvalues λ with e^{λt}. The constant eigenvalues (22)-(23) are at best local/asymptotic because the linearized coefficient α+2ψ_i(ξ) depends on ξ; they are not a global stability criterion. This is a load-bearing flaw because the abstract's claims rest on it.\n\nThe no-viscosity conclusion is even thinner. The system (27) is integrated for four values of α and a few initial conditions, and from that the authors state 'no wave is stable in the absence of viscosity.' That is a universal claim from three or four phase portraits. No Lyapunov function, no spectral argument.\n\nMinor: the conclusions mention 'MTs' (microtubules) in a DNA paper, which is a copy-paste error.\n\nBottom line: this is not ready for review. The stability analysis would need to be redone with the proper temporal linearization. The paper's value is limited to the re-derivation of known solutions. I wouldn't send it to a referee; desk reject.\n\nWho would get value? Someone working on HPB solitons might want to see the explicit kink forms, but they are already in Ref. [10]. The paper teaches more about how not to do stability analysis.","headline":"The paper's stability conclusion is not supported: the analysis is of the traveling-wave ODE, not the time-dependent PDE, and the no-viscosity claim is overgeneralized from a few numerics.","tokens_in":8884,"tokens_out":7010,"would_cite":false,"duration_ms":63278,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34D20","35Q51","37K40","92C40"],"pacs":[],"model":"deepseek-v4-flash","headline":"The helicoidal Peyrard-Bishop model of DNA has stable solitary waves only as subsonic kinks in a viscous medium; supersonic kinks and all inviscid waves are unstable.","keywords":["helicoidal Peyrard-Bishop model","DNA nonlinear dynamics","kink solitons","stability analysis","viscosity","subsonic solitons","supersonic solitons","continuum approximation"],"falsifier":"Numerically integrate system (27) with $\\rho=0$, $\\alpha=0.3$, and initial condition $(X_0,Y_0)=(10^{-3},10^{-3})$ over a much longer interval than those shown in the paper. If $(X,Y)$ stays bounded and its phase curve closes around $(0,0)$ instead of drifting off or diverging at finite $\\xi$, then the inviscid solution $\\psi_{10}$ would be stable, contradicting the paper's central claim.","tokens_in":7949,"feed_emoji":"🧬","tokens_out":9021,"duration_ms":81772,"temperature":0.7,"pith_summary":"The paper studies the helicoidal Peyrard-Bishop model of DNA, takes its continuum limit, and asks which traveling kink-soliton solutions are stable against small perturbations. Its central claims are that with viscosity the subsonic kink ($\\alpha<0$) is stable while the supersonic kink ($\\alpha>0$) is unstable, and that without viscosity no solution of the model is stable. The viscous result follows from the eigenvalues of the linearized perturbation system; the inviscid result comes from numerical integration of the nonlinear perturbation equations, since the linear eigenvalues are purely imaginary. If these claims hold, an effective dissipation from the surrounding medium is not a nuisance but an essential condition for solitary-wave transport in DNA.","feed_headline":"Only subsonic DNA kinks survive, and only with viscosity","feed_subtitle":"In the helicoidal Peyrard-Bishop model, faster-than-sound kinks and all inviscid waves decay under perturbation.","key_machinery":"The central object is the reduced traveling-wave equation (6), a damped anharmonic oscillator in the moving-coordinate variable $\\xi$, whose coefficients $\\alpha$ and $\\rho$ encode soliton speed, sound speed, and viscosity. Stability is decided by writing $\\psi=\\psi_i+f$, linearizing the resulting system (19) about each kink, and computing eigenvalues of the associated Jacobian: a positive real part means instability, so the supersonic kink fails while the subsonic kink survives. In the inviscid case the eigenvalues from this linearization are purely imaginary, and the argument switches to the nonlinear perturbation system (27), whose numerically integrated orbits are used to judge stability.","core_discovery":"In the continuum approximation the HPB equation of motion reduces to a damped nonlinear oscillator for the traveling-wave profile. The modified extended tanh-function method yields two kink solutions when viscosity is present: one with $\\alpha>0$ that travels faster than the linear sound speed and one with $\\alpha<0$ that is subsonic. Linearizing the perturbation equation around these kinks gives eigenvalues whose real parts show the supersonic kink to be unstable and the subsonic kink to be stable. In the inviscid limit the same linearization gives purely imaginary eigenvalues, so stability has to be decided at nonlinear order; numerical solutions of the perturbation system for a few values of $\\alpha$ and several initial conditions show orbits that are either non-recurrent around the origin or blow up at finite $\\xi$. The authors therefore conclude that viscosity selects the subsonic kink as the only stable solitary wave in this model.","pith_inferences":["The no-viscosity instability is inferred from a handful of numerical runs rather than from a Lyapunov argument; a natural next step would be to test system (27) at many $\\alpha$ values and random initial conditions to see whether the 'no stable wave' conclusion is generic.","By the same eigenvalue logic, the threshold $\\alpha=0$ ($V=c$) should mark a sharp transition; the paper does not discuss what happens exactly at the sound speed, where the linearized eigenvalues become degenerate.","If confirmed in the discrete HPB chain, the result suggests that thermal fluctuations plus viscosity could act as a filter favoring subsonic kinks, which may matter for DNA-RNA transcription bubbles and other local-opening phenomena.","The viscous stabilization mechanism may be generic for damped anharmonic chains, so an analogous stability split between subsonic and supersonic solitons may appear in other Peyrard-Bishop-type models."],"forward_implications":["In the viscous HPB model, a supersonic kink launched into the chain will not persist; small perturbations grow, so only slower-than-sound solitary waves can carry energy coherently.","A purely conservative (inviscid) HPB chain cannot support a stable dark or bright soliton of the form (16)-(17); some damping mechanism is required for stability.","The stability boundary can be rephrased as a speed cutoff at the linear sound speed $c$: stable propagation requires $V<c$ when $\\rho>0$.","The result gives a concrete role to viscosity in DNA biophysics: local openings modeled by kinks are stable only if the environment dissipates energy."],"supporting_citations":[{"why":"Supplies the continuum equation of motion, the parameter definitions, and the kink and bright/dark soliton solutions that are the subjects of the stability analysis.","marker":"[10]"},{"why":"Introduces the helicoidal Peyrard-Bishop Hamiltonian whose helicoidal coupling and Morse potential give the equation of motion.","marker":"[5]"},{"why":"The original Peyrard-Bishop model on which the HPB model is built; it supplies the harmonic and Morse potential structure.","marker":"[11]"},{"why":"Adds the viscosity force to the DNA soliton equation, producing the damping parameter $\\rho$ used throughout the stability analysis.","marker":"[17]"},{"why":"Another reference for the viscosity term in the same class of DNA and microtubule soliton models.","marker":"[18]"},{"why":"The modified extended tanh-function method used to solve the reduced ordinary differential equation and obtain the two kink solutions.","marker":"[36]"}],"fun_headline_variants":["Viscosity is the gatekeeper: only slow kinks survive","Supersonic kinks unstable, inviscid waves doomed","DNA kinks: only subsonic with viscosity are stable","No viscosity, no stability: subsonic kinks win","Viscosity decides DNA kink fate: slow wins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The sweeping conclusion that no wave is stable in the absence of viscosity rests on numerical integration of the nonlinear perturbation equations for only a few values of the speed parameter and a few small initial perturbations, with no Lyapunov function or spectral proof to back the generalization.","fun_headline_variants_meta":{"raw":{"variants":["Viscosity is the gatekeeper: only slow kinks survive","Supersonic kinks unstable, inviscid waves doomed","DNA kinks: only subsonic with viscosity are stable","No viscosity, no stability: subsonic kinks win","Viscosity decides DNA kink fate: slow wins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":2849,"prompt_tokens":796,"completion_tokens":2053,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":412,"completion_tokens_details":{"reasoning_tokens":1967}},"tokens_in":412,"tokens_out":2053,"duration_ms":14024,"temperature":1.0,"reasoning_tokens":1967,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:50:06.698554+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate system (27) with $\\rho=0$, $\\alpha=0.3$, and initial condition $(X_0,Y_0)=(10^{-3},10^{-3})$ over a much longer interval than those shown in the paper. If $(X,Y)$ stays bounded and its phase curve closes around $(0,0)$ instead of drifting off or diverging at finite $\\xi$, then the inviscid solution $\\psi_{10}$ would be stable, contradicting the paper's central claim.","supporting_citations":[{"cited_title":"Zdravković, D","cited_arxiv_id":null,"evidence_quote":"Supplies the continuum equation of motion, the parameter definitions, and the kink and bright/dark soliton solutions that are the subjects of the stability analysis."},{"cited_title":"Dauxois, Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the helicoidal Peyrard-Bishop Hamiltonian whose helicoidal coupling and Morse potential give the equation of motion."},{"cited_title":"Peyrard and A","cited_arxiv_id":null,"evidence_quote":"The original Peyrard-Bishop model on which the HPB model is built; it supplies the harmonic and Morse potential structure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Adds the viscosity force to the DNA soliton equation, producing the damping parameter $\\rho$ used throughout the stability analysis."},{"cited_title":"Vasumathi and M","cited_arxiv_id":null,"evidence_quote":"Another reference for the viscosity term in the same class of DNA and microtubule soliton models."},{"cited_title":"Fan, Phys","cited_arxiv_id":null,"evidence_quote":"The modified extended tanh-function method used to solve the reduced ordinary differential equation and obtain the two kink solutions."}],"review_version":1}