Solvation of Ions and Polar Solutes

Ion solvation is not a historical curiosity confined to physical chemistry textbooks — it is operational machinery inside every living cell and every modern drug and material built to interact with aqueous ions. The selectivity of K⁺ channels over Na⁺, first explained atomically by Roderick MacKinnon's KcsA structure, depends on the selectivity filter's backbone carbonyl oxygens presenting a cage of dipoles geometrically matched to dehydrated K⁺ (ionic radius 1.38 Å) but mismatched to smaller Na⁺ (1.02 Å), which would rattle in an oversized cage and fail to be fully stabilized — precisely the ion-dipole geometry problem this article develops. Understanding this quantitatively is central to modern electrophysiology, to the design of antiarrhythmic and anesthetic drugs that target specific ion channels, and to engineering biomimetic nanopores for water desalination and DNA sequencing. In pharmacology, the aqueous solvation free energy of a candidate drug molecule — how strongly it is stabilized in water versus how strongly it would be stabilized in a lipid bilayer — directly predicts its membrane permeability and oral bioavailability, and is estimated computationally using descendants of the Born model (generalized Born, Poisson-Boltzmann continuum solvation) as a fast approximation before committing to expensive explicit-solvent molecular dynamics. In battery and materials chemistry, the solvation shell structure of Li⁺ in organic carbonate electrolytes governs the rate and reversibility of charge-discharge cycles in lithium-ion batteries, and engineers deliberately tune solvent dielectric properties and additive molecules to control how tightly Li⁺ is held versus how easily it desolvates at the electrode surface — a direct industrial descendant of the Born-equation logic. In geochemistry and environmental science, the solvation energetics of Mg²⁺, Ca²⁺, and heavy metal ions determine mineral solubility, groundwater ion transport, and the bioavailability of toxic metals. And in structural biology, the observation that Mg²⁺ almost always appears in enzyme active sites still partially hydrated — coordinating substrate phosphate oxygens with only some of its six coordination positions, the rest still occupied by water — is a direct, visible consequence of the enormous hydration free energy computed in Example 2, which no enzyme can afford to strip away for free.

Ask the tutor: How does the size of an ion's hydrated radius compare to its bare crystallographic radius, and why does this reverse the usual periodic trend for ionic mobility in solution?

Water is uniquely suited to solvation because it is small, has a large permanent dipole moment (1.85 D in the gas phase, effectively larger in the liquid due to cooperative polarization), and carries both a hydrogen-bond donor and acceptor capacity on the same molecule. When an ion enters water, the solvent does not merely "surround" it passively; individual water molecules reorient so that their electrostatically favorable pole faces the ion. For a cation such as Na⁺ or K⁺, the oxygen atom — bearing the partial negative charge and two lone pairs — turns to face the ion directly, so the first hydration shell is a cage of oxygen atoms pointing lone pairs inward. For an anion such as Cl⁻ or the carboxylate of glutamate, the hydrogens face inward instead, each contributing an O–H···X⁻ hydrogen bond directed at the negative center. This reorientation is an ion-dipole interaction, and its strength falls off much more gently with distance than the dipole-dipole interactions between two neutral water molecules. A rough electrostatic estimate for the interaction energy between a point charge $q$ and a dipole $\mu$ at distance $r$, oriented for maximum stabilization, is

which decays as $1/r^2$, compared to the $1/r^3$ decay of a dipole-dipole interaction and the much gentler $1/r^6$ decay of a van der Waals dispersion interaction. This is why ion-dipole forces dominate the local structure of water for several ångströms around any charged species, and why hydration shells are physically real, energetically deep, and kinetically long-lived compared to the fleeting hydrogen-bond network of bulk water. The first coordination shell is the most tightly bound and most extensively characterized. For small, highly charged cations, this shell is essentially locked in place: water molecules exchange with the bulk on a timescale of nanoseconds to microseconds for ions like Mg²⁺, versus picoseconds for bulk water hydrogen bonds. Beyond the first shell lies a second, looser shell, oriented by the residual field of the ion and by hydrogen bonding to the first shell, and beyond that the solvent gradually relaxes into the normal, fluctuating hydrogen-bond network of bulk water. This shell structure is why the "effective size" of an ion in solution — its hydrodynamic radius — can be substantially larger than its crystallographic (bare) ionic radius, especially for small, highly charged ions that drag a heavy hydration shell with them as they move.

The central quantitative tool for estimating ionic solvation energy is the Born equation, derived from a purely electrostatic, continuum-dielectric argument. It says that the free energy of transferring a charged sphere of radius $r_i$ and charge $z_ie$ from vacuum into a medium of dielectric constant $\varepsilon_r$ is

where:

