Balance detallado

Principio en sistemas cinéticos

El principio de equilibrio detallado se puede utilizar en sistemas cinéticos que se descomponen en procesos elementales (colisiones, pasos o reacciones elementales). Establece que en el equilibrio , cada proceso elemental está en equilibrio con su proceso inverso.

Historia

El principio de equilibrio detallado fue introducido explícitamente para las colisiones por Ludwig Boltzmann . En 1872, demostró su teorema H utilizando este principio. [1] Los argumentos a favor de esta propiedad se basan en la reversibilidad microscópica . [2]

Cinco años antes de Boltzmann, James Clerk Maxwell utilizó el principio de equilibrio detallado para la cinética de los gases con referencia al principio de razón suficiente . [3] Comparó la idea del equilibrio detallado con otros tipos de equilibrio (como el equilibrio cíclico) y encontró que "Ahora es imposible asignar una razón" por la cual el equilibrio detallado debería ser rechazado (pág. 64).

En 1901, Rudolf Wegscheider introdujo el principio de equilibrio de detalle para la cinética química. [4] En particular, demostró que los ciclos irreversibles son imposibles y encontró explícitamente las relaciones entre las constantes cinéticas que se desprenden del principio de equilibrio de detalle. En 1931, Lars Onsager utilizó estas relaciones en sus trabajos, [5] por lo que recibió el Premio Nobel de Química en 1968 . A 1 A 2 A norte A 1 {\displaystyle {\ce {A1->A2->\cdots ->A_{\mathit {n}}->A1}}}

Albert Einstein utilizó en 1916 el principio del equilibrio detallado como base para su teoría cuántica de emisión y absorción de radiación. [6]

El principio de equilibrio detallado se ha utilizado en los métodos de Monte Carlo de cadenas de Markov desde su invención en 1953. [7] En particular, en el algoritmo Metropolis-Hastings y en su importante caso particular, el muestreo de Gibbs , se utiliza como una condición simple y confiable para proporcionar el estado de equilibrio deseable.

En la actualidad, el principio de equilibrio detallado es una parte estándar de los cursos universitarios de mecánica estadística, química física , cinética química y física. [8] [9] [10]

Fondo microscópico

La microscópica "inversión del tiempo" se convierte, a nivel cinético, en una "inversión de las flechas": los procesos elementales se transforman en sus procesos inversos. Por ejemplo, la reacción

