Loading [MathJax]/jax/output/HTML-CSS/jax.js
Saltar al contenido principal
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
LibreTexts Español

3.3: Conjunto canónico

( \newcommand{\kernel}{\mathrm{null}\,}\)

La termodinámica de equilibrio describe sistemas que están en equilibrio térmico. En una imagen de conjunto, esto se puede considerar asumiendo que el sistema está en contacto con un baño de calor muy grande, para fines matemáticos infinitamente grande. Debido a esto, los sistemas individuales en el conjunto pueden diferir en energía. Sin embargo, la distribución de densidad de probabilidad en el espacio de fase o espacio de estado debe ser consistente con la temperatura constanteT, que es la temperatura del baño de calor. En experimentos, es la temperatura del ambiente.

Concepto3.3.1: Canonical Ensemble

Un conjunto con un número constanteN de partículas en un volumen constanteV y en equilibrio térmico con un baño de calor a temperatura constanteT puede considerarse como un conjunto de subconjuntos microcanónicos con diferentes energíasϵi. La dependencia energética de la densidad de probabilidad se ajusta a la distribución de Boltzmann. Tal conjunto se llama conjunto canónico.

Nota

Debido a que cada sistema puede intercambiar calor con el baño y así cambiar su energía, los sistemas se transferirán entre subconjuntos durante la evolución. Esto no invalida la idea de subconjuntos microcanónicos con números de partículas constantesNi. Para un conjunto suficientemente grande en equilibrio térmicoNi son constantes de movimiento.

Existen diferentes formas de derivar la distribución de Boltzmann. La mayoría de ellos son bastante abstractos y se basan en un gran aparato matemático. La derivación se alarga si se quiere crear la ilusión de que sabemos por qué la constanteβ introducida a continuación siempre es igual1/kBT, dondekB=R/NAv está la constante de Boltzmann, que a su vez es la relación de la constante de gas universalR y la constante de Avogadro NAv. Aquí seguimos una derivación que es físicamente transparente y se basa en un mínimo de aparato matemático que ya hemos introducido.

Distribución Boltzmann

Aquí nos desviamos de la imagen de conjunto y utilizamos un sistema deN partículas que pueden existir enr diferentes estados con energíasϵi coni=0r1. El número de partículas con energíaϵi esNi. Las partículas no interactúan, son completamente independientes entre sí. Por lo tanto, podríamos asociar estas partículas con subconjuntos microcanónicos de un conjunto canónico, pero la situación es más fácil de imaginar con partículas. La probabilidad dePi=Ni/N encontrar una partícula con energíaϵi puede asociarse con la densidad de probabilidad para el subconjunto microcanónico en energíaϵi. La diferencia entre esta derivación simple y la derivación más elaborada para un conjunto canónico es, por lo tanto, esencialmente la diferencia entre la teoría de probabilidad discreta y continua. Asumimos además que las partículas son partículas clásicas y, por lo tanto, distinguibles.

Para calcular la distribución de probabilidadPi=Ni/N, observamos que

r10Ni=N

y

r10Niϵi=E ,

dondeE es una energía total constante del sistema. Hay que tener cuidado al interpretar esta última ecuación en la imagen de conjunto. La cantidadE corresponde a la energía de todo el conjunto canónico, que en efecto es una constante de movimiento, si consideramos un número suficientemente grande de sistemas en contacto con un baño termal. Así podemos usar nuestro modelo simple deN partículas para adivinar la distribución de densidad de probabilidad en el conjunto canónico.

Lo que estamos buscando es la distribución más probable de lasN partículas en los niveles der energía. Esto equivale a poner bolasN distinguibles enr cajas. Ya resolvimos el problema de distribuirN objetos a 2 estados al considerar la distribución binomial en Sección [binomial_distribución]. El peso estadístico de una configuración conn objetos en el primer estado yNn objetos en el segundo estado fue(Nn). Con esta información ya podríamos resolver el problema de un conjunto canónico deN espinesS=1/2 en contacto térmico con el entorno, sin tener en cuenta por el momento las diferencias entre la estadística clásica y la cuántica (ver Sección [sección:quantum_statistics]).

