Electrochemical response of biological membranes to localized currents and external electric fields

Hyeongjoo Row; Joshua B. Fernandes; Karthik Shekhar; Kranthi K. Mandadapu

REVIEW 3 major objections 4 minor 23 references

Electrochemical response of biological membranes to localized currents and external electric fields

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

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

Pith's one-line read This paper shows that under physiological conditions all free charge in an electrode-bounded membrane-electrolyte system is confined to nanometer-scale interfacial layers, so the transmembrane potential follows closed-form expressions in di

desk verdict Genuinely useful multiscale membrane-electrolyte reduction with honest caveats about cable theory; the missing derivation and numeric details make it a conditional accept rather than a clear pass. read the letter →

arxiv 2508.14001 v1 pith:R2IR6H2T submitted 2025-08-19 cond-mat.soft physics.bio-phq-bio.SC

Joshua B. Fernandes , Hyeongjoo Row , Kranthi K. Mandadapu , Karthik Shekhar This is my paper · ORCID

classification cond-mat.softphysics.bio-phq-bio.SC

keywords Poisson-Nernst-Planckmatchedasymptoticexpansionsdiffusechargelayerstransmembranepotentialelectrodeconfinementcabletheoryelectrophysiologylipidmembranes

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 analyzes a flat lipid membrane sitting between two parallel electrodes, with an ion-channel-like current crossing the membrane while the electrodes impose a step voltage. Using matched asymptotic expansions of the Poisson–Nernst–Planck equations, it argues that under physiological conditions the nanometer-scale diffuse charge layers at the membrane and electrode surfaces reach a quasi-steady state quickly, while the bulk electrolyte stays electroneutral. As a result, all free charge is confined to the interfacial screening layers, the bulk potential satisfies Laplace's equation with time-dependent boundary conditions, and the response separates into distinct space-time regimes. In each regime the paper derives closed-form expressions for the transmembrane potential, and at long times shows the system maps onto an equivalent circuit of the cable-theory type; the formulas are checked against nonlinear numerical simulations. These results matter because electrode-bounded measurements are common in electrophysiology, and this provides a predictive picture of how an externally applied field and a local ionic current jointly set membrane voltage.

What carries the argument

The central machinery is a matched asymptotic expansion of the Poisson–Nernst–Planck equations (the standard continuum description of ion transport coupled to electrostatics) in the limit where the Debye screening length is far smaller than the electrode gap. The expansion splits the domain into quasi-steady diffuse charge layers at the membrane and electrode interfaces plus an electroneutral bulk governed by Laplace's equation, with time-dependent flux boundary conditions coupling the layers to the bulk. This asymptotic splitting is what generates the regime partition, the closed-form transmembrane potentials, and the long-time equivalent-circuit reduction.

What would settle it

Solve the full nonlinear Poisson–Nernst–Planck equations for the same geometry but with an electrode gap comparable to the Debye length (or with a concentrated electrolyte where the Debye relaxation time is not much shorter than the bulk charging time). If the solution does not show rapidly equilibrated diffuse layers and an electroneutral bulk, the paper's regime boundaries and closed-form transmembrane potentials will not match, falsifying the central claim for regimes outside its stated physiological range.

Watch

Extended reading notes

Core claim

In a membrane-electrolyte cell closed by parallel electrodes, the full Poisson–Nernst–Planck system reduces via matched asymptotic expansions under physiological scale separation. The diffuse layers equilibrate quasi-statically while the bulk stays electroneutral, so all free charge sits in nanometer-scale interfacial layers; the bulk potential obeys Laplace's equation, with time-dependent boundary conditions coupling the layers. This partitions the response into capacitive, transitional, and diffusive regimes for blocking electrodes, and a screened steady state for Faradaic electrodes. Closed-form transmembrane potentials follow in each regime and match nonlinear simulations. At long times

Load-bearing premise