• $z_i$ is the ionic charge number (with sign; the equation depends only on $z_i^2$, so cations and anions of equal charge and radius solvate equally well in this model),
• $e$ is the elementary charge, $\varepsilon_0$ the vacuum permittivity,
• $r_i$ is the ionic (Born) radius, essentially the crystallographic radius plus a small correction,
• $\varepsilon_r$ is the solvent's relative dielectric constant (about 78 for water at 25 °C),
• $N_A$ converts the single-ion energy into a molar quantity.

Reading this equation: solvation free energy is always negative (favorable) and grows in magnitude as the square of the charge and inversely with the ionic radius — small, highly charged ions are solvated far more strongly than large, weakly charged ones. The factor $(1 - 1/\varepsilon_r)$ approaches 1 for water (since $\varepsilon_r \approx 78 \gg 1$), meaning nearly the full theoretical maximum of solvation stabilization is realized in water; the same ion in a low-dielectric solvent like hexane ($\varepsilon_r \approx 2$) would be solvated far more weakly. A second essential relationship links the overall enthalpy of dissolving an ionic solid to two much larger, partially cancelling quantities:

where:

• $\Delta H_{\text{lattice}}$ is the lattice energy of the solid, defined as the energy released when gaseous ions come together to form the crystal (a positive number by convention here, so breaking the lattice costs $+\Delta H_{\text{lattice}}$),
• $\Delta H_{\text{hydration}}$ is the (negative) enthalpy released when the gaseous ions are surrounded by water,
• $\Delta H_{\text{soln}}$ is the net, experimentally measurable heat of solution — often only a few kJ/mol despite the individual terms being hundreds of kJ/mol.

Reading this: dissolution is a Born-Haber-type cycle. Two enormous, opposing numbers determine a small residual, so any solvation theory must get both terms right to a few percent to predict even the sign of $\Delta H_{\text{soln}}$ — this is the central numerical challenge of solvation thermodynamics and the reason the Born model, though crude, remains pedagogically central: it gives the correct trend even when it misses the absolute value by 10–20%. Finally, the overall spontaneity of dissolution is governed by the Gibbs free energy,

where $\Delta S_{\text{soln}}$ combines a large, usually negative contribution from ordering water molecules into hydration shells (loss of solvent entropy) with a positive contribution from the translational entropy of the dispersed ions (mixing entropy, always favorable and roughly $R\ln(1/x)$ for dilute solutions). Whether a salt dissolves therefore depends on a genuine three-way competition between lattice energy, hydration enthalpy, and the entropy cost of ordering solvent — which is why some salts (like NaCl) dissolve readily and endothermically (cooling the solution) while others (like Al₂O₃) essentially never dissolve in water at room temperature.

Example 1 — Foundational: Born hydration free energy of Na⁺. Estimate the Born hydration free energy of Na⁺ using an effective ionic radius of $r_i = 1.02\ \text{Å} = 1.02\times10^{-10}$ m (the Pauling crystallographic radius for six-coordinate Na⁺), $z = 1$, and $\varepsilon_r = 78.4$ for water at 25 °C. Constants: $e = 1.602\times10^{-19}$ C, $\varepsilon_0 = 8.854\times10^{-12}$ F/m, $N_A = 6.022\times10^{23}\ \text{mol}^{-1}$. First compute the prefactor:

Numerator: $(6.022\times10^{23})(2.566\times10^{-38}) = 1.546\times10^{-14}$. Denominator: $8\pi(8.854\times10^{-12}) = 2.225\times10^{-10}$.

Divide by $r_i$:

Multiply by $(1 - 1/78.4) = 0.9872$ and apply the negative sign:

This compares reasonably with experimental single-ion hydration free energy estimates for Na⁺ (roughly $-375$ to $-410$ kJ/mol on the most common extrathermodynamic conventions), overestimating in magnitude by nearly a factor of 1.6–1.8 — the well-known systematic deviation of the naive Born model, mostly because the true "Born radius" of an ion in solution is larger than its bare crystallographic radius (it should include roughly half a water-molecule's radius to account for the position of the water dipole's effective center, and the dielectric near the ion is not the bulk value of 78.4 due to dielectric saturation). The trend, however — highly negative, dominated by $z^2/r$ — is exactly right, and is the reason the Born model remains the standard qualitative teaching tool despite its known ~40–60% overestimate in magnitude for small monovalent ions.

Example 2 — Applied: comparing hydration strength of Mg²⁺ and Ca²�sup>+ using charge density. Mg²⁺ has an ionic radius of about 0.72 Å (six-coordinate) and Ca²⁺ about 1.00 Å. Using the Born model's scaling $\Delta G_{\text{Born}} \propto z^2/r$, compare their relative hydration free energies.

