The interplay of cation/anion and monovalent/divalent selectivity in negatively charged nanopores: local charge inversion and anion leakage

This study demonstrates that the anomalous mole fraction effect and anion leakage in negatively charged wide nanopores are governed by a delicate interplay between charge inversion, anion leakage, and ionic mobility, which can be accurately reproduced by matching the distance of closest approach between ions and surface charges regardless of the specific microscopic model used for surface groups.

The Gist

Imagine a tiny, microscopic tunnel bored through a plastic membrane. This tunnel is lined with tiny, negatively charged "hooks" (like Velcro). In the water flowing through this tunnel, there are different types of swimmers: heavy, double-charged swimmers (Calcium, or Ca²⁺) and lighter, single-charged swimmers (Potassium, or K⁺) and their opposites, the negative swimmers (Chloride, or Cl⁻).

Usually, you'd expect the tunnel to be picky, letting only the light swimmers pass or only the heavy ones. But in this specific tunnel, something weird happens called the Anomalous Mole Fraction Effect (AMFE).

The Mystery: The "Traffic Jam" Paradox

Imagine you are mixing two types of cars on a highway: fast sports cars (Calcium) and slow sedans (Potassium).

• If the road is full of only sedans, traffic moves at a steady pace.
• If the road is full of only sports cars, traffic also moves at a steady pace.
• But, if you mix them 20% sports cars and 80% sedans, the traffic grinds to a complete halt.

This is the "anomalous" part. The paper tries to figure out why adding a little bit of the "fast" Calcium actually slows down the whole system more than having just Calcium or just Potassium.

The Detective Work: How the Tunnel is Built

The scientists in this paper are like engineers trying to build a digital model of this tunnel to see what's happening inside. They know the tunnel is lined with "hooks" (carboxyl groups) that grab onto the Calcium swimmers.

They realized that how they drew these hooks in their computer model mattered immensely. They tested two main ways of drawing the hooks:

1. The "Sticky Note" Model: They drew the hooks as fixed, frozen dots stuck right on the wall.
2. The "Wiggly Arm" Model: They drew the hooks as actual atoms (Oxygen) that could wiggle and move slightly, like a person reaching out with a flexible arm.

The Big Discovery: It's All About the "Elbow Room"

The most surprising finding is that it doesn't matter if the hooks are frozen or wiggly, as long as you get one specific measurement right: The Distance of Closest Approach (DCA).

Think of it like a dance floor:

• If the Calcium swimmers can get too close to the hooks (like dancers bumping into each other), they get stuck. They grab onto the hooks so tightly they can't move. They become "frozen" in place.
• When the Calcium is frozen, the Potassium swimmers are pushed away by the electric charge.
• Meanwhile, the negative Chloride swimmers (who usually get pushed out) find a weird loophole. Because the Calcium is stuck, the Chloride swimmers can actually slip through the middle of the tunnel and start carrying the current!

The paper found that if you set the "elbow room" (the distance between the swimmer and the hook) correctly in the model, you get the same result whether you use frozen hooks or wiggly arms. The distance is the master key, not the shape of the hook.

The "Leaky" Tunnel

In the past, scientists thought these tunnels were perfect filters that blocked all negative swimmers (Chloride). But this paper shows that in wide tunnels, the tunnel is actually leaky.

• When the Calcium gets stuck on the walls, the tunnel becomes so crowded with stuck Calcium that the negative Chloride swimmers take over the job of carrying the electricity.
• In fact, in pure Calcium water, the negative swimmers might actually carry more current than the Calcium swimmers themselves!

The Conclusion: Why the Traffic Stops

So, why does the traffic jam happen at that specific mix (20% Calcium)?

1. The Trap: The Calcium swimmers are attracted to the hooks and get stuck, creating a "traffic jam" of stuck Calcium ions.
2. The Blockade: This jam blocks the Potassium swimmers from passing.
3. The Leak: The negative Chloride swimmers try to help, but the jam is so bad that the total flow drops to a minimum.
4. The Recovery: If you add more Calcium, the jam eventually clears out because there are so many Calcium ions that they start pushing each other through, and the flow picks up again.

The Takeaway for Everyone

This research teaches us that when we try to understand how tiny filters work (like in water purification or even in our own cells), we don't need to model every single atom perfectly. We just need to get the distance right between the particles and the wall.

It's like trying to predict how people move through a crowded hallway. You don't need to know if everyone is wearing a hat or a scarf (the microscopic details); you just need to know how much space they have to walk (the distance). If you get the space right, your prediction will be spot on, even if your model is much simpler than reality.

Technical Summary

1. Problem Statement

The Anomalous Mole Fraction Effect (AMFE) is a hallmark phenomenon in ion transport where the total ionic current through a negatively charged pore exhibits a non-monotonic minimum at an intermediate mole fraction of divalent cations (e.g., $Ca^{2+}$) in a mixture with monovalent cations (e.g., $K^+$).

• Context: While AMFE is well-understood in narrow biological ion channels (where anions are excluded and competition is strictly between cations), its microscopic origin in wide synthetic nanopores (radius $\sim$2.7 nm) remains unclear.
• The Challenge: In wide pores, anion leakage ($Cl^-$) is significant and can contribute substantially to the total current. Previous models struggled to reproduce the full experimental AMFE curve and the specific behavior of anion leakage in pure $CaCl_2$ solutions, where recent Molecular Dynamics (MD) simulations and experiments suggest $Cl^-$ current may exceed $Ca^{2+}$ current.
• Key Question: How does the microscopic modeling of surface carboxyl ($COO^-$) groups influence charge inversion, ionic currents, and the resulting AMFE in wide nanopores?