The load-bearing premise is that the Debye screening length and its relaxation time are far smaller than the electrode gap and the bulk charging time; if that scale separation fails, the quasi-steady diffuse layers and the electroneutral-bulk partition—and with them the closed-form regimes—no longer hold.

Editorial extensions

If this is right

With blocking electrodes, the long-time response reduces to a radial diffusion equation and an equivalent circuit; with Faradaic electrodes, it reaches a screened steady state, so the electrode kinetics qualitatively select the regime structure.
The equivalent circuit applies only after the bulk charging time and only globally over the electrode area; it does not describe local membrane dynamics, and the externally measured current deviates slightly from the transmembrane current while charge accumulates in the membrane-side boundary layer.
Electrode confinement radially screens the electric field that would be long-ranged in an unbounded geometry, changing where injected charge is deposited and how the transmembrane potential builds up.
Within each regime, closed-form expressions for the transmembrane potential can replace full nonlinear simulation, and the paper verifies them against such simulations.

Reading between the lines

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

If the scale-separation assumption carries over to curved geometries with local curvature radii large compared with the Debye length, the same reduction would give fast, analytically tractable predictions for field stimulation of cells and tissues in confined setups. (Editorial inference.)
The global-versus-local limitation of the equivalent circuit implies that standard cable-theory estimates of local membrane voltage from recorded external current carry a small systematic error during membrane-side charge accumulation; a combined electrode-current and optical voltage measurement could detect it. (Editorial inference.)
The regime boundaries suggest a testable scaling collapse: plotting transmembrane potential against time scaled by the bulk charging time should overlay curves from different electrode gaps while the asymptotic assumption holds, and deviations would mark where higher-order corrections or full PNP are needed. (Editorial inference.)

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.

Desk Editor's Note

read the letter

Short version: this is a genuinely useful paper. It takes the authors' previous unbounded-domain membrane work and adds a realistic electrode-bounded geometry, with a clean blocking/Faradaic distinction and an honest statement about when cable theory works. The main reason for a conditional verdict is not the physics but the missing derivation and numerical validation details.

The new content is real. The regime partition - capacitive, transitional, diffusive for blocking electrodes, capacitive and screened steady state for Faradaic - is a concrete advance over Ref [3], and the demonstration that the equivalent circuit is only global and only valid for t > τ_B is exactly the kind of caveat that makes a theory paper useful. The discussion is candid: the circuit model does not describe local membrane dynamics. Credit where due.

The matched-asymptotics reduction follows the standard diffuse-charge program. If the scale separation holds - Debye length much smaller than the electrode gap, and Debye relaxation much faster than bulk charging - then the bulk Laplace description and the quasi-steady boundary layers are the right leading-order picture. The abstract says these are physiological conditions and points to Table I, but the text we have doesn't show the actual parameter values, the derivation, or a quantitative validity bound. That is a real gap. The stress-test note is fair: without explicit criteria, a user doesn't know when the closed forms start to fail. But this is a missing-support issue, not a sign that the argument is wrong.

Likewise, the verification is only described, not shown. No solver details, no convergence study. For a paper whose whole claim rests on an asymptotic limit, the numerical section needs to demonstrate it actually probes that limit. I'd want that in the supplement.

The heavy self-citation to Refs [3] and [9] is not a problem by itself; the prior results appear to be the foundation. But it means errors there would propagate, so the referee should check the chain.

Bottom line: the physics is plausible, the caveats are honest, and the toolkit could be useful to anyone modeling electrostimulation or transporters under electrodes. It deserves peer review, not a desk reject. I'd send it back for a revision that shows the asymptotic derivation and the numerical setup. If those check out, this is a solid reference.

Referee Report

3 major / 4 minor

