{"id":"a6273732-497c-41b2-a72a-c83fbb9898ce","arxiv_id":"1908.05086","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":15,"one_line_summary":"A four-variable electrodiffusion-based reformulation of the Hodgkin-Huxley squid axon model predicts action potential speed, rebound spiking, and low-calcium repetitive firing without HH's empirical power-law gating.","lead":"This preprint offers a simpler, electrodiffusion-based replacement for the classic Hodgkin-Huxley squid axon model, keeping four gating variables but using constant relaxation times. It reports action potential speeds matching experiment and rebound spiking at temperatures where the original Hodgkin-Huxley model fails.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central no-spike-train prediction for prolonged constant current is asserted without bifurcation analysis; a unique stable rest state does not rule out limit cycles, especially since HH's tonic firing arises from a Hopf bifurcation of a single fixed point.","rationale":"The reader's nominated weakest assumption is the Goldman constant-field approximation combined with a spatially constant PMF. That is a real idealization, and the authors acknowledge it explicitly ('We adopt the simplest possible model for the PMF'), though they do not quantify the resulting error. However, this approximation would mainly shift quantitative predictions, and the model's qualitative phenomena might survive a more structured PMF; it is therefore not the most load-bearing point. The more serious gap is internal to the paper's central claim of superiority over HH: the claim that prolonged constant-current stimulation produces no spike trains 'of any size' is a universal negative, yet the only dynamical evidence offered is a unique, linearly stable fixed point for zero current. That evidence cannot exclude limit cycles or Hopf bifurcations for nonzero injected current, and indeed the classic HH mechanism for tonic firing is a Hopf bifurcation of a unique fixed point. The paper reports no bifurcation diagram in the stimulus current, no Lyapunov or global-stability argument, and no long-time constant-current simulations. This is directly checkable, and the outcome could either strengthen or falsify a headline prediction. The correct verdict remains CONDITIONAL, since the concern is unresolved rather than demonstrably fatal, but the condition should explicitly include a constant-current bifurcation analysis.","tokens_in":22319,"tokens_out":5177,"duration_ms":54645,"concrete_test":"Fix all Table 1 parameters and add a constant term I_stim to Eq. (9b): C_m dV/dt = I_stim - (I_Na + I_K + I_Cl). For I_stim ranging over at least -200 to +200 μA/cm^2 in fine steps (e.g., 401 values), integrate the four ODEs (9b, 11a-c) from the resting initial condition and from several perturbed initial conditions for t ≥ 500 ms, and separately compute fixed points and Jacobian eigenvalues to locate Hopf bifurcations. If no oscillatory asymptotic states and no eigenvalue crossings are found for any I_stim, the no-spike-train claim is supported; if a limit cycle or Hopf bifurcation appears for any I_stim, the central claim fails as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.1 establishes only that the zero-stimulus steady state is the unique fixed point of the autonomous system and that it is linearly stable (Fig. 3). The concluding claim that 'unlike its HH ... counterpart, our new ... model ... does not predict spike trains ... during prolonged constant-current stimulation of any size' (Section 3) does not follow from this analysis. A system can have a unique, linearly stable fixed point and still possess a stable limit cycle for nonzero injected current; in the HH model itself, maintained current produces repetitive firing precisely through a Hopf bifurcation that leaves the number of fixed points unchanged. Equations (9b) and (11a-c) are nonlinear, so the existence of oscillations for I_stim ≠ 0 is a global question requiring either a bifurcation diagram in I_stim, a Lyapunov argument, or direct long-time integration. The paper reports none of these. Because the assertion is a universal negative over 'any size' of current and is one of the four headline experimental agreements claimed, the missing analysis is load-bearing rather than cosmetic.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a four-variable model for squid giant axon action potentials, replacing the HH phenomenological gating expressions with a dynamic GHK electrodiffusion formulation. The model assumes spatially constant potentials of mean force inside channels, first-order gating kinetics with constant time constants, and sigmoidal steady-state gating functions. It reports that the model reproduces the stationary action potential waveform, propagates a pulse at 22.3 m/s at 20C (vs HH's 21.2 m/s at 18.5C), predicts rebound spiking at 20C (where HH fails), predicts repetitive firing under low external calcium after brief depolarizing current, and predicts no spike trains under prolonged constant current in perfused axons. The paper also discusses temperature dependencies qualitatively. Parameter values are explicitly tabulated.","tokens_in":22719,"tokens_out":8684,"duration_ms":79062,"significance":"The model is a serious attempt to simplify the HH framework while grounding it in electrodiffusion. It avoids the arbitrary m^3h and n^4 powers, uses only four variables, and produces several falsifiable predictions, including rebound spiking at 20C and low-calcium repetitive firing. The manuscript is commendably transparent: it explicitly states that the seven PMF barriers, three relaxation times, steepness parameters, and area fractions were selected to reproduce the HH action potential shape and voltage-clamp transients, so those agreements are not independent validations. The genuinely predictive elements are the propagation speed, rebound spiking, and low-calcium firing. The paper also acknowledges discrepancies with experimental permeability ratios. However, the universal negative claim about constant-current spike trains is not justified by the local stability analysis, and the speed result lacks any sensitivity or error assessment. If the missing global bifurcation analysis and sensitivity studies are added, the model would be a valuable contribution.","major_comments":[{"comment":"The concluding claim that the model 'does not predict spike trains of repetitive action potentials during prolonged constant-current stimulation of any size' is not supported by the analysis. Section 2.1 establishes only that the zero-stimulus steady state is the unique fixed point and is linearly stable (Fig. 3). A dynamical system can have a unique, linearly stable fixed point and still possess stable limit cycles for nonzero parameter values; in the HH model, tonic firing arises through a Hopf bifurcation of a single fixed point. Equations (9b) and (11a-c) are nonlinear, so the existence of oscillations for I_stim ≠ 0 is a global question. A universal negative over all stimulus amplitudes requires either a bifurcation diagram in I_stim, a Lyapunov argument, or direct long-time simulations for a range of amplitudes. Because this claim is one of the four headline agreements, the missing analysis is load-bearing.","section":"Section 3; Section 2.1"},{"comment":"The similarity of the model's stationary action potential to the HH recordings is partly built in by construction. The authors state: 'Other parameters were selected to ensure that stationary membrane action potentials (see Fig. 4) have temporal profiles similar to those recorded by HH,' and that the seven PMF barriers, three gating relaxation times, and channel area fractions were chosen with this criterion and with the additional requirement that voltage-clamp transients 'at least qualitatively' resemble those of HH. Consequently, the agreement in Fig. 4 should not be presented as an independent success. The independent tests are the propagation speed (Section 2.4), rebound spiking at 20C (Figs. 1B and 7), and low-calcium repetitive firing (Section 2.3); these should be foregrounded.","section":"Section 2.2 (paragraph before Table 1)"},{"comment":"The propagation speed of 22.3 m/s at 20C is compared with the HH experimental estimate of 21.2 m/s at 18.5C, but the result is a single deterministic simulation with no error bars, no sensitivity analysis, and no discussion of how the speed depends on key parameters such as axon radius a, internal resistivity R, gating time constants, and PMF barriers. Given the 1.5C temperature difference and the acknowledged fitting of the parameters to stationary AP properties, the agreement could be fortuitous. A parameter scan or at least a statement of the expected sensitivity of v to the main parameters is needed to support the 'good agreement' claim.","section":"Section 2.4 (Fig. 6)"},{"comment":"The derivation of the current-voltage relation assumes the potential of mean force w(x,t) inside a channel is spatially constant, so the only x-dependence in the flux integral is the linear electrostatic term. The paper gives no evidence that this approximation is adequate during fast transients, when the gating variables are changing rapidly. If the true PMF is spatially structured, the predicted currents, resting potential, and action potential shape could shift materially. The authors should either justify this approximation quantitatively (e.g., by comparing with a spatially varying PMF or with known free-energy profiles from molecular dynamics) or state clearly that the model is a coarse-graining whose quantitative predictions are approximate.","section":"Section 1, Eq. (2)"}],"minor_comments":[{"comment":"The GHK equation is garbled in the typeset; the subscripts and superscripts for internal and external concentrations are misaligned, obscuring the numerator and denominator structure. Please re-set the equation carefully.","section":"Eq. (7)"},{"comment":"The steepness parameter s_h is listed as '11' without units; since h_ss depends on the dimensionless variable m, s_h is dimensionless, but this should be stated explicitly. The notation for τ_m, τ_h and s_m, s_h is also cramped and should be clarified.","section":"Table 1"},{"comment":"The stimulus current amplitude switches from -65 and -69 μA cm^-2 in the first part of the paragraph to -269 μA cm^-2 later; the relationship between these values and the bifurcation description should be clarified.","section":"Section 2.3"},{"comment":"The labels 'A action potential V(t) (mV)' and 'B action potential V(t) (mV)' are malformed; they should read 'A. Action potential...' and 'B. Action potential...'.","section":"Fig. 1 caption"},{"comment":"The reference to 'Stiles and Gray, 2019' is the manuscript under review itself (arXiv:1908.05086); citing the submitted manuscript as a reference is inappropriate and should be removed.","section":"References"},{"comment":"The full text contains only placeholder text for Fig. 6 and Fig. 7 with no graphical content; the final version must include these plots so the propagation and rebound waveforms can be inspected.","section":"Figures 6 and 7"},{"comment":"The reference list contains an editorial note appended to the Perram and Stiles (2010) entry ('This model incorrectly describes...'); such commentary belongs in the main text, not in the reference list.","section":"Reference list entry for Perram and Stiles (2010)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its fitting procedure, but the no-spike-train claim is the most serious gap; I would request a bifurcation analysis before publication. Also, the self-citation to the same preprint and the editorial note in the reference list should be removed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a real reformulation, and the physics story is attractive. But the paper's most advertised prediction—no spike trains under prolonged constant current—is not actually established by anything in the paper. That's a load-bearing gap.\n\nWhat's new: They take GHK electrodiffusion and let the permeabilities become time-dependent through gate-dependent PMF barriers. That removes HH's m^3h and n^4 factors and gives a four-variable resonator model with constant relaxation times. The derivation is explicit. They match the HH action potential shape, get a propagation speed of 22.3 m/s versus the experimental 21.2 at 18.5 C, and predict rebound spiking at temperatures where HH fails. They also show low-calcium repetitive firing by lowering the Na activation barrier. If the model works, it is simpler and more physical than HH for perfused squid axon.\n\nCredit where due: the paper is honest about its own history—it explicitly says their 2010 model was wrong about integrator dynamics and had three fixed points. The parameters are all listed in a table, and the fitting criteria are stated. That's more than many modeling papers do.\n\nThe soft spots, in order of size:\n\nFirst, the no-spike-trains claim. The stress-test note is exactly right. Section 2.1 shows a unique, linearly stable fixed point for the autonomous zero-stimulus system. That says nothing about whether a stable limit cycle exists for constant I_stim ≠ 0. HH itself has a Hopf bifurcation from a single fixed point. Equations (9b) and (11) are nonlinear; the claim \"does not predict spike trains ... during prolonged constant-current stimulation of any size\" is a universal negative over all current amplitudes. That requires either a bifurcation diagram in I_stim, a Lyapunov function, or direct long-time integration. None is provided. This is not a cosmetic omission; it is one of the four headline results.\n\nSecond, parameter circularity. Many parameters are chosen so that the space-clamped AP matches HH's recorded shapes. So the agreement in Fig. 4 is partially built in. The propagation speed is a single number with no sensitivity or error analysis. The low-calcium spiking is a genuine prediction, since those parameters weren't fit to it, but it needs a bifurcation analysis too.\n\nThird, no code or detailed numerical procedures. The paper says \"Mathematica permits ... routine\" solution, but without code or step sizes/error tolerances, reproduction is harder than it should be.\n\nThe PMF assumption (spatially constant w inside the channel) is clearly flagged as the simplest possible model. That's a limitation, not a flaw.\n\nBottom line: this paper deserves a serious referee. The missing bifurcation analysis should be requested in revision, not treated as a desk-reject issue. The model is plausible, the writing is clear, and the physiological claims are concrete and testable. I would send it out.","headline":"A genuine GHK-based reformulation of HH with a clear head, but the central no-spike-train claim lacks the required bifurcation analysis.","tokens_in":23173,"tokens_out":2290,"would_cite":true,"duration_ms":22721,"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":"Electrodiffusion model predicts squid rebound spikes at 20 °C","keywords":["action potential","squid giant axon","Goldman-Hodgkin-Katz equation","electrodiffusion","potential of mean force","gating kinetics","rebound spiking","Hodgkin-Huxley model"],"falsifier":"Measure the spatial profile of the potential of mean force for sodium and potassium ions inside squid voltage-gated channels under voltage-clamp conditions; if the profile is not approximately flat or its gating dependence cannot be captured by the three PMF barriers, Eq. (3) fails. A simpler check: voltage-clamp the fully open sodium conductance and test whether the current follows the GHK form across the full range of membrane potentials; a systematic deviation would break the predicted I-V curve and threshold.","tokens_in":2032,"feed_emoji":"⚡","tokens_out":2886,"duration_ms":72832,"temperature":0.7,"pith_summary":"The paper argues that the classic Hodgkin-Huxley (HH) model's empirical conductance rules can be replaced by a simpler four-variable model built from the Goldman-Hodgkin-Katz (GHK) electrodiffusion current equation. It claims this model reproduces the squid giant axon's main behaviors: action potentials, propagation speed around 22 m/s, rebound spiking at 18.5 to 20 °C where the HH model fails, and repetitive firing after external calcium is lowered. It also predicts that perfused axons do not fire spike trains under sustained constant-current stimulation, in line with experiment but contrary to HH. The result matters because it ties action-potential generation to electrodiffusion quantities, such as permeabilities, free-energy barriers, and gating relaxation times, rather than to fitted gating powers, and it removes several known HH discrepancies.","feed_headline":"Electrodiffusion model predicts squid rebound spikes at 20 °C","feed_subtitle":"A four-variable GHK-based model matches propagation speed and low-calcium firing where Hodgkin-Huxley falls short.","key_machinery":"The central object is the generalized GHK current equation, Eq. (3), in which the resting-state Goldman-Hodgkin-Katz formula is extended by replacing the constant membrane potential and permeabilities with time-dependent ones: $i = -e^2 P(t) V_m(t)\\,[c^{\\mathrm{int}}-c^{\\mathrm{ext}}\\exp(e\\beta V_m(t))]/[1-\\exp(e\\beta V_m(t))]$. The permeabilities are exponential functions of potentials of mean force, $P(t)=(fD/L)\\exp(-\\beta w(t))$, and the potentials of mean force for sodium and potassium channels depend linearly on the gating variables $m$, $h$, and $n$. The dynamical system closes with three first-order relaxation equations using constant time constants $\\tau_m$, $\\tau_h$, $\\tau_n$ and sigmoidal steady-state gating functions. This machinery carries the argument by turning gating from empirical conductance multipliers into modulated free-energy barriers, so that the same electrodiffusion formalism produces both the resting potential and the action potential.","core_discovery":"The central claim is that a four-variable dynamical system, consisting of membrane depolarization $V$ plus gating variables $m$, $h$, and $n$, with ionic currents from the GHK constant-field equation, reproduces squid axon firing without HH's $m^3 h$ and $n^4$ dependencies. Gating enters through time-dependent permeabilities $P(t)=(fD/L)\\exp(-\\beta w(t))$, where the potential of mean force $w(t)$ for a sodium or potassium channel is a linear combination of open- and closed-gate free-energy barriers weighted by the gating variables; chloride is ungated and static. The model is a resonator with a single stable resting fixed point, and it predicts stationary action potentials peaking near 120 mV, a propagating speed of 22.3 m/s at 20 °C versus the experimental 21.2 m/s at 18.5 °C, rebound spiking after prolonged hyperpolarization at both 6.3 °C and 20 °C, and repetitive firing when the sodium activation barrier is lowered from $\\beta w_\\mathrm{Na}^O=3.0$ to $1.48$ to mimic low external calcium. Unlike HH, it does not produce spike trains under prolonged constant current in perfused, pumps-off axons.","pith_inferences":["If the effective free-energy barrier interpretation holds, the same equations could connect extracellular calcium concentration quantitatively to sodium-channel permeability, yielding a testable prediction for how threshold current shifts with $[\\mathrm{Ca}^{2+}]_\\mathrm{out}$.","The model's success would suggest that atomic-level free-energy profiles from molecular dynamics could be inserted into Eq. (3) in place of the single constant PMF, producing a more mechanistic bridge from channel structure to whole-axon firing.","Because the chloride concentration choice is load-bearing for the resting potential, the model could be tested by switching internal chloride to the higher values reported by later experiments and checking whether the predicted action-potential shape survives.","A related testable extension would be to apply the same GHK-based gating to other axon types, such as myelinated fibers, replacing the empirical gating powers still used in Frankenhaeuser-Huxley-type models."],"forward_implications":["HH's $m^3 h$ and $n^4$ gating terms are not necessary to reproduce squid action potentials; gating can be represented as exponential changes in channel free-energy barriers.","Rebound spiking after prolonged hyperpolarization is available near physiological temperature, so inhibitory rebound excitation may be common in squid axons, not limited to the cold temperatures at which HH can produce it.","The low-calcium repetitive firing regime is tied to a single control parameter, the sodium activation barrier $\\beta w_\\mathrm{Na}^O$, lowered from 3.0 to 1.48, giving a concrete bifurcation that could be checked against calcium-binding measurements.","Perfused squid axons with pumps off should not sustain spike trains under constant current, matching Clay's experiments and distinguishing perfused from freshly dissected fibers.","The model is a resonator with one stable fixed point, so subthreshold oscillations and resonate-and-fire behavior arise naturally from gating delays rather than from any explicit inductive element."],"supporting_citations":[{"why":"Supplies the constant-field approximation that the model extends from resting-state electrodiffusion to time-dependent action potentials.","marker":"Goldman, 1943"},{"why":"The original quantitative model and experimental dataset whose gating terms, voltage-clamp curves, propagation speed, and rebound-spiking observations the paper compares against.","marker":"Hodgkin and Huxley, 1952"},{"why":"Gives the GHK resting-potential formula that the model generalizes to time-dependent permeabilities.","marker":"Hodgkin and Katz, 1949a"},{"why":"Earlier electrodiffusion dynamical treatment of squid action potentials whose integrator behavior and three fixed points the present resonator model corrects.","marker":"Perram and Stiles, 2010"},{"why":"Experimental evidence that lowering external calcium enhances excitability and drives repetitive firing in squid axons.","marker":"Frankenhaeuser and Hodgkin, 1957"},{"why":"Experimental claim that perfused squid axons do not fire repetitive action potentials under prolonged constant-current stimulation, which the new model is designed to match.","marker":"Clay, 1998"},{"why":"Defines the perfused, pumps-off squid axon preparation whose resting potential and ionic currents the paper models.","marker":"Baker et al., 1961, 1962"},{"why":"Voltage-clamp current-voltage data used to validate that the model's steady-state I-V curve has a single fixed point.","marker":"Hodgkin et al., 1952"}],"fun_headline_variants":["Squid axon model reboots Hodgkin-Huxley with electrodiffusion","GHK-based model sparks rebound spikes in squid simulation","Four-variable electrodiffusion model beats Hodgkin-Huxley on squid","Simpler squid model: electrodiffusion predicts rebound firing","Squid action potentials: electrodiffusion model outperforms HH"],"cache_read_input_tokens":25216,"weakest_assumption_plain":"The whole current-voltage relation rests on the Goldman constant-field approximation together with the assumption that each channel's potential of mean force is spatially constant inside the pore; if the real free-energy profile is strongly structured or changes during fast transients, the predicted currents, resting potential, and action potentials would shift materially.","fun_headline_variants_meta":{"raw":{"variants":["Squid axon model reboots Hodgkin-Huxley with electrodiffusion","GHK-based model sparks rebound spikes in squid simulation","Four-variable electrodiffusion model beats Hodgkin-Huxley on squid","Simpler squid model: electrodiffusion predicts rebound firing","Squid action potentials: electrodiffusion model outperforms HH"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000783,"raw_usage":{"total_tokens":3524,"prompt_tokens":1076,"completion_tokens":2448,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":2357}},"tokens_in":692,"tokens_out":2448,"duration_ms":13806,"temperature":1.0,"reasoning_tokens":2357,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:24:07.991163+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spatial profile of the potential of mean force for sodium and potassium ions inside squid voltage-gated channels under voltage-clamp conditions; if the profile is not approximately flat or its gating dependence cannot be captured by the three PMF barriers, Eq. (3) fails. A simpler check: voltage-clamp the fully open sodium conductance and test whether the current follows the GHK form across the full range of membrane potentials; a systematic deviation would break the predicted I-V curve and threshold.","supporting_citations":[],"review_version":1}