Viscosity as a Smoking Gun for Complex Formation in Solution: Fe$^{2+}$ and Mg$^{2+}$ Chlorides as Examples

This paper demonstrates that viscosity serves as a reliable indicator for determining the extent of complex formation in concentrated electrolyte solutions, using comparative simulations and experiments on FeCl$_2$ and MgCl$_2$ to show that higher viscosity correlates with greater speciation, a finding validated by recent X-ray and neutron scattering data.

The Gist

Imagine you are stirring a pot of soup. If you add a little salt, the water flows easily. But if you add a huge amount of salt, the soup gets thick and sluggish, like honey. Scientists have long known that high concentrations of salt (electrolytes) make water viscous, but they've been arguing for decades about why it happens, especially with certain metals like Iron (Fe) and Magnesium (Mg).

The big mystery is speciation: Do the salt ions just float around freely, or do they stick together to form little "clumps" or "complexes"?

The Problem: A Heated Argument

Scientists have been trying to count these clumps using expensive microscopes and chemical tests. The results? A mess. Some say Iron forms big clumps; others say it stays lonely. It's like trying to count the number of people holding hands in a crowded, dark room by listening to the noise. Everyone hears something different.

The New Idea: The "Traffic Jam" Analogy

This paper proposes a clever new way to solve the mystery. Instead of trying to see the clumps directly, they decided to look at how fast the soup flows (viscosity).

Think of the water molecules as cars on a highway and the salt ions as giant trucks.

• Free Ions: When a truck (ion) is alone, it drags a huge entourage of water cars behind it, creating a massive traffic jam. The water moves slowly.
• Complexes (Clumps): When two trucks stick together to form a "double-truck" (a complex), they actually drag fewer water cars with them. The traffic jam clears up, and the water flows faster (lower viscosity).

The Twist: Usually, we think adding more stuff makes things thicker. But here, forming clumps actually makes the solution thinner (less viscous) because the clumps are more efficient at moving through the water than the lonely, dragging ions.

The Experiment: Iron vs. Magnesium

The researchers compared two salts: Iron Chloride (FeCl₂) and Magnesium Chloride (MgCl₂).

• At low concentrations: Both act the same. The soup flows similarly.
• At high concentrations: They behave very differently. The Iron solution flows much faster (is less thick) than the Magnesium solution.

The Deduction:
Since the Iron solution flows faster, the "Traffic Jam" theory says the Iron ions must be forming more clumps (complexes) than the Magnesium ions. If they were all free, the Iron solution would be just as thick as the Magnesium one. The fact that it's thinner is the "smoking gun" that Iron is busy forming partnerships.

How They Proved It (The Simulation)

Since they couldn't wait for the ions to naturally form clumps in a computer simulation (it would take too long, like waiting for a snowflake to form in a blizzard), they played a game of "What If."

1. They built a virtual world of water and ions.
2. They manually "glued" some Iron and Magnesium ions together to force them into clumps.
3. They watched how the "soup" flowed.

The Result:

• When they forced the ions to clump, the viscosity dropped.
• To match the real-world experiment where Iron flowed faster than Magnesium, they had to assume Iron formed significantly more clumps than Magnesium.

The Takeaway

This paper is a detective story. Instead of trying to find the suspect (the clumps) directly, the detectives looked at the footprints (the viscosity).

• Old Way: "I think Iron forms clumps because my microscope says so." (But the microscope was blurry).
• New Way: "Iron flows faster than Magnesium. Therefore, Iron must be forming more clumps to reduce the drag."

They concluded that viscosity is a reliable detective. By measuring how thick a solution is, we can infer how many molecular "hand-holding" partnerships are happening inside, even if we can't see them directly. This helps us understand everything from how batteries work to how our cells function, without needing to guess in the dark.

Technical Summary

Here is a detailed technical summary of the paper "Viscosity as a Smoking Gun for Complex Formation in Solution: Fe2+ and Mg2+ Chlorides as Examples."

1. Problem Statement

Concentrated electrolyte solutions are ubiquitous in nature and industry (e.g., seawater, biological cells, energy storage), yet their theoretical understanding remains limited, particularly regarding speciation (the formation of ion pairs and complexes).

• The Challenge: Determining the extent of complex formation (e.g., $M^{2+} + Cl^- \rightleftharpoons MCl^+$ or $MCl_2^0$) in concentrated solutions is difficult. Experimental estimates of equilibrium constants ($K$) for species like $FeCl^+$ vary wildly (log $K$ from 0.74 to -0.89), leading to contradictory conclusions about dominant species.
• The Limitation of Current Methods:
  • Experiments: Different techniques (spectrophotometry, X-ray absorption, THz/FIR) yield conflicting results.
  • Theory/Models: Chemical equilibrium models (e.g., SIT, Davies equation) rely on activity coefficients and association constants that are often inaccurate at high concentrations.
  • Simulations: Standard Molecular Dynamics (MD) with empirical force fields often fail to spontaneously form long-lived complexes within feasible simulation timescales (microseconds are required for equilibrium, but simulations are typically nanoseconds). Furthermore, simple point-charge models struggle to capture the specific interaction dynamics of transition metals like $Fe^{2+}$.

2. Methodology

The authors propose a novel approach: using viscosity as a proxy to infer the extent of complexation, rather than trying to calculate complex populations directly.

A. Experimental Work