Summary. The paper studies a flat lipid membrane flanked by electrolyte-filled gaps and bounded by parallel electrodes, subject to a localized channel-like current and a step voltage applied across the electrodes. The authors claim that matched asymptotic expansions of the Poisson-Nernst-Planck (PNP) equations show that, under physiological conditions, diffuse charge layers rapidly quasi-steady while the bulk remains electroneutral, so that all free charge is confined to nanoscale interfacial layers. In this limit the bulk potential satisfies Laplace's equation with time-dependent boundary conditions, and the response is partitioned into capacitive, transitional, and diffusive regimes. Closed-form expressions for the transmembrane potential are derived for each regime and are stated to be verified against nonlinear numerical simulations. At long times an equivalent-circuit description is recovered, with the explicit caveat that this circuit model applies globally and only for t > tau_B, and that for blocking electrodes the external current deviates slightly from the transmembrane current.

Significance. If the central claim holds, the paper provides a useful multiscale reduction of a complex electrochemical problem, showing how electrode-induced screening and confinement alter the transmembrane-potential response. The approach is commendable for being parameter-free in the sense that the asymptotic derivation starts from the PNP equations with physiological inputs rather than fitting circuit elements, and for testing the closed forms against nonlinear numerical simulations. The discussion is also candid about the limits of the equivalent-circuit representation, acknowledging that it holds only in a global sense and only for times longer than the electrolyte charging time. These strengths make the work potentially valuable for membrane biophysics and for connecting continuum electrokinetics to cable-theory models. However, the quantitative validity range of the central asymptotic reduction is not demonstrated, and the numerical verification is described only qualitatively. The paper builds heavily on the authors' prior work (Ref. [3]), so the novel contribution must be clearly delineated.

major comments (3)

[Abstract / Sec. II C (Table I) / Discussion (Fig. 10)] The central reduction to quasi-steady diffuse layers and an electroneutral bulk presumes a separation of scales: the Debye length and its relaxation time must be much smaller than the electrode gap and the bulk charging time. The manuscript states this under 'physiological conditions' but does not derive a quantitative criterion or an error estimate. The discussion explicitly provides a time bound for the circuit model (t > tau_B), but no analogous bound is given for the quasi-steady diffuse-layer assumption. If the separation fails (small electrode gaps, concentrated electrolytes, or fast voltage steps), the bulk is not electroneutral to leading order and the closed-form regime expressions break down. Please state the asymptotic ordering, the parameter ranges in Table I, and the expected accuracy of the closed forms across that range.
[Abstract / 'verified against nonlinear numerical simulations'] The paper reports numerical verification but provides no solver details, mesh resolution, time-step control, convergence studies, or quantitative error metrics. Without these details the reader cannot judge whether the simulations actually probe the asymptotic limit or whether the comparison is meaningful. Please report the numerical method, the parameter sweep, convergence checks, and the error between the asymptotic and simulated transmembrane potentials for each regime.
[Sec. I and Discussion (Eqs. (26)-(27), (36))] The matched asymptotic expansion itself is not presented in the supplied text. The regime partition, the reduction to Laplace's equation, and the equivalent-circuit equations are described verbally and via schematic diagrams. Since these are the main theoretical results, the manuscript should include the expansions, the boundary-layer scalings, the matching conditions, and the derivation leading to Eqs. (26)-(27) and (36). The relationship to Ref. [3] should also be explicit so the reader can identify which steps are new and which are inherited.

minor comments (4)

[Fig. 10] The axes are described as schematic. Please label the regime boundaries with the relevant dimensionless parameters (e.g., t/tau_B, r/lambda_D) or provide representative scales, so the schematic can be connected to the asymptotic results.
[Introduction] The phrase 'an ion channel-like current flows across the membrane' is vague. Please specify the current-injection condition used in the model (e.g., fixed current density, current source at a point, or a prescribed flux boundary condition).
[Throughout] The term 'physiological conditions' is used without definition. Please provide the parameter ranges intended (ion concentration, electrode gap, membrane thickness, voltage step, current magnitude) and cite typical values, since these determine whether the asymptotic scale separation holds.
[Discussion, blocking electrodes] The statement that the external current 'deviates slightly (but consistently)' should be quantified, ideally with an expression or a bound showing the order of the deviation in the asymptotic parameter.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the PNP asymptotics, regime analysis, and equivalent-circuit result are derived from the governing equations and checked against nonlinear simulations; self-citations are background, not load-bearing.

