Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Capítulo 5: Termodinámica y Cinética del Transporte Iónico

Cinvestav Monterrey

En este capítulo, estudiaremos cómo las fuerzas químicas y eléctricas se integran para gobernar el movimiento de partículas cargadas a través de interfaces, un proceso fundamental tanto en la fisiología celular como en la electroquímica industrial.

El Trabajo Eléctrico y la Energía Interna

Para comprender el transporte iónico, debemos reconocer que los iones no solo están sujetos a gradientes de concentración, sino también a campos eléctricos. En termodinámica, el cambio en la energía interna UU de un sistema depende del calor y de todas las formas de trabajo realizado. Si el sistema contiene partículas cargadas, debemos incluir el trabajo eléctrico (WeW_e) necesario para traer una carga qq desde un punto de referencia (infinito) hasta un potencial eléctrico ϕ\phi:

We=qϕ.W_e = q \cdot \phi.

Dado que trabajamos con sustancias químicas, expresamos la carga en términos de moles (NiN_i). La carga total de una especie ii es el producto del número de moles por su valencia (ziz_i) y la constante de Faraday (FF), que representa la carga de un mol de electrones (F96,485F \approx 96,485 C/mol):

qi=ziFNi.q_i = z_i F N_i.

Al integrar este trabajo en la relación fundamental de Gibbs para un sistema abierto, obtenemos la expresión de la energía interna que incluye la contribución de la energía potencial eléctrica:

dU=TdSPdV+iμiquimdNi+iϕd(ziFNi).dU = TdS - PdV + \sum_i \mu_i^{quim} dN_i + \sum_i \phi d(z_i F N_i).

Aquí, el término ϕd(ziFNi)\sum \phi d(z_i F N_i) representa la energía necesaria para añadir iones a un sistema que ya posee un potencial eléctrico determinado.

Potencial Electroquímico (μ~i\tilde{\mu}_i)

Para la mayoría de los sistemas biológicos y químicos, el control experimental se ejerce sobre la presión (PP) y la temperatura (TT). Por ello, utilizamos la energía libre de Gibbs (GG). Al aplicar la transformación de Legendre (G=U+PVTSG = U + PV - TS) a la ecuación anterior, la diferencial se convierte en:

dG=SdT+VdP+i(μiquim+ziFϕ)dNi.dG = -SdT + VdP + \sum_i \left( \mu_i^{quim} + z_i F \phi \right) dN_i.

El término entre paréntesis representa la energía total aportada por un mol de la especie ii, considerando tanto sus propiedades químicas como su interacción con el campo eléctrico. Definimos formalmente este valor como el potencial electroquímico (μ~i\tilde{\mu}_i):

μ~i=(GNi)T,P,Nji=μ~i=μiquim+ziFϕ.\tilde{\mu}_i = \left( \frac{\partial G}{\partial N_i} \right)_{T, P, N_{j \neq i}} = \tilde{\mu}_i = \mu_i^{quim} + z_i F \phi.

Sustituyendo la definición del potencial químico para una solución ideal (μiquim=μi0+RTlnci\mu_i^{quim} = \mu_i^0 + RT \ln c_i), obtenemos la ecuación fundamental del transporte iónico:

μ~i=μi0+RTlnci+ziFϕ††.\tilde{\mu}_i = \underbrace{\mu_i^0 + RT \ln c_i}_{†} + \underbrace{z_i F \phi}_{††}.

Condición de Equilibrio y Ecuación de Nernst

El equilibrio termodinámico entre dos compartimentos (1 y 2) separados por una membrana se define como el estado donde no existe una fuerza impulsora neta; es decir, mover un ion del lado 1 al lado 2 no produce un cambio en la energía libre del sistema. Esto ocurre cuando los potenciales electroquímicos se igualan:

μ~i(1)=μ~i(2),μi0+RTlnci(1)+ziFϕ(1)=μi0+RTlnci(2)+ziFϕ(2).\begin{align*} \tilde{\mu}_i^{(1)} &= \tilde{\mu}_i^{(2)}, \\ \mu_i^0 + RT \ln c_i^{(1)} + z_i F \phi^{(1)} &= \mu_i^0 + RT \ln c_i^{(2)} + z_i F \phi^{(2)}. \end{align*}