• System: Aqueous solutions of Ferrous Chloride ($FeCl_2$) and Magnesium Chloride ($MgCl_2$).
• Conditions: Concentrations up to 4.1 mol/kg at 298.15 K.
• Techniques:
  • Density: Measured using a vibrating tube densimeter (Anton Paar DMA5000).
  • Viscosity: Measured using a falling ball viscometer (AMVn).
  • Precautions: Experiments were conducted in a $N_2$ atmosphere to prevent $Fe^{2+}$ oxidation.

B. Computational Work (Molecular Dynamics)

• Force Field: Extended the Madrid-2019 force field (based on TIP4P/2005 water and scaled ion charges) to model $Fe^{2+}$.
  • Justification: $Fe^{2+}$ and $Mg^{2+}$ have similar hydration free energies and cation-oxygen distances. The authors used the same Lennard-Jones parameters for $Fe^{2+}$ as for $Mg^{2+}$ but adjusted the charge scaling ($q=0.85e$ and $q=0.80e$).
• The "Seeding" Strategy: Since complexes do not form spontaneously in standard MD runs within accessible timescales, the authors introduced a static complex population strategy:
  • They artificially "froze" the distance between cations and chloride ions to create monomers ($MCl^+$) and dimers ($MCl_2^0$).
  • They varied the fraction of free ions ($\alpha$), monomers ($\alpha'$), and dimers ($\alpha''$) in the simulation box.
  • This bypasses the need for microsecond-scale simulations to reach equilibrium, allowing the study of how specific complex distributions affect transport properties.

3. Key Contributions

1. Viscosity as a Speciation Probe: The paper establishes that viscosity is highly sensitive to the presence of complexes. It demonstrates that viscosity can be used inversely: if the viscosity is known, one can infer the degree of complexation.
2. Force Field Extension: Successfully extended the Madrid-2019 force field to a transition metal ($Fe^{2+}$) by leveraging the similarity with alkaline earth metals ($Mg^{2+}$), validating this approach through density predictions.
3. Resolution of Contradictions: The study reconciles conflicting literature on $FeCl_2$ speciation by showing that high complexation is required to match experimental viscosity data, supporting recent X-ray absorption studies over older spectrophotometric claims.
4. Toy Model for Viscosity: Developed a heuristic linear model relating viscosity to the percentage of free ions, monomers, and dimers, allowing for qualitative prediction of speciation based on transport data.

4. Key Results

A. Experimental Findings

• Low Concentration (< 2 m): $FeCl_2$ and $MgCl_2$ exhibit very similar viscosities. The Jones-Dole $B$ coefficients (representing ion-water interactions) are nearly identical ($0.405$vs$0.375$ kg/mol).
• High Concentration (> 2 m): A stark divergence occurs. $FeCl_2$ viscosity is significantly lower than $MgCl_2$ viscosity at 4 m (4.17 mPa·s for $FeCl_2$ vs 4.73 mPa·s for $MgCl_2$).
• Interpretation: Lower viscosity in $FeCl_2$ suggests a reduction in the number of "encumbered" water molecules, implying a higher degree of ion pairing/complexation in $FeCl_2$ compared to $MgCl_2$.

B. Simulation Findings

• Without Complexes: Simulations assuming only free ions overestimate the viscosity for both salts at high concentrations and fail to reproduce the experimental difference between $FeCl_2$ and $MgCl_2$.
• With Complexes:
  • Introducing monomers and dimers significantly reduces the simulated viscosity.
  • Mechanism: Complexes have a lower net charge than free ions, exerting a weaker electric field on surrounding water. Additionally, water molecules in the solvation shell of a complex diffuse faster than those around a free ion.
  • Quantitative Match:
    • Using the $q=0.80e$ model, the experimental viscosity of $MgCl_2$ is matched with ~40% free ions.
    • The experimental viscosity of $FeCl_2$ is matched only when assuming 0% free ions (100% complexed, primarily monomers).
  • Conclusion: $FeCl_2$ exhibits a significantly higher degree of complexation (35–40% more) than $MgCl_2$ at 4 m.

C. The "Toy Model"

The authors proposed a linear equation: $\eta = \eta_0 - \delta'\alpha' - \delta''\alpha''$.

• By inputting equilibrium constants ($K_1, K_2$) from NIST and literature, and solving for the complex distribution that fits the experimental viscosity, they confirmed that a high population of monomers and dimers in $FeCl_2$ is consistent with the data.
• This model predicts a viscosity of 4.46 mPa·s for $FeCl_2$ (close to the experimental 4.17 mPa·s) and 4.96 mPa·s for $MgCl_2$ (close to 4.73 mPa·s).

5. Significance

• Methodological Shift: The paper flips the traditional problem-solving approach. Instead of using theory to predict speciation (which is error-prone), it uses a macroscopic transport property (viscosity) to constrain and infer the microscopic speciation.
• Validation of Recent Experiments: The results strongly support recent X-ray absorption spectroscopy findings that $FeCl_2$ is dominated by neutral dichloro complexes and monochloro complexes at high concentrations, contradicting older theories that suggested minimal complexation.
• Practical Application: This approach offers a pathway to characterize complex speciation in other concentrated electrolyte systems where direct structural measurement is difficult or where experimental data is contradictory.
• Future Directions: The authors suggest that combining this viscosity-based inference with Neutron Diffraction with Isotopic Substitution (NDIS) could provide a definitive quantitative map of complex populations in concentrated solutions.

In summary, the paper demonstrates that viscosity is a "smoking gun" for complex formation. The significant drop in viscosity for concentrated $FeCl_2$ compared to $MgCl_2$ is direct evidence of extensive complexation in the former, a conclusion that resolves long-standing discrepancies in the literature.