full rationale

The paper's central derivation is a matched asymptotic expansion of the Poisson-Nernst-Planck equations in the limit of thin Debye layers relative to the electrode gap, with the diffuse-layer relaxation assumed fast compared with bulk charging. This is a stated scale-separation hypothesis, not an output secretly defined as an input; the paper itself acknowledges the limits of the reduction, noting in the discussion that the equivalent-circuit model 'only emerges on timescales longer than the electrolyte charging time τB' and 'is not universally valid.' The bulk-Laplace/quasi-steady diffuse-layer structure and the circuit representation are obtained from the governing electrodiffusion equations rather than from fitting parameters or from importing an external uniqueness theorem. No free parameters are fitted to quantities later called predictions; the closed-form expressions are verified against nonlinear numerical simulations, an external benchmark. The authors' prior work (Refs. [3] and [9]) is used for background and for comparison—e.g., Fig. 10 contrasts the unbounded case of Ref. [3] with the present electrode-bounded case—but the core asymptotic reduction and circuit derivation are not justified solely by those citations. The supplied excerpt omits the detailed asymptotic calculation and the numerical-methods section, so the derivation cannot be fully audited from the text provided; that is an evidentiary limitation, not circularity. No step in the visible text reduces an output to an input by construction, and no self-citation chain forces the central claims.

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

No parameters are fitted to data. All inputs (membrane thickness, gap size, permittivities, salt concentration) are stated physiological values from Table I, treated as inputs from prior literature; the closed-form expressions are verified against independent nonlinear simulations, so the central claim contains no tuned constants. No new particles, forces, dimensions, or conserved quantities are introduced; the equivalent circuit and the screened steady state are coarse-grained representations of existing electrostatics, not new entities.

assumptions (3)

domain assumption Poisson-Nernst-Planck continuum equations with constant electrolyte and membrane permittivities describe the ionic dynamics.
Invoked in the abstract as the starting model. Standard for dilute electrolytes but neglects finite ion size, steric effects, and permittivity variation near interfaces.
domain assumption Strong scale separation: Debye screening length and boundary-layer relaxation time are much smaller than the electrode gap L and the bulk charging time tau_B (the "physiological conditions").
This separation justifies the matched asymptotic expansion and the quasi-steady diffuse-charge layers. It is stated, not derived, and the discussion's circuit limit explicitly requires t > tau_B.
domain assumption The membrane is a thin dielectric slab with a prescribed channel-like current source, and electrode kinetics are idealized as purely blocking or purely Faradaic.
Defines the setup in Fig. 1 and the abstract. The discussion shows the blocking/Faradaic choice changes the steady state, so the conclusions depend on these idealizations.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Electrochemical response of biological membranes to localized currents and external electric fields." pith.science (2026). https://pith.science/paper/R2IR6H2T

@misc{pith2026250814001,
  author       = {Pith},
  title        = {Pith review of: Electrochemical response of biological membranes to localized currents and external electric fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R2IR6H2T}},
  note         = {Machine review of arXiv:2508.14001}
}

read the original abstract