Volviendo aN partículas y niveles der energía, todavía tenemosN! permutaciones. Si asignamos las primerasN0 partículas al estado con energíaϵ0, las siguientesN1 partículas aϵ1 y así sucesivamente, necesitamos dividir cada vez por el número de permutacionesNi! en el mismo estado energético, porque la secuencia de partículas con la misma energía no importa. Llamamos al vector de los números de ocupaciónNi una configuración. La configuración especifica un macroestado particular del sistema y la probabilidad relativa de los macroestados para partículas distinguibles y estados no degenerados viene dada por sus pesos estadísticos,

Ω=N!N0!N1!Nr1! .

El caso con niveles de energía degenerada se trata en la Sección [sec:Maxwell-Boltzmann].

El macroestado más probable es el que tiene el peso estadístico máximoΩ. Debido al pico de distribuciones de probabilidad para grandesN, necesitamos calcular solo este macroestado más probable; es representativo para todo el conjunto. En lugar de maximizar tambiénΩ podemos maximizarlnΩ, ya que el logaritmo natural es una función estrictamente monótona. Esto nos permite aplicar la fórmula de Stirling,

lnΩ=lnN!r1i=0lnNi!NlnNN+1r1i=0NilnNi+r10Nir .

Al insertar Ecuación\ ref {EQ:Conservation_n} encontramos

lnΩNlnNr1i=0NilnNi+1r .

Obsérvese que el segundo término en el lado derecho de la Ecuación\ ref {EQ:LN_Omega} tiene cierta similitud con la entropía de la mezcla, lo que sugiere quelnΩ está relacionado con la entropía.

Al máximo delnΩ la derivada delnΩ respecto alNi mosto de desaparecer,

0=δiNilnNi=i(NiδlnNi+δNilnNi)=iδNi+ilnNiδNi .

Además, debemos considerar las condiciones límite del número constante de partículas, Ecuación\ ref {EQ:Conservation_n},

δN=iδNi=0

y energía total constante, Ecuación\ ref {EQ:Conservación_e},

δE=iϵiδNi=0 .

Podría parecer que la Ecuación\ ref {EQ:Conservation_N_diff} podría usarse para cancelar un término en la Ecuación\ ref {EQ:Max_LN_Omega}, pero esto estaría mal ya que la Ecuación\ ref {EQ:Conservation_N_diff} es una restricción que debe cumplirse por separado. Para la maximización restringida podemos utilizar el método de multiplicadores Lagrange.

El máximo o mínimo de una funciónf(x1,xn) den variables es un punto estacionario que se alcanza en

δf=ni=1(fxi)xkxixi=0 .

Consideramos ahora el caso en el que los posibles conjuntos de lasn variables están restringidos por ecuacionesc adicionales

gj(x1,x2,,xn)=0 ,

donde el índicej recorre lasc restricciones (j=1c). Cada restricción introduce otra ecuación de la misma forma que la de Ecuación\ ref {eq:extremum_multi},

δgj=ni=1(gjxi)xkxixi=0 .

Las restricciones se pueden introducir multiplicando cada una de lasc ecuaciones por un multiplicadorλj y restándola de la ecuación para el punto estacionario sin las restricciones,

δL=ni=1[(fxi)xkxicj=1λj(gjxi)xkxi]xi .

Si un conjunto de variables{x0,1,x0,n} resuelve el problema restringido entonces existe un conjunto{λ0,1λ0,r} para el cual{x0,1,x0,2,,x0,n} también corresponde a un punto estacionario de la función LagrangianaL(x1,,xn,λ1,λr). Tenga en cuenta que no todos los puntos estacionarios de la función lagrangiana son necesariamente soluciones del problema restringido. Esto debe verificarse por separado. [Concepto:Lagrangian_multiplicadores]

Con este método, podemos escribir

0=iδNi+ilnNiδNi+αiδNi+βiϵiδNi=iδNi(1+lnNi+α+βϵi) .