2. Methodology

The authors employed a reduced modeling framework combining the Nernst-Planck (NP) equation with Local Equilibrium Monte Carlo (LEMC) simulations (NP+LEMC).

• System Geometry: A biconical PET nanopore (length 12 $\mu$m, varying radius 2.7–79 nm) was simplified. The simulation explicitly modeled the central narrowest section ($H=12$ nm, $R=2.7$ nm), treating the rest of the pore via adjusted diffusion coefficients to match experimental conductance.
• Implicit Solvent: Water was treated as a continuum with dielectric screening and position-dependent diffusion coefficients, rather than explicit molecules, to allow efficient sampling of low concentrations.
• Surface Charge Modeling: The study systematically varied the representation of the negatively charged $COO^-$ groups on the pore wall to identify controlling parameters:
  1. Fixed Point Charges: Charges placed on a grid with spacing $\Delta x$.
  2. Distance of Closest Approach (DCA): The distance ($r_0$) between the surface charge and the pore wall (or the center of the ion).
  3. Explicit Particles: Modeling surface oxygens as charged hard spheres (either single or double oxygen atoms per group) confined by harmonic potentials, allowing for some vibrational mobility.
• Fitting Strategy: Instead of fitting only to the endpoints ($\eta=0$ and $\eta=1$), the authors fitted the pore diffusion coefficients ($D^p_i$) to three experimental data points:
  1. Pure $KCl$($\eta=0$).
  2. Pure $CaCl_2$ ($\eta=1$).
  3. The minimum of the experimental conductance curve ($\eta=0.2$).
     This allowed for the reproduction of the entire AMFE curve and the specific behavior of anion leakage.

3. Key Contributions

• Identification of Control Parameters: The study identifies two critical parameters governing charge inversion and transport:
  1. Localization ($\Delta x$): The grid spacing of surface charges.
  2. Distance of Closest Approach (DCA, $r_0$): The physical distance between the ionic charge center and the pore surface charge.
• Anion Leakage Mechanism: The paper demonstrates that in wide nanopores, AMFE is not solely driven by cation competition. Instead, it arises from a delicate interplay between charge inversion (strong $Ca^{2+}$ adsorption), anion leakage ($Cl^-$ transport), and ionic mobility.
• Device-Level vs. Microscopic Equivalence: A major finding is that vastly different microscopic models of surface groups (fixed point charges vs. flexible explicit oxygen atoms) yield indistinguishable device-level conductance curves if the DCA is matched. This suggests that accurate representation of the effective ion-surface distance is more critical than the detailed atomistic structure of the surface groups.

4. Key Results

• Reproduction of AMFE and Anion Leakage: By tuning the DCA ($r_0$) to match the radius of oxygen atoms ($r_0 \approx 0.14$ nm), the model successfully reproduced the entire experimental AMFE curve. Crucially, it reproduced the finding that in pure $CaCl_2$, the $Cl^-$ conductance exceeds the $Ca^{2+}$ conductance ($G_{Cl^-} > G_{Ca^{2+}}$), consistent with recent MD simulations and experiments.
• Role of Charge Inversion:
  • Strong charge inversion (small $r_0$ or small $\Delta x$) leads to deep depletion zones between adsorption sites for $Ca^{2+}$, severely limiting $Ca^{2+}$ mobility.
  • As $r_0$ increases (weakening charge inversion), $Ca^{2+}$ mobility improves, but the fitting procedure requires a reduction in the fitted diffusion coefficient ($D^p_{Ca}$) to match experimental conductance.
  • Conversely, $Cl^-$ conductance increases with $r_0$ because the "overcharged" pore wall attracts fewer co-ions, but the fitted diffusion coefficient ($D^p_{Cl}$) must increase significantly to match the experimental total current, reflecting bulk-like behavior where $Cl^-$ is more mobile than $Ca^{2+}$.
• Selectivity Mechanism: In wide pores, monovalent vs. divalent cation selectivity is modulated by cation vs. anion selectivity. The strong adsorption of $Ca^{2+}$ to surface groups creates a local positive charge layer that repels $K^+$, while allowing $Cl^-$ to pass more freely in the central region of the pore.
• Limitations with Small Ions: The model overestimated the competition of small monovalent ions ($Na^+$, $Li^+$) against $Ca^{2+}$ compared to $K^+$. This is attributed to the use of Pauling radii which neglect the larger effective hydrated radii of $Na^+$ and $Li^+$, a factor known to hinder their ability to displace $Ca^{2+}$ from the surface.

5. Significance

• Theoretical Advancement: The work clarifies that AMFE in wide nanopores is a complex phenomenon involving anion leakage and charge inversion, distinct from the single-file cation competition seen in narrow biological channels.
• Modeling Efficiency: It validates the use of reduced models for nanofluidic systems. The study proves that complex atomistic details of surface groups are not strictly necessary to predict macroscopic transport properties, provided the effective ion-surface distance (DCA) is accurately represented.
• Experimental Alignment: The model bridges the gap between experimental conductance data, MD simulations, and theoretical predictions, specifically resolving the paradox of high anion leakage in negatively charged pores.
• Design Guidelines: The findings provide clear guidelines for designing synthetic nanopores and interpreting ion selectivity: controlling the physical distance between ions and surface charges is the primary lever for tuning selectivity and conductance.