i alfa i A i yo β yo B yo {\displaystyle \sum _{i}\alpha _{i}{\ce {A}}_{i}{\ce {->}}\sum _{j}\beta _{j}{\ce { B}}_ {j}} se transforma en yo β yo B yo i alfa i A i {\displaystyle \sum _{j}\beta _{j}{\ce {B}}_{j}{\ce {->}}\sum _{i}\alpha _{i}{\ce { Ai}}

y viceversa. (Aquí, son símbolos de componentes o estados, son coeficientes). El conjunto de equilibrio debería ser invariante con respecto a esta transformación debido a la microreversibilidad y la unicidad del equilibrio termodinámico. Esto nos lleva inmediatamente al concepto de equilibrio detallado: cada proceso se equilibra mediante su proceso inverso. A i , B yo {\displaystyle {\ce {A}}_{i},{\ce {B}}_{j}} alfa i , β yo 0 {\displaystyle \alpha _{i},\beta _{j}\geq 0}

Este razonamiento se basa en tres supuestos:

  1. A i {\displaystyle {\ce {A}}_{i}} no cambia con la inversión del tiempo;
  2. El equilibrio es invariante bajo la inversión del tiempo;
  3. Los procesos elementales macroscópicos son microscópicamente distinguibles, es decir, representan conjuntos disjuntos de eventos microscópicos.

Cualquiera de estos supuestos puede ser violado. [11] Por ejemplo, la colisión de Boltzmann puede ser representada como , donde es una partícula con velocidad v . Bajo la inversión del tiempo se transforma en . Por lo tanto, la colisión se transforma en la colisión inversa por la transformación PT , donde P es la inversión espacial y T es la inversión temporal. El balance detallado para la ecuación de Boltzmann requiere PT -invariancia de la dinámica de las colisiones, no solo T -invariancia. De hecho, después de la inversión temporal la colisión , se transforma en . Para el balance detallado necesitamos la transformación en . Para este propósito, necesitamos aplicar adicionalmente la inversión espacial P . Por lo tanto, para el balance detallado en la ecuación de Boltzmann no T -invariancia es necesaria sino PT -invariancia. A v + A w A v + A w {\displaystyle {\ce {{A_{\mathit {v}}}+A_{\mathit {w}}->{A_{\mathit {v'}}}+A_{\mathit {w'}}}}} A v {\displaystyle {\ce {A}}_{v}} A v {\displaystyle {\ce {A}}_{v}} A v {\displaystyle {\ce {A}}_{-v}} A v + A w A v + A w {\displaystyle {\ce {{A_{\mathit {v}}}+A_{\mathit {w}}->{A_{\mathit {v'}}}+A_{\mathit {w'}}}}} A v + A w A v + A w {\displaystyle {\ce {{A_{\mathit {-v'}}}+A_{\mathit {-w'}}->{A_{\mathit {-v}}}+A_{\mathit {-w}}}}} A v + A w A v + A w {\displaystyle {\ce {{A_{\mathit {v'}}}+A_{\mathit {w'}}->{A_{\mathit {v}}}+A_{\mathit {w}}}}}

El equilibrio puede no ser invariante respecto de T o PT incluso si las leyes del movimiento son invariantes. Esta no invariancia puede ser causada por la ruptura espontánea de la simetría . Existen medios no recíprocos (por ejemplo, algunos materiales biisotrópicos ) sin invariancia respecto de T y PT . [11]

Si se toman muestras de diferentes procesos macroscópicos a partir de los mismos eventos microscópicos elementales, entonces el equilibrio detallado macroscópico [ aclaración necesaria ] puede violarse incluso cuando se mantiene el equilibrio detallado microscópico. [11] [12]

Ahora, después de casi 150 años de desarrollo, el alcance de la validez y las violaciones del equilibrio detallado en la cinética parecen estar claras.

Balance detallado

Reversibilidad

Un proceso de Markov se denomina proceso de Markov reversible o cadena de Markov reversible si existe una distribución estacionaria positiva π que satisface las ecuaciones de equilibrio detalladas [13] donde P ij es la probabilidad de transición de Markov del estado i al estado j , es decir, P ij  =  P ( X t  =  j  |  X t  − 1  =  i ) , y π i y π j son las probabilidades de equilibrio de estar en los estados i y j , respectivamente. [13] Cuando Pr( X t −1  =  i ) = π i para todo i , esto es equivalente a la matriz de probabilidad conjunta, Pr( X t −1  =  iX t  =  j ) siendo simétrica en i y j ; o simétrica en t  − 1t . π i P i j = π j P j i , {\displaystyle \pi _{i}P_{ij}=\pi _{j}P_{ji}\,,}

La definición se aplica directamente a las variables continuas, donde π se convierte en una densidad de probabilidad y P ( s ′,  s ) en una densidad de probabilidad del núcleo de transición del estado s ′ al estado  s : La condición de equilibrio detallado es más fuerte que la requerida simplemente para una distribución estacionaria , porque hay procesos de Markov con distribuciones estacionarias que no tienen equilibrio detallado. π ( s ) P ( s , s ) = π ( s ) P ( s , s ) . {\displaystyle \pi (s')P(s',s)=\pi (s)P(s,s')\,.}

Las matrices de transición que son simétricas ( P ij  =  P ji o P ( s ′,  s ) =  P ( ss ′)) siempre tienen un equilibrio detallado. En estos casos, una distribución uniforme sobre los estados es una distribución de equilibrio.

Criterio de Kolmogorov

La reversibilidad es equivalente al criterio de Kolmogorov : el producto de las tasas de transición en cualquier bucle cerrado de estados es el mismo en ambas direcciones.

Por ejemplo, implica que, para todos a , b y c , Por ejemplo, si tenemos una cadena de Markov con tres estados tales que solo son posibles estas transiciones: , entonces violan el criterio de Kolmogorov. P ( a , b ) P ( b , c ) P ( c , a ) = P ( a , c ) P ( c , b ) P ( b , a ) . {\displaystyle P(a,b)P(b,c)P(c,a)=P(a,c)P(c,b)P(b,a)\,.} A B , B C , C A , B A {\displaystyle A\to B,B\to C,C\to A,B\to A}

En el caso de sistemas continuos con un equilibrio detallado, puede ser posible transformar continuamente las coordenadas hasta que la distribución de equilibrio sea uniforme, con un núcleo de transición que entonces sea simétrico. En el caso de estados discretos, puede ser posible [ aclaración necesaria ] lograr algo similar descomponiendo los estados de Markov en subestados degenerados de tamaño adecuado.

En el caso de una matriz de transición de Markov y una distribución estacionaria, las ecuaciones de equilibrio detalladas pueden no ser válidas. Sin embargo, se puede demostrar que existe una única matriz de transición de Markov que es la más cercana según la distribución estacionaria y una norma dada. La matriz más cercana se puede calcular resolviendo un problema de optimización cuadrático-convexa.

Balance detallado y aumento de entropía.

Para muchos sistemas de cinética física y química, el equilibrio detallado proporciona condiciones suficientes para el aumento estricto de la entropía en sistemas aislados. Por ejemplo, el famoso teorema H de Boltzmann [1] establece que, según la ecuación de Boltzmann, el principio de equilibrio detallado implica positividad de la producción de entropía. La fórmula de Boltzmann (1872) para la producción de entropía en la cinética de gases enrarecidos con equilibrio detallado [1] [2] sirvió como prototipo de muchas fórmulas similares para la disipación en la cinética de acción de masas [14] y la cinética de acción de masas generalizada [15] con equilibrio detallado.

Sin embargo, el principio de equilibrio detallado no es necesario para el crecimiento de la entropía. Por ejemplo, en el ciclo lineal irreversible , la producción de entropía es positiva, pero el principio de equilibrio detallado no se cumple. A 1 A 2 A 3 A 1 {\displaystyle {\ce {A1 -> A2 -> A3 -> A1}}}

Por lo tanto, el principio de equilibrio detallado es una condición suficiente pero no necesaria para el aumento de la entropía en la cinética de Boltzmann. Estas relaciones entre el principio de equilibrio detallado y la segunda ley de la termodinámica se aclararon en 1887 cuando Hendrik Lorentz objetó el teorema H de Boltzmann para gases poliatómicos. [16] Lorentz afirmó que el principio de equilibrio detallado no es aplicable a las colisiones de moléculas poliatómicas.

Boltzmann inventó inmediatamente una nueva condición, más general, suficiente para el crecimiento de la entropía. [17] La ​​condición de Boltzmann es válida para todos los procesos de Markov, independientemente de la reversibilidad temporal. Más tarde, el aumento de la entropía se demostró para todos los procesos de Markov mediante un método directo. [18] [19] Estos teoremas pueden considerarse como simplificaciones del resultado de Boltzmann. Más tarde, esta condición se denominó condición de "equilibrio cíclico" (porque se cumple para ciclos irreversibles) o "equilibrio semidetallado" o "equilibrio complejo". En 1981, Carlo Cercignani y Maria Lampis demostraron que los argumentos de Lorentz eran erróneos y que el principio de equilibrio detallado es válido para moléculas poliatómicas. [20] Sin embargo, las condiciones de equilibrio semidetallado extendidas inventadas por Boltzmann en esta discusión siguen siendo la notable generalización del equilibrio detallado.

Condiciones de Wegscheider para la ley de acción de masas generalizada

En cinética química , las reacciones elementales se representan mediante ecuaciones estequiométricas donde son las componentes y son los coeficientes estequiométricos. Aquí, las reacciones inversas con constantes positivas se incluyen en la lista por separado. Necesitamos esta separación de reacciones directas e inversas para aplicar más adelante el formalismo general a los sistemas con algunas reacciones irreversibles. El sistema de ecuaciones estequiométricas de reacciones elementales es el mecanismo de reacción . i α r i A i j β r j A j ( r = 1 , , m ) , {\displaystyle \sum _{i}\alpha _{ri}{\ce {A}}_{i}{\ce {->}}\sum _{j}\beta _{rj}{\ce {A}}_{j}\;\;(r=1,\ldots ,m)\,,} A i {\displaystyle {\ce {A}}_{i}} α r i , β r j 0 {\displaystyle \alpha _{ri},\beta _{rj}\geq 0}

La matriz estequiométrica es , (ganancia menos pérdida). Esta matriz no necesita ser cuadrada. El vector estequiométrico es la fila r de con coordenadas . Γ = ( γ r i ) {\displaystyle {\boldsymbol {\Gamma }}=(\gamma _{ri})} γ r i = β r i α r i {\displaystyle \gamma _{ri}=\beta _{ri}-\alpha _{ri}} γ r {\displaystyle \gamma _{r}} Γ {\displaystyle {\boldsymbol {\Gamma }}} γ r i = β r i α r i {\displaystyle \gamma _{ri}=\beta _{ri}-\alpha _{ri}}

Según la ley de acción de masas generalizada , la velocidad de reacción para una reacción elemental es donde es la actividad (la "concentración efectiva") de . w r = k r i = 1 n a i α r i , {\displaystyle w_{r}=k_{r}\prod _{i=1}^{n}a_{i}^{\alpha _{ri}}\,,} a i 0 {\displaystyle a_{i}\geq 0} A i {\displaystyle A_{i}}

El mecanismo de reacción incluye reacciones con constantes de velocidad de reacción . Para cada r se utilizan las siguientes notaciones: ; ; es la constante de velocidad de reacción para la reacción inversa si está en el mecanismo de reacción y 0 si no lo está; es la velocidad de reacción para la reacción inversa si está en el mecanismo de reacción y 0 si no lo está. Para una reacción reversible, es la constante de equilibrio . k r > 0 {\displaystyle k_{r}>0} k r + = k r {\displaystyle k_{r}^{+}=k_{r}} w r + = w r {\displaystyle w_{r}^{+}=w_{r}} k r {\displaystyle k_{r}^{-}} w r {\displaystyle w_{r}^{-}} K r = k r + / k r {\displaystyle K_{r}=k_{r}^{+}/k_{r}^{-}}

El principio de equilibrio detallado para la ley generalizada de acción de masas es: Para valores dados existe un equilibrio positivo que satisface el equilibrio detallado, es decir, . Esto significa que el sistema de ecuaciones de equilibrio detallado lineal es solucionable ( ). El siguiente resultado clásico da las condiciones necesarias y suficientes para la existencia de un equilibrio positivo con equilibrio detallado (véase, por ejemplo, el libro de texto [9] ). k r {\displaystyle k_{r}} a i e q > 0 {\displaystyle a_{i}^{\rm {eq}}>0} w r + = w r {\displaystyle w_{r}^{+}=w_{r}^{-}} i γ r i x i = ln k r + ln k r = ln K r {\displaystyle \sum _{i}\gamma _{ri}x_{i}=\ln k_{r}^{+}-\ln k_{r}^{-}=\ln K_{r}} x i = ln a i e q {\displaystyle x_{i}=\ln a_{i}^{\rm {eq}}} a i e q > 0 {\displaystyle a_{i}^{\rm {eq}}>0}

Dos condiciones son suficientes y necesarias para la solubilidad del sistema de ecuaciones de equilibrio detalladas:

  1. Si entonces y, a la inversa, si entonces (reversibilidad); k r + > 0 {\displaystyle k_{r}^{+}>0} k r > 0 {\displaystyle k_{r}^{-}>0} k r > 0 {\displaystyle k_{r}^{-}>0} k r + > 0 {\displaystyle k_{r}^{+}>0}
  2. Para cualquier solución del sistema λ = ( λ r ) {\displaystyle {\boldsymbol {\lambda }}=(\lambda _{r})} λ Γ = 0 ( i.e. r λ r γ r i = 0 for all i ) {\displaystyle {\boldsymbol {\lambda \Gamma }}=0\;\;\left({\mbox{i.e.}}\;\;\sum _{r}\lambda _{r}\gamma _{ri}=0\;\;{\mbox{for all}}\;\;i\right)}

La identidad de Wegscheider [21] es válida: r = 1 m ( k r + ) λ r = r = 1 m ( k r ) λ r . {\displaystyle \prod _{r=1}^{m}(k_{r}^{+})^{\lambda _{r}}=\prod _{r=1}^{m}(k_{r}^{-})^{\lambda _{r}}\,.}

Observación. Es suficiente utilizar en las condiciones de Wegscheider una base de soluciones del sistema . λ Γ = 0 {\displaystyle {\boldsymbol {\lambda \Gamma }}=0}

En particular, para cualquier ciclo de reacciones monomoleculares (lineales), el producto de las constantes de velocidad de reacción en el sentido horario es igual al producto de las constantes de velocidad de reacción en el sentido antihorario. La misma condición es válida para los procesos Markov reversibles (es equivalente a la condición de "sin flujo neto").

Un ejemplo no lineal simple nos da un ciclo lineal complementado con un paso no lineal: [21]

  1. A 1 A 2 {\displaystyle {\ce {A1 <=> A2}}}
  2. A 2 A 3 {\displaystyle {\ce {A2 <=> A3}}}
  3. A 3 A 1 {\displaystyle {\ce {A3 <=> A1}}}
  4. A 1 + A 2 2 A 3 {\displaystyle {\ce {{A1}+A2 <=> 2A3}}}

Hay dos identidades de Wegscheider independientes no triviales para este sistema: y Corresponden a las siguientes relaciones lineales entre los vectores estequiométricos: y k 1 + k 2 + k 3 + = k 1 k 2 k 3 {\displaystyle k_{1}^{+}k_{2}^{+}k_{3}^{+}=k_{1}^{-}k_{2}^{-}k_{3}^{-}} k 3 + k 4 + / k 2 + = k 3 k 4 / k 2 {\displaystyle k_{3}^{+}k_{4}^{+}/k_{2}^{+}=k_{3}^{-}k_{4}^{-}/k_{2}^{-}} γ 1 + γ 2 + γ 3 = 0 {\displaystyle \gamma _{1}+\gamma _{2}+\gamma _{3}=0} γ 3 + γ 4 γ 2 = 0. {\displaystyle \gamma _{3}+\gamma _{4}-\gamma _{2}=0.}

El aspecto computacional de las condiciones de Wegscheider fue estudiado por D. Colquhoun con coautores. [22]

Las condiciones de Wegscheider demuestran que, si bien el principio de equilibrio detallado enuncia una propiedad local del equilibrio, implica las relaciones entre las constantes cinéticas que son válidas para todos los estados alejados del equilibrio. Esto es posible porque se conoce una ley cinética y las relaciones entre las velocidades de los procesos elementales en equilibrio se pueden transformar en relaciones entre constantes cinéticas que se utilizan globalmente. Para las condiciones de Wegscheider, esta ley cinética es la ley de acción de masas (o la ley generalizada de acción de masas).

Disipación en sistemas con balance detallado

Para describir la dinámica de los sistemas que obedecen a la ley generalizada de acción de masas, hay que representar las actividades como funciones de las concentraciones c j y la temperatura . Para ello, se utiliza la representación de la actividad a través del potencial químico: donde μ i es el potencial químico de la especie en las condiciones de interés, es el potencial químico de esa especie en el estado estándar elegido , R es la constante de los gases y T es la temperatura termodinámica . El potencial químico puede representarse como una función de c y T , donde c es el vector de concentraciones con componentes c j . Para los sistemas ideales, y : la actividad es la concentración y la ley generalizada de acción de masas es la ley usual de acción de masas . a i = exp ( μ i μ i R T ) {\displaystyle a_{i}=\exp \left({\frac {\mu _{i}-\mu _{i}^{\ominus }}{RT}}\right)} μ i {\displaystyle \mu _{i}^{\ominus }} μ i = R T ln c i + μ i {\displaystyle \mu _{i}=RT\ln c_{i}+\mu _{i}^{\ominus }} a j = c j {\displaystyle a_{j}=c_{j}}

Consideremos un sistema en condición isotérmica ( T = constante) isocórica (el volumen V = constante). Para estas condiciones, la energía libre de Helmholtz ⁠ ⁠ F ( T , V , N ) {\displaystyle F(T,V,N)} mide el trabajo “útil” obtenible de un sistema. Es una función de la temperatura T , el volumen V y las cantidades de componentes químicos N j (generalmente medidos en moles ), N es el vector con componentes N j . Para los sistemas ideales, F = R T i N i ( ln ( N i V ) 1 + μ i ( T ) R T ) . {\displaystyle F=RT\sum _{i}N_{i}\left(\ln \left({\frac {N_{i}}{V}}\right)-1+{\frac {\mu _{i}^{\ominus }(T)}{RT}}\right).}

El potencial químico es una derivada parcial: . μ i = F ( T , V , N ) / N i {\displaystyle \mu _{i}=\partial F(T,V,N)/\partial N_{i}}

Las ecuaciones cinéticas químicas son d N i d t = V r γ r i ( w r + w r ) . {\displaystyle {\frac {dN_{i}}{dt}}=V\sum _{r}\gamma _{ri}(w_{r}^{+}-w_{r}^{-}).}

Si el principio de equilibrio detallado es válido entonces para cualquier valor de T existe un punto positivo de equilibrio detallado c eq : El álgebra elemental da donde w r + ( c e q , T ) = w r ( c e q , T ) = w r e q {\displaystyle w_{r}^{+}(c^{\rm {eq}},T)=w_{r}^{-}(c^{\rm {eq}},T)=w_{r}^{\rm {eq}}} w r + = w r e q exp ( i α r i ( μ i μ i e q ) R T ) ; w r = w r e q exp ( i β r i ( μ i μ i e q ) R T ) ; {\displaystyle w_{r}^{+}=w_{r}^{\rm {eq}}\exp \left(\sum _{i}{\frac {\alpha _{ri}(\mu _{i}-\mu _{i}^{\rm {eq}})}{RT}}\right);\;\;w_{r}^{-}=w_{r}^{\rm {eq}}\exp \left(\sum _{i}{\frac {\beta _{ri}(\mu _{i}-\mu _{i}^{\rm {eq}})}{RT}}\right);} μ i e q = μ i ( c e q , T ) {\displaystyle \mu _{i}^{\rm {eq}}=\mu _{i}(c^{\rm {eq}},T)}

Para la disipación obtenemos de estas fórmulas: La desigualdad se cumple porque ln es una función monótona y, por tanto, las expresiones y tienen siempre el mismo signo. d F d t = i F ( T , V , N ) N i d N i d t = i μ i d N i d t = V R T r ( ln w r + ln w r ) ( w r + w r ) 0 {\displaystyle {\frac {dF}{dt}}=\sum _{i}{\frac {\partial F(T,V,N)}{\partial N_{i}}}{\frac {dN_{i}}{dt}}=\sum _{i}\mu _{i}{\frac {dN_{i}}{dt}}=-VRT\sum _{r}(\ln w_{r}^{+}-\ln w_{r}^{-})(w_{r}^{+}-w_{r}^{-})\leq 0} ln w r + ln w r {\displaystyle \ln w_{r}^{+}-\ln w_{r}^{-}} w r + w r {\displaystyle w_{r}^{+}-w_{r}^{-}}

Desigualdades similares [9] son ​​válidas para otras condiciones clásicas de los sistemas cerrados y las funciones características correspondientes: para condiciones isobáricas isotérmicas la energía libre de Gibbs disminuye, para los sistemas isocóricos con energía interna constante ( sistemas aislados ) la entropía aumenta, así como para los sistemas isobáricos con entalpía constante .

Relaciones recíprocas de Onsager y balance detallado

Supongamos que el principio de equilibrio detallado es válido. Entonces, para pequeñas desviaciones del equilibrio, la respuesta cinética del sistema puede aproximarse como linealmente relacionada con su desviación del equilibrio químico, dando las velocidades de reacción para la ley de acción de masas generalizada como: w r + = w r e q ( 1 + i α r i ( μ i μ i e q ) R T ) ; w r = w r e q ( 1 + i β r i ( μ i μ i e q ) R T ) ; {\displaystyle w_{r}^{+}=w_{r}^{\rm {eq}}\left(1+\sum _{i}{\frac {\alpha _{ri}(\mu _{i}-\mu _{i}^{\rm {eq}})}{RT}}\right);\;\;w_{r}^{-}=w_{r}^{\rm {eq}}\left(1+\sum _{i}{\frac {\beta _{ri}(\mu _{i}-\mu _{i}^{\rm {eq}})}{RT}}\right);}

Por lo tanto, nuevamente en el régimen de respuesta lineal cerca del equilibrio, las ecuaciones cinéticas son ( ): γ r i = β r i α r i {\displaystyle \gamma _{ri}=\beta _{ri}-\alpha _{ri}} d N i d t = V j [ r w r e q γ r i γ r j ] μ j μ j e q R T . {\displaystyle {\frac {dN_{i}}{dt}}=-V\sum _{j}\left[\sum _{r}w_{r}^{\rm {eq}}\gamma _{ri}\gamma _{rj}\right]{\frac {\mu _{j}-\mu _{j}^{\rm {eq}}}{RT}}.}

Esta es exactamente la forma de Onsager: siguiendo el trabajo original de Onsager, [5] deberíamos introducir las fuerzas termodinámicas y la matriz de coeficientes en la forma X j {\displaystyle X_{j}} L i j {\displaystyle L_{ij}} X j = μ j μ j e q T ; d N i d t = j L i j X j {\displaystyle X_{j}={\frac {\mu _{j}-\mu _{j}^{\rm {eq}}}{T}};\;\;{\frac {dN_{i}}{dt}}=\sum _{j}L_{ij}X_{j}}

La matriz de coeficientes es simétrica: L i j {\displaystyle L_{ij}} L i j = V R r w r e q γ r i γ r j {\displaystyle L_{ij}=-{\frac {V}{R}}\sum _{r}w_{r}^{\rm {eq}}\gamma _{ri}\gamma _{rj}}

Estas relaciones de simetría, , son exactamente las relaciones recíprocas de Onsager . La matriz de coeficientes no es positiva. Es negativa en el intervalo lineal de los vectores estequiométricos . L i j = L j i {\displaystyle L_{ij}=L_{ji}} L {\displaystyle L} γ r {\displaystyle \gamma _{r}}

Así pues, las relaciones de Onsager se derivan del principio de equilibrio detallado en la aproximación lineal cerca del equilibrio.

Balance semidetallado

Para formular el principio de balance semidetallado, es conveniente contar por separado las reacciones elementales directas e inversas. En este caso, las ecuaciones cinéticas tienen la forma: Usemos las notaciones , para los vectores de entrada y salida de los coeficientes estequiométricos de la r -ésima reacción elemental. Sea el conjunto de todos estos vectores . d N i d t = V r γ r i w r = V r ( β r i α r i ) w r {\displaystyle {\frac {dN_{i}}{dt}}=V\sum _{r}\gamma _{ri}w_{r}=V\sum _{r}(\beta _{ri}-\alpha _{ri})w_{r}} α r = α r i {\displaystyle \alpha _{r}=\alpha _{ri}} β r = β r i {\displaystyle \beta _{r}=\beta _{ri}} Y {\displaystyle Y} α r , β r {\displaystyle \alpha _{r},\beta _{r}}

Para cada , definamos dos conjuntos de números: ν Y {\displaystyle \nu \in Y} R ν + = { r | α r = ν } ; R ν = { r | β r = ν } {\displaystyle R_{\nu }^{+}=\{r|\alpha _{r}=\nu \};\;\;\;R_{\nu }^{-}=\{r|\beta _{r}=\nu \}}

r R ν + {\displaystyle r\in R_{\nu }^{+}} si y solo si es el vector de los coeficientes estequiométricos de entrada para la r -ésima reacción elemental; si y solo si es el vector de los coeficientes estequiométricos de salida para la r -ésima reacción elemental. ν {\displaystyle \nu } α r {\displaystyle \alpha _{r}} r R ν {\displaystyle r\in R_{\nu }^{-}} ν {\displaystyle \nu } β r {\displaystyle \beta _{r}}

El principio de equilibrio semidetallado significa que en equilibrio se cumple la condición de equilibrio semidetallado: para cada ν Y {\displaystyle \nu \in Y} r R ν w r = r R ν + w r {\displaystyle \sum _{r\in R_{\nu }^{-}}w_{r}=\sum _{r\in R_{\nu }^{+}}w_{r}}

La condición de equilibrio semidetallado es suficiente para la estacionariedad: implica que d N d t = V r γ r w r = 0. {\displaystyle {\frac {dN}{dt}}=V\sum _{r}\gamma _{r}w_{r}=0.}

Para la cinética de Markov, la condición de equilibrio semidetallada es simplemente la ecuación de equilibrio elemental y se cumple para cualquier estado estacionario. Para la ley de acción de masas no lineal, es, en general, una condición suficiente pero no necesaria para la estacionariedad.

La condición de equilibrio semidetallado es más débil que la de equilibrio detallado: si se cumple el principio de equilibrio detallado, entonces también se cumple la condición de equilibrio semidetallado.

Para los sistemas que obedecen la ley de acción de masas generalizada, la condición de equilibrio semidetallado es suficiente para la desigualdad de disipación (para la energía libre de Helmholtz en condiciones isocóricas isotérmicas y para las desigualdades de disipación en otras condiciones clásicas para los potenciales termodinámicos correspondientes). d F / d t 0 {\displaystyle dF/dt\geq 0}

Boltzmann introdujo la condición de equilibrio semidetallado para colisiones en 1887 [17] y demostró que garantiza la positividad de la producción de entropía. Para la cinética química, esta condición (como condición de equilibrio complejo ) fue introducida por Horn y Jackson en 1972. [23]

Los antecedentes microscópicos para el equilibrio semidetallado se encontraron en la microcinética de Markov de los compuestos intermedios que están presentes en pequeñas cantidades y cuyas concentraciones están en cuasiequilibrio con los componentes principales. [24] Bajo estos supuestos microscópicos, la condición de equilibrio semidetallado es simplemente la ecuación de equilibrio para la microcinética de Markov según el teorema de MichaelisMentenStueckelberg . [25]

Disipación en sistemas con balance semidetallado

Representemos la ley de acción de masas generalizada en la forma equivalente: la velocidad del proceso elemental es donde es el potencial químico y es la energía libre de Helmholtz . El término exponencial se llama factor de Boltzmann y el multiplicador es el factor cinético. [25] Contemos la reacción directa e inversa en la ecuación cinética por separado: Una función auxiliar de una variable es conveniente para la representación de la disipación para la ley de acción de masas Esta función puede considerarse como la suma de las velocidades de reacción para coeficientes estequiométricos de entrada deformados . Porque es solo la suma de las velocidades de reacción. La función es convexa porque . i α r i A i i β r i A i {\displaystyle \sum _{i}\alpha _{ri}{\ce {A}}_{i}{\ce {->}}\sum _{i}\beta _{ri}{\ce {A}}_{i}} w r = φ r exp ( i α r i μ i R T ) {\displaystyle w_{r}=\varphi _{r}\exp \left(\sum _{i}{\frac {\alpha _{ri}\mu _{i}}{RT}}\right)} μ i = F ( T , V , N ) / N i {\displaystyle \mu _{i}=\partial F(T,V,N)/\partial N_{i}} F ( T , V , N ) {\displaystyle F(T,V,N)} φ r 0 {\displaystyle \varphi _{r}\geq 0} d N i d t = V r γ r i w r {\displaystyle {\frac {dN_{i}}{dt}}=V\sum _{r}\gamma _{ri}w_{r}} θ ( λ ) {\displaystyle \theta (\lambda )} λ [ 0 , 1 ] {\displaystyle \lambda \in [0,1]} θ ( λ ) = r φ r exp ( i ( λ α r i + ( 1 λ ) β r i ) ) μ i R T ) {\displaystyle \theta (\lambda )=\sum _{r}\varphi _{r}\exp \left(\sum _{i}{\frac {(\lambda \alpha _{ri}+(1-\lambda )\beta _{ri}))\mu _{i}}{RT}}\right)} θ ( λ ) {\displaystyle \theta (\lambda )} α ~ ρ ( λ ) = λ α ρ + ( 1 λ ) β ρ {\displaystyle {\tilde {\alpha }}_{\rho }(\lambda )=\lambda \alpha _{\rho }+(1-\lambda )\beta _{\rho }} λ = 1 {\displaystyle \lambda =1} θ ( λ ) {\displaystyle \theta (\lambda )} θ ( λ ) 0 {\displaystyle \theta ''(\lambda )\geq 0}

El cálculo directo da como resultado que, según las ecuaciones cinéticas, esta es la fórmula general de disipación para la ley de acción de masas generalizada . [25] d F d t = V R T d θ ( λ ) d λ | λ = 1 {\displaystyle {\frac {dF}{dt}}=-VRT\left.{\frac {d\theta (\lambda )}{d\lambda }}\right|_{\lambda =1}}

La convexidad de da las condiciones suficientes y necesarias para la desigualdad de disipación adecuada: θ ( λ ) {\displaystyle \theta (\lambda )} d F d t < 0  if and only if  θ ( λ ) < θ ( 1 )  for some  λ < 1 ; {\displaystyle {\frac {dF}{dt}}<0{\text{ if and only if }}\theta (\lambda )<\theta (1){\text{ for some }}\lambda <1;} d F d t 0  if and only if  θ ( λ ) θ ( 1 )  for some  λ < 1. {\displaystyle {\frac {dF}{dt}}\leq 0{\text{ if and only if }}\theta (\lambda )\leq \theta (1){\text{ for some }}\lambda <1.}

La condición de equilibrio semidetallado se puede transformar en identidad . Por lo tanto, para los sistemas con equilibrio semidetallado . [23] θ ( 0 ) θ ( 1 ) {\displaystyle \theta (0)\equiv \theta (1)} d F / d t 0 {\displaystyle {dF}/{dt}\leq 0}

Teorema del cono y equivalencia local del equilibrio detallado y complejo

Para cualquier mecanismo de reacción y un equilibrio positivo dado, se define un cono de velocidades posibles para los sistemas con equilibrio detallado para cualquier estado de no equilibrio N , donde el cono representa la envoltura cónica y las funciones constantes por partes no dependen de valores (positivos) de velocidades de reacción de equilibrio y están definidas por cantidades termodinámicas bajo el supuesto de equilibrio detallado. Q D B ( N ) = c o n e { γ r s g n ( w r + ( N ) w r ( N ) )   |   r = 1 , , m } , {\displaystyle \mathbf {Q} _{\rm {DB}}(N)={\rm {cone}}\{\gamma _{r}{\rm {sgn}}(w_{r}^{+}(N)-w_{r}^{-}(N))\ |\ r=1,\ldots ,m\},} s g n ( w r + ( N ) w r ( N ) ) {\displaystyle {\rm {sgn}}(w_{r}^{+}(N)-w_{r}^{-}(N))} w r e q {\displaystyle w_{r}^{\rm {eq}}}

El teorema del cono establece que para el mecanismo de reacción dado y el equilibrio positivo dado, la velocidad ( dN/dt ) en un estado N para un sistema con equilibrio complejo pertenece al cono . Es decir, existe un sistema con equilibrio detallado, el mismo mecanismo de reacción, el mismo equilibrio positivo, que da la misma velocidad en el estado N . [26] Según el teorema del cono, para un estado dado N , el conjunto de velocidades de los sistemas de equilibrio semidetallado coincide con el conjunto de velocidades de los sistemas de equilibrio detallado si sus mecanismos de reacción y equilibrios coinciden. Esto significa equivalencia local de equilibrio detallado y complejo. Q D B ( N ) {\displaystyle \mathbf {Q} _{\rm {DB}}(N)}

Balance detallado para sistemas con reacciones irreversibles

El balance detallado establece que en el equilibrio cada proceso elemental se equilibra mediante su proceso inverso y requiere la reversibilidad de todos los procesos elementales. Para muchos sistemas físicos y químicos complejos reales (por ejemplo, combustión homogénea, oxidación catalítica heterogénea, la mayoría de las reacciones enzimáticas, etc.), los mecanismos detallados incluyen reacciones reversibles e irreversibles. Si uno representa las reacciones irreversibles como límites de pasos reversibles, entonces resulta obvio que no todos los mecanismos de reacción con reacciones irreversibles pueden obtenerse como límites de sistemas o reacciones reversibles con balance detallado. Por ejemplo, el ciclo irreversible no puede obtenerse como tal límite, pero el mecanismo de reacción sí. [27] A 1 A 2 A 3 A 1 {\displaystyle {\ce {A1 -> A2 -> A3 -> A1}}} A 1 A 2 A 3 A 1 {\displaystyle {\ce {A1 -> A2 -> A3 <- A1}}}

Teorema de Gorban - Yablonsky . Un sistema de reacciones con algunas reacciones irreversibles es un límite de sistemas con equilibrio detallado cuando algunas constantes tienden a cero si y solo si (i) la parte reversible de este sistema satisface el principio de equilibrio detallado y (ii) la envoltura convexa de los vectores estequiométricos de las reacciones irreversibles tiene intersección vacía con el espacio lineal de los vectores estequiométricos de las reacciones reversibles. [21] Físicamente, la última condición significa que las reacciones irreversibles no pueden incluirse en vías cíclicas orientadas.

Véase también

Referencias

  1. ^ abc Boltzmann, L. (1964), Conferencias sobre teoría de los gases, Berkeley, CA, EE. UU.: U. of California Press.
  2. ^ ab Tolman, RC (1938). Los principios de la mecánica estadística . Oxford University Press, Londres, Reino Unido.
  3. ^ Maxwell, JC (1867), Sobre la teoría dinámica de los gases, Philosl Trans R Soc London, 157, págs. 49–88
  4. ^ Wegscheider, R. (1901) Über simultane Gleichgewichte und die Beziehungen zwischen Thermodynamik und Reactionskinetik homogener Systeme, Monatshefte für Chemie / Chemical Monthly 32(8), 849–906.
  5. ^ ab Onsager, L. (1931), Relaciones recíprocas en procesos irreversibles. Archivado el 26 de octubre de 2011 en Wayback Machine. I, Phys. Rev. 37, 405–426; II 38, 2265–2279
  6. ^ Einstein, A. (1916). Strahlungs-Emission und -Absortion nach der Quantentheorie [=Emisión y absorción de radiación en la teoría cuántica], Verhandlungen der Deutschen Physikalischen Gesellschaft 18 (13/14). Braunschweig: Vieweg, 318–323. Véase también: A. Einstein (1917). Zur Quantentheorie der Strahlung [=Sobre la teoría cuántica de la radiación], Physikalische Zeitschrift 18 (1917), 121-128. Traducción al inglés: D. ter Haar (1967): La antigua teoría cuántica. Prensa de Pérgamo, págs. 167–183.
  7. ^ Metropolis, N. ; Rosenbluth, AW; Rosenbluth, MN ; Teller, AH; Teller, E. (1953). "Cálculos de ecuaciones de estado mediante máquinas de computación rápida" . Journal of Chemical Physics . 21 (6): 1087–1092. Bibcode :1953JChPh..21.1087M. doi :10.1063/1.1699114. OSTI  4390578. S2CID  1046577.
  8. ^ van Kampen, NG "Procesos estocásticos en física y química", Elsevier Science (1992).
  9. ^ abc Yablonskii, GS, Bykov, VI, Gorban, AN , Elokhin, VI (1991), Modelos cinéticos de reacciones catalíticas, Ámsterdam, Países Bajos: Elsevier.
  10. ^ Lifshitz, EM; Pitaevskii, LP (1981). Cinética física . Londres: Pergamon. ISBN 978-0-08-026480-6.Vol. 10 del Curso de Física Teórica (3ª Ed).
  11. ^ abc Gorban, AN (2014), Balance detallado en micro y macrocinética y microdistinguibilidad de macroprocesos, Results in Physics 4, 142–147
  12. ^ Joshi, B. (2013), El equilibrio detallado determinista en redes de reacciones químicas es suficiente pero no necesario para el equilibrio detallado estocástico, arXiv:1312.4196 [math.PR].
  13. ^ ab O'Hagan, Anthony; Forster, Jonathan (2004). "Sección 10.3". Kendall's Advanced Theory of Statistics, Volumen 2B: Inferencia bayesiana . Nueva York: Oxford University Press. pág. 263. ISBN 978-0-340-80752-1.
  14. ^ Volpert, AI , Khudyaev, SI (1985), Análisis en clases de funciones discontinuas y ecuaciones de física matemática. Dordrecht, Países Bajos: Nijoff. (Traducción de la 1.ª edición rusa, Moscú, Nauka publ., 1975.)
  15. ^ Schuster, S., Schuster R. (1989). Una generalización de la condición de Wegscheider. Implicaciones para las propiedades de los estados estacionarios y para la aproximación al estado cuasi estacionario. J. Math. Chem, 3 (1), 25–42.
  16. ^ Lorentz H.-A. (1887) Über das Gleichgewicht der lebendigen Kraft unter Gasmolekülen. Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften en Viena. 95 (2), 115-152.
  17. ^ ab Boltzmann L. (1887) Neuer Beweis zweier Sätze über das Wärmegleichgewicht unter mehratomigen Gasmolekülen. Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften en Viena. 95 (2), 153-164.
  18. ^ Shannon, CE (1948) Una teoría matemática de la comunicación, Bell System Technical Journal , vol. 27, págs. 379–423, 623–656. [1] [2]
  19. ^ Hugh Everett, Teoría de la función de onda universal, Tesis, Universidad de Princeton, (1956, 1973), Apéndice I, págs. 121 y siguientes. En su tesis, Everett utilizó el término "balance detallado" de manera poco convencional, en lugar de ecuación de balance.
  20. ^ Cercignani, C. y Lampis, M. (1981). Sobre el teorema H para gases poliatómicos, Journal of Statistical Physics, V. 26 (4), 795–801.
  21. ^ abc Gorban, AN , Yablonsky, GS (2011) Balance detallado extendido para sistemas con reacciones irreversibles, Chemical Engineering Science 66, 5388–5399.
  22. ^ Colquhoun, D., Dowsland, KA, Beato, M. y Plested, AJR (2004) Cómo imponer reversibilidad microscópica en mecanismos de reacción complejos, Biophysical Journal 86, junio de 2004, 3510–3518
  23. ^ ab Horn, F., Jackson, R. (1972) Cinética general de acción de masas. Arch. Ration. Mech. Anal. 47, 87–116.
  24. ^ Stueckelberg , ECG (1952) Teorema de H et unitarité de S. Helv. Phys. Acta 25, 577–-580
  25. ^ abc Gorban, AN, Shahzad, M. (2011) El teorema de Michaelis–Menten–Stueckelberg. Entropy 13, n.º 5, 966–1019.
  26. ^ Mirkes, Evgeny M. (2020). "Entropías universales de Gorban: estudio de caso geométrico". Entropía . 22 (3): 264. arXiv : 2004.14249 . Código Bibliográfico :2020Entrp..22..264M. doi : 10.3390/e22030264 . PMC 7516716 . PMID  33286038. 
  27. ^ Chu, Ch. (1971), Absorción de gas acompañada de un sistema de reacciones de primer orden, Chem. Eng. Sci. 26(3), 305–312.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Detailed_balance&oldid=1243791118"