Pith. sign in

REVIEW 4 major objections 4 minor 40 references

Stability analysis of solutions in the helicoidal Peyrard-Bishop model of DNA molecule

T0 review · 4 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2501.06967 v1 pith:76CRLFCU submitted 2025-01-12 physics.bio-ph

Anna Batova , Dragana Ranković , Slobodan Zdravković This is my paper · ORCID
classification physics.bio-ph MSC 34D2035Q5137K4092C40
keywords helicoidalPeyrard-BishopmodelDNAnonlineardynamicskinksolitonsstabilityanalysisviscositysubsonicsupersoniccontinuumapproximation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

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.

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 (4)
  1. [§3, Eqs. (19)–(23)] 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.
  2. [§3, Eq. (20)] 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.
  3. [§3, Eqs. (24)–(27), Figs. 1–7] 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.
  4. [§3] 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.
minor comments (4)
  1. [Eq. (6)] 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.
  2. [§3 heading] The word "soluitons" should be "solitons."
  3. [§4] The conclusions refer to "MTs" (microtubules), although the paper concerns DNA; this appears to be a leftover from another manuscript and should be corrected.
  4. [Figs. 1–7] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation: stability conclusions are computed from the linearized ODE, not assumed in the ansatz.

full rationale

The paper derives kink solutions from the traveling-wave ODE (6) via the METHF ansatz and cites prior work [10] for these profiles; citing one's own earlier derivation of the solutions is not circular because the stability analysis performed here linearizes the ODE around those profiles and computes eigenvalues (22)-(23) whose sign depends on α. The supersonic/subsonic labels are linked to the sign of α through Eq. (14), but the stability verdict is not put into that linkage; it follows from the algebra of the linearized system. The no-viscosity claim is obtained by numerical integration of system (27) for selected α values; while this is an inductive generalization that may be too strong as a scientific conclusion, it is not a reduction of the conclusion to the input data. No fitted parameter is renamed as a prediction, no uniqueness claim is imported from the authors' prior work, and the central stability classification retains independent content beyond the cited solutions. Any concern about ODE-in-ξ stability versus PDE temporal stability is a correctness issue, not a circularity issue.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The paper adds no new entities or forces. It relies on the continuum limit, a solution ansatz, and a simplified constant-coefficient linearization for the viscous case. The no-viscosity conclusion depends on a small set of numerical experiments.

free parameters (3)
  • a2^(1), a2^(2) (solution amplitudes in no-viscosity case) = free (Eq. 15)
    In the no-viscosity case, the parameters a2^(1) and a2^(2) are left free in the METHF ansatz; the specific solutions (16) and (17) used for stability analysis correspond to particular values, but the stability conclusion is claimed for the family.
  • alpha values in phase portraits = 0.8, 1.5, 0.5, -0.8
    The no-viscosity numerical stability analysis is carried out at a few hand-picked values of alpha (Figs. 1-7), not over a range or analytically.
  • Initial perturbation amplitudes X0, Y0 = 10^-4, 10^-2, -10^-4
    The conclusion that perturbations grow is based on selected initial conditions; one sign leads to periodic orbits, the other to blow-up.
assumptions (4)
  • domain assumption Continuum approximation of the discrete lattice (Eq. 4)
    The discrete DNA chain is replaced by a continuous field y(x,t); validity for the short lattice spacing is assumed.
  • ad hoc to paper METHF ansatz: solution is a finite series in tanh (Eq. 8)
    The method assumes the traveling wave is a polynomial in tanh; this restricts the class of solutions and is not derived from the equation.
  • ad hoc to paper Frozen-coefficient linearization (Eq. 20)
    The linearized perturbation equation for psi1 is treated with a constant coefficient f/4, ignoring the xi-dependence of the background kink; eigenvalues (22) follow from this simplification.
  • domain assumption Numerical phase portraits of the nonlinear perturbation ODE (27) determine stability
    The paper assumes that boundedness/escape of solutions in the phase plane for a few parameter values establishes stability/instability of the original wave.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stability analysis of solutions in the helicoidal Peyrard-Bishop model of DNA molecule." pith.science (2026). https://pith.science/paper/76CRLFCU

