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

Using scaled-charge molecular dynamics simulations, this study reveals that while negatively charged silica nanopores exhibit conventional cation selectivity for NaCl, they lose this selectivity for CaCl$_2$ due to calcium ion immobilization and charge inversion, which shifts dominant conduction to chloride ions in the pore interior.

The Gist

Imagine a tiny, microscopic tunnel made of glass (silica) that is so narrow it's only a few atoms wide. The walls of this tunnel are negatively charged, like a magnet with a negative pole. Usually, when you push water with salt through such a tunnel, the negative walls act like a bouncer, letting the positive ions (cations) pass through easily while blocking the negative ones (anions). This is called "cation selectivity."

However, this paper investigates what happens when you change the type of salt. Specifically, the researchers looked at two scenarios:

1. Sodium Chloride (NaCl): The common table salt.
2. Calcium Chloride (CaCl₂): A salt containing calcium, which has a stronger electrical charge (it's "multivalent").

Here is the simple breakdown of what they found, using everyday analogies:

The "Bouncer" vs. The "Sticky Trap"

In the Sodium (NaCl) case, the negative walls act like a standard bouncer. They attract the positive sodium ions, creating a crowd of them right next to the wall. These sodium ions are still free to move around, so they zip through the tunnel easily. The tunnel works as expected: it lets the positive ions through and blocks the negative ones.

In the Calcium (CaCl₂) case, things get weird. Calcium ions are like "super-sticky" magnets. When they hit the negative wall, they don't just hang out nearby; they stick to the wall so tightly that they become frozen in place.

• The Analogy: Imagine a hallway where the walls are covered in super-strong Velcro. If you throw a regular ball (Sodium) at the wall, it bounces off or rolls along. But if you throw a heavy, sticky ball (Calcium), it slams into the wall and sticks there, unable to move.

The "Traffic Jam" and the "Middle Lane"

Because the Calcium ions are stuck to the walls, they stop contributing to the flow of electricity. They are available, but they aren't moving.

• The Result: The layer of water right next to the wall (the "surface layer") stops conducting electricity effectively because the ions are immobilized.
• The Twist: Since the Calcium ions are stuck to the negative wall, they actually overcompensate for the wall's negative charge. They make the wall effectively positive.
• The Consequence: Now that the wall acts positive, it repels the negative Chloride ions, pushing them away from the wall and into the center of the tunnel.

So, the flow of electricity in the Calcium solution doesn't happen near the walls (where the ions are stuck); it happens in the middle of the tunnel. In this middle section, the negative Chloride ions are actually moving faster than the Calcium ions. This causes the tunnel to lose its "positive-ion-only" rule and start behaving more like a normal, open pipe where both types of ions can pass, or even favoring the negative ones slightly.

The "Driver" of the Story: Force Fields

The researchers used computer simulations to watch this happen. They had to be very careful about the "rules" they programmed into the computer (called "force fields").

• The Metaphor: Think of the force field as the rulebook for how the atoms interact. If the rulebook says Calcium is too sticky, the simulation shows the ions getting stuck forever. If the rulebook says they are too slippery, they don't stick enough.
• The Finding: The researchers found that the general story (Calcium sticks, Chloride moves to the middle, selectivity is lost) is true no matter which rulebook they used. However, the exact details (how fast they move, exactly how much current flows) changed significantly depending on which rulebook they picked. This means that while we understand the big picture, getting the numbers right requires very precise modeling.

The "Water Flow" Surprise

The study also looked at the water itself. When ions move, they drag water molecules along with them (like a crowd of people moving through a hallway, bumping into the air).

• The Finding: Because the Calcium ions are stuck and the Chloride ions are moving in the middle, the water flow is a messy mix. Sometimes the water flows one way, sometimes the other, depending on exactly which "rulebook" was used in the simulation. It's a delicate balance where a tiny change in the rules can flip the direction of the water flow.

Summary

In short, this paper explains why a negatively charged nanopore acts like a one-way gate for simple salt (Sodium) but acts like a confused, mixed-traffic zone for calcium salt.

• Sodium: Stays mobile near the walls; the tunnel selects for positive ions.
• Calcium: Gets stuck to the walls; the tunnel loses its selectivity because the "traffic" moves in the middle of the pipe instead of the walls.

The researchers emphasize that while this mechanism is robust, the exact numbers depend heavily on how accurately we model the interactions between the ions, the water, and the glass walls.

Technical Summary

1. Problem Statement

Ion transport through negatively charged nanopores is traditionally expected to exhibit cation selectivity due to the formation of an Electrical Double Layer (EDL) where counter-ions (cations) accumulate near the wall. However, experimental and simulation studies have shown that this behavior breaks down in the presence of multivalent ions (e.g., $Ca^{2+}$).

• The Paradox: In multivalent electrolytes like $CaCl_2$, strong ion-surface correlations can lead to charge inversion (where the surface effectively becomes positive due to over-adsorption of $Ca^{2+}$), potentially allowing anions to conduct or reversing selectivity.
• The Gap: While static adsorption (structure) is well-understood, the quantitative link between static interfacial structure and dynamic perm-selectivity (actual ion transport) remains unclear. Furthermore, it is uncertain how sensitive these transport properties are to the specific Force Field (FF) choices (specifically ion charge scaling and water models) used in molecular dynamics (MD) simulations.

2. Methodology

The authors employed atomistic Molecular Dynamics (MD) simulations to investigate $NaCl$ and $CaCl_2$ transport through negatively charged silica nanopores.

• System Setup:

  • Geometry: Cylindrical silica nanopores with a radius ($R_P$) ranging from $\approx 1$ nm to $3$ nm and a length of $\approx 7$ nm.
  • Surface: The pore walls contain deprotonated silanol groups (OS) creating a negative surface charge density ($\sigma \approx -0.85$ to $-1.54$ $e/nm^2$).
  • Electrolyte: 1 M $NaCl$ or $CaCl_2$ solutions.
  • Driving Force: An external electric field ($E = 0.066$ V/nm) applied along the pore axis to induce ionic current.

• Force Fields (FF) & Models:

  • Ions: The study compares Full-Charge (FULL) models against Scaled-Charge (ECC) models. In ECC models, ion charges are scaled by $1/\sqrt{\epsilon_\infty}$ (approx. 0.75) to account for electronic polarization and charge transfer.
  • Specific Models Tested:
    • $NaCl$: TIP4P/2005 water.
    • $CaCl_2$: Various models including FULL, ECC, ECCR2 (scaled charge + reduced diameter), and ECCR (scaled charge + further reduced diameter).
  • Water Models: TIP4P/2005 and SPC/E were compared to assess hydration shell effects.
  • Surface Charges: Both full-charge ($-0.74e$) and scaled-charge ($-0.555e$) silanol oxygen models were used to ensure consistency with ion models.

• Analysis Framework:

  • The core analytical approach decomposes the radial particle current density ($j_i(r)$) into two components:
    $$j_i(r) = c_i(r) \cdot v_i(r)$$
    Where $c_i(r)$ is the concentration (availability of charge carriers) and $v_i(r)$ is the local velocity (mobility under the electric field).
  • Selectivity Metrics:
    • Pore Selectivity ($S^P_+$): Ratio of cation to anion charge currents in the pore.
    • Bulk Selectivity ($S^B_+$): Ratio in bulk solution (used as a baseline).
    • Reduced Selectivity ($S^*_+$): $S^P_+ / S^B_+$, isolating the pore's effect from intrinsic bulk mobility differences.

3. Key Contributions

1. Decoupling Structure and Dynamics: The paper explicitly separates the contributions of ion availability (concentration profile) and mobility (velocity profile) to total current, providing a mechanistic explanation for selectivity loss.
2. Mechanism of Selectivity Loss in $CaCl_2$: It identifies that the loss of cation selectivity in $CaCl_2$ is not due to a lack of cations near the wall, but rather their immobilization.
3. Force Field Sensitivity Analysis: It systematically evaluates how different FF parameterizations (charge scaling, ion size, water models) affect electrokinetic predictions, highlighting that while qualitative mechanisms are robust, quantitative results (especially Electroosmotic Flow) are highly sensitive.

4. Key Results

A. $NaCl$ vs. $CaCl_2$ Transport Mechanisms

• $NaCl$ (Monovalent): Exhibits classical cation selectivity. $Na^+$ ions form a diffuse layer near the wall with high availability and moderate mobility. The surface region contributes significantly to the total current.
• $CaCl_2$ (Multivalent):
  • Surface Layer: $Ca^{2+}$ ions adsorb strongly to the silanol groups, leading to charge inversion (the surface layer becomes effectively positive). However, these adsorbed $Ca^{2+}$ ions are immobilized (residence time is long, hopping between sites is rare). Consequently, their contribution to surface conduction is negligible.
  • Anion Behavior: Due to charge inversion, the surface layer repels $Cl^-$, making them unavailable in the surface region.
  • Dominant Mechanism: The total current is dominated by the bulk-like region in the center of the pore. Here, $Cl^-$ ions have higher mobility than $Ca^{2+}$ (similar to bulk behavior).
  • Result: The pore exhibits bulk-like or even slight anion selectivity ($S^*_+ \approx 0.6 - 0.9$), contrasting sharply with the cation selectivity seen in $NaCl$($S^*_+ \approx 2.4 - 2.9$).

B. Impact of Force Field Parameters

• Ion Models: The ECCR2 model (scaled charge + reduced diameter) provided the best agreement with experimental diffusion and conductivity data.
  • FULL models showed "practically non-ergodic" behavior (ions stuck to the surface), leading to unphysically slow dynamics.
  • ECCR (smallest diameter) showed stronger binding to the surface and reduced $Cl^-$ mobility, slightly altering selectivity trends.
• Water Models:
  • SPC/E water is "stickier" (stronger ion-water interaction) than TIP4P/2005. This hinders $Ca^{2+}$ binding to the surface, keeping $Ca^{2+}$ more mobile in the near-wall region, which slightly increases cation selectivity compared to TIP4P/2005.
• Pore Radius: As the pore radius increases ($>1.3$ nm), the bulk-like central region dominates, and selectivity converges to bulk-like behavior. In very narrow pores ($\approx 1$ nm), double layers overlap, and selectivity increases slightly due to confinement effects.

C. Electroosmotic Flow (EOF)

• EOF is highly sensitive to model choices.
• In most simulations, water flows in the direction of $Ca^{2+}$ migration (positive EOF).
• However, for specific combinations (Full-charge OS + TIP4P/2005), the flow reverses sign (negative EOF), consistent with experimental observations of charge inversion.
• The direction of EOF depends on a delicate balance between the dominance of $Ca^{2+}$ vs. $Cl^-$ in different radial regions and their coupling to water molecules.

5. Significance and Conclusion

• Qualitative Robustness: The fundamental mechanism—that strong adsorption immobilizes multivalent cations, shifting transport dominance to the bulk-like core where anions may dominate—is robust across different force fields.
• Quantitative Sensitivity: Quantitative predictions of selectivity and EOF direction are highly sensitive to the balance of Ion-Ion, Ion-Water, and Ion-Surface interactions encoded in the force field.
• Implication for Modeling: The study validates the use of scaled-charge models (like ECCR2) for simulating confined electrolytes, as they correct the unphysical slow dynamics seen in full-charge models.
• Broader Impact: The findings suggest that simplified models (implicit water, hard walls) can capture essential physics for device-level behavior, but careful validation of FF parameters is critical for predicting specific electrokinetic phenomena like EOF reversal in nanofluidic devices.

In summary, the paper demonstrates that the "loss" of cation selectivity in $CaCl_2$ nanopores is a dynamic phenomenon driven by the immobilization of adsorbed cations, rendering the surface layer non-conductive and allowing the bulk-like interior (often anion-selective) to dictate the overall transport behavior.