pith. sign in

arxiv: 2602.10600 · v2 · submitted 2026-02-11 · ❄️ cond-mat.stat-mech

Interplay of ion availability and mobility in the loss of cation selectivity for CaCltextsubscript{2} in negatively charged nanopores: molecular dynamics using scaled-charge models

Pith reviewed 2026-05-16 03:45 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords nanopore ion transportcharge inversionmultivalent electrolytesmolecular dynamicsperm-selectivityCaCl2 conductionelectrical double layerscaled-charge models
0
0 comments X

The pith

Negatively charged nanopores lose expected cation selectivity for CaCl2 because charge inversion immobilizes calcium ions near the walls and lets chloride ions dominate conduction in the interior.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses molecular dynamics simulations with scaled-charge models to examine NaCl and CaCl2 flow through negatively charged silica nanopores. It separates the concentration and velocity parts of the radial particle current density to link how ions adsorb statically to how they move under an applied field. NaCl follows the usual picture of cation preference, but CaCl2 does not. Instead, strong surface correlations invert the wall charge, pinning Ca2+ near the surface while Cl- carries most of the current through the pore center, producing bulk-like or even anion-favored transport.

Core claim

CaCl2 conduction through negatively charged nanopores is nearly bulk-like or anion-favored because Ca2+ ions become immobilized near the surface after charge inversion, leaving dominant Cl- conduction in the pore interior; this is shown by decomposing the radial current density into concentration and velocity contributions, while the same pores remain cation-selective for NaCl.

What carries the argument

Decomposition of the radial particle current density into separate concentration and velocity profiles, which directly connects equilibrium adsorption layers to nonequilibrium perm-selectivity.

If this is right

  • Monovalent salts retain conventional cation selectivity under the same surface charge and pore size.
  • Multivalent salts can produce charge-inverted transport whose sign depends on the relative mobility of the counter-ion after surface pinning.
  • The qualitative reversal is robust across reasonable force-field choices, but quantitative selectivity ratios shift with small changes in ion-surface interaction strength.
  • Pore design for selective transport must account for both static double-layer structure and the velocity distribution inside the pore rather than surface charge alone.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar immobilization effects could appear in biological channels or synthetic membranes exposed to divalent ions, altering expected rectification or filtering behavior.
  • The same current-density decomposition offers a general diagnostic for any electrolyte where adsorption and mobility compete.
  • Refining surface interaction parameters in scaled-charge models would allow quantitative prediction of the crossover concentration at which selectivity inverts.

Load-bearing premise

The scaled-charge force field correctly balances ion-surface attraction against ion-water solvation so that immobilization and charge inversion occur as observed.

What would settle it

Experimental or simulated radial profiles that show Ca2+ velocity dropping to near zero within one nanometer of the wall while Cl- velocity remains high in the pore center, together with measured reversal of the sign of the total ionic current for CaCl2.

Figures

Figures reproduced from arXiv: 2602.10600 by Dezs\H{o} Boda, Salman Shabbir, Zolt\'an Hat\'o.

