The transmembrane potential across a charged nanochannel subjected to asymmetric electrolytes
Pith reviewed 2026-05-16 17:54 UTC · model grok-4.3
The pith
Two analytical expressions for the zero-current transmembrane voltage in charged nanochannels with asymmetric electrolytes are derived from the Poisson-Nernst-Planck equations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We derive two expressions for V_{i=0}. First, we consider the much simpler scenario of two species. Then, we can consider an electrolyte comprised of an arbitrary number of species. The difference in the models is that the latter solution utilizes an ad-hoc assumption of a linear concentration profile, while the former solution does not require such an ad-hoc assumption. We analyze both models and show how to reduce them to several known models. We also verify both models with numerical simulations of the one-dimensional Poisson-Nernst-Planck equations. We show how the interplay between diffusion coefficients and ionic valencies significantly varies the system response and why it isessential
What carries the argument
The zero-current transmembrane voltage V_{i=0} obtained by setting the total ionic current to zero in the Nernst-Planck flux equations under Poisson electrostatics.
If this is right
- The expressions reduce to several previously published models in appropriate limits.
- The voltage depends on the full set of diffusion coefficients and ionic valences rather than on a single effective parameter.
- The models can be used to interpret experimental ion-transport measurements in nanoporous membranes.
- The two-species exact solution provides a benchmark against which the linear-profile approximation can be tested.
Where Pith is reading between the lines
- The exact two-species solution may serve as a test case for developing higher-order approximations that avoid the linear-profile assumption entirely.
- Because the voltage enters energy-harvesting and desalination efficiency calculations, small changes in predicted V_{i=0} translate directly into changes in predicted power density.
- The same framework could be applied to physiological ion channels once surface-charge and geometry parameters are supplied.
Load-bearing premise
The multi-species expression relies on assuming a linear concentration profile across the channel without detailed justification or bounds on its validity.
What would settle it
A direct numerical solution of the one-dimensional Poisson-Nernst-Planck equations for a multi-species electrolyte with non-linear concentration profiles that deviates from the analytical V_{i=0} prediction by more than numerical discretization error.
read the original abstract
The transmembrane voltage, $V$, which is the potential drop required to nullify the electrical current ($i=0$), is a key characteristic of water desalination and energy harvesting systems that utilize macroscopically large nanoporous membranes, as well as for physiological ion channels subjected to asymmetric salt concentrations. To date, existing analytical expressions for $V_{i=0}$ have been limited to simple scenarios or under simplifying assumptions. In this work, we derive two expressions for $V_{i=0}$. First, we consider the much simpler scenario of two species. Then, we can consider an electrolyte comprised of an arbitrary number of species. The difference in the models is that the latter solution utilizes an ad-hoc assumption of a linear concentration profile, while the former solution does not require such an ad-hoc assumption. We analyze both models and show how to reduce them to several known models. We also verify both models with numerical simulations of the one-dimensional Poisson-Nernst-Planck equations. We show how the interplay between diffusion coefficients and ionic valencies significantly varies the system response and why it is essential to account for all system parameters. Ultimately, this model can be used to improve experimental interpretation of ion transport measurements
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives closed-form expressions for the zero-current transmembrane voltage V_{i=0} in a charged nanochannel with asymmetric electrolytes. For binary electrolytes the expression follows directly from the steady-state PNP current-null condition. For an arbitrary number of species an ad-hoc linear concentration profile is imposed to close the integrated flux equations, yielding an explicit formula. Both expressions are reduced to known limits and compared to numerical solutions of the one-dimensional Poisson-Nernst-Planck equations. The work stresses the interplay between diffusion coefficients and ionic valences.
Significance. If the linear-profile approximation holds in the regimes of interest, the results supply practical analytical tools for interpreting transmembrane potentials in nanoporous membranes used for desalination and energy harvesting as well as in physiological ion channels. The binary result is free of extra assumptions and directly verifiable; the multi-species result extends generality but its accuracy is conditional on the fidelity of the imposed profile. Numerical checks are provided, yet the absence of quantified error bounds restricts immediate use for arbitrary electrolytes.
major comments (2)
- [§3.2] §3.2 (multi-species derivation): the linear concentration profile c_i(x) is introduced without derivation, asymptotic analysis (e.g., small-Debye-length or weak-charge limits), or error estimates as a function of surface charge, concentration ratio, or valence. Because this assumption closes the integrated Nernst-Planck fluxes for the arbitrary-species case, the headline claim that the model works for “an electrolyte comprised of an arbitrary number of species” rests entirely on the unquantified accuracy of this profile.
- [Numerical verification] Numerical verification section: comparisons with 1D PNP simulations are shown, but the deviation between the linear-profile formula and the full numerical solution is not mapped versus Debye length, surface charge density, or valence combinations. Consequently the practical range of validity of the multi-species expression remains uncharted.
minor comments (2)
- [Abstract] Abstract: the statement that both models “reduce to several known models” is not accompanied by explicit citations or section references; the reductions should be identified.
- [§3.2] Notation: the definition of the effective diffusion coefficient or averaged mobility used in the multi-species formula should be stated explicitly rather than left implicit.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to strengthen the presentation of the multi-species approximation and its validation.
read point-by-point responses
-
Referee: [§3.2] §3.2 (multi-species derivation): the linear concentration profile c_i(x) is introduced without derivation, asymptotic analysis (e.g., small-Debye-length or weak-charge limits), or error estimates as a function of surface charge, concentration ratio, or valence. Because this assumption closes the integrated Nernst-Planck fluxes for the arbitrary-species case, the headline claim that the model works for “an electrolyte comprised of an arbitrary number of species” rests entirely on the unquantified accuracy of this profile.
Authors: We agree that the linear profile is an ad-hoc closure assumption rather than an exact solution of the PNP system. It is motivated by the expectation that, for channels much longer than the Debye length, the concentration variation between reservoirs is approximately linear when the net current is zero. In the revised manuscript we add a dedicated paragraph in §3.2 that (i) derives the linear profile as the leading-order solution in the thin-double-layer limit, (ii) provides the corresponding asymptotic expansion of V_{i=0}, and (iii) reports relative-error bounds obtained from the existing PNP simulations for representative values of surface charge, concentration ratio, and valence. The abstract and introduction are also updated to state explicitly that the multi-species formula is approximate. revision: yes
-
Referee: [Numerical verification] Numerical verification section: comparisons with 1D PNP simulations are shown, but the deviation between the linear-profile formula and the full numerical solution is not mapped versus Debye length, surface charge density, or valence combinations. Consequently the practical range of validity of the multi-species expression remains uncharted.
Authors: We accept this criticism. The revised numerical section now includes two additional figures: one showing the relative error of the multi-species formula versus normalized Debye length (spanning two orders of magnitude in bulk concentration) at fixed surface charge, and a second contour plot of the same error over the plane of surface charge density versus valence asymmetry. These maps are generated from the same 1D PNP solver already used in the manuscript and are accompanied by a short table of maximum errors for the parameter ranges typical of desalination and biological channels. revision: yes
Circularity Check
No significant circularity; derivations start from PNP equations and remain independent of target result
full rationale
The paper derives the two V_{i=0} expressions directly from the steady-state Poisson-Nernst-Planck current-null condition. The two-species case closes without extra assumptions. The arbitrary-species case introduces an explicit ad-hoc linear c_i(x) profile to integrate the flux equations, but this is stated as an assumption rather than obtained by re-arranging the target voltage or by fitting. Verification uses separate numerical PNP solutions. No self-citations, fitted parameters, or self-definitional steps appear; the expressions are not equivalent to their inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- ad hoc to paper Linear concentration profile across the nanochannel for the multi-species model
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We derive two expressions for V_{i=0}... the latter solution utilizes an ad-hoc assumption of a linear concentration profile
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the transport of the two species is governed by the Poisson-Nernst-Planck (PNP) equations
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
S. J. Kim, S. H. Ko, K. H. Kang, J. Han, Nature Nanotech 2010, 5, 297
work page 2010
-
[2]
V. V. Nikonenko, A. V. Kovalenko, M. K. Urtenov, N. D. Pismenskaya, J. Han, P. Sistat, G. Pourcelly, Desalination 2014, 342, 85
work page 2014
- [3]
-
[4]
Q. Wen, P. Jia, L. Cao, J. Li, D. Quan, L. Wang, Y. Zhang, D. Lu, L. Jiang, W. Guo, Advanced Materials 2020, 32, 1903954
work page 2020
-
[5]
F. H. J. van der Heyden, D. J. Bonthuis, D. Stein, C. Meyer, C. Dekker, Nano Lett. 2007, 7, 1022
work page 2007
- [6]
-
[7]
Q. Wang, X. Zhang, H. Zhu, X. Zhang, Q. Liu, M. Fu, J. Zhu, J. Pu, Z. Qu, Advanced Physics Research 2023, 2, 2300016
work page 2023
-
[8]
Z. Ding, T. Gu, M. Zhang, K. Wang, D. Sun, J. Li, Small 2024, 20, 2403593
work page 2024
-
[9]
A. Alcaraz, M. L. López, M. Queralt -Martín, V. M. Aguilella, ACS Nano 2017, 11, 10392
work page 2017
- [10]
-
[11]
M. J. Gallenito, M. T. Clabbers, J. Lin, J. Hattne, T. Gonen, Advanced Science 2025, 12, e04881
work page 2025
-
[12]
I. Vlassiouk, T. R. Kozel, Z. S. Siwy, J. Am. Chem. Soc. 2009, 131, 8211
work page 2009
- [13]
- [14]
-
[15]
M. Aramesh, C. Forró, L. Dorwling -Carter, I. Lüchtefeld, T. Schlotter, S. J. Ihle, I. Shorubalko, V. Hosseini, D. Momotenko, T. Zambelli, E. Klotzsch, J. Vörös, Nat. Nanotechnol. 2019, 14, 791
work page 2019
- [16]
- [17]
-
[18]
Y. Shu, K. Yuan, Z. Xiang, P. Chen, H. Wang, D. Ye, Advanced Materials n.d., n/a, e10095
-
[19]
D. Branton, D. W. Deamer, A. Marziali, H. Bayley, S. A. Benner, T. Butler, M. Di Ventra, S. Garaj, A. Hibbs, X. Huang, S. B. Jovanovich, P. S. Krstic, S. Lindsay, X. S. Ling, C. H. Mastrangelo, A. Meller, J. S. Oliver, Y. V. Pershin, J. M. Ramsey, R. Riehn, G. V. Soni, V. Tabard -Cossa, M. Wanunu, M. Wiggin, J. A. Schloss, Nat Biotech 2008, 26, 1146
work page 2008
- [20]
- [21]
- [22]
-
[23]
J. G. Hedley, K. Coshic, A. Aksimentiev, A. A. Kornyshev, Phys. Rev. X 2024, 14, 031042
work page 2024
-
[24]
T. Emmerich, N. Ronceray, K. V. Agrawal, S. Garaj, M. Kumar, A. Noy, A. Radenovic, Nat Rev Methods Primers 2024, 4, 69
work page 2024
-
[25]
N. Soni, Z. Rosenstock, N. C. Verma, K. Siddharth, N. Talor, J. Liu, B. Marom, A. B. Kolomeisky, A. Aksimentiev, A. Meller, Nat. Nanotechnol. 2025, 20, 1482
work page 2025
-
[26]
J. Feng, M. Graf, K. Liu, D. Ovchinnikov, D. Dumcenco, M. Heiranian, V. Nandigana, N. R. Aluru, A. Kis, A. Radenovic, Nature 2016, 536, 197
work page 2016
-
[27]
C. Zhu, P. Liu, B. Niu, Y. Liu, W. Xin, W. Chen, X.- Y. Kong, Z. Zhang, L. Jiang, L. Wen, J. Am. Chem. Soc. 2021, 143, 1932
work page 2021
-
[28]
Y.-C. Yao, A. Taqieddin, M. A. Alibakhshi, M. Wanunu, N. R. Aluru, A. Noy, ACS Nano 2019, 13, 12851
work page 2019
-
[29]
Hille, Ion Channels of Excitable Membranes , Sinauer Associate, Sunderland, Mass, 2001
B. Hille, Ion Channels of Excitable Membranes , Sinauer Associate, Sunderland, Mass, 2001
work page 2001
-
[30]
Teorell, Proceedings of the Society for Experimental Biology and Medicine 1935, 33, 282
T. Teorell, Proceedings of the Society for Experimental Biology and Medicine 1935, 33, 282
work page 1935
-
[31]
K. H. Meyer, J. -F. Sievers, Helvetica Chimica Acta 1936, 19, 649
work page 1936
-
[32]
D. E. Goldman, J Gen Physiol 1943, 27, 37
work page 1943
-
[33]
A. L. Hodgkin, B. Katz, The Journal of Physiology 1949, 108, 37
work page 1949
-
[34]
A. H. Galama, J. W. Post, H. V. M. Hamelers, V. V. Nikonenko, P. M. Biesheuvel, Journal of Membrane Science and Research 2016, 2, 128
work page 2016
-
[35]
O. Lavi, Y. Green, Commun Mater 2024, 5, 1
work page 2024
- [36]
-
[37]
R. Abu -Rjal, V. Chinaryan, M. Z. Bazant, I. Rubinstein, B. Zaltzman, Phys. Rev. E 2014, 89, 012302
work page 2014
- [38]
- [39]
-
[40]
Rubinstein, Electro-Diffusion of Ions, Society For Industrial And Applied Mathematics, 1990
I. Rubinstein, Electro-Diffusion of Ions, Society For Industrial And Applied Mathematics, 1990
work page 1990
-
[41]
W. Chen, T. Dong, Y. Xiang, Y. Qian, X. Zhao, W. Xin, X. -Y. Kong, L. Jiang, L. Wen, Advanced Materials 2022, 34, 2108410
work page 2022
-
[42]
S. N. Sankar, J. Fernandes, M. B. dos Santos, B. Espiña, P. Alpuim, A. G. Díez, S. Lanceros-Mendez, L. Saini, S. Kaushik, G. Kalon, A. Capasso, Advanced Functional Materials 2023, 33, 2214889
work page 2023
-
[43]
Z. Liu, L. Li, L. Qin, S. Guo, G. Fang, Z. Luo, S. Liang, Advanced Materials 2022, 34, 2204681
work page 2022
- [44]
-
[45]
P. M. Biesheuvel, M. Z. Bazant, Phys. Rev. E 2016, 94, 050601
work page 2016
- [46]
-
[47]
Z. Li, A. Noy, Chem. Soc. Rev. 2025, 54, 8582
work page 2025
-
[48]
Green, 2025, DOI 10.48550/arXiv.2512.15284
Y. Green, 2025, DOI 10.48550/arXiv.2512.15284
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.