Las dos condiciones límite fijan sólo dos de los números de poblaciónNi. Podemos elegir los multiplicadoresα y de unaβ manera que(1+lnNi+α+βϵi)=0 para estos dosNi, lo que asegura que las derivadas parciales delnΩ con respecto a estos dosNi desaparezcan. Los otros números der2 población pueden, en principio, ser elegidos libremente, pero nuevamente debemos tener

1+lnNi+α+βϵi=0

para que todosi nos aseguremos de que encontremos un máximo con respecto a la variación de cualquiera de los números der población. Esto da

Ni=γeβϵi

conγ=e(1+α). Podemos eliminarγ usando la Ecuación\ ref {EQ:Conservation_n},

iNi=γieβϵi=N ,

dando

γ=Nieβϵi ,

y finalmente llevando a

Pi=NiN=eβϵiieβϵi .

Para muchos problemas en la termodinámica estadística, el multiplicador de Lagrangeα está relacionado con el potencial químico porα=μ/(kBT). El multiplicador Lagrangeβ debe tener la dimensión recíproca de una energía, ya que el exponente debe ser adimensional. Como se indicó anteriormente, no podemos en esta etapa probar queβ es la misma energía para todos los problemas del tipo que aquí nos hemos planteado, y mucho menos para todos los problemas análogos de los conjuntos canónicos. Todo el formalismo se puede conectar a la termodinámica fenomenológica a través de la teoría cinética del gas de Maxwell (ver también Sección [subsección:equipartición]). Para este problema uno encuentra

β=1kBT .

Concepto3.3.2: Boltzmann Distribution

Para un conjunto canónico clásico con niveles de energía,ϵi la distribución de probabilidad para las poblaciones de nivel viene dada por la distribución de Boltzmann

Pi=NiN=eϵi/kBTieϵi/kBT .

La suma sobre los estados

Z(N,V,T)=ieϵi/kBT

requerido para la normalización se llama función de partición canónica. 10 La función de partición es una función de estado termodinámico.

Para la función de partición, utilizamos el símboloZ relativo al término alemán Zustandssumme (“sum over states”), que es una descripción más lúcida de esta cantidad.

Teorema de Equipartición

La comparación de la teoría cinética de gases de Maxwell con la ecuación de estado del gas ideal de la termodinámica fenomenológica proporciona una energía cinética media de una partícula puntual deϵkin=3kBT/2. Esta energía corresponde a

ϵtrans=12mv2=12mp2 ,

es decir, es cuadrática en las coordenadas de velocidad del espacio dinámico o las coordenadas de momento del espacio de fase. La energía traslacional se distribuye a través de tres grados de libertad, ya que las velocidades o momentos tienen componentes a lo largo de tres direcciones ortogonales por pares en el espacio. Cada grado cuadrático de libertad aporta así una energía media dekBT/2.

Si aceptamos que el multiplicador Lagrangeβ asume un valor1/kBT, encontramos una energía mediakBT de un oscilador armónico en el límite de alta temperatura. Dicho oscilador tiene dos grados de libertad que contribuyen cuadráticamente en los grados de libertad a la energía,

ϵvib=12μv2+12fx2 ,

dondeμ está la masa reducida yf la constante de fuerza. El primer término contribuye a la energía cinética, el segundo a la energía potencial. En el promedio de tiempo, cada término aporta la misma energía y asumiendo ergodicidad esto significa que cada uno de los dos grados de libertadkBT/2 contribuye con la energía promedio de un sistema en equilibrio térmico.

El mismo ejercicio se puede realizar para grados de libertad rotacional con energía

ϵrot=12Iω2 ,

dondeI es el momento angular y la frecuenciaω angular. Cada grado de libertad rotacional, siendo cuadrático en deω nuevo aporta una energía media dekBT/2.

Basado en la ecuación\ ref {eq:boltzmann_distribution_0} se puede demostrar que para una energía

ϵi=η0+η1+η2+=fk=1ηk ,

donde índicek recorre los grados individuales de libertad, el número de moléculas que aportan energíaηk no depende de los términosηj conjk. Se puede demostrar además que

ηk=12β

