Cuando se coloca un cristal en un disolvente, las moléculas/iones del cristal se disocian y son disueltos por el disolvente. El proceso se detiene cuando hay un equilibrio entre el cristal y la solución de moléculas/iones más el disolvente. En este momento la solución está saturada, pero sólo si todavía hay algún sólido en contacto con la solución. En este punto la energía libre es mínima.
Algunas moléculas/iones siempre se disociarán del cristal debido a la ganancia de entropía al diluirse o mezclarse, con el disolvente. La solubilidad depende entonces de si estas moléculas/iones son más estables en el disolvente que en el cristal.
(Estos procesos no indican la rapidez con la que se disuelve un cristal. Eso depende de la barrera de energía de activación para salir del cristal y entrar en la solución, no de las energías relativas del cristal y de las moléculas/iones solvatados).
En un cristal iónico los iones se mantienen unidos por fuerzas culombianas, que son de largo alcance y tienen un $1/r$ dependencia de la distancia. Para determinar la energía de enlace exacta en un cristal hay que tener en cuenta todas las interacciones con otros iones, no sólo con los vecinos más cercanos. Se trata de un cálculo estándar cuya parte más difícil es encontrar la constante de Madelung para el cristal en cuestión.
En la superficie del cristal, las interacciones son mucho menores que en la masa, ya que estas moléculas/iones están en contacto con sólo la mitad de los iones del cristal y también interactúan con el disolvente. Para estimar la solubilidad, la energía de Coulomb de separar dos iones a una distancia r en un disolvente de constante dieléctrica (permitividad relativa) $\epsilon$ se calcula. Esta energía es $$w =\frac{z_1z_2e^2}{4\pi\epsilon_0\epsilon r}$$ donde $\epsilon _0$ es la permitividad del espacio libre, e la carga electrónica y $z_i$ la carga del ion i y la energía en $r=\infty$ se toma como cero.
La ley de Coulomb no es realmente válida a pequeñas separaciones, ya que la molecularidad del medio hace que la hipótesis del continuo, en términos de constante dieléctrica, se rompa. No obstante, ofrece una visión, aunque sea aproximada, de la disolución.
La energía libre para separar dos iones con radios $r_+$ y $r_-$ es por lo tanto $$\Delta\mu \approx \frac{+z_1z_2e^2}{4\pi\epsilon_0\epsilon (r_+ + r_-)}$$ que es positiva, ya que la interacción coulombiana siempre será negativa para combinar cargas opuestas.
La distribución de Boltzmann para la concentración $X_1$ y $X_2$ en dos fases coexistentes es $$ X_1=X_2exp(-\frac{(\mu_1-\mu_2)}{k_BT})=X_2exp(-\frac{\Delta\mu}{k_BT})$$ Si asociamos la concentración $X$ como una fracción molar, entonces esto puede ser una medida de la solubilidad, $X_s$ Por lo tanto $$ X_s \approx exp(-\frac{(\Delta\mu}{k_BT})=exp[-\frac{e^2}{4\pi\epsilon_0\epsilon k_BT}\frac{1}{(r_1+r_2)}]$$
En el agua la fracción molar de solubilidad del NaCl es $0.11$ que es mayor que el calculado $X_s =0.08 $ pero no por mucho teniendo en cuenta la simplicidad del modelo.
La última ecuación tiene la forma general $$X_s = exp(-C/\epsilon)$$ donde $C$ es una constante que depende de los pares de iones y de la temperatura. Por lo tanto, un gráfico del logaritmo de la solubilidad (como fracción molar) frente al recíproco de la constante dieléctrica para los mismos pares de iones en diferentes disolventes debería dar una línea recta. El siguiente gráfico muestra los datos del NaCl en diferentes disolventes, la relación aproximada se confirma. Debido a la naturaleza recíproca de los radios iónicos en la ecuación, los iones más grandes deberían ser generalmente más solubles que los más pequeños.
Los disolventes( $\epsilon$ ) de izquierda a derecha son metilformamida(182,4), formamida(109,5), agua(78,5), etilenglicol(40,7), etanolamina(37,72), metanol(32,6), etanol(24,3), propanol(20,2),butanol(17,8), pentanol(13,9). (Datos adaptados de J. Israelachvilli, Intermolecular and Surface Forces).
El enorme rango de solubilidad es bastante dramático, y muestra que el principal efecto sobre la solubilidad se debe a la constante dieléctrica del disolvente $\epsilon$ . Esto se puede apreciar en la ley de Coulomb. El campo eléctrico alrededor de un ion se "apaga" en un dieléctrico alto, lo que significa que la energía de interacción disminuye rápidamente con la distancia. De este modo, un ion queda aislado de todos los demás y está rodeado únicamente por los dipolos del disolvente. En un dieléctrico bajo, el campo eléctrico del ion se extiende a una distancia mucho mayor (en la relación de constantes dieléctricas, que puede ser 20 veces mayor), por lo que un ion puede "sentir" el efecto de otros cercanos cuando se difunden cerca de otro y, si se atraen, pueden unirse para formar una molécula iónica, por lo que dejan de ser solubles.
En la figura todos estos disolventes son bastante similares y disolventes de enlace de hidrógeno en mayor o menor medida. Parece que los disolventes sin enlace de hidrógeno son menos eficaces para disolver compuestos iónicos, por ejemplo, la solubilidad del NaCl en acetona ( $\epsilon = 20.7$ ) es $X_s=10^{-7} $ mucho menor que para el propanol con una $\epsilon$