Mg²⁺ is predicted to be hydrated about 39% more strongly than Ca²⁺, purely from the smaller ionic radius at equal charge. Using the Example 1 prefactor scaled to $z=2$ (multiply the monovalent prefactor by $z^2=4$) and $r_{Mg} = 0.72\times10^{-10}$ m:

The experimental hydration free energy of Mg²⁺ is around $-1900$ to $-2000$ kJ/mol on common conventions — again the Born model overestimates by roughly a factor of two, but correctly predicts that divalent, small-radius Mg²⁺ is hydrated several times more strongly than monovalent Na⁺ (compare to the $-672$ kJ/mol computed above), which is precisely why Mg²⁺ is never found "naked" in biological systems: it moves through the cell essentially always as a hydrated or partially hydrated species, and enzymes that use Mg²⁺ as a cofactor (kinases, polymerases, many phosphatases) must actively strip or partially strip this shell, an energetically expensive step that is often coupled to substrate binding to make it favorable.

Example 3 — Exam/competition level: heat of solution of a hypothetical salt from lattice and hydration data, testing the near-cancellation. A hypothetical 1:1 salt $MX$ has lattice energy $\Delta H_{\text{lattice}} = +742\ \text{kJ/mol}$. Individual (single-ion, extrathermodynamic) hydration enthalpies are $\Delta H_{\text{hyd}}(M^+) = -418\ \text{kJ/mol}$ and $\Delta H_{\text{hyd}}(X^-) = -335\ \text{kJ/mol}$. (a) Compute $\Delta H_{\text{soln}}$. (b) If the entropy of solution is $\Delta S_{\text{soln}} = +45\ \text{J/(mol}\,\cdot\,\text{K)}$, find $\Delta G_{\text{soln}}$ at 298 K and determine whether the salt dissolves spontaneously. (c) Estimate the solubility product order of magnitude implied, roughly, by this $\Delta G$.

(a) Using the Hess's law relation from the second derivation:

Note the near-total cancellation: two terms near 750 kJ/mol summing to leave only $-11$ kJ/mol, illustrating exactly the fragility discussed in the second derivation — a 1–2% error in either lattice or hydration energy could flip the sign of $\Delta H_{\text{soln}}$ entirely.

(b) Compute $\Delta G_{\text{soln}}$:

Since $\Delta G_{\text{soln}} < 0$, dissolution is spontaneous under standard conditions — and notice that here the entropy term (favorable mixing of ions into solution, plus possibly some release of ordered water if this salt disrupts more solvent structure than it creates) is actually larger in magnitude than the enthalpy term, driving the process forward even though $\Delta H_{\text{soln}}$ alone is only mildly favorable.

(c) Relating $\Delta G_{\text{soln}}^\circ$ to an equilibrium constant (here treated loosely as a solubility-type equilibrium constant $K_{sp}$ for illustration, using $\Delta G^\circ = -RT\ln K$):

A solubility product this large (≫1) is consistent with a highly soluble salt, in line with the negative $\Delta G_{\text{soln}}$ found — this is the qualitative regime that real, highly soluble alkali halides like NaCl or NaBr actually occupy, in contrast to sparingly soluble salts like AgCl, whose $K_{sp} \sim 10^{-10}$ corresponds to a positive $\Delta G_{\text{soln}}$ of comparable magnitude (about $+55$ kJ/mol), reflecting a lattice energy that is not adequately compensated by hydration for that particular ion pair.

Max Born (1882–1970) is remembered overwhelmingly for his foundational contributions to quantum mechanics — the probabilistic (Born) interpretation of the wavefunction, for which he received the 1954 Nobel Prize in Physics — but his 1920 paper on the energetics of ion solvation, developed with collaborators while at Göttingen, gave chemistry one of its most durable order-of-magnitude tools. Born approached the solvation problem exactly as a physicist trained in classical electrostatics would: rather than trying to model the messy molecular detail of hydrogen-bonded water, he asked what a continuum dielectric theory alone could predict, and discovered that a startlingly simple electrostatic charging argument reproduced the correct qualitative trends in ionic solvation across the periodic table. The Born model's genius is not its precision — it is known to overestimate hydration free energies of small ions by 30–60% because it ignores the discrete, saturating, and partially ordered structure of the first few solvent shells — but its scope: with only two physical parameters (an ionic radius and the solvent's bulk dielectric constant), it explains why divalent ions are solvated more strongly than monovalent ones, why solvation weakens as ionic radius increases down a periodic group, and why polar solvents dissolve salts that nonpolar solvents cannot touch. Refinements by Bernal and Fowler, Latimer, Pitzer and Slansky, and later generations of computational chemists added discrete first-shell water molecules, adjusted "effective" Born radii, and ultimately full molecular dynamics simulations — but the qualitative scaffolding, $z^2/r$ times a dielectric factor, remains the first equation every chemistry and biochemistry student reaches for when asked to rank the solvation strength of two ions.