Electrochemical phenomena in biology often unfold in confined geometries where micrometer- to millimeter-scale domains coexist with nanometer-scale interfacial diffuse charge layers. We analyze a model lipid membrane-electrolyte system where an ion channel-like current flows across the membrane while parallel electrodes simultaneously apply a step voltage, emulating an extrinsic electric field. Matched asymptotic expansions of the Poisson-Nernst-Planck equations show that, under physiological conditions, the diffuse charge layers rapidly reach a quasi-steady state, and the bulk electrolyte remains electroneutral. As a result, all free charge is confined to the nanometer-scale screening layers at the membrane and electrode interfaces. The bulk electric potential satisfies Laplace's equation, and is dynamically coupled to the interfacial layers through time-dependent boundary conditions. This multiscale coupling partitions the space-time response into distinct regimes. At sufficiently long times, we show that the system can be represented by an equivalent circuit analogous to those used in classical cable theory. We derive closed-form expressions of the transmembrane potential within each regime, and verify them against nonlinear numerical simulations. Our results show how electrode-induced screening and confinement effects influence the electrochemical response over multiple length and time scales in biological systems.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

[3]

H. Row, J. B. Fernandes, K. K. Mandadapu, and K. Shekhar, Spatiotemporal dynamics of ionic reorga- nization near biological membrane interfaces, Physical Review Research 7, 013185 (2025)

work page 2025
[9]

Farhadi, J

J. Farhadi, J. B. Fernandes, K. Shekhar, and K. K. Mandadapu, Capacitive response of biological membranes, Physical Review E 111, 064412 (2025)

work page 2025
[1]

Lodish, Molecular Cell Biology (W.H

H. Lodish, Molecular Cell Biology (W.H. Freeman, New York, NY, 2000)

work page 2000
[2]

Hille, Ionic Channels of Excitable Membranes (Oxford University Press, Oxford, 1992)

B. Hille, Ionic Channels of Excitable Membranes (Oxford University Press, Oxford, 1992)

work page 1992
[4]

Bikson, M

M. Bikson, M. Inoue, H. Akiyama, J. K. Deans, J. E. Fox, H. Miyakawa, and J. G. Jefferys, Effects of uniform extracellular DC electric fields on excitability in rat hippocampal slices in vitro, The Journal of Physiology 557, 175 (2004)

work page 2004
[5]

Radman, R

T. Radman, R. L. Ramos, J. C. Brumberg, and M. Bikson, Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation in vitro, Brain Stimulation 2, 215 (2009)

work page 2009
[6]

J. T. Francis, B. J. Gluckman, and S. J. Schiff, Sensitivity of neurons to weak electric fields, Journal of Neuroscience 23, 7255 (2003)

work page 2003
[7]

Neher and B

E. Neher and B. Sakmann, The patch clamp technique, Scientific American 266, 44 (1992)

work page 1992

Show all 23 references

[8]

J. W. Moore and K. S. Cole, Voltage clamp techniques, in Physical Techniques in Biological Research , edited by W. L. Nastuk (Academic Press, Cambridge, MA, 1963) Chap. 6, pp. 263–321

1963
[10]

J. R. Macdonald, Binary electrolyte small-signal frequency response, Journal of Electroanalytical Chem- istry and Interfacial Electrochemistry 53, 1 (1974)

1974
[11]

A. A. Kornyshev and M. A. Vorotyntsev, Conductivity and space charge phenomena in solid electrolytes with one mobile charge carrier species, a review with original material, Electrochimica Acta 26, 303 (1981)

1981
[12]

M. Z. Bazant, K. Thornton, and A. Ajdari, Diffuse-charge dynamics in electrochemical systems, Physical Review E 70, 021506 (2004)

2004
[13]

Janssen and M

M. Janssen and M. Bier, Transient dynamics of elec- tric double-layer capacitors: Exact expressions within the debye-falkenhagen approximation, Physical Review E97, 052616 (2018)

2018
[14]

Dayan and L

P. Dayan and L. F. Abbott, Theoretical Neuroscience, Computational Neuroscience (MIT Press, Cambridge, MA, 2001)

2001
[15]

K. D. Fong, H. K. Bergstrom, B. D. McCloskey, and K. K. Mandadapu, Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statis- tical mechanics, AIChE Journal 66, e17091 (2020)

2020
[16]

Nernst, Zur kinetik der in l¨ osung befindlichen k¨ orper, Zeitschrift f¨ ur Physikalische Chemie2U, 613 (1888)

W. Nernst, Zur kinetik der in l¨ osung befindlichen k¨ orper, Zeitschrift f¨ ur Physikalische Chemie2U, 613 (1888)
[17]

Nernst, Die elektromotorische wirksamkeit der jonen, Zeitschrift f¨ ur Physikalische Chemie4U, 129 (1889)

W. Nernst, Die elektromotorische wirksamkeit der jonen, Zeitschrift f¨ ur Physikalische Chemie4U, 129 (1889)
[18]

Planck, Ueber die potentialdifferenz zwischen zwei verd¨ unnten l¨ osungen bin¨ arer electrolyte, Annalen der Physik 276, 561 (1890)

M. Planck, Ueber die potentialdifferenz zwischen zwei verd¨ unnten l¨ osungen bin¨ arer electrolyte, Annalen der Physik 276, 561 (1890)
[19]

A. G. Volkov and T. Hampton, Energetics of membrane permeability, Advances in Planar Lipid Bilayers and Liposomes 8, 155 (2008)

2008
[20]

Nymeyer and H.-X

H. Nymeyer and H.-X. Zhou, A method to determine dielectric constants in nonhomogeneous systems: appli- cation to biological membranes, Biophysical Journal 94, 1185 (2008)

2008
[21]

Gouy, Sur la constitution de la charge ´ electrique ` a la surface d’un ´ electrolyte, Journal de Physique Th´ eorique et Appliqu´ ee9, 457–468 (1910)

M. Gouy, Sur la constitution de la charge ´ electrique ` a la surface d’un ´ electrolyte, Journal de Physique Th´ eorique et Appliqu´ ee9, 457–468 (1910)

1910
[22]

D. L. Chapman, LI. A contribution to the theory of electrocapillarity, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 25, 475 (1913)

1913
[23]

Moldenhauer, I

H. Moldenhauer, I. D ´ ıaz-Franulic, F. Gonz´ alez-Nilo, and

Pith tools

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

[1] [3]

H. Row, J. B. Fernandes, K. K. Mandadapu, and K. Shekhar, Spatiotemporal dynamics of ionic reorga- nization near biological membrane interfaces, Physical Review Research 7, 013185 (2025)

work page 2025

[2] [9]

Farhadi, J

J. Farhadi, J. B. Fernandes, K. Shekhar, and K. K. Mandadapu, Capacitive response of biological membranes, Physical Review E 111, 064412 (2025)

work page 2025

[3] [1]

Lodish, Molecular Cell Biology (W.H

H. Lodish, Molecular Cell Biology (W.H. Freeman, New York, NY, 2000)

work page 2000

[4] [2]

Hille, Ionic Channels of Excitable Membranes (Oxford University Press, Oxford, 1992)

B. Hille, Ionic Channels of Excitable Membranes (Oxford University Press, Oxford, 1992)

work page 1992

[5] [4]

Bikson, M

M. Bikson, M. Inoue, H. Akiyama, J. K. Deans, J. E. Fox, H. Miyakawa, and J. G. Jefferys, Effects of uniform extracellular DC electric fields on excitability in rat hippocampal slices in vitro, The Journal of Physiology 557, 175 (2004)

work page 2004

[6] [5]

Radman, R

T. Radman, R. L. Ramos, J. C. Brumberg, and M. Bikson, Role of cortical cell type and morphology in subthreshold and suprathreshold uniform electric field stimulation in vitro, Brain Stimulation 2, 215 (2009)

work page 2009

[7] [6]

J. T. Francis, B. J. Gluckman, and S. J. Schiff, Sensitivity of neurons to weak electric fields, Journal of Neuroscience 23, 7255 (2003)

work page 2003

[8] [7]

Neher and B

E. Neher and B. Sakmann, The patch clamp technique, Scientific American 266, 44 (1992)

work page 1992