para todos los términos que contribuyen cuadráticamente a la energía. 11

Este resultado tiene dos consecuencias. Primero, podemos generalizarβ=1/kBT, que conocíamos estrictamente solo por grados de libertad traslacionales, a cualquier conjunto canónico para el que todas las contribuciones individuales de energía sean cuadráticas a lo largo de una dimensión en el espacio de fases. Segundo, podemos formular el

Cada grado de libertad, cuya energía se escala cuadráticamente con una de las coordenadas del espacio estatal, aporta una energía media dekBT/2.

El teorema de equipartición se aplica a todos los grados de libertad que se activan. Los grados de libertad traslacionales siempre se activan y los grados de libertad rotacional se activan a temperatura ambiente, lo que corresponde al límite de alta temperatura de la dinámica rotacional. A los grados de libertad vibratoria el teorema de equipartición se aplica sólo en el límite de alta temperatura. En general, el teorema de equipartición falla para grados de libertad cuantificados si el espaciamiento de energía cuántica es comparablekBT/2 o excede este valor. Volveremos a este punto al discutir la función de partición vibracional.

Energía interna y capacidad calorífica del conjunto canónico

La energía internau de un sistema que consiste enN partículas que se distribuyen a niveles der energía puede identificarse como la energía totalE del sistema considerado en la Sección ([Subsección:Boltzmann]). Usando Eqs. \ ref {eq:conservation_e} y\ ref {eq:boltzmann_distribution} encontramos

u=Niϵieϵi/kBTieϵi/kBT=Niϵieϵi/kBTZ .

La suma en el numerador puede ser expresada por la función de partición, ya que

dZdT=1kBT2iϵieϵi/kBT .

Así obtenemos

u=NkBT21ZdZdT=NkBT2dlnZdT .

Nuevamente se mantiene la analogía de nuestro sistema sencillo con el conjunto canónico. En este punto hemos calculado una de las funciones estatales de la termodinámica fenomenológica a partir del conjunto de niveles de energía. La derivación de la distribución de Boltzmann también ha indicado quelnΩ, y así las funciones de partición probablementeZ estén relacionadas con la entropía. Veremos en la Sección [section:state_fct_partition_fct] que este es efectivamente el caso y que podemos calcular todas las funciones de estado termodinámico a partir deZ.

Aquí todavía podemos derivar la capacidad caloríficacV a volumen constante, que es la derivada parcial de la energía interna con respecto a la temperatura. Para ello observamos que la función de partición para el conjunto canónico se relaciona con el volumen constante y el número constante de partículas.

cV=(uT)V=NT(kBT2lnZT)V=NT(kBlnZ1/T)V=NkB([lnZ/1/T]T)V=NkBT2([lnZ/1/T]1/T)V=kBT2(2lnz(1/T)2)V .

En la última línea de Ecuación\ ref {eq:cv} hemos sustituido la función de partición molecularZ por la función de partición para todo el sistema,lnz=NlnZ. Tenga en cuenta que esto implica una generalización. Antes, estábamos considerando un sistema de partículasN idénticas. Ahora asumimos implícitamente que la Ecuación\ ref {eq:cv}, así comou=kBT2dlnzdT se mantendrá para cualquier sistema, siempre y cuando derivemos correctamente la función de partición del sistemaz.

Observamos aquí que el conjunto canónico describe un sistema cerrado que puede intercambiar calor con su entorno, pero por definición no puede intercambiar trabajo, porque su volumenV es constante. Esto no presenta ningún problema, ya que las funciones de estado se pueden computar en diferentesV. En particular, la presión también sep puede calcular a partir de la función de partición (ver Sección [section:state_fct_partition_fct]). Sin embargo, debido a que el conjunto canónico está cerrado, no se puede aplicar fácilmente a todos los problemas que involucran reacciones químicas. Para ello necesitamos eliminar la restricción de un número constante de partículas en los sistemas que conforman el conjunto.


This page titled 3.3: Conjunto canónico is shared under a CC BY-NC 3.0 license and was authored, remixed, and/or curated by Gunnar Jeschke via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?