@misc{pith2026250106967,
  author       = {Pith},
  title        = {Pith review of: Stability analysis of solutions in the helicoidal Peyrard-Bishop model of DNA molecule},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/76CRLFCU}},
  note         = {Machine review of arXiv:2501.06967}
}
read the original abstract

We use the helicoidal Peyrard-Bishop model of DNA in the current work. We solve a dynamical equation of motion using a continuum approximation, resulting in kink-solitary waves that travel along the chain. We demonstrate that, whereas supersonic kink solitons are not stable, subsonic ones are. Moreover, we demonstrate the importance of viscosity by showing that no wave is stable in the absence of viscosity.

Figures

Figures reproduced from arXiv: 2501.06967 by the authors.

Figure 1
Figure 1. The functions ( ) X1  (a), ( ) Y1  (b), and ( ) Y1 X1 (c), for   0.8 and 4 0 0 10 X  Y  [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. The functions ( ) X1  (a), ( ) Y1  (b), and ( ) Y1 X1 (c), for   0.8 and X0  Y0  0.01 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [10]

    Zdravković, D

    S. Zdravković, D. Chevizovich, A. N. Bugay, and A. Maluckov, Chaos 29 (2019) 053118

    work page 2019

  2. [1]

    L. V. Yakushevich, Nonlinear Physics of DNA , Wiley Series in Nonlinear Science, John Wiley, Chichester 1998

    work page 1998

  3. [2]

    Gaeta, C

    G. Gaeta, C. Reiss, M. Peyrard, and T. Dauxois, Riv. Nuovo Cimento 17 (1994) 1

    work page 1994

  4. [3]

    Zdravković, Nonlinea r Dynamics of DNA Chain, In Nonlinear Dynamics of Nanobiophysics , Edited by S

    S. Zdravković, Nonlinea r Dynamics of DNA Chain, In Nonlinear Dynamics of Nanobiophysics , Edited by S. Zdravković and D. Chevizovich, Springer, 2022, 29-66

    work page 2022

  5. [4]

    S. W. Englander, N. R. Kallenbach, A. J. Heeger, J. A. Krumhansl, and S. Litwin, Proc. Natl. Acad. Sci. (USA) 777 (1980) 7222

    work page 1980

  6. [5]

    Dauxois, Phys

    T. Dauxois, Phys. Lett. A 159 (1991) 390

    work page 1991

  7. [6]

    Zdravković, J

    S. Zdravković, J. Nonlin Math. Phys. 18, Suppl. 2 (2011) 463

    work page 2011

  8. [7]

    N. K. Voulgarakis, A. Redondo, A. R. Bishop, and K. Ø. Rasmussen, Phys. Rev. Lett. 96 (2006) 248101

    work page 2006

Show all 40 references
  1. [8]

    C. L. Gninzanlong, F. Th. Ndjomatchoua, and C. Tchawoua, Phys. Rev. E 102 (2020) 052212

    2020

  2. [9]

    R. P. A. Lima and A. V. Malyshev, Phys. Rev. E 106 (2022) 024414

    2022

  3. [11]

    Peyrard and A

    M. Peyrard and A. R. Bishop, Phys. Rev. Lett. 62 (1989) 2755

    1989

  4. [12]

    Dauxois, M

    T. Dauxois, M. Peyrard , and A. R. Bishop, Phys. Rev. E 47 (1993) R44

    1993

  5. [13]

    B. S. Alexandrov, V. G. Stanev, A. R. Bishop, and K. Ø. Rasmussen, Phys. Rev. E 86 (2012) 061913

    2012

  6. [14]

    Hillebrand, G

    M. Hillebrand, G. Kalosakas, A. Schwellnus, and Ch. Skokos, Phys. Rev. E 99 (2019) 022213

    2019

  7. [15]

    T. R. Strick, M. N. Dessinges, G. Charvin, N. H. Dekker, J. F. Allemand, D. Bensimon, and V. Croquette, Rep. Prog. Phys. 66 (2003) 1

    2003

  8. [16]

    Zdravković , M

    S. Zdravković , M. V. Satarić , A. Yu. Parkhomenko , and A. N. Bugay, Chaos 28 (2018) 113103

    2018

  9. [17]

    M. V. Satarić, J. A. Tuszyńsky, and R. B. Žakula, Phys. Rev. E 48 (1993) 589

    1993

  10. [18]

    Vasumathi and M

    V. Vasumathi and M. Daniel, Phys. Rev. E 80 (2009) 061904

    2009

  11. [19]

    J. B. Okaly, A. Mvogo, R.L. Woulaché, and T. C. Kofané, Commun. Nonlinear Sci. Numer. Simulat. 55 (2018) 183

    2018

  12. [20]

    Chevizovich, D

    D. Chevizovich, D. Michieletto, A. Mvogo, F. Zakiryanov, and S. Zdravković, R. Soc. Open Sci. 7 (2020) 200774

    2020

  13. [21]

    J. A. D. Wattis, Phys. Rev. E 105 (2022) 054212

    2022

  14. [22]

    Zdravković, M

    S. Zdravković, M. V. Satarić, and Lj. Hadžievski, Chaos 20 (2010) 043141

    2010

  15. [23]

    Zdravković and M

    S. Zdravković and M. V. Satarić, Europhys. Lett. 78 (2007) 38004

    2007

  16. [24]

    Zdravković and M

    S. Zdravković and M. V. Satarić, Europhys. Lett. 80 (2007) 38003

    2007

  17. [25]

    N. A. Kudryashov, Phys. Lett. A 342 (2005) 99

    2005

  18. [26]

    N. A. Kudryashov, Chaos Solitons Fract. 24 (2005) 1217

    2005

  19. [27]

    Zdravković and G

    S. Zdravković and G. Gligorić, Chaos 26 (2016) 063101

    2016

  20. [28]

    N. A. Kudryashov and N. B. Loguinova, Appl. Math. Comput . 205 (2008) 396

    2008

  21. [29]

    Dai and J

    C. Dai and J. Zhang, Chaos Solitons Fract. 27 (2006) 1042

    2006

  22. [30]

    Zeković, S

    S. Zeković, S. Zdravković, L. Kavitha, and A. Muniyappan, Chin. Phys. B 23 (2014) 020504

    2014

  23. [31]

    Cornejo -Pérez, J

    O. Cornejo -Pérez, J. Negro, L. M. Nieto, and H. C. Rosu, Found. Phys. 36 (2006) 1587

    2006

  24. [32]

    W. Alka, A. Goyal, and C. N. Kumar, Phys. Lett. A 375 (2011) 480

    2011

  25. [33]

    M. N. Alam, M. G. Hafez, M. A. Akbar, and H. -O. Roshid, J. Sci. Res. 7 (2015) 1

    2015

  26. [34]

    M. N. Alam and F. B. M. Belgacem, Mathematics 4 (2016) 1

    2016

  27. [35]

    Gordon, Physica B 146 (1987) 373

    A. Gordon, Physica B 146 (1987) 373

    1987

  28. [36]

    Fan, Phys

    E. Fan, Phys. Lett. A 277 (2000) 212

    2000

  29. [37]

    S. A. El -Wakil and M. A. Abdou, Chaos Solitons Fract. 31 (2007) 840

    2007

  30. [38]

    A. H. A. Ali, Phys. Lett. A 363 (2007) 420

    2007

  31. [39]

    A. J. M. Jawad, M. D. Petković, and A. Biswas, Appl. Math. Comput. 217 (2010) 869

    2010

  32. [40]

    Zdravković and S

    S. Zdravković and S. Zeković, Chin. J. Phys. 55 (2017) 2400

    2017

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.