Figure 1
Figure 1. Figure 1: FIG. 1. From left to right: radial dependence of the axial ( [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Radial distribution functions (RDF) for pairs of O [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. From top to bottom: radial dependence of the axial ( [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. From left to right: radial dependence of the axial ( [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. From left to right: radial dependence of the axial ( [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: The 1st peak in the OS-Ca2+ RDF is higher for the TIP4P/2005 model, while the reverse behavior is observed for the OS-Ow RDF. The OS-Ca2+ curve for SPC/E (left panel, blue color) is not something that we expect for the RDF between two oppositely charged par￾ticles. The 1st peak is too small, indicating a weak bind￾ing of Ca2+ ions to the OS atoms hindered by the SPC/E [PITH_FULL_IMAGE:figures/full_fig_p01… view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. From left to right: radial dependence of the axial ( [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Radial distribution functions for pairs of O [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Radial electroosmotic current density profiles [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
read the original abstract

Ion transport through charged nanopores is commonly interpreted in terms of electrical double layer structure, leading to the expectation of cation-selective conduction in negatively charged pores. This picture can break down for multivalent electrolytes, where strong ion-urface correlations and charge inversion modify transport behavior. Here, we study NaCl and CaCl$_2$ conduction through negatively charged silica nanopores using atomistic molecular dynamics simulations with scaled-charge ion models. By separating concentration and velocity contributions to the radial particle current density, we connect static adsorption to dynamic perm-selectivity. While NaCl exhibits conventional cation selectivity, CaCl$_2$ shows nearly bulk-like or even anion-favored transport due to Ca$^{2+}$ immobilization near the surface and dominant Cl$^-$ conduction in the pore interior following charge inversion. Although this qualitative mechanism is robust, its detailed manifestation depends sensitively on the balance of ion-surface and ion-water interactions encoded in the force field.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper uses atomistic molecular dynamics simulations with scaled-charge ion models to investigate NaCl and CaCl2 transport through negatively charged silica nanopores. It separates concentration and velocity contributions to the radial particle current density to link static ion adsorption profiles to dynamic perm-selectivity. While NaCl exhibits conventional cation selectivity consistent with the electrical double layer picture, CaCl2 shows nearly bulk-like or anion-favored conduction because Ca2+ ions immobilize near the pore wall following charge inversion, leaving Cl- to dominate transport in the pore interior. The abstract notes that the qualitative mechanism is robust but its quantitative details depend on the balance of ion-surface and ion-water interactions in the force field.

Significance. If the central mechanism holds, the work provides a concrete dynamical explanation for selectivity reversal in multivalent electrolytes inside charged nanopores, moving beyond static EDL interpretations. The separation of radial current into concentration and velocity fields is a useful technical contribution that directly connects adsorption structure to conductance. The study also underscores the practical importance of force-field choices in MD modeling of confined ion transport.

major comments (2)
  1. [Abstract and Results] Abstract and Results: The central claim that CaCl2 transport becomes bulk-like or anion-favored rests on scaled-charge models, yet the manuscript provides no cross-validation against full-charge or polarizable force fields. Given the abstract's own statement that detailed manifestation depends sensitively on ion-surface and ion-water balance, the absence of any quantification of how the radial current decomposition (concentration vs. velocity) shifts with alternative parameters leaves the load-bearing mechanism unanchored.
  2. [Methods and Results] Methods and Results: No direct comparison is made to experimental pore conductance or selectivity data for the same silica geometry and electrolyte conditions. Without such anchors, it is difficult to assess whether the reported immobilization of Ca2+ and interior Cl- conduction are physical or specific to the chosen scaling factors.
minor comments (1)
  1. [Results] Notation for the radial current density decomposition should be defined explicitly in the main text rather than only in supplementary material to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and positive assessment of our manuscript. We address each major comment below, providing clarifications on our methodological choices and the scope of the study.

read point-by-point responses
  1. Referee: [Abstract and Results] Abstract and Results: The central claim that CaCl2 transport becomes bulk-like or anion-favored rests on scaled-charge models, yet the manuscript provides no cross-validation against full-charge or polarizable force fields. Given the abstract's own statement that detailed manifestation depends sensitively on ion-surface and ion-water balance, the absence of any quantification of how the radial current decomposition (concentration vs. velocity) shifts with alternative parameters leaves the load-bearing mechanism unanchored.

    Authors: We appreciate the referee's emphasis on force-field sensitivity. The scaled-charge models were selected specifically because they reproduce experimental bulk ion activity coefficients and pairing behavior more accurately than full-charge models, which tend to overestimate electrostatic attractions and produce unphysical over-adsorption even in monovalent cases. While we agree that a systematic parameter sweep would be ideal, performing additional simulations with full-charge or polarizable models lies outside the scope of the present work. The radial current decomposition is presented to illustrate the mechanism (immobilized Ca2+ near the wall enabling interior Cl- flow) rather than to claim quantitative universality. We will add a short paragraph in the revised Methods and Discussion sections justifying the model choice with references to bulk validation studies and reiterating the abstract's caveat on quantitative details. revision: partial

  2. Referee: [Methods and Results] Methods and Results: No direct comparison is made to experimental pore conductance or selectivity data for the same silica geometry and electrolyte conditions. Without such anchors, it is difficult to assess whether the reported immobilization of Ca2+ and interior Cl- conduction are physical or specific to the chosen scaling factors.

    Authors: We concur that experimental anchoring would be valuable. However, to the best of our knowledge, no published conductance or selectivity measurements exist for silica nanopores with the precise diameter (~2 nm), surface charge density, and CaCl2 concentrations simulated here. Existing experimental work on multivalent ion transport in charged nanopores typically employs different pore materials (e.g., PET or alumina), larger diameters, or different surface chemistries, precluding direct quantitative comparison. We will expand the Discussion section to reference relevant experimental literature on selectivity reversal in multivalent electrolytes and clarify the limitations of direct comparison, while emphasizing that the qualitative mechanism (charge-inversion-induced Ca2+ immobilization) aligns with broader experimental observations of reduced cation selectivity in divalent salts. revision: partial

Circularity Check

0 steps flagged

No significant circularity: results from explicit MD trajectories

full rationale

The paper reports atomistic molecular dynamics simulations of NaCl and CaCl2 transport in charged silica nanopores using scaled-charge ion models. Central claims about cation vs. anion selectivity, charge inversion, and Ca2+ immobilization are obtained by decomposing radial particle current density into concentration and velocity profiles directly from simulation trajectories. No mathematical derivation chain exists that reduces predictions to inputs by construction, no parameters are fitted to a subset of data and then relabeled as predictions, and no self-citations serve as load-bearing uniqueness theorems. The work explicitly notes force-field sensitivity rather than claiming model-independent results. This is standard computational evidence generation, not circular reasoning.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard molecular dynamics assumptions plus the validity of scaled-charge ion models whose parameters are taken from prior literature.

free parameters (1)
  • ion charge scaling factors
    Scaled-charge models introduce empirical scaling of ion charges to account for polarization; these factors are fitted in the force field and directly affect immobilization and selectivity.
axioms (2)
  • domain assumption Force fields accurately capture ion-surface and ion-water interactions for the chosen silica model
    The abstract states that detailed manifestation depends sensitively on this balance.
  • standard math Periodic boundary conditions and finite pore length do not alter the radial current decomposition
    Implicit in any MD nanopore study.

pith-pipeline@v0.9.0 · 5488 in / 1275 out tokens · 48892 ms · 2026-05-16T03:45:08.886127+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

81 extracted references · 81 canonical work pages

  1. [1]

    Zheng, D

    Z. Zheng, D. J. Hansford, and A. T. Conlisk. Effect of multivalent ions on electroosmotic flow in micro- and nanochannels.Electrophoresis, 24(17):3006–3017, 2003

  2. [2]

    F. H. J. van der Heyden, D. Stein, K. Besteman, S. G. Lemay, and C. Dekker. Charge inversion at high ionic strength studied by streaming currents.Phys. Rev. Lett., 96(22):224502, 2006

  3. [3]

    J. Li, Q. Liang, L. Wen, Z. Jiang, and L. Jiang. Effects of ion size, ion valence and pH of electrolyte on electroos- motic flow reversal in fused silica capillaries.Analytica Chimica Acta, 1074:99—-108, 2019

  4. [4]

    S. X. Li, W. Guan, B. Weiner, and M. A. Reed. Direct observation of charge inversion in divalent nanofluidic de- vices.Nano Lett., 15(8):5046–5051, 2015

  5. [5]

    Lin, C.-Y

    K. Lin, C.-Y. Lin, J. W. Polster, Y. Chen, and Z. S. Siwy. Charge inversion and calcium gating in mixtures of ions in nanopores.J. Am. Chem. Soc., 142(6):2925– 2934, 2020

  6. [6]

    Siwy and A

    Z. Siwy and A. Fulinski. Fabrication of a synthetic nanopore ion pump.Phys. Rev. Lett., 89(19):198103, 2002

  7. [7]

    Siwy and S

    Z. Siwy and S. Howorka. Nanofluidics: Ion–current recti- fication in nanopores waits for transport studies.Chem- istry – A European Journal, 12:1317–1324, 2006

  8. [8]

    Y. He, D. Gillespie, D. Boda, I. Vlassiouk, R. S. Eisen- berg, and Z. S. Siwy. Tuning transport properties of nanofluidic devices with local charge inversion.JACS, 131(14):5194–5202, 2009

  9. [9]

    Gillespie, D

    D. Gillespie, D. Boda, Y. He, P. Apel, and Z.S. Siwy. Synthetic nanopores as a test case for ion channel theo- ries: The anomalous mole fraction effect without single filing.Biophys. J., 95(2):609–619, 2008

  10. [10]

    C. D. Lorenz and A. Travesset. Charge inversion of di- valent ionic solutions in silica channels.Phys. Rev. E, 75(6):061202, 2007

  11. [11]

    I. C. Bourg and G. Sposito. Molecular dynamics simu- lations of the electrical double layer on smectite surfaces contacting concentrated mixed electrolyte (nacl–cacl 2) solutions.J. Colloid Interf. Sci., 360(2):701–715, 2011

  12. [12]

    Hartkamp, B

    R. Hartkamp, B. Siboulet, J.-F. Dufrˆ eche, and B. Coasne. Ion-specific adsorption and electroosmosis in charged amorphous porous silica.Phys. Chem. Chem. Phys., 17(38):24683–24695, 2015

  13. [13]

    M. F. D¨ opke, J. L¨ utzenkirchen, O. A. Moultos, B. Si- boulet, J.-F. Dufrˆ eche, J. T. Padding, and R. Hartkamp. Preferential adsorption in mixed electrolytes confined by charged amorphous silica.J. Phys. Chem. C, 123(27):16711–16720, 2019

  14. [14]

    K. Wang, B. Siboulet, D. R´ ebiscoul, and J.-F. Dufrˆ eche. How ion pair formation drives adsorption in the electrical double layer: Molecular dynamics of charged silica–water interfaces in the presence of divalent alkaline earth ions. 15 J. Phys. Chem. C, 125(37):20551–20569, 2021

  15. [15]

    F´ abi´ an, Zs

    H. F´ abi´ an, Zs. Sarkadi, M. Valisk´ o, D. Gillespie, and D. Boda. Calcium versus potassium selectivity in a nanopore: The effect of charge inversion at localized pore charges.J. Mol. Liq., 368:120715, 2022

  16. [16]

    Rojano, D

    A. Rojano, D. Becerra, J. H. Walther, S. Prakash, and H. A. Zambrano. Effect of charge inversion on the elec- trokinetic transport of nanoconfined multivalent ionic so- lutions.Physics of Fluids, 36(10):102025, 2024

  17. [17]

    Hocine, R

    S. Hocine, R. Hartkamp, B. Siboulet, M. Duvail, B. Coasne, P. Turq, and J.-F. Dufrˆ eche. How ion conden- sation occurs at a charged surface: A molecular dynamics investigation of the stern layer for water–silica interfaces. J. Phys. Chem. C, 120(2):963–973, 2016

  18. [18]

    Malloggi, S

    F. Malloggi, S. ben Jabrallah, L. Girard, B. Siboulet, K. Wang, P. Fontaine, and J. Daillant. X-ray standing waves and molecular dynamics studies of ion surface in- teractions in water at a charged silica interface.J. Phys. Chem. C, 123(50):30294–30304, 2019

  19. [19]

    Selmani, B

    A. Selmani, B. Siboulet, M. ˇSpadina, Y. Foucaud, G. Draˇ zi´ c, B. Radatovi´ c, K. Korade, I. Nemet, D. Kovaˇ cevi´ c, J.-F. Dufrˆ eche, and K. Bohinc. Cation adsorption in tio 2 nanotubes: Implication for water de- contamination.ACS Appl. Nano Mater., 6(14):12711– 12725, 2023

  20. [20]

    Hartkamp, A.-L

    R. Hartkamp, A.-L. Biance, L. Fu, J.-F. Dufrˆ eche, O. Bonhomme, and L. Joly. Measuring surface charge: Why experimental characterization and molecular mod- eling should be coupled.Curr. Opin. Colloid Interface Sci., 37:101–114, 2018

  21. [21]

    Siboulet, S

    B. Siboulet, S. Hocine, R. Hartkamp, and J.-F. Dufrˆ eche. Scrutinizing electro-osmosis and surface con- ductivity with molecular dynamics.J. Phys. Chem. C, 121(12):6756–6769, 2017

  22. [22]

    Cervera, B

    J. Cervera, B. Schiedt, and P. Ram´ ırez. A Poisson/Nernst-Planck model for ionic transport through synthetic conical nanopores.Europhys. Lett., 71(1):35– 41, 2005

  23. [23]

    Matejczyk, M

    B. Matejczyk, M. Valisk´ o, M.-T. Wolfram, J.-F. Pietschmann, and D. Boda. Multiscale modeling of a rectifying bipolar nanopore: Comparing Poisson-Nernst- Planck to Monte Carlo.J. Chem. Phys., 146(12):124125, 2017

  24. [24]

    Boda and D

    D. Boda and D. Gillespie. Steady state electrodiffusion from the Nernst-Planck equation coupled to Local Equi- librium Monte Carlo simulations.J. Chem. Theor. Com- put., 8(3):824–829, 2012

  25. [25]

    Lorenz, Paul Crozier, Christian D

    Christian D. Lorenz, Paul Crozier, Christian D. Lorenz, Paul S. Crozier, Paul S. Crozier, Paul S. Crozier, Paul S. Crozier, Jon Anderson, Joshua A. Anderson, Alex Trav- esset, and Alex Travesset. Molecular dynamics of ionic transport and electrokinetic effects in realistic silica chan- nels.The Journal of Physical Chemistry, 2008

  26. [26]

    Ma and F

    E. Ma and F. M. Geiger. Divalent ion specific outcomes on stern layer structure and total surface potential at the silica:water interface.J. Phys. Chem. A, 125(46):10079– 10088, 2021

  27. [27]

    K. Wang, B. Siboulet, and J.-F. Dufrˆ eche. Collective ion adsorption on silica surfaces driven by ion pairs.J. Phys. Chem. C, 127(45):22315–22335, 2023

  28. [28]

    Shabbir, D

    S. Shabbir, D. Boda, and Z. Hat´ o. Transport properties of cacl 2 are determined by the balance of ion-ion and ion-water interactions as revealed by molecular dynamics simulations.J. Mol. Liq., 426:127308, 2025

  29. [29]

    Martinek, E

    T. Martinek, E. Dubou´ e-Dijon, ˇS. Timr, P. E. Mason, K. Baxov´ a, H. E. Fischer, B. Schmidt, E. Pluhaˇ rov´ a, and P. Jungwirth. Calcium ions in aqueous solutions: Accurate force field description aided by ab initio molec- ular dynamics and neutron scattering.J. Chem. Phys., 148(22):222813, 2018

  30. [30]

    I. V. Leontyev and A. A. Stuchebrukhov. Electronic con- tinuum model for molecular dynamics simulations.J. Chem. Phys., 130(8):085102, 2009

  31. [31]

    Leontyev and A

    I. Leontyev and A. Stuchebrukhov. Accounting for elec- tronic polarization in non-polarizable force fields.Phys. Chem. Chem. Phys., 13(7):2613–2626, 2011

  32. [32]

    Z. R. Kann and J. L. Skinner. A scaled-ionic-charge sim- ulation model that reproduces enhanced and suppressed water diffusion in aqueous salt solutions.J. Chem. Phys., 141(10):104507, 2014

  33. [33]

    Biriukov, O

    D. Biriukov, O. Kroutil, and M. Pˇ redota. Modeling of solid–liquid interfaces using scaled charges: rutile (110) surfaces.Phys. Chem. Chem. Phys., 20(37):23954–23966, 2018

  34. [34]

    Yue and A

    S. Yue and A. Z. Panagiotopoulos. Dynamic properties of aqueous electrolyte solutions from non-polarisable, po- larisable, and scaled-charge models.Mol. Phys., 117(23– 24):3538–3549, 2019

  35. [35]

    Kohagen, P

    M. Kohagen, P. E. Mason, and P. Jungwirth. Accurate description of calcium solvation in concentrated aqueous solutions.J. Phys. Chem. B, 118(28):7902–7909, 2014

  36. [36]

    I. M. Zeron, J. L. F. Abascal, and C. Vega. A force field of Li +, Na +, K +, Mg 2+, Ca 2+, Cl −, and SO 2− 4 in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions.J. Chem. Phys., 151(13):134504, 2019

  37. [37]

    Duboue-Dijon, M

    E. Duboue-Dijon, M. Javanainen, P. Delcroix, P. Jung- wirth, and H. Martinez-Seara. A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization.J. Chem. Phys., 153(5):050901, 2020

  38. [38]

    Pˇ redota and D

    M. Pˇ redota and D. Biriukov. Electronic continuum cor- rection without scaled charges.J. Mol. Liq., 314:113571, 2020

  39. [39]

    Megyes, I

    T. Megyes, I. Bak´ o, Sz. B´ alint, T. Gr´ osz, and T. Rad- nai. Ion pairing in aqueous calcium chloride solution: Molecular dynamics simulation and diffraction studies. J. Mol. Liq., 129(1–2):63–74, 2006

  40. [40]

    good”, the “bad

    D. Biriukov, H.-W. Wang, N. Rampal, C. Tempra, P. Kula, J. C. Neuefeind, A. G. Stack, and M. Pˇ redota. The “good”, the “bad”, and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solu- tions.J. Chem. Phys., 156(19):194505, 2022

  41. [41]

    J. M. Young, C. Tietz, and A. Z. Panagiotopoulos. Ac- tivity coefficients and solubility of CaCl 2 from molecular simulations.J. Chem. Eng. Data, 65(2):337–348, 2019

  42. [42]

    Oostenbrink, A

    C. Oostenbrink, A. Villa, A. E Mark, and W. F. Van Gunsteren. A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6.J. Comp. Chem., 25(13):1656–1676, 2004

  43. [43]

    Mean force potential for the calcium–chloride ion pair in wa- ter

    L. X. Dang and D. E. Smith. Comment on “Mean force potential for the calcium–chloride ion pair in wa- ter”[J. Chem. Phys. 99, 4229 (1993)].J. Chem. Phys., 102(8):3483–3484, 1995

  44. [44]

    Kohagen, P

    M. Kohagen, P. E. Mason, and P. Jungwirth. Accounting for electronic polarization effects in aqueous sodium chlo- ride via molecular dynamics aided by neutron scattering. 16 J. Phys. Chem. B, 120(8):1454–1460, 2015

  45. [45]

    A. D. MacKerell Jr, D. Bashford, M.L.D.R. Bellott, R. L. Dunbrack Jr, J. D. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B, 102(18):3586–3616, 1998

  46. [46]

    Deublein, S

    S. Deublein, S. Reiser, J. Vrabec, and H. Hasse. A set of molecular models for alkaline-earth cations in aqueous solution.J. Phys. Chem. B, 116(18):5448–5457, 2012

  47. [47]

    Mamatkulov, M

    S. Mamatkulov, M. Fyta, and R. R. Netz. Force fields for divalent cations based on single-ion and ion-pair proper- ties.J. Chem. Phys., 138(2):024505, 2013

  48. [48]

    Vazdar, P

    M. Vazdar, P. Jungwirth, and P. E. Mason. Aque- ous guanidinium–carbonate interactions by molecular dy- namics and neutron scattering: Relevance to ion–protein interactions.J. Phys. Chem. B, 117(6):1844–1848, 2013

  49. [49]

    Kohagen, M

    M. Kohagen, M. Lepˇ s´ ık, and P. Jungwirth. Calcium binding to calmodulin by molecular dynamics with effec- tive polarization.J. Phys. Chem. Lett., 5(22):3964–3969, 2014

  50. [50]

    Melcrov´ a, S

    A. Melcrov´ a, S. Pokorna, S. Pullanchery, M. Kohagen, P. Jurkiewicz, M. Hof, P. Jungwirth, P. S. Cremer, and L. Cwiklik. The complex nature of calcium cation inter- actions with phospholipid bilayers.Sci. Rep., 6(1):38035, 2016

  51. [51]

    Magarkar, P

    A. Magarkar, P. Jurkiewicz, C. Allolio, M. Hof, and P. Jungwirth. Increased binding of calcium ions at pos- itively curved phospholipid membranes.J. Phys. Chem. Lett., 8(2):518–523, 2017

  52. [52]

    Biriukov, P

    D. Biriukov, P. Fibich, and M. Pˇ redota. Zeta poten- tial determination from molecular simulations.J. Phys. Chem. C, 124(5):3159–3170, 2020

  53. [53]

    Le Breton and L

    G. Le Breton and L. Joly. Molecular modeling of aque- ous electrolytes at interfaces: Effects of long-range dis- persion forces and of ionic charge rescaling.J. Chem. Phys., 152(24):241102, 2020

  54. [54]

    Kraus, J

    H. Kraus, J. Rybka, A. H¨ oltzel, N. Trebel, U. Tallarek, and N. Hansen. Porems: a software tool for generating silica pore models with user-defined surface functionali- sation and pore dimensions.Mol. Sim., 47(4):306–316, 2021

  55. [55]

    T. S. Gulmen and W. H. Thompson. Testing a two-state model of nanoconfined liquids: Conformational equilib- rium of ethylene glycol in amorphous silica pores.Lang- muir, 22(26):10919–10923, 2006

  56. [56]

    Br´ odka and T

    A. Br´ odka and T. W. Zerda. Properties of liquid acetone in silica pores: Molecular dynamics simulation.J. Chem. Phys., 104(16):6319–6326, 1996

  57. [57]

    L. T. Zhuravlev. Concentration of hydroxyl groups on the surface of amorphous silicas.Langmuir, 3(3):316–318, 1987

  58. [58]

    Luo and J

    R.-S. Luo and J. Jonas. Raman scattering study of liquid ethylene glycol confined to nanoporous silica glasses.J. Raman Spectrosc., 32(11):975–978, 2001

  59. [59]

    Pˇ redota, M

    M. Pˇ redota, M. L. Machesky, and D. J. Wesolowski. Molecular origins of the zeta potential.Langmuir, 32(40):10189–10198, 2016

  60. [60]

    J. Yang, H. Su, C. Lian, Y. Shang, H. Liu, and J. Wu. Un- derstanding surface charge regulation in silica nanopores. Phys. Chem. Chem. Phys., 22:15373–15380, 2020

  61. [61]

    B. Hess, C. Kutzner, D. Van Der Spoel, and E. Lin- dahl. GROMACS 4: algorithms for highly efficient, load- balanced, and scalable molecular simulation.J. Chem. Theor. Comput., 4(3):435–447, 2008

  62. [62]

    B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije. Lincs: A linear constraint solver for molecular simulations.J. Comput. Chem., 18(12):1463–1472, 1997

  63. [63]

    Essmann, L

    U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. A smooth particle mesh Ewald method.J. Chem. Phys., 103(19):8577–8593, 1995

  64. [64]

    Bussi, D

    G. Bussi, D. Donadio, and M. Parrinello. Canonical sampling through velocity rescaling.J. Chem. Phys., 126(1):014101, 2007

  65. [65]

    Parrinello and A

    M. Parrinello and A. Rahman. Polymorphic transitions in single crystals: A new molecular dynamics method.J. Appl. Phys., 52(12):7182–7190, 1981

  66. [66]

    D. Boda, W. R. Fawcett, D. Henderson, and S. Soko lowski. Monte Carlo, density functional theory, and Poisson-Boltzmann theory study of the structure of an electrolyte near an electrode.J. Chem. Phys., 116(16):7170–7176, 2002

  67. [67]

    Sarkadi, D

    Zs. Sarkadi, D. Fertig, Z. Hat´ o, M. Valisk´ o, and D. Boda. From nanotubes to nanoholes: Scaling of selectivity in uniformly charged nanopores through the Dukhin num- ber for 1:1 electrolytes.J. Chem. Phys., 154(15):154704, 2021

  68. [68]

    Sarkadi, D

    Zs. Sarkadi, D. Fertig, M. Valisk´ o, and D. Boda. The Dukhin number as a scaling parameter for selectivity in the infinitely long nanopore limit: extension to multiva- lent electrolytes.J. Mol. Liq., 357:119072, 2022

  69. [69]

    Sarkadi, Z

    Zs. Sarkadi, Z. Hat´ o, M. Valisk´ o, and D. Boda. Scaling for selectivity in finite nanopores for 1:1 electrolytes: the dependence of predictability of device behavior on system parameters.J. Mol. Liq., 387:122571, 2023

  70. [70]

    Almers, E

    W. Almers, E. W. McCleskey, and P. T. Palade. A non- selective cation conductance in frog muscle membrane blocked by micromolar external calcium ions.J. Physiol., 353(1):565–583, 1984

  71. [71]

    D. Boda, M. Valisk´ o, B. Eisenberg, W. Nonner, D. Hen- derson, and D. Gillespie. The effect of protein dielectric coefficient on the ionic selectivity of a calcium channel. J. Chem. Phys., 125(3):034901, 2006

  72. [72]

    Rezaei, A

    M. Rezaei, A. R. Azimian, A. R. Pishevar, and D. J. Bonthuis. Viscous interfacial layer formation causes elec- troosmotic mobility reversal in monovalent electrolytes. Phys. Chem. Chem. Phys., 20(35):22517–22524, 2018

  73. [73]

    Vincze, M

    J. Vincze, M. Valisk´ o, and D. Boda. The nonmonotonic concentration dependence of the mean activity coefficient of electrolytes is a result of a balance between solvation and ion-ion correlations.J. Chem. Phys., 133(15):154507, 2010

  74. [74]

    Valisk´ o and D

    M. Valisk´ o and D. Boda. Resurrection of H¨ uckel’s idea: Decoupling ion–ion and ion–water terms in activity coef- ficients via the state-dependent dielectric constant.Fluid Phase Equilib., 572:113826, 2023

  75. [75]

    Henderson and D

    D. Henderson and D. Boda. Insights from theory and simulation on the electrical double layer.Phys. Chem. Chem. Phys., 11(20):3822–3830, 2009

  76. [76]

    D. Boda, W. Nonner, M. Valisk´ o, D. Henderson, B. Eisenberg, and D. Gillespie. Steric selectivity in Na channels arising from protein polarization and mobile side chains.Biophys. J, 93(6):1960–1980, 2007

  77. [77]

    Gillespie and D

    D. Gillespie and D. Boda. The anomalous mole fraction effect in calcium channels: A measure of preferential se- lectivity.Biophys. J., 95(6):2658–2672, 2008

  78. [78]

    D. Boda, M. Valisk´ o, D. Henderson, B. Eisenberg, D. Gillespie, and W. Nonner. Ion selectivity in L-type 17 calcium channels by electrostatics and hard-core repul- sion.J. Gen. Physiol., 133(5):497–509, 2009

  79. [79]

    Malasics, D

    M. Malasics, D. Boda, M. Valisk´ o, D. Henderson, and D. Gillespie. Simulations of calcium channel block by trivalent ions: Gd3+ competes with permeant ions for the selectivity filter.Biochim. et Biophys. Acta - Biomem- branes, 1798(11):2013–2021, 2010

  80. [80]

    Valisk´ o, B

    M. Valisk´ o, B. Matejczyk, Z. Hat´ o, T. Krist´ of, E. M´ adai, D. Fertig, D. Gillespie, and D. Boda. Multiscale analysis of the effect of surface charge pattern on a nanopore’s rectification and selectivity properties: from all-atom model to poisson-nernst-planck.J. Chem. Phys., 150(14):144703, 2019

Showing first 80 references.