Effects of Surface Trapping and Contact Ion Pairing on Ion Transport in Nanopores
Pith reviewed 2026-05-24 15:10 UTC · model grok-4.3
The pith
Contact ion pairs reduce ion mobility in neutral nanopores while surface trapping affects charged nanopores, with a modified first-passage time model accounting for both.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In nanopores without surface charges the formation of contact ion pairs plays a critical role in reducing ion mobility, with mobility for both cations and anions decreasing as pore size shrinks because pair formation becomes easier in smaller neutral pores. Inside charged nanopores, surface charges reduce counterion mobility through surface trapping in addition to pair formation; Na+ mobility rises then falls with increasing surface charge density while Cl- shows the opposite trend. A modified first-passage time model that incorporates ion pair formation and trapped ions describes ion transport through a nanopore.
What carries the argument
Modified first-passage time model that accounts for both contact ion pair formation and surface-trapped ions.
If this is right
- Mobility decreases with reduced pore diameter in neutral nanopores because contact ion pairs form more readily.
- In charged nanopores Na+ mobility first increases then decreases with surface charge density while Cl- mobility follows the opposite trend.
- The modified first-passage time model supplies a quantitative description of transport that includes both pairing and trapping.
Where Pith is reading between the lines
- Device designs for selective ion passage could target specific charge densities to exploit the non-monotonic mobility curves.
- The same pairing-plus-trapping picture may apply to other monovalent ion pairs or to pores with different wall chemistries.
- Testing the model against measured transit times in pores of controlled charge would directly check its predictive range.
Load-bearing premise
The all-atom molecular dynamics force fields and pore models accurately reproduce real ion-ion and ion-surface interactions without significant artifacts from finite-size effects or incomplete sampling.
What would settle it
Experimental measurement of Na+ and Cl- mobility versus pore diameter in neutral nanopores and versus surface charge density in charged nanopores, testing whether mobility decreases monotonically with size in neutral pores and shows the predicted non-monotonic charge-density trends in charged pores.
Figures
read the original abstract
Ion transport in highly-confined space is important to various applications, such as biosensing and seawater desalination with nanopores. All-atom molecular dynamics simulations are conducted to investigate the transport of Na$^+$ and Cl$^-$ ions through nanopores with the diameter below 6 nm. It is found that the formation of the contact ion pair plays a critical role in reducing the ion mobility inside a nanopore without surface charges. The mobility for both cations and anions decreases with the reduced pore size because it is easier to form the contact ion pairs inside the neutral nanopore with smaller diameter. Inside a charged nanopores, besides the contact ion pair formation, the surface charges also play a significant role in reducing the counterion mobility through surface trapping. It is uncovered that the mobility of Na$^+$ ions increases first and then decreases with the surface charge density, while Cl$^-$ ions have the opposite trend. A modified first-passage time model is proposed to take into account the ion pair formation and the trapped ions inside a nanopore, which provides a clear picture in describing ion transport through a nanopore.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports all-atom molecular dynamics simulations of Na+ and Cl- transport through nanopores with diameters below 6 nm. It claims that contact ion pair formation is the dominant mechanism reducing ion mobility in neutral pores (with the effect strengthening at smaller diameters), while in charged pores surface trapping by the wall charges additionally reduces counterion mobility, producing non-monotonic mobility versus surface-charge-density curves that differ in sign for cations and anions. A modified first-passage-time model is introduced to incorporate both ion-pairing and surface-trapping effects.
Significance. If the simulation trends prove robust, the work supplies a concrete mechanistic picture of two distinct confinement-induced slowing mechanisms that are directly relevant to nanopore-based desalination and biosensing. The modified first-passage model offers a compact interpretive tool that could be used for rapid parameter exploration once calibrated. The significance is limited by the absence of force-field validation or quantitative error bars on the reported mobilities.
major comments (2)
- [Methods and Results] The central attribution of mobility reduction to contact-ion-pair formation (neutral pores) and surface trapping (charged pores) rests on the assumption that the chosen all-atom force fields correctly capture ion-ion and ion-surface energetics. No comparison of simulated contact-pair association constants or surface adsorption free energies to experimental benchmarks is provided, leaving open the possibility that the observed trends are force-field artifacts rather than general physical effects.
- [Abstract and Results] The abstract states that the reported trends 'come directly from the simulations,' yet the manuscript supplies neither tabulated mobility values with statistical uncertainties nor convergence tests with respect to simulation length or system size. Without these data it is impossible to judge whether the non-monotonic mobility-versus-charge-density curves or the pore-size dependence are statistically significant.
minor comments (2)
- [Abstract] The phrase 'Inside a charged nanopores' in the abstract is grammatically incorrect and should read 'Inside charged nanopores.'
- [Figures] Figure captions should explicitly state the number of independent trajectories and the total sampling time used to compute each mobility datum.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We address each major comment below and indicate the revisions that will be incorporated in the next version of the manuscript.
read point-by-point responses
-
Referee: [Methods and Results] The central attribution of mobility reduction to contact-ion-pair formation (neutral pores) and surface trapping (charged pores) rests on the assumption that the chosen all-atom force fields correctly capture ion-ion and ion-surface energetics. No comparison of simulated contact-pair association constants or surface adsorption free energies to experimental benchmarks is provided, leaving open the possibility that the observed trends are force-field artifacts rather than general physical effects.
Authors: We agree that explicit validation against experimental ion-pairing constants and surface adsorption energies would strengthen the manuscript. The simulations use standard CHARMM parameters for ions together with TIP3P water, which have been employed in multiple prior studies of confined electrolytes. While we did not recompute bulk association constants here, the contact-pairing behavior we observe is consistent with the known increase in pairing under confinement. In the revised manuscript we will add a dedicated paragraph in the Methods section that cites literature benchmarks for these force fields in nanopore and bulk settings and explicitly states that the reported trends are qualitative. No new simulations are required for this addition. revision: partial
-
Referee: [Abstract and Results] The abstract states that the reported trends 'come directly from the simulations,' yet the manuscript supplies neither tabulated mobility values with statistical uncertainties nor convergence tests with respect to simulation length or system size. Without these data it is impossible to judge whether the non-monotonic mobility-versus-charge-density curves or the pore-size dependence are statistically significant.
Authors: The referee correctly identifies a presentational gap. In the revised manuscript we will (i) add a table in the main text or SI that lists the computed mobilities together with standard errors obtained from block averaging over independent trajectory segments, and (ii) include a short convergence subsection (or SI note) that reports mobility values as a function of simulation length and system size for representative pore diameters and charge densities. These additions will allow readers to assess the statistical robustness of the non-monotonic trends and the pore-size dependence. revision: yes
Circularity Check
No significant circularity; claims rest on MD simulation outputs and interpretive model
full rationale
The paper derives its claims from all-atom MD trajectories examining ion mobility, contact pair formation, and surface trapping effects in nanopores of varying size and charge. The modified first-passage time model is introduced post-simulation as an interpretive framework to account for those observed mechanisms, without any equation reducing by construction to a fitted parameter or self-citation chain. No load-bearing self-definitional steps, uniqueness theorems, or ansatz smuggling appear. The central results are externally falsifiable via independent simulations or experiments and do not presuppose their own outputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption All-atom force fields and periodic boundary conditions accurately capture ion-ion and ion-wall interactions in sub-6 nm pores.
Reference graph
Works this paper leans on
-
[1]
Deamer, D.; Akeson, M.; Branton, D., Three Decades of Nanopore Sequencing. Nat. Biotechnol. 2016, 34, 518-524
work page 2016
-
[2]
Rollings, R. C.; Kuan, A. T.; Golovchenko, J. A., Ion Selectivity of Graphene Nanopores. Nat. Commun. 2016, 7, 11408
work page 2016
-
[3]
Zhang, Z.; Sui, X.; Li, P.; Xie, G.; Kong, X. Y .; Xiao, K.; Gao, L.; Wen, L.; Jiang, L., 19 Ultrathin and Ion -Selective Janus Membranes for High -Performance Osmotic Energy Conversion. J. Am. Chem. Soc. 2017, 139, 8905-8914
work page 2017
-
[4]
J.; Dekker, C., Graphene Nanodevices for DNA Sequencing
Heerema, S. J.; Dekker, C., Graphene Nanodevices for DNA Sequencing. Nat. Nanotechnol. 2016, 11, 127-136
work page 2016
-
[5]
Liu, K.; Feng, J.; Kis, A.; Radenovic, A., Atomically Thin Molybdenum Disulfide Nanopores with High Sensitivity for DNA Translocation. ACS Nano 2014, 8, 2504-2511
work page 2014
-
[6]
O'Hern, S. C .; Boutilier, M. S.; Idrobo, J. C.; Song, Y .; Kong, J.; Laoui, T.; Atieh, M.; Karnik, R., Selective Ionic Transport through Tunable Subnanometer Pores in Single -Layer Graphene Membranes. Nano Lett. 2014, 14, 1234-1241
work page 2014
-
[7]
Surwade, S. P.; Smirnov, S. N.; Vl assiouk, I. V .; Unocic, R. R.; Veith, G. M.; Dai, S.; Mahurin, S. M., Water Desalination Using Nanoporous Single -Layer Graphene. Nat. Nanotechnol. 2015, 10, 459-464
work page 2015
-
[8]
D.; Graf, M.; Liu, K.; Ovchinnikov, D.; Dumcenco, D.; Heiranian, M.; Nandigana, V .; Aluru, N
Feng, J. D.; Graf, M.; Liu, K.; Ovchinnikov, D.; Dumcenco, D.; Heiranian, M.; Nandigana, V .; Aluru, N. R.; Kis, A.; Radenovic, A., Single -Layer Mos2 Nanopores as Nanopower Generators. Nature 2016, 536, 197-200
work page 2016
-
[9]
N., Chapter 14 - Electrostatic Forces between Surfaces in Liquids
Israelachvili, J. N., Chapter 14 - Electrostatic Forces between Surfaces in Liquids. In Intermolecular and Surface Forces (Third Edition), Israelachvili, J. N., Ed. Academic Press: San Diego, 2011; pp 291-340
work page 2011
-
[10]
Stein, D.; Kruithof, M.; Dekker, C., Surface -Charge-Governed Ion Transport in Nanofluidic Channels. Phys. Rev. Lett. 2004, 93, 035901
work page 2004
-
[11]
S.; Luijten, E., Dielectric Modulation of Ion Transport near Interfaces
Antila, H. S.; Luijten, E., Dielectric Modulation of Ion Transport near Interfaces. Phys. Rev. Lett. 2018, 120, 135501
work page 2018
-
[12]
Pang, P.; He, J.; Park, J. H.; Krstić, P. S.; Lindsay, S., Origin of Giant Ionic Currents in Carbon Nanotube Channels. ACS Nano 2011, 5, 7277-7283
work page 2011
-
[13]
Cheng, L. J.; Guo, L. J., Nanofluidic Diodes. Chem. Soc. Rev. 2010, 39, 923-938
work page 2010
-
[14]
S., Ion-Current Rectification in Nanopores and Nanotubes with Broken Symmetry
Siwy, Z. S., Ion-Current Rectification in Nanopores and Nanotubes with Broken Symmetry. Adv. Funct. Mater. 2006, 16, 735-746
work page 2006
-
[15]
Kis, A.; Di Ventra, M.; Radenovic, A., Observation of Ionic Coulomb Blockade in Nanopores
Feng, J.; Liu, K.; Graf, M.; Dumcenco, D. ; Kis, A.; Di Ventra, M.; Radenovic, A., Observation of Ionic Coulomb Blockade in Nanopores. Nat. Mater. 2016, 15, 850-855
work page 2016
-
[16]
Qi, H.; Li, Z.; Tao, Y .; Zhao, W.; Lin, K.; Ni, Z.; Jin, C.; Zhang, Y .; Bi, K.; Chen, Y ., 20 Fabrication of Sub -Nanometer Pores on Graphene Membrane for Ion Selective Transport. Nanoscale 2018, 10, 5350-5357
work page 2018
-
[17]
Jain, T.; Rasera, B. C.; Guerrero, R. J.; Boutilier, M. S.; O'Hern, S. C.; Idrobo, J. C.; Karnik, R., Heterogeneous Sub -Continuum Ionic Transport in Statistically Isolated Graphene Nanopores. Nat. Nanotechnol. 2015, 10, 1053-1057
work page 2015
-
[18]
Duan, C.; Majumdar, A., Anomalous Ion Transport in 2 -Nm Hydrophilic Nanochannels. Nat. Nanotechnol. 2010, 5, 848-852
work page 2010
-
[19]
Esfandiar, A.; Radha, B.; Wang, F.; Yang, Q.; Hu, S.; Garaj, S.; Nair, R.; Geim, A.; Gopinadhan, K., Size Effect in Ion Transport through Angstrom-Scale Slits. Science 2017, 358, 511-513
work page 2017
-
[20]
Ma, J., et al., Drastically Reduced Ion Mobility in a Nanopore Due to Enhanced Pairing and Collisions between Dehydrated Ions. J. Am. Chem. Soc. 2019, 141, 4264-4272
work page 2019
-
[21]
Redner, S., A Guide to First-Passage Processes; Cambridge University Press, 2001
work page 2001
-
[22]
Shankla, M.; Aksimentiev, A., Conformational Transitions and Stop -and-Go Nanopore Transport of Single-Stranded DNA on Charged Graphene. Nat. Commun. 2014, 5, 5171
work page 2014
-
[23]
Humphrey, W.; Dalke, A.; Schulten, K., Vmd: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33-38
work page 1996
-
[24]
C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E
Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E. ; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K., Scalable Molecular Dynamics with Namd. J. Comput. Chem. 2005, 26, 1781-1802
work page 2005
-
[25]
L.; Chandrasekhar, J.; Madura, J
Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L., Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926-935
work page 1983
-
[26]
D.; Bashford, D.; Bellott, M.; Dunbrack Jr, R
MacKerell Jr, A. D.; Bashford, D.; Bellott, M.; Dunbrack Jr, R. L.; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S., All -Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586-3616
work page 1998
-
[27]
B.; Belkin, M.; Comer, J.; Aksimentiev, A., Assessing Graphene Nanopores for Sequencing DNA
Wells, D. B.; Belkin, M.; Comer, J.; Aksimentiev, A., Assessing Graphene Nanopores for Sequencing DNA. Nano Lett. 2012, 12, 4117-4123
work page 2012
-
[28]
L.; Darden, T.; Lee, H.; Pedersen, L
Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G., A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577-8593. 21
work page 1995
-
[29]
A., Settle - an Analytical Version of the Shake and Rattle Algorithm for Rigid Water Models
Miyamoto, S.; Kollman, P. A., Settle - an Analytical Version of the Shake and Rattle Algorithm for Rigid Water Models. J. Chem. Phys. 1992, 13, 952-962
work page 1992
-
[30]
Martyna, G. J.; Tobias, D. J.; Klein, M. L., Constant -Pressure Molecular -Dynamics Algorithms. J. Chem. Phys. 1994, 101, 4177-4189
work page 1994
-
[31]
B.; Timp, G.; Schulten, K., Microscopic Kinetics of DNA Translocation through Synthetic Nanopores
Aksimentiev, A.; Heng, J. B.; Timp, G.; Schulten, K., Microscopic Kinetics of DNA Translocation through Synthetic Nanopores. Biophys. J. 2004, 87, 2086-2097
work page 2004
-
[32]
Si, W.; Aksimentiev, A., Nanopore Sensing of Protein Folding. ACS Nano 2017, 11, 7091- 7100
work page 2017
- [33]
-
[34]
R., Scaling of Electrokinetic Transport in Nanometer Channels
Qiao, R.; Aluru, N. R., Scaling of Electrokinetic Transport in Nanometer Channels. Langmuir 2005, 21, 8972-8977
work page 2005
-
[35]
Aksimentiev, A.; Schulten, K., Imaging Alpha-Hemolysin with Molecular Dynamics: Ionic Conductance, Osmotic Permeability, and the Electrostatic Potential Map. Biophys. J. 2005, 88, 3745-3761
work page 2005
-
[36]
E., Access Resistance of a Small Circular Pore
Hall, J. E., Access Resistance of a Small Circular Pore. J Gen Physiol 1975, 66, 531-2
work page 1975
-
[37]
Wanunu, M.; Dadosh, T.; Ray, V .; Jin, J. M.; McReynolds, L.; Drndic, M., Rapid Electronic Detection of Probe-Specific Micrornas Using Thin Nanopore Sensors. Nat. Nanotechnol. 2010, 5, 807-814
work page 2010
-
[38]
Smeets, R. M. M.; Keyser, U. F.; Krapf, D.; Wu, M. Y .; Dekker, N. H.; Dekker, C., Salt Dependence of Ion Transport and DNA Translocation through Solid -State Nanopores. Nano Lett. 2006, 6, 89-95
work page 2006
-
[39]
Lee, C.; Joly, L.; Siria, A.; Biance, A. L.; Fulcrand, R.; Bocquet, L., Large Apparent Electric Size of Solid-State Nanopores Due to Spatially Extended Surface Conduction. Nano Lett. 2012, 12, 4037-4044
work page 2012
-
[40]
Bocquet, L.; Charlaix, E., Nanofluidics, from Bulk to Interfaces. Chem. Soc. Rev. 2010, 39, 1073-1095
work page 2010
-
[41]
R., Charge Inversion and Flow Reversal in a Nanochannel Electro - Osmotic Flow
Qiao, R.; Aluru, N. R., Charge Inversion and Flow Reversal in a Nanochannel Electro - Osmotic Flow. Phys. Rev. Lett. 2004, 92, 198301
work page 2004
-
[42]
Grosberg, A. Y .; Nguyen, T. T.; Shklovskii, B. I., Colloquium: The Physics of Charge 22 Inversion in Chemical and Biological Systems. Reviews of Modern Physics 2002, 74, 329-345
work page 2002
-
[43]
Li, Z.; Qiu, Y .; Li, K.; Sha, J.; Li, T.; Chen, Y ., Optimal Design of Graphene Nanopores for Seawater Desalination. J. Chem. Phys. 2018, 148, 014703
work page 2018
-
[44]
Report on Conductivity of Strong Electrolytes in Dilute Solutions
Debye, P., Introductory Paper—Part I: Mobilities of Ions. Report on Conductivity of Strong Electrolytes in Dilute Solutions. Trans. Faraday Soc. 1927, 23, 334-340
work page 1927
-
[45]
Physics of Fluids 2011, 23, 072005
Wu, P.; Qiao, R., Physical Origins of Apparently Enhanced Viscosity of Interfacial Fluids in Electrokinetic Transport. Physics of Fluids 2011, 23, 072005
work page 2011
-
[46]
Marcus, Y .; Hefter, G., Ion Pairing. Chem. Rev. 2006, 106, 4585-4621
work page 2006
-
[47]
Schoch, R. B.; Han, J. Y .; Renaud, P., Transport Phenomena in Nanofluidics. Reviews of Modern Physics 2008, 80, 839-883
work page 2008
-
[48]
Suk, M. E.; Aluru, N. R., Ion Transport in Sub-5-Nm Graphene Nanopores. J. Chem. Phys. 2014, 140, 084707
work page 2014
-
[49]
Meller, A.; Nivon, L.; Branton, D., V oltage -Driven DNA Translocations through a Nanopore. Phys. Rev. Lett. 2001, 86, 3435-3438
work page 2001
-
[50]
Yeh, I.-C.; Hummer, G., Nucleic Acid Transport through Carbon Nanotube Membranes. P . Natl. Acad. Sci. U.S.A. 2004, 101, 12177-12182
work page 2004
-
[51]
Mathé, J.; Aksimentiev, A.; Nelson, D. R.; Schulten, K.; Meller, A., Orientation Discrimination of Single-Stranded DNA inside the Α-Hemolysin Membrane Channel. P . Natl. Acad. Sci. U.S.A. 2005, 102, 12377-12382
work page 2005
-
[52]
-C.; Bhattacharya, A., Sequence Dependence of DNA Translocation through a Nanopore
Luo, K.; Ala -Nissila, T.; Ying, S. -C.; Bhattacharya, A., Sequence Dependence of DNA Translocation through a Nanopore. Phys. Rev. Lett. 2008, 100, 058101
work page 2008
-
[53]
S.; Li, J., Single -Molecule Protein Unfolding in Solid State Nanopores
Talaga, D. S.; Li, J., Single -Molecule Protein Unfolding in Solid State Nanopores. J. Am. Chem. Soc. 2009, 131, 9287-9297
work page 2009
-
[54]
Chen, Z., et al., DNA Translocation through an Array of Kinked Nanopores. Nat. Mater. 2010, 9, 667-675
work page 2010
-
[55]
Single-Molecule Protein Unfolding in Solid State Nanopores
Talaga, D. S.; Li, J., Correction to “Single-Molecule Protein Unfolding in Solid State Nanopores”. J. Am. Chem. Soc. 2013, 135, 13220-13220
work page 2013
-
[56]
Y .; Muthukumar, M.; Rosenstein, J
Larkin, J.; Henley, R. Y .; Muthukumar, M.; Rosenstein, J. K.; Wanunu, M., High - Bandwidth Protein Analysis Using Solid-State Nanopores. Biophys. J. 2014, 106, 696-704
work page 2014
-
[57]
McMullen, A.; de Haan, H. W.; Tang, J. X.; Stein, D., Stiff Filamentous Virus 23 Translocations through Solid-State Nanopores. Nat. Commun. 2014, 5, 4171
work page 2014
-
[58]
Ling, D. Y .; Ling, X. S., On the Distribution of DNA Translocat ion Times in Solid -State Nanopores: An Analysis Using Schrödinger ’s First-Passage-Time Theory. J. Phys.: Condens. Mat. 2013, 25, 375102
work page 2013
-
[59]
Risken, B., The Fokker-Planck Equation : Methods of Solution and Applications; Springer- Verlag, 1984, p 1206-1207
work page 1984
-
[60]
Reimann, P.; Meyer, A.; Getfert, S., On the Lubensky -Nelson Model of Polymer Translocation through Nanopores. Biophys. J. 2012, 103, 889-897
work page 2012
-
[61]
Metzler, R.; Redner, S.; Oshanin, G., First-Passage Phenomena and Their Applications ; World Scientific, 2014
work page 2014
-
[62]
Mittal, J.; Truskett, T. M.; Errington, J. R.; Hummer, G., Layering and Position-Dependent Diffusive Dynamics of Confined Fluids. Phys. Rev. Lett. 2008, 100, 145901
work page 2008
-
[63]
Qiu, Y . H.; Ma, J.; Chen, Y . F., Ionic Behavior in Highly Concentrated Aqueous Solutions Nanoconfined between Discretely Charged Silicon Surfaces. Langmuir 2016, 32, 4806-4814. 24 Figure for Table of Contents Table of Contents graphic showing the ion mobility as a function of surface charge density along with the schematic of surface trapping and contac...
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.