Sustituyendo los términos correspondientes y asumiendo que el potencial estándar μi0\mu_i^0 es el mismo en ambos lados de la membrana:

ziF(ϕ(2)ϕ(1))=RTlnci(1)RTlnci(2),ziFΔϕ=RTlnci(1)ci(2).\begin{align*} z_i F (\phi^{(2)} - \phi^{(1)}) &= RT \ln c_i^{(1)} - RT \ln c_i^{(2)}, \\ z_i F \Delta \phi &= RT \ln \frac{c_i^{(1)}}{c_i^{(2)}}. \end{align*}

Para hallar la diferencia de potencial eléctrico (Δϕ=ϕ(2)ϕ(1)\Delta \phi = \phi^{(2)} - \phi^{(1)}) necesaria para mantener esta diferencia de concentraciones en equilibrio, reordenamos los términos:

ziF(ϕ(2)ϕ(1))=RTlnci(1)RTlnci(2).z_i F (\phi^{(2)} - \phi^{(1)}) = RT \ln c_i^{(1)} - RT \ln c_i^{(2)}.

Utilizando la propiedad de los logaritmos (lnAlnB=lnAB\ln A - \ln B = \ln \frac{A}{B}):

ziFΔϕ=RTlnci(1)ci(2).z_i F \Delta \phi = RT \ln \frac{c_i^{(1)}}{c_i^{(2)}}.

Finalmente, despejando el potencial, llegamos a la ecuación de Nernst:

Δϕ=RTziFlnci(1)ci(2).\Delta \phi = \frac{RT}{z_i F} \ln \frac{c_i^{(1)}}{c_i^{(2)}}.

Interpretación Física: El potencial de Nernst es el voltaje específico que debe existir a través de una membrana para que un ion esté en equilibrio pese a tener concentraciones desiguales en ambos lados. En este estado, el flujo neto es cero porque la fuerza eléctrica empuja a los iones con la misma intensidad, pero en sentido opuesto, que la fuerza difusiva (gradiente químico).

Cinética del Transporte: La Ecuación de Nernst-Planck

Cuando los potenciales electroquímicos no son iguales (μ~i(1)μ~i(2)\tilde{\mu}_i^{(1)} \neq \tilde{\mu}_i^{(2)}), el sistema se encuentra fuera del equilibrio y se genera un flujo neto (JiJ_i). Según la termodinámica de procesos irreversibles, el flujo es proporcional a la fuerza impulsora, la cual se define como el negativo del gradiente del potencial electroquímico (μ~i-\nabla \tilde{\mu}_i). La relación fenomenológica es:

Ji=uiciμ~i,J_i = -u_i c_i \nabla \tilde{\mu}_i,

donde uiu_i es la movilidad del ion y cic_i es la concentración local. Procedemos a derivar el gradiente de μ~i\tilde{\mu}_i:

μ~i=(μi0+RTlnci+ziFϕ).\nabla \tilde{\mu}_i = \nabla (\mu_i^0 + RT \ln c_i + z_i F \phi).

Considerando que μi0\mu_i^0 es constante en el espacio (μi0=0\nabla \mu_i^0 = 0) y aplicando la regla de la cadena para el término logarítmico (lnc=1cc\nabla \ln c = \frac{1}{c} \nabla c):

μ~i=RTcici+ziFϕ.\nabla \tilde{\mu}_i = \frac{RT}{c_i} \nabla c_i + z_i F \nabla \phi.

Al sustituir este gradiente en la relación fenomenológica inicial:

Ji=uici(RTcici+ziFϕ)=uiRTciuiziFciϕ\begin{align*} J_i &= -u_i c_i \left( \frac{RT}{c_i} \nabla c_i + z_i F \nabla \phi \right) &= -u_i RT \nabla c_i - u_i z_i F c_i \nabla \phi \end{align*}

Finalmente, aplicamos la relación de Einstein, que vincula la movilidad iónica con el coeficiente de difusión macroscópico (Di=uiRTD_i = u_i RT). Obtenemos así la ecuación de Nernst-Planck:

Ji=Di(ci+ziFciRTϕ).J_i = -D_i \left( \nabla c_i + \frac{z_i F c_i}{RT} \nabla \phi \right).

Análisis de los términos:

  1. Difusión: Representa el flujo provocado por el movimiento aleatorio térmico (Ley de Fick). Los iones tienden a desplazarse de zonas de mayor a menor concentración para maximizar la entropía del sistema.

  2. Migración Eléctrica: Representa el flujo provocado por la fuerza electrostática sobre las partículas cargadas. Los cationes migran hacia potenciales negativos y los aniones hacia potenciales positivos.

Esta ecuación es la piedra angular para modelar procesos como el potencial de acción neuronal o la electrodiálisis, permitiendo predecir cómo varían las concentraciones y el voltaje en sistemas dinámicos fuera del equilibrio.

Recuperación de la Ecuación de Nernst

La ecuación de Nernst-Planck describe el flujo iónico como una respuesta dinámica a gradientes químicos y eléctricos. Sin embargo, en el equilibrio termodinámico no hay transporte neto de materia (Ji=0J_i = 0), a pesar de que el movimiento térmico individual de los iones persiste.

Para demostrar cómo la cinética converge hacia la termodinámica clásica, consideremos el transporte en una sola dimensión (xx) a través del espesor de una membrana. Si igualamos el flujo a cero en la ecuación de Nernst-Planck:

0=Di(dcidx+ziFciRTdϕdx).0 = -D_i \left( \frac{dc_i}{dx} + \frac{z_i F c_i}{RT} \frac{d\phi}{dx} \right).

Dado que el coeficiente de difusión DiD_i es una propiedad física del sistema distinta de cero, la expresión dentro del paréntesis debe anularse necesariamente para satisfacer la igualdad. Reordenando los términos para separar las variables de concentración (cic_i) y de potencial eléctrico (ϕ\phi), obtenemos:

dcidx=ziFciRTdϕdx.\frac{dc_i}{dx} = -\frac{z_i F c_i}{RT} \frac{d\phi}{dx}.

Dividiendo ambos lados por cic_i y eliminando dxdx, eliminamos la dependencia espacial explícita, transformando la ecuación en una relación de diferenciales exactas:

dcici=ziFRTdϕ.\frac{dc_i}{c_i} = -\frac{z_i F}{RT} d\phi.

Para encontrar la diferencia de potencial entre dos compartimentos, integramos la ecuación desde el compartimento 1 (estado inicial) hasta el compartimento 2 (estado final):

ci(1)ci(2)1cidci=ziFRTϕ(1)ϕ(2)dϕ.\int_{c_i^{(1)}}^{c_i^{(2)}} \frac{1}{c_i} dc_i = -\frac{z_i F}{RT} \int_{\phi^{(1)}}^{\phi^{(2)}} d\phi.

Al resolver las integrales, obtenemos:

lnci(2)lnci(1)=ziFRT(ϕ(2)ϕ(1)),lnci(2)ci(1)=ziFRTΔϕ.\begin{align*} \ln c_i^{(2)} - \ln c_i^{(1)} &= -\frac{z_i F}{RT} (\phi^{(2)} - \phi^{(1)}), \\ \ln \frac{c_i^{(2)}}{c_i^{(1)}} &= -\frac{z_i F}{RT} \Delta \phi. \\ \end{align*}

Finalmente, despejando Δϕ\Delta \phi e invirtiendo el argumento del logaritmo para eliminar el signo negativo recuperamos la ecuación de Nernst:

Δϕ=RTziFlnci(1)ci(2).\Delta \phi = \frac{RT}{z_i F} \ln \frac{c_i^{(1)}}{c_i^{(2)}}.

Esto demuestra que el equilibrio de Nernst es el caso límite de la ecuación de Nernst-Planck cuando las fuerzas impulsoras se cancelan exactamente. Físicamente, esto significa que el potencial eléctrico de la membrana se ajusta hasta compensar la tendencia del ion a difundirse a favor de su gradiente, logrando que el potencial electroquímico sea uniforme (μ~i=0\nabla \tilde{\mu}_i = 0) y, por ende, el flujo neto sea nulo.