{"id":"8b97ac2a-8116-4108-b196-58a0f79b8cb5","arxiv_id":"2607.20509","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Exact Lambert W kink solutions of an extended nerve-membrane wave equation admit supersymmetric partner kinks via reversed factorization order.","lead":"This paper adds a supersymmetric quantum mechanics style pairing to exact Lambert W kink solutions of an extended Heimburg Jackson model of nerve pulses. The result is a way to generate partner solitons with modified width and amplitude, proposed as a mathematical description of mechanically perturbed nerve membranes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (6)-(7) do not map Eq. (4) to Eq. (5): the stated rescaling yields +y^5 and an extra s-tilde^{1/4} factor in the linear term, so the Lambert W kinks are not verified as solutions of Eq. (4) for the plotted parameters.","rationale":"The reader identified the zero-integration-constant step as the weakest assumption. That concern is legitimate but secondary: for bounded kink solutions the xi-linear constant must vanish, and a constant remainder can often be absorbed by shifting the background density. The more load-bearing defect is the rescaling in Eqs. (6)-(7): a direct two-fold integration of Eq. (4) gives a quintic coefficient of +1 under the stated s-tilde, while the entire factorization construction is written for -y^5. The Lambert W solutions and their SUSY partners do appear to satisfy the abstract Liénard equation Eq. (17) and the displayed partner equations, which is real evidence for the abstract ODE part of the paper. However, the claimed bridge to the membrane model is unsubstantiated unless Eqs. (5)/(7) are re-derived and the figures are recomputed with consistent physical parameter signs. I would keep the CONDITIONAL verdict, but the condition should be a corrected traveling-wave reduction and a residual check of the exact solutions against Eq. (4), rather than only a discussion of integration constants. There is no machine-checked proof and no reproducible code in the submission, so this direct algebraic check is the appropriate arbiter.","tokens_in":16652,"tokens_out":46697,"duration_ms":383101,"concrete_test":"Use a CAS to integrate Eq. (4) twice with xi = k zeta - v t and zero integration constants, then apply z = y/s-tilde^{1/4} with s-tilde = s/(5 Lambda); compare the resulting ODE with Eq. (5). Then substitute Eq. (27) for the Fig. 4 parameters p=q=300, s=27, r=-141, k=2, delta=v=1 into both the actual reduced ODE and Eq. (5) and print the residuals. If the residual against Eq. (5) is zero while the residual against the actual reduced ODE is not, the reduction is the failure point. Also repeat with s-tilde = -s/(5 Lambda) and s<0 to test whether a sign-corrected mapping exists.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central reduction is algebraically inconsistent as printed. Integrating Eq. (4) twice with the traveling-wave ansatz and zero constants gives z'' + (gamma v / Lambda) z' + [(k^2 - v^2)/(k^2 Lambda)] z + (p/(2 Lambda)) z^2 + (q/(3 Lambda)) z^3 + (r/(4 Lambda)) z^4 + (s/(5 Lambda)) z^5 = 0. With z = y / s-tilde^{1/4} and s-tilde = s/(5 Lambda) from Eq. (7), the y^5 coefficient becomes s/(5 Lambda s-tilde) = 1, while Eq. (5) and the factorization Eq. (17) both use -y^5. The linear coefficient also acquires a factor s-tilde^{1/4}, so it does not match a_1 = (k^2 - v^2)/(Lambda k^2) in Eq. (7). For the plotted parameters p=q=300, s=27, k=2, delta=v=1, Lambda=3 > 0 and s>0, so no real choice of s-tilde makes the quintic term -1. Thus Eq. (27) and its SUSY partner Eq. (35c) are verified only for the abstract Liénard equation Eq. (17), not for the stated reduction of the extended Heimburg-Jackson model Eq. (4). This is a checkable algebraic mismatch, not a matter of biological interpretation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies an extended Heimburg–Jackson (HJ) model with higher-order polynomial nonlinearities, Eq. (1)/(4). Under a traveling-wave reduction, the authors claim to obtain the Liénard-type equation Eq. (5), which they factorize to construct exact Lambert W-kink solitons, Eq. (27), and supersymmetric partner solutions, Eq. (35c). The paper also performs a Painlevé analysis, derives effective potentials, and discusses biological interpretations related to traumatic brain injury.","tokens_in":17141,"tokens_out":26997,"duration_ms":202308,"significance":"If the claimed reduction were correct, the paper would establish a novel connection between nonlinear electromechanical waves in biological membranes and supersymmetric quantum mechanics, providing exact solutions for a non-integrable extension of a physically relevant model. The factorization technique is applied carefully, and the internal algebra of the Liénard equation appears consistent, including the same-damping partner construction behind Eq. (36). However, the central link to the physical model is broken by an algebraic error in the reduction, so the advertised physical significance is not established.","major_comments":[{"comment":"The traveling-wave reduction of Eq. (4) is algebraically incorrect. Substituting z(ξ), ξ = kζ - v t̃, integrating twice with zero constants, and rescaling z = y/s̃^{1/4} with s̃ = s/(5Λ) yields y'' + (γv/Λ)y' - a1 y - a2 y² - a3 y³ - a4 y⁴ - y⁵ = 0, where a1 = (k² - v²)/(Λk²), a2 = p/(2Λs̃^{1/4}), a3 = q/(3Λs̃^{1/2}), a4 = r/(4Λs̃^{3/4}). Eq. (5) instead has +a2 y². Thus the sign of the quadratic term is opposite to that obtained from Eq. (4). Consequently, the Lambert W-kink (27) and its SUSY partner (35c) are verified only for the abstract Liénard equation (17), not for the stated reduction of the extended HJ model. This invalidates the central claim that these are exact traveling-wave solutions of the physical model.","section":"§2.1, Eqs. (5)–(7)"},{"comment":"The assertion that setting both integration constants to zero is 'without loss of generality' is not justified. The linear-in-ξ constant must vanish for any bounded solution, but the additive constant from the second integration cannot generally be set to zero: it selects the background state and restricts the admissible asymptotic values of the kink. As a result, even if the sign issue were corrected, the derived exact solutions would cover only a special subfamily of traveling waves, not the full physical family.","section":"§2.1, paragraph after Eq. (5)"},{"comment":"Because the reduction error changes the sign in front of a2, the coefficient-matching conditions (18a)–(18d), the cubic (19), and the subsequent explicit forms of α, A, and B are all tied to the incorrect Eq. (5). The figures and parameter choices (e.g., p=q=300, s=27, r=-141, k=2, δ=v=1) therefore illustrate solutions of a Liénard equation that is not the reduction of Eq. (4). The paper should either re-derive the analysis for the correct sign or demonstrate explicitly that the error does not affect the existence of Lambert W kinks for physically admissible parameters; as presented, the connection to the Heimburg–Jackson model is unsupported.","section":"§3.2–§3.3, Eqs. (18)–(27), (35a)–(35c)"}],"minor_comments":[{"comment":"The phrase 'third and fourth order nonlinearities' is misleading: the model (1) contains nonlinearities up to fourth power in the density, not third and fourth order derivatives. Please clarify the terminology.","section":"Abstract and §1"},{"comment":"The notation y^{(1)}_{2=}(ξ) and the switching between subscripts >, <, and = is difficult to follow; please add a table or a clearer explanation of the notation in the captions.","section":"§3.2.2, Eq. (32c) and Fig. 7"},{"comment":"The text says the partner equation has 'the same damping coefficient' but Eq. (36) shows ∓γ̃1; the sign convention should be explained, and Figure 4 should be cross-referenced with the sign choices in Eqs. (27) and (35c).","section":"§3.3.1, Eq. (36)"},{"comment":"There are several typographical errors, e.g., 'explicity' in §3.2.2 and inconsistent spacing of 'Lambert W' and 'LambertW'. A careful proofreading would improve readability.","section":"General"}],"recommendation":"reject","confidential_remarks":"The paper presents a potentially interesting factorization analysis of a Liénard equation, but the advertised physical connection to the extended Heimburg–Jackson model is invalidated by a sign error in the traveling-wave reduction, Eq. (5). This is not a local typo: the sign propagates through the matching conditions, the cubic for α, and the explicit solutions. The 'without loss of generality' claim about the integration constants is also problematic. In my view, the manuscript cannot be accepted without re-deriving the core results for the correct reduced equation, which may change the parameter regimes and the existence of the Lambert W solutions. I recommend rejection, though a substantially revised version addressing these points could be reconsidered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nShort version: the Lambert W kinks and their supersymmetric partners are legitimate solutions of the abstract Liénard equation (17), and the partner construction is new, but the claimed reduction from the extended Heimburg–Jackson model (4) to Eq. (5) does not survive a sign check. The stress-test note's specific algebra is partly wrong, but its conclusion is right for a different reason.\n\nWhat is good: the factorization method is applied carefully, the consistency conditions between phi1 and phi2 check out, and the explicit partner potentials and Lambert W solutions are not in the earlier literature. The Painlevé analysis is correct and motivates the multivalued nature of the solutions. The paper is clearly written and the authors credit refs. [42] and [47] properly.\n\nThe real problem: integrating Eq. (4) twice with zero constants and rescaling z = y / s-tilde^{1/4} gives\n\ny'' + (gamma v / Lambda) y' - a1 y - a2 y^2 - a3 y^3 - a4 y^4 - y^5 = 0,\n\nwith a2 = p / (2 Lambda s-tilde^{1/4}). Eq. (5) instead has + a2 y^2. For the plotted parameters (p = q = 300 > 0), that sign difference flips the shape of the nonlinearity, so the solutions shown are not solutions of the stated reduction of the HJ model. The stress-test note claims the y^5 term comes out positive and that s-tilde appears in the linear term; neither is correct. The actual flaw is the y^2 sign.\n\nThere is also the 'without loss of generality' setting of integration constants to zero in Section 2.1. The second constant C2 shifts the equilibrium levels; setting it to zero restricts to kinks connecting zeros of the polynomial, which is not the full physical family. This is a real loss of generality, though a milder issue than the sign error.\n\nNet: the mathematical construction is sound on its own terms, but the bridge to the nerve model is broken as written. The fix is a sign correction (or a redefinition of a2) plus a discussion of C2. If fixed, this is a publishable contribution to nonlinear biophysics and SUSY-inspired factorization. As is, I would ask for major revision, not desk rejection.\n\nRecommendation: send it to peer review. The algebra is checkable, the error is fixable, and the explicit partner solutions are worth preserving. No experimental claims are made, so the bar is internal consistency—which the paper currently misses in one place.\n\nBest,\n[Your name]","headline":"The SUSY pairing machinery is real, but the paper's central reduction to Eq. (5) has a sign error that severs the link to the nerve model for the plotted parameters.","tokens_in":17533,"tokens_out":17935,"would_cite":false,"duration_ms":126884,"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":"A nonlinear nerve-membrane model yields exact Lambert W-kink solitons with supersymmetric partner pulses.","keywords":["Lambert W kink solitons","Heimburg-Jackson model","nerve impulse","supersymmetric pairing","factorization method","Liénard equation","Boussinesq equation","membrane phase transition"],"falsifier":"Set the dimensionless parameters to one of the paper's example sets, e.g. $p=q=300$, $s=27$, $r=-141$, $k=2$, $\\delta=v=1$, keep a small nonzero constant $C$ in the second integration that led to Eq. (5), and integrate the travelling-wave ODE numerically; if the resulting front deviates from the Lambert W-kink profile (27) in a way that does not vanish as $C\\to 0$, the 'without loss of generality' step fails and the exact family is incomplete.","tokens_in":16474,"feed_emoji":"🧠","tokens_out":8877,"duration_ms":70864,"temperature":0.7,"pith_summary":"The paper claims that an improved Heimburg-Jackson model of nerve impulses, extended with third- and fourth-order nonlinear terms, possesses exact travelling-wave solutions whose profile is governed by the Lambert W function, and that these come in supersymmetric partner pairs. The argument runs through a travelling-wave reduction to a Liénard-type equation, a Painlevé test showing the equation is not integrable in the usual sense, and a factorization of the Liénard operator into two first-order factors. Reversing the factors while keeping the damping coefficient fixed produces a partner equation and a partner soliton with the same damping but different amplitude, width, and profile. The paper interprets the partner solutions as mechanically modulated membrane states, with possible relevance to abnormal neural activity and traumatic brain injury.","feed_headline":"Nerve impulses pair up as supersymmetric Lambert W kinks","feed_subtitle":"Exact solutions connect electromechanical membrane waves to supersymmetric quantum mechanics.","key_machinery":"The load-bearing object is the factorization of the Liénard-type equation (5) into a product of two first-order differential operators, $(d/d\\xi-\\phi_2(y))(d/d\\xi-\\phi_1(y))y=0$, with consistency conditions $\\phi_1\\phi_2=f(y)/y$ and $\\phi_1+\\phi_2+y\\,d\\phi_1/dy=-\\tilde{\\gamma}$. Imposing the first-order compatibility condition $(d/d\\xi-\\phi_1)y=0$ turns the second-order equation into a separable first-order one; integrating $dy/[(y-\\alpha)^2 y]=d\\xi/\\sqrt{3}$ produces the Lambert W function, whose multivalued branch structure matches the movable branch point found in the Painlevé test. Reversing $\\phi_1$ and $\\phi_2$ while holding $\\tilde{\\gamma}$ fixed generates the partner Liénard equation and a new potential $V_2(y)$. The cubic equation (19) for $\\alpha$ and its discriminant (22) control how many admissible factorization roots exist; the region $D<0$, with three real roots, is what allows heteroclinic kink connections.","core_discovery":"On the paper's own terms, the central discovery is that the generalized Heimburg-Jackson density-wave equation (1), after the rescaling and travelling-wave reduction leading to Eq. (5), admits an exact Lambert W-kink solution, Eq. (27) for one factorization order and Eq. (32c) or Eq. (41) for the other, and that reversing the order of the two first-order factors yields a supersymmetric partner equation whose exact solutions, Eqs. (35c) and (41), are partner solitons. The two partner solitons share the same damping coefficient $\\tilde{\\gamma}$, which the model associates with the axoplasmic fluid, and the same front speed, differing in amplitude, width, and pulse shape. The associated effective sextic potentials $V_1(y)$ and $V_2(y)$ are explicit, and the paper shows they do not satisfy shape invariance in the strict supersymmetric quantum mechanics sense. This establishes a link between nonlinear electromechanical wave propagation in nerve membranes and supersymmetric quantum mechanics.","pith_inferences":["An extension not pursued in the paper: if the dropped integration constants are restored, the exact Lambert W form may survive only as a restricted subclass, and the supersymmetric partner construction would need to be re-derived for the full family; checking this is the most direct stress test of the claim.","The $D<0$ existence condition ties the kinks to multistability of the effective potential, so the paper implicitly predicts that varying the elastic coefficients $p,q,r,s$ can switch membrane pulse propagation between soliton-supporting and monostable regimes; this is testable in principle by tuning lipid composition or compression.","Because the factorization machinery is generic for Liénard-type equations, the same supersymmetric pairing could be applied to other Boussinesq-type biophysical models, such as DNA soliton equations or mechanical metamaterial domain walls, to produce partner waveforms without any membrane-specific assumption.","The paper's proposal that partner pulses model mechanically induced perturbations could be made quantitative by comparing fitted amplitude and width changes from the partner formulas against experimentally recorded deformed action potentials under compression, but the paper does not provide such data."],"forward_implications":["Strongly nonlinear extensions of the Heimburg-Jackson model can support coherent localized density waves even though the reduced equation fails the Painlevé test.","Every Lambert W-kink soliton constructed here has a partner soliton with the same damping coefficient and front speed but modified amplitude and width, so supersymmetry offers a systematic way to generate new pulse morphologies from a known one.","The two partner potentials $V_1$ and $V_2$ are explicit and not shape invariant, meaning the pairing does not reduce to the usual shape-invariant supersymmetric quantum mechanics families.","Reversing the factorization order reverses the phase-transition front, so the same gel-liquid transition can propagate in either direction with different pulse geometry."],"supporting_citations":[{"why":"Foundation: Heimburg-Jackson model of the action potential as a propagating density pulse, which Eq. (1) extends.","marker":"[15]"},{"why":"Foundation: thermodynamic rationale for why the nerve impulse is a reversible compression wave; supplies the sound-speed analogy.","marker":"[16]"},{"why":"Source of the extended density-wave equation (1) with third- and fourth-order nonlinearities and double dispersion.","marker":"[48]"},{"why":"Earlier derivation of Lambert W-kink solitons for higher-order lipid-membrane nonlinearities; the current parameter ranges and kink form build on it.","marker":"[42]"},{"why":"Introduces supersymmetric pairing of kinks through reversed factorization order; the paper extends this to Lambert W kinks.","marker":"[47]"},{"why":"Provides the factorization method for nonlinear ODEs that reduces the Liénard equation to first-order compatibility conditions.","marker":"[45]"},{"why":"Supplies the general factorization conditions for nonlinear second-order differential equations, used to match coefficients in Eqs. (15a)-(15b).","marker":"[46]"},{"why":"Gives the phi-six kink with a long-range tail; its asymptotics match the power-law and exponential tails of the Lambert W kink.","marker":"[44]"},{"why":"Sets out the Painlevé test steps used to establish non-integrability of the reduced equation.","marker":"[51]"},{"why":"Supports the claim that the $D<0$ multistability region is needed for heteroclinic kink connections.","marker":"[60]"}],"fun_headline_variants":["Nerve impulses reveal supersymmetric partner kinks","Lambert W-kinks pair up as supersymmetric solitons","Exact kink solutions link nerve waves to quantum symmetry","Supersymmetric pairing found in nerve impulse kinks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two integration constants dropped in the reduction from Eq. (4) to Eq. (5) can be set to zero without losing any travelling-wave solutions.","fun_headline_variants_meta":{"raw":{"variants":["Nerve impulses reveal supersymmetric partner kinks","Lambert W-kinks pair up as supersymmetric solitons","Exact kink solutions link nerve waves to quantum symmetry","Supersymmetric pairing found in nerve impulse kinks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1284,"prompt_tokens":910,"completion_tokens":374,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":308}},"tokens_in":526,"tokens_out":374,"duration_ms":3645,"temperature":1.0,"reasoning_tokens":308,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:35:23.762115+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Set the dimensionless parameters to one of the paper's example sets, e.g. $p=q=300$, $s=27$, $r=-141$, $k=2$, $\\delta=v=1$, keep a small nonzero constant $C$ in the second integration that led to Eq. (5), and integrate the travelling-wave ODE numerically; if the resulting front deviates from the Lambert W-kink profile (27) in a way that does not vanish as $C\\to 0$, the 'without loss of generality' step fails and the exact family is incomplete.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundation: Heimburg-Jackson model of the action potential as a propagating density pulse, which Eq. (1) extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Foundation: thermodynamic rationale for why the nerve impulse is a reversible compression wave; supplies the sound-speed analogy."},{"cited_title":"Mendoza-Mill´ an, J.L","cited_arxiv_id":null,"evidence_quote":"Earlier derivation of Lambert W-kink solitons for higher-order lipid-membrane nonlinearities; the current parameter ranges and kink form build on it."},{"cited_title":"Cornejo-P´ erez and H","cited_arxiv_id":null,"evidence_quote":"Provides the factorization method for nonlinear ODEs that reduces the Liénard equation to first-order compatibility conditions."},{"cited_title":"Aϕ 6 soliton with a long-range tail","cited_arxiv_id":null,"evidence_quote":"Gives the phi-six kink with a long-range tail; its asymptotics match the power-law and exponential tails of the Lambert W kink."},{"cited_title":"Kudryashov","cited_arxiv_id":null,"evidence_quote":"Sets out the Painlevé test steps used to establish non-integrability of the reduced equation."},{"cited_title":"Solutions of the generalized Heim- burg–Jackson model for membrane pulses","cited_arxiv_id":null,"evidence_quote":"Supports the claim that the $D<0$ multistability region is needed for heteroclinic kink connections."}],"review_version":2}