Distribución de Maxwell-Boltzmann

Función de distribución de probabilidad específica, importante en física
Maxwell-Boltzmann
Función de densidad de probabilidad
Función de distribución acumulativa
Parámetros a > 0 {\displaystyle a>0}
Apoyo incógnita ( 0 ; ) {\displaystyle x\in (0;\infty )}
PDF

2 π incógnita 2 a 3 exp ( incógnita 2 2 a 2 ) {\displaystyle {\sqrt {\frac {2}{\pi }}}\,{\frac {x^{2}}{a^{3}}}\,\exp \left({\frac {-x^{2}}{2a^{2}}}\right)}

(donde exp es la función exponencial )
CDF

campo ( incógnita 2 a ) 2 π incógnita a exp ( incógnita 2 2 a 2 ) {\displaystyle \operatorname {erf} \left({\frac {x}{{\sqrt {2}}a}}\right)-{\sqrt {\frac {2}{\pi }}}\,{\frac {x}{a}}\,\exp \left({\frac {-x^{2}}{2a^{2}}}\right)}

(donde erf es la función de error )
Significar micras = 2 a 2 π {\displaystyle \mu = 2a{\sqrt {\frac {2}{\pi }}}}
Modo 2 a {\displaystyle {\sqrt {2}}a}
Diferencia σ 2 = a 2 ( 3 π 8 ) π {\displaystyle \sigma ^{2}={\frac {a^{2}(3\pi -8)}{\pi }}}
Oblicuidad gamma 1 = 2 2 ( 16 5 π ) ( 3 π 8 ) 3 / 2 {\displaystyle \gamma _{1}={\frac {2{\sqrt {2}}(16-5\pi )}{(3\pi -8)^{3/2}}}}
Exceso de curtosis gamma 2 = 4 ( 96 + 40 π 3 π 2 ) ( 3 π 8 ) 2 {\displaystyle \gamma _{2}={\frac {4(-96+40\pi -3\pi ^{2})}{(3\pi -8)^{2}}}}
Entropía En ( a 2 π ) + gamma 1 2 {\displaystyle \ln \left(a{\sqrt {2\pi }}\right)+\gamma -{\frac {1}{2}}}

En física (en particular en mecánica estadística ), la distribución de Maxwell-Boltzmann , o distribución Maxwell(iana) , es una distribución de probabilidad particular que lleva el nombre de James Clerk Maxwell y Ludwig Boltzmann .

Se definió y utilizó por primera vez para describir las velocidades de partículas en gases idealizados , donde las partículas se mueven libremente dentro de un recipiente estacionario sin interactuar entre sí, excepto por colisiones muy breves en las que intercambian energía y momento entre sí o con su entorno térmico. El término "partícula" en este contexto se refiere solo a partículas gaseosas ( átomos o moléculas ), y se supone que el sistema de partículas ha alcanzado el equilibrio termodinámico . [1] Las energías de tales partículas siguen lo que se conoce como estadísticas de Maxwell-Boltzmann , y la distribución estadística de velocidades se deriva al igualar las energías de las partículas con la energía cinética .

Matemáticamente, la distribución de Maxwell-Boltzmann es la distribución chi con tres grados de libertad (los componentes del vector de velocidad en el espacio euclidiano ), con un parámetro de escala que mide las velocidades en unidades proporcionales a la raíz cuadrada de (la relación entre la temperatura y la masa de la partícula). [2] yo / metro {\estilo de visualización T/m}

La distribución de Maxwell-Boltzmann es un resultado de la teoría cinética de los gases , que proporciona una explicación simplificada de muchas propiedades gaseosas fundamentales, incluyendo la presión y la difusión . [3] La distribución de Maxwell-Boltzmann se aplica fundamentalmente a las velocidades de partículas en tres dimensiones, pero resulta que depende solo de la velocidad (la magnitud de la velocidad) de las partículas. Una distribución de probabilidad de velocidad de partículas indica qué velocidades son más probables: una partícula elegida al azar tendrá una velocidad seleccionada aleatoriamente de la distribución, y es más probable que esté dentro de un rango de velocidades que otro. La teoría cinética de los gases se aplica al gas ideal clásico , que es una idealización de los gases reales. En los gases reales, hay varios efectos (por ejemplo, interacciones de van der Waals , flujo vorticial , límites de velocidad relativistas e interacciones de intercambio cuántico ) que pueden hacer que su distribución de velocidad sea diferente de la forma de Maxwell-Boltzmann. Sin embargo, los gases enrarecidos a temperaturas ordinarias se comportan casi como un gas ideal y la distribución de velocidad de Maxwell es una excelente aproximación para tales gases. Esto también es cierto para los plasmas ideales , que son gases ionizados de densidad suficientemente baja. [4]

La distribución fue derivada por primera vez por Maxwell en 1860 sobre bases heurísticas . [5] [6] Más tarde, en la década de 1870, Boltzmann llevó a cabo importantes investigaciones sobre los orígenes físicos de esta distribución. La distribución puede derivarse sobre la base de que maximiza la entropía del sistema. A continuación, se presenta una lista de derivaciones:

  1. Distribución de probabilidad de máxima entropía en el espacio de fases, con la restricción de conservación de la energía media yo = mi ; {\displaystyle \langle H\rangle =E;}
  2. Conjunto canónico .

Función de distribución

Para un sistema que contiene una gran cantidad de partículas clásicas idénticas, no interactuantes y no relativistas, en equilibrio termodinámico, la fracción de las partículas dentro de un elemento infinitesimal del espacio de velocidad tridimensional d  3 v , centrado en un vector de velocidad de magnitud , está dada por donde: en {\displaystyle \mathbf {v}} en {\estilo de visualización v} F ( en )   d 3 en = [ metro 2 π a B yo ] 3 / 2 exp ( metro en 2 2 a B yo )   d 3 en , {\displaystyle f(\mathbf {v} )~d^{3}\mathbf {v} ={\biggl [}{\frac {m}{2\pi k_{\text{B}}T}}{\biggr ]}^{{3}/{2}}\,\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right)~d^{3}\mathbf {v} ,}

  • m es la masa de la partícula;
  • k B es la constante de Boltzmann ;
  • T es la temperatura termodinámica ;
  • F ( en ) {\displaystyle f(\mathbf {v} )} es una función de distribución de probabilidad, adecuadamente normalizada de modo que todas las velocidades sean la unidad. F ( en ) d 3 en {\textstyle \int f(\mathbf {v} )\,d^{3}\mathbf {v} }
Funciones de densidad de probabilidad de velocidad de las velocidades de algunos gases nobles a una temperatura de 298,15 K (25 °C). El eje y está en s/m, de modo que el área bajo cualquier sección de la curva (que representa la probabilidad de que la velocidad se encuentre en ese rango) es adimensional.

Se puede escribir el elemento del espacio de velocidad como , para velocidades en un sistema de coordenadas cartesianas estándar, o como en un sistema de coordenadas esféricas estándar, donde es un elemento de ángulo sólido y . d 3 en = d en incógnita d en y d en el {\displaystyle d^{3}\mathbf {v} =dv_{x}\,dv_{y}\,dv_{z}} d 3 en = en 2 d en d Ohmio {\displaystyle d^{3}\mathbf {v} =v^{2}\,dv\,d\Omega } d Ohmio = pecado en θ   d en ϕ   d en θ {\displaystyle d\Omega =\sin {v_{\theta }}~dv_{\phi }~dv_{\theta }} en 2 = | en | 2 = en incógnita 2 + en y 2 + en el 2 {\textstyle v^{2}=|\mathbf {v} |^{2}=v_{x}^{2}+v_{y}^{2}+v_{z}^{2}}

La función de distribución maxwelliana para partículas que se mueven en una sola dirección, si esta dirección es x , es la que se puede obtener integrando la forma tridimensional dada anteriormente sobre v y y v z . F ( en incógnita )   d en incógnita = metro 2 π a B yo exp ( metro en incógnita 2 2 a B yo )   d en incógnita , {\displaystyle f(v_{x})~dv_{x}={\sqrt {\frac {m}{2\pi k_{\text{B}}T}}}\,\exp \left(-{\frac {mv_{x}^{2}}{2k_{\text{B}}T}}\right)~dv_{x},}

Reconociendo la simetría de , se puede integrar sobre un ángulo sólido y escribir una distribución de probabilidad de velocidades como la función [7] F ( en ) {\estilo de visualización f(v)}

F ( en ) = [ metro 2 π a B yo ] 3 / 2 4 π en 2 exp ( metro en 2 2 a B yo ) . {\displaystyle f(v)={\biggl [}{\frac {m}{2\pi k_{\text{B}}T}}{\biggr ]}^{{3}/{2}}\,4\pi v^{2}\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right).}

Esta función de densidad de probabilidad proporciona la probabilidad, por unidad de velocidad, de encontrar la partícula con una velocidad cercana a v . Esta ecuación es simplemente la distribución de Maxwell-Boltzmann (dada en el cuadro de información) con parámetro de distribución La distribución de Maxwell-Boltzmann es equivalente a la distribución chi con tres grados de libertad y parámetro de escala a = a B yo / metro . {\textstyle a={\sqrt {k_{\text{B}}T/m}}\,.} a = a B yo / metro . {\textstyle a={\sqrt {k_{\text{B}}T/m}}\,.}

La ecuación diferencial ordinaria más simple que satisface la distribución es: 0 = a B yo en F " ( en ) + F ( en ) ( metro en 2 2 a B yo ) , F ( 1 ) = 2 π [ metro a B yo ] 3 / 2 exp ( metro 2 a B yo ) ; {\displaystyle {\begin{aligned}0&=k_{\text{B}}Tvf'(v)+f(v)\left(mv^{2}-2k_{\text{B}}T\right),\\[4pt]f(1)&={\sqrt {\frac {2}{\pi }}}\,{\biggl [}{\frac {m}{k_{\text{B}}T}}{\biggr ]}^{3/2}\exp \left(-{\frac {m}{2k_{\text{B}}T}}\right);\end{aligned}}}

o en presentación sin unidades : Con el método de valores medios de Darwin-Fowler , se obtiene como resultado exacto la distribución de Maxwell-Boltzmann. 0 = a 2 incógnita F " ( incógnita ) + ( incógnita 2 2 a 2 ) F ( incógnita ) , F ( 1 ) = 1 a 3 2 π exp ( 1 2 a 2 ) . {\displaystyle {\begin{aligned}0&=a^{2}xf'(x)+\left(x^{2}-2a^{2}\right)f(x),\\[4pt]f(1)&={\frac {1}{a^{3}}}{\sqrt {\frac {2}{\pi }}}\exp \left(-{\frac {1}{2a^{2}}}\right).\end{aligned}}}

Simulación de un gas 2D que se relaja hacia una distribución de velocidad de Maxwell-Boltzmann

Relajación de la distribución Maxwell-Boltzmann 2D

Para partículas confinadas a moverse en un plano, la distribución de velocidad está dada por

PAG ( s < | en | < s + d s ) = metro s a B yo exp ( metro s 2 2 a B yo ) d s {\displaystyle P(s<|\mathbf {v} |<s+ds)={\frac {ms}{k_{\text{B}}T}}\exp \left(-{\frac {ms^{2}}{2k_{\text{B}}T}}\right)ds}

Esta distribución se utiliza para describir sistemas en equilibrio. Sin embargo, la mayoría de los sistemas no comienzan en su estado de equilibrio. La evolución de un sistema hacia su estado de equilibrio está regida por la ecuación de Boltzmann . La ecuación predice que para interacciones de corto alcance, la distribución de velocidad de equilibrio seguirá una distribución de Maxwell-Boltzmann. A la derecha hay una simulación de dinámica molecular (MD) en la que 900  partículas esféricas duras están restringidas a moverse en un rectángulo. Interactúan a través de colisiones perfectamente elásticas . El sistema se inicializa fuera del equilibrio, pero la distribución de velocidad (en azul) converge rápidamente a la distribución de Maxwell-Boltzmann 2D (en naranja).

Velocidades típicas

Distribución de Maxwell-Boltzmann en la atmósfera solar.
Distribución de Maxwell-Boltzmann correspondiente a la atmósfera solar. Las masas de las partículas son la masa de un protón , m p =1,67 × 10 −27  kgDa , y la temperatura es la temperatura efectiva de la fotosfera del Sol , T = 5800 K. ,,y V rms marcan las velocidades más probables, medias y de raíz cuadrada media, respectivamente. Sus valores son≈ V ~ {\displaystyle {\tilde {V}}} V ¯ {\estilo de visualización {\bar {V}}} V ~ {\displaystyle {\tilde {V}}} 9,79 km/s , ≈ V ¯ {\estilo de visualización {\bar {V}}} 11,05 km/s y V rms12,00 kilómetros por segundo .

La velocidad media , la velocidad más probable ( moda ) v p y la velocidad cuadrática media se pueden obtener a partir de las propiedades de la distribución de Maxwell. en {\displaystyle \langle v\rangle } en 2 {\textstyle {\sqrt {\langle v^{2}\rangle }}}

Esto funciona bien para gases monoatómicos casi ideales como el helio , pero también para gases moleculares como el oxígeno diatómico . Esto se debe a que a pesar de la mayor capacidad térmica (mayor energía interna a la misma temperatura) debido a su mayor número de grados de libertad , su energía cinética traslacional (y, por lo tanto, su velocidad) no cambia. [8]

  • La velocidad más probable, v p , es la velocidad que es más probable que posea cualquier molécula (de la misma masa m ) en el sistema y corresponde al valor máximo o la moda de f ( v ) . Para encontrarla, calculamos la derivada ⁠ ⁠ d F d en , {\displaystyle {\frac {df}{dv}},} la fijamos a cero y resolvemos para v : con la solución: donde: d F ( en ) d en = 8 π [ metro 2 π a B yo ] 3 / 2 en [ metro en 2 2 a B yo 1 ] exp ( metro en 2 2 a B yo ) = 0 {\displaystyle {\frac {df(v)}{dv}}=-8\pi {\biggl [}{\frac {m}{2\pi k_{\text{B}}T}}{\biggr ]}^{3/2}\,v\,\left[{\frac {mv^{2}}{2k_{\text{B}}T}}-1\right]\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right)=0} m v p 2 2 k B T = 1 ; v p = 2 k B T m = 2 R T M {\displaystyle {\frac {mv_{\text{p}}^{2}}{2k_{\text{B}}T}}=1;\quad v_{\text{p}}={\sqrt {\frac {2k_{\text{B}}T}{m}}}={\sqrt {\frac {2RT}{M}}}}
    • R es la constante del gas ;
    • M es la masa molar de la sustancia y, por lo tanto, puede calcularse como un producto de la masa de la partícula, m , y la constante de Avogadro , N A : M = m N A . {\displaystyle M=mN_{\mathrm {A} }.}

Para el nitrógeno diatómico ( N 2 , el componente principal del aire ) [nota 1] a temperatura ambiente (300 K ), esto da v p 2 8.31   J mol 1 K 1   300   K 0.028   kg mol 1 422   m/s . {\displaystyle v_{\text{p}}\approx {\sqrt {\frac {2\cdot 8.31\ {\text{J}}\cdot {\text{mol}}^{-1}{\text{K}}^{-1}\ 300\ {\text{K}}}{0.028\ {\text{kg}}\cdot {\text{mol}}^{-1}}}}\approx 422\ {\text{m/s}}.}

  • La velocidad media es el valor esperado de la distribución de velocidad, estableciendo : b = 1 2 a 2 = m 2 k B T {\textstyle b={\frac {1}{2a^{2}}}={\frac {m}{2k_{\text{B}}T}}} v = 0 v f ( v ) d v = 4 π [ b π ] 3 / 2 0 v 3 e b v 2 d v = 4 π [ b π ] 3 / 2 1 2 b 2 = 4 π b = 8 k B T π m = 8 R T π M = 2 π v p {\displaystyle {\begin{aligned}\langle v\rangle &=\int _{0}^{\infty }v\,f(v)\,dv\\[1ex]&=4\pi \left[{\frac {b}{\pi }}\right]^{3/2}\int _{0}^{\infty }v^{3}e^{-bv^{2}}dv=4\pi \left[{\frac {b}{\pi }}\right]^{3/2}{\frac {1}{2b^{2}}}={\sqrt {\frac {4}{\pi b}}}\\[1.4ex]&={\sqrt {\frac {8k_{\text{B}}T}{\pi m}}}={\sqrt {\frac {8RT}{\pi M}}}={\frac {2}{\sqrt {\pi }}}v_{\text{p}}\end{aligned}}}
  • La velocidad cuadrática media es el momento bruto de segundo orden de la distribución de velocidad. La "raíz cuadrática media de la velocidad" es la raíz cuadrada de la velocidad cuadrática media, que corresponde a la velocidad de una partícula con energía cinética media , configurando : v 2 {\displaystyle \langle v^{2}\rangle } v r m s {\displaystyle v_{\mathrm {rms} }} b = 1 2 a 2 = m 2 k B T {\textstyle b={\frac {1}{2a^{2}}}={\frac {m}{2k_{\text{B}}T}}} v r m s = v 2 = [ 0 v 2 f ( v ) d v ] 1 / 2 = [ 4 π ( b π ) 3 / 2 0 v 4 e b v 2 d v ] 1 / 2 = [ 4 π ( b π ) 3 / 2 3 8 ( π b 5 ) 1 / 2 ] 1 / 2 = 3 2 b = 3 k B T m = 3 R T M = 3 2 v p {\displaystyle {\begin{aligned}v_{\mathrm {rms} }&={\sqrt {\langle v^{2}\rangle }}=\left[\int _{0}^{\infty }v^{2}\,f(v)\,dv\right]^{1/2}\\[1ex]&=\left[4\pi \left({\frac {b}{\pi }}\right)^{3/2}\int _{0}^{\infty }v^{4}e^{-bv^{2}}dv\right]^{1/2}\\[1ex]&=\left[4\pi \left({\frac {b}{\pi }}\right)^{3/2}{\frac {3}{8}}\left({\frac {\pi }{b^{5}}}\right)^{1/2}\right]^{1/2}={\sqrt {\frac {3}{2b}}}\\[1ex]&={\sqrt {\frac {3k_{\text{B}}T}{m}}}={\sqrt {\frac {3RT}{M}}}={\sqrt {\frac {3}{2}}}v_{\text{p}}\end{aligned}}}

En resumen, las velocidades típicas se relacionan de la siguiente manera: v p 88.6 %   v < v < 108.5 %   v v r m s . {\displaystyle v_{\text{p}}\approx 88.6\%\ \langle v\rangle <\langle v\rangle <108.5\%\ \langle v\rangle \approx v_{\mathrm {rms} }.}

La velocidad cuadrática media está directamente relacionada con la velocidad del sonido c en el gas, por donde es el índice adiabático , f es el número de grados de libertad de la molécula de gas individual. Para el ejemplo anterior, el nitrógeno diatómico (que se aproxima al aire ) a c = γ 3   v r m s = f + 2 3 f   v r m s = f + 2 2 f   v p , {\displaystyle c={\sqrt {\frac {\gamma }{3}}}\ v_{\mathrm {rms} }={\sqrt {\frac {f+2}{3f}}}\ v_{\mathrm {rms} }={\sqrt {\frac {f+2}{2f}}}\ v_{\text{p}},} γ = 1 + 2 f {\textstyle \gamma =1+{\frac {2}{f}}} 300 K , [nota 2] y el valor real del aire se puede aproximar utilizando el peso molar promedio del aire ( f = 5 {\displaystyle f=5} c = 7 15 v r m s 68 %   v r m s 84 %   v p 353   m / s , {\displaystyle c={\sqrt {\frac {7}{15}}}v_{\mathrm {rms} }\approx 68\%\ v_{\mathrm {rms} }\approx 84\%\ v_{\text{p}}\approx 353\ \mathrm {m/s} ,} 29 g/mol ), obteniéndose347 m/s a300 K (las correcciones por humedad variable son del orden de 0,1% a 0,6%).

La velocidad relativa promedio donde la distribución de velocidad tridimensional es v rel | v 1 v 2 | = d 3 v 1 d 3 v 2 | v 1 v 2 | f ( v 1 ) f ( v 2 ) = 4 π k B T m = 2 v {\displaystyle {\begin{aligned}v_{\text{rel}}\equiv \langle |\mathbf {v} _{1}-\mathbf {v} _{2}|\rangle &=\int \!d^{3}\mathbf {v} _{1}\,d^{3}\mathbf {v} _{2}\left|\mathbf {v} _{1}-\mathbf {v} _{2}\right|f(\mathbf {v} _{1})f(\mathbf {v} _{2})\\[2pt]&={\frac {4}{\sqrt {\pi }}}{\sqrt {\frac {k_{\text{B}}T}{m}}}={\sqrt {2}}\langle v\rangle \end{aligned}}} f ( v ) [ 2 π k B T m ] 3 / 2 exp ( 1 2 m v 2 k B T ) . {\displaystyle f(\mathbf {v} )\equiv \left[{\frac {2\pi k_{\text{B}}T}{m}}\right]^{-3/2}\exp \left(-{\frac {1}{2}}{\frac {m\mathbf {v} ^{2}}{k_{\text{B}}T}}\right).}

La integral se puede realizar fácilmente cambiando las coordenadas y u = v 1 v 2 {\displaystyle \mathbf {u} =\mathbf {v} _{1}-\mathbf {v} _{2}} U = v 1 + v 2 2 . {\displaystyle \mathbf {U} ={\tfrac {\mathbf {v} _{1}\,+\,\mathbf {v} _{2}}{2}}.}

Limitaciones

La distribución de Maxwell-Boltzmann supone que las velocidades de las partículas individuales son mucho menores que la velocidad de la luz, es decir, que . Para los electrones, la temperatura de los electrones debe ser . T m c 2 k B {\displaystyle T\ll {\frac {mc^{2}}{k_{\text{B}}}}} T e 5.93 × 10 9   K {\displaystyle T_{e}\ll 5.93\times 10^{9}~\mathrm {K} }

Estadísticas de Maxwell-Boltzmann

La derivación original en 1860 por James Clerk Maxwell fue un argumento basado en colisiones moleculares de la teoría cinética de los gases , así como ciertas simetrías en la función de distribución de velocidad; Maxwell también dio un argumento temprano de que estas colisiones moleculares implican una tendencia hacia el equilibrio. [5] [6] [9] Después de Maxwell, Ludwig Boltzmann en 1872 [10] también derivó la distribución sobre bases mecánicas y argumentó que los gases deberían con el tiempo tender hacia esta distribución, debido a las colisiones (ver Teorema H ). Más tarde (1877) [11] derivó la distribución nuevamente bajo el marco de la termodinámica estadística . Las derivaciones en esta sección siguen las líneas de la derivación de Boltzmann de 1877, comenzando con el resultado conocido como estadísticas de Maxwell-Boltzmann (de termodinámica estadística). Las estadísticas de Maxwell-Boltzmann dan el número promedio de partículas encontradas en un microestado de una sola partícula dado . Bajo ciertas suposiciones, el logaritmo de la fracción de partículas en un microestado dado es lineal en la relación entre la energía de ese estado y la temperatura del sistema: existen constantes y tales que, para todo , Las suposiciones de esta ecuación son que las partículas no interactúan y que son clásicas; esto significa que el estado de cada partícula puede considerarse independientemente de los estados de las otras partículas. Además, se supone que las partículas están en equilibrio térmico. [1] [12] k {\displaystyle k} C {\displaystyle C} i {\displaystyle i} log ( N i N ) = 1 k E i T + C . {\displaystyle -\log \left({\frac {N_{i}}{N}}\right)={\frac {1}{k}}\cdot {\frac {E_{i}}{T}}+C.}

Esta relación se puede escribir como una ecuación introduciendo un factor normalizador:

N i N = exp ( E i k B T ) j exp ( E j k B T ) {\displaystyle {\frac {N_{i}}{N}}={\frac {\exp \left(-{\frac {E_{i}}{k_{\text{B}}T}}\right)}{\displaystyle \sum _{j}\exp \left(-{\tfrac {E_{j}}{k_{\text{B}}T}}\right)}}} ( 1 )

dónde:

  • N i es el número esperado de partículas en el microestado de partícula única i ,
  • N es el número total de partículas en el sistema,
  • E i es la energía del microestado i ,
  • La suma sobre el índice j tiene en cuenta todos los microestados,
  • T es la temperatura de equilibrio del sistema,
  • k B es la constante de Boltzmann .

El denominador en la ecuación 1 es un factor normalizador de modo que las proporciones suman la unidad; en otras palabras, es una especie de función de partición (para el sistema de una sola partícula, no la función de partición habitual de todo el sistema). N i : N {\displaystyle N_{i}:N}

Como la velocidad y la rapidez están relacionadas con la energía, la ecuación ( 1 ) se puede utilizar para derivar relaciones entre la temperatura y las velocidades de las partículas de gas. Todo lo que se necesita es descubrir la densidad de microestados en energía, que se determina dividiendo el espacio de momento en regiones de igual tamaño.

Distribución del vector de momento

La energía potencial se considera igual a cero, de modo que toda la energía se encuentra en forma de energía cinética. La relación entre la energía cinética y el momento para partículas masivas no relativistas es

E = p 2 2 m {\displaystyle E={\frac {p^{2}}{2m}}} ( 2 )

donde p 2 es el cuadrado del vector de momento p = [ p x , p y , p z ] . Por lo tanto, podemos reescribir la ecuación ( 1 ) como:

N i N = 1 Z exp ( p i , x 2 + p i , y 2 + p i , z 2 2 m k B T ) {\displaystyle {\frac {N_{i}}{N}}={\frac {1}{Z}}\exp \left(-{\frac {p_{i,x}^{2}+p_{i,y}^{2}+p_{i,z}^{2}}{2mk_{\text{B}}T}}\right)} ( 3 )

dónde:

Esta distribución de N i  : N es proporcional a la función de densidad de probabilidad f p para encontrar una molécula con estos valores de componentes de momento, por lo que:

f p ( p x , p y , p z ) exp ( p x 2 + p y 2 + p z 2 2 m k B T ) {\displaystyle f_{\mathbf {p} }(p_{x},p_{y},p_{z})\propto \exp \left(-{\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{2mk_{\text{B}}T}}\right)} ( 4 )

La constante de normalización se puede determinar reconociendo que la probabilidad de que una molécula tenga cierto momento debe ser 1. Integrando la exponencial en la ecuación 4 sobre todos p x , p y y p z se obtiene un factor de + exp ( p x 2 + p y 2 + p z 2 2 m k B T ) d p x d p y d p z = [ π 2 m k B T ] 3 {\displaystyle \iiint _{-\infty }^{+\infty }\exp \left(-{\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{2mk_{\text{B}}T}}\right)dp_{x}\,dp_{y}\,dp_{z}={\Bigl [}{\sqrt {\pi }}{\sqrt {2mk_{\text{B}}T}}{\Bigr ]}^{3}}

De modo que la función de distribución normalizada es:

f p ( p x , p y , p z ) = [ 1 2 π m k B T ] 3 / 2 exp ( p x 2 + p y 2 + p z 2 2 m k B T ) {\displaystyle f_{\mathbf {p} }(p_{x},p_{y},p_{z})=\left[{\frac {1}{2\pi mk_{\text{B}}T}}\right]^{3/2}\exp \left(-{\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{2mk_{\text{B}}T}}\right)}    ( 6 )

Se observa que la distribución es el producto de tres variables independientes distribuidas normalmente , , y , con varianza . Además, se puede ver que la magnitud del momento se distribuirá como una distribución de Maxwell-Boltzmann, con . La distribución de Maxwell-Boltzmann para el momento (o igualmente para las velocidades) se puede obtener de manera más fundamental utilizando el teorema H en equilibrio dentro del marco de la teoría cinética de los gases . p x {\displaystyle p_{x}} p y {\displaystyle p_{y}} p z {\displaystyle p_{z}} m k B T {\displaystyle mk_{\text{B}}T} a = m k B T {\textstyle a={\sqrt {mk_{\text{B}}T}}}

Distribución de energía

La distribución de energía se encuentra imponente

f E ( E ) d E = f p ( p ) d 3 p , {\displaystyle f_{E}(E)\,dE=f_{p}(\mathbf {p} )\,d^{3}\mathbf {p} ,} ( 7 )

donde es el volumen infinitesimal del espacio de fases de los momentos correspondientes al intervalo de energía dE . Haciendo uso de la simetría esférica de la relación de dispersión de energía-momento esto se puede expresar en términos de dE como d 3 p {\displaystyle d^{3}\mathbf {p} } E = | p | 2 2 m , {\displaystyle E={\tfrac {|\mathbf {p} |^{2}}{2m}},}

d 3 p = 4 π | p | 2 d | p | = 4 π m 2 m E   d E . {\displaystyle d^{3}\mathbf {p} =4\pi |\mathbf {p} |^{2}d|\mathbf {p} |=4\pi m{\sqrt {2mE}}\ dE.} ( 8 )

Usando entonces ( 8 ) en ( 7 ), y expresando todo en términos de la energía E , obtenemos y finalmente f E ( E ) d E = [ 1 2 π m k B T ] 3 / 2 exp ( E k B T ) 4 π m 2 m E   d E = 2 E π [ 1 k B T ] 3 / 2 exp ( E k B T ) d E {\displaystyle {\begin{aligned}f_{E}(E)dE&=\left[{\frac {1}{2\pi mk_{\text{B}}T}}\right]^{3/2}\exp \left(-{\frac {E}{k_{\text{B}}T}}\right)4\pi m{\sqrt {2mE}}\ dE\\[1ex]&=2{\sqrt {\frac {E}{\pi }}}\,\left[{\frac {1}{k_{\text{B}}T}}\right]^{3/2}\exp \left(-{\frac {E}{k_{\text{B}}T}}\right)\,dE\end{aligned}}}

f E ( E ) = 2 E π [ 1 k B T ] 3 / 2 exp ( E k B T ) {\displaystyle f_{E}(E)=2{\sqrt {\frac {E}{\pi }}}\,\left[{\frac {1}{k_{\text{B}}T}}\right]^{3/2}\exp \left(-{\frac {E}{k_{\text{B}}T}}\right)}    ( 9 )

Dado que la energía es proporcional a la suma de los cuadrados de los tres componentes del momento distribuidos normalmente, esta distribución de energía se puede escribir de manera equivalente como una distribución gamma , utilizando un parámetro de forma y un parámetro de escala. k shape = 3 / 2 {\displaystyle k_{\text{shape}}=3/2} θ scale = k B T . {\displaystyle \theta _{\text{scale}}=k_{\text{B}}T.}

Utilizando el teorema de equipartición , dado que la energía se distribuye uniformemente entre los tres grados de libertad en equilibrio, también podemos dividir en un conjunto de distribuciones de chi-cuadrado , donde la energía por grado de libertad, ε se distribuye como una distribución de chi-cuadrado con un grado de libertad, [13] f E ( E ) d E {\displaystyle f_{E}(E)dE} f ε ( ε ) d ε = 1 π ε k B T   exp ( ε k B T ) d ε {\displaystyle f_{\varepsilon }(\varepsilon )\,d\varepsilon ={\sqrt {\frac {1}{\pi \varepsilon k_{\text{B}}T}}}~\exp \left(-{\frac {\varepsilon }{k_{\text{B}}T}}\right)\,d\varepsilon }

En equilibrio, esta distribución se mantendrá para cualquier número de grados de libertad. Por ejemplo, si las partículas son dipolos de masa rígidos con un momento dipolar fijo, tendrán tres grados de libertad de traslación y dos grados de libertad de rotación adicionales. La energía en cada grado de libertad se describirá de acuerdo con la distribución de chi-cuadrado anterior con un grado de libertad, y la energía total se distribuirá de acuerdo con una distribución de chi-cuadrado con cinco grados de libertad. Esto tiene implicaciones en la teoría del calor específico de un gas.

Distribución del vector de velocidad

Reconociendo que la densidad de probabilidad de velocidad f v es proporcional a la función de densidad de probabilidad de momento por

f v d 3 v = f p ( d p d v ) 3 d 3 v {\displaystyle f_{\mathbf {v} }d^{3}\mathbf {v} =f_{\mathbf {p} }\left({\frac {dp}{dv}}\right)^{3}d^{3}\mathbf {v} }

y usando p = m v obtenemos

f v ( v x , v y , v z ) = [ m 2 π k B T ] 3 / 2 exp ( m ( v x 2 + v y 2 + v z 2 ) 2 k B T ) {\displaystyle f_{\mathbf {v} }(v_{x},v_{y},v_{z})={\biggl [}{\frac {m}{2\pi k_{\text{B}}T}}{\biggr ]}^{3/2}\exp \left(-{\frac {m\left(v_{x}^{2}+v_{y}^{2}+v_{z}^{2}\right)}{2k_{\text{B}}T}}\right)}

que es la distribución de velocidad de Maxwell-Boltzmann. La probabilidad de encontrar una partícula con velocidad en el elemento infinitesimal [ dv x , dv y , dv z ] con una velocidad v = [ v x , vy , v z ] es

f v ( v x , v y , v z ) d v x d v y d v z . {\displaystyle f_{\mathbf {v} }{\left(v_{x},v_{y},v_{z}\right)}\,dv_{x}\,dv_{y}\,dv_{z}.}

Al igual que el momento, se observa que esta distribución es el producto de tres variables independientes distribuidas normalmente , , y , pero con varianza . También se puede observar que la distribución de velocidad de Maxwell-Boltzmann para el vector velocidad [ v x , vy , v z ] es el producto de las distribuciones para cada una de las tres direcciones: donde la distribución para una sola dirección es v x {\displaystyle v_{x}} v y {\displaystyle v_{y}} v z {\displaystyle v_{z}} k B T / m {\textstyle k_{\text{B}}T/m} f v ( v x , v y , v z ) = f v ( v x ) f v ( v y ) f v ( v z ) {\displaystyle f_{\mathbf {v} }{\left(v_{x},v_{y},v_{z}\right)}=f_{v}(v_{x})f_{v}(v_{y})f_{v}(v_{z})} f v ( v i ) = m 2 π k B T exp ( m v i 2 2 k B T ) . {\displaystyle f_{v}(v_{i})={\sqrt {\frac {m}{2\pi k_{\text{B}}T}}}\exp \left(-{\frac {mv_{i}^{2}}{2k_{\text{B}}T}}\right).}

Cada componente del vector de velocidad tiene una distribución normal con media y desviación estándar , por lo que el vector tiene una distribución normal tridimensional, un tipo particular de distribución normal multivariada , con media y covarianza , donde es la matriz identidad de 3 × 3 . μ v x = μ v y = μ v z = 0 {\displaystyle \mu _{v_{x}}=\mu _{v_{y}}=\mu _{v_{z}}=0} σ v x = σ v y = σ v z = k B T / m {\textstyle \sigma _{v_{x}}=\sigma _{v_{y}}=\sigma _{v_{z}}={\sqrt {k_{\text{B}}T/m}}} μ v = 0 {\displaystyle \mu _{\mathbf {v} }=\mathbf {0} } Σ v = ( k B T m ) I {\textstyle \Sigma _{\mathbf {v} }=\left({\frac {k_{\text{B}}T}{m}}\right)I} I {\displaystyle I}

Distribución de la velocidad

La distribución de Maxwell-Boltzmann para la velocidad se deduce inmediatamente de la distribución del vector de velocidad, que se muestra arriba. Nótese que la velocidad es y el elemento de volumen en coordenadas esféricas, donde y son los ángulos de coordenadas esféricas del vector de velocidad. La integración de la función de densidad de probabilidad de la velocidad sobre los ángulos sólidos produce un factor adicional de . La distribución de velocidad con la sustitución de la velocidad por la suma de los cuadrados de los componentes del vector: v = v x 2 + v y 2 + v z 2 {\displaystyle v={\sqrt {v_{x}^{2}+v_{y}^{2}+v_{z}^{2}}}} d v x d v y d v z = v 2 sin θ d v d θ d ϕ = v 2 d v d Ω {\displaystyle dv_{x}\,dv_{y}\,dv_{z}=v^{2}\sin \theta \,dv\,d\theta \,d\phi =v^{2}\,dv\,d\Omega } ϕ {\displaystyle \phi } θ {\displaystyle \theta } d Ω {\displaystyle d\Omega } 4 π {\displaystyle 4\pi }

f ( v ) = 2 π [ m k B T ] 3 / 2 v 2 exp ( m v 2 2 k B T ) . {\displaystyle f(v)={\sqrt {\frac {2}{\pi }}}\,{\biggl [}{\frac {m}{k_{\text{B}}T}}{\biggr ]}^{3/2}v^{2}\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right).}

Ennorte-espacio dimensional

En el espacio n -dimensional, la distribución de Maxwell-Boltzmann se convierte en: f ( v )   d n v = [ m 2 π k B T ] n / 2 exp ( m | v | 2 2 k B T )   d n v {\displaystyle f(\mathbf {v} )~d^{n}\mathbf {v} ={\biggl [}{\frac {m}{2\pi k_{\text{B}}T}}{\biggr ]}^{n/2}\exp \left(-{\frac {m|\mathbf {v} |^{2}}{2k_{\text{B}}T}}\right)~d^{n}\mathbf {v} }

La distribución de velocidad se convierte en: donde es una constante normalizadora. f ( v )   d v = A exp ( m v 2 2 k B T ) v n 1   d v {\displaystyle f(v)~dv=A\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right)v^{n-1}~dv} A {\displaystyle A}

El siguiente resultado integral es útil: donde es la función Gamma . Este resultado se puede utilizar para calcular los momentos de la función de distribución de velocidad: que es la velocidad media en sí misma 0 v a exp ( m v 2 2 k B T ) d v = [ 2 k B T m ] a + 1 2 0 e x x a / 2 d x 1 / 2 = [ 2 k B T m ] a + 1 2 0 e x x a / 2 x 1 / 2 2 d x = [ 2 k B T m ] a + 1 2 Γ ( a + 1 2 ) 2 {\displaystyle {\begin{aligned}\int _{0}^{\infty }v^{a}\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right)dv&=\left[{\frac {2k_{\text{B}}T}{m}}\right]^{\frac {a+1}{2}}\int _{0}^{\infty }e^{-x}x^{a/2}\,dx^{1/2}\\[2pt]&=\left[{\frac {2k_{\text{B}}T}{m}}\right]^{\frac {a+1}{2}}\int _{0}^{\infty }e^{-x}x^{a/2}{\frac {x^{-1/2}}{2}}\,dx\\[2pt]&=\left[{\frac {2k_{\text{B}}T}{m}}\right]^{\frac {a+1}{2}}{\frac {\Gamma {\left({\frac {a+1}{2}}\right)}}{2}}\end{aligned}}} Γ ( z ) {\displaystyle \Gamma (z)} v = 0 v v n 1 exp ( m v 2 2 k B T ) d v 0 v n 1 exp ( m v 2 2 k B T ) d v = 2 k B T m     Γ ( n + 1 2 ) Γ ( n 2 ) {\displaystyle \langle v\rangle ={\frac {\displaystyle \int _{0}^{\infty }v\cdot v^{n-1}\exp \left(-{\tfrac {mv^{2}}{2k_{\text{B}}T}}\right)\,dv}{\displaystyle \int _{0}^{\infty }v^{n-1}\exp \left(-{\tfrac {mv^{2}}{2k_{\text{B}}T}}\right)\,dv}}={\sqrt {\frac {2k_{\text{B}}T}{m}}}~~{\frac {\Gamma {\left({\frac {n+1}{2}}\right)}}{\Gamma {\left({\frac {n}{2}}\right)}}}} v a v g = v = 2 k B T m   Γ ( n + 1 2 ) Γ ( n 2 ) . {\textstyle v_{\mathrm {avg} }=\langle v\rangle ={\sqrt {\frac {2k_{\text{B}}T}{m}}}\ {\frac {\Gamma \left({\frac {n+1}{2}}\right)}{\Gamma \left({\frac {n}{2}}\right)}}.}

v 2 = 0 v 2 v n 1 exp ( m v 2 2 k B T ) d v 0 v n 1 exp ( m v 2 2 k B T ) d v = [ 2 k B T m ] Γ ( n + 2 2 ) Γ ( n 2 ) = [ 2 k B T m ] n 2 = n k B T m {\displaystyle {\begin{aligned}\langle v^{2}\rangle &={\frac {\displaystyle \int _{0}^{\infty }v^{2}\cdot v^{n-1}\exp \left(-{\tfrac {mv^{2}}{2k_{\text{B}}T}}\right)\,dv}{\displaystyle \int _{0}^{\infty }v^{n-1}\exp \left(-{\tfrac {mv^{2}}{2k_{\text{B}}T}}\right)\,dv}}\\[1ex]&=\left[{\frac {2k_{\text{B}}T}{m}}\right]{\frac {\Gamma {\left({\frac {n+2}{2}}\right)}}{\Gamma {\left({\frac {n}{2}}\right)}}}\\[1.2ex]&=\left[{\frac {2k_{\text{B}}T}{m}}\right]{\frac {n}{2}}={\frac {nk_{\text{B}}T}{m}}\end{aligned}}} que da la velocidad cuadrática media v rms = v 2 = n k B T m . {\textstyle v_{\text{rms}}={\sqrt {\langle v^{2}\rangle }}={\sqrt {\frac {nk_{\text{B}}T}{m}}}.}

La derivada de la función de distribución de velocidad: d f ( v ) d v = A exp ( m v 2 2 k B T ) [ m v k B T v n 1 + ( n 1 ) v n 2 ] = 0 {\displaystyle {\frac {df(v)}{dv}}=A\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right){\biggl [}-{\frac {mv}{k_{\text{B}}T}}v^{n-1}+(n-1)v^{n-2}{\biggr ]}=0}

Esto produce la velocidad más probable ( modo ) v p = ( n 1 ) k B T / m . {\textstyle v_{\text{p}}={\sqrt {\left(n-1\right)k_{\text{B}}T/m}}.}

Véase también

Notas

  1. ^ El cálculo no se ve afectado por el hecho de que el nitrógeno sea diatómico. A pesar de la mayor capacidad calorífica (mayor energía interna a la misma temperatura) de los gases diatómicos en relación con los gases monoatómicos, debido a su mayor número de grados de libertad , sigue siendo la energía cinética traslacional media . El hecho de que el nitrógeno sea diatómico solo afecta el valor de la masa molar M = 3 R T M m {\displaystyle {\frac {3RT}{M_{\text{m}}}}} 28 g/mol . Véase, por ejemplo, K. Prakashan, Engineering Physics (2001), 2.278.
  2. ^ El nitrógeno a temperatura ambiente se considera un gas diatómico "rígido", con dos grados de libertad rotacionales adicionales a los tres traslacionales, y el grado de libertad vibracional no accesible.

Referencias

  1. ^ ab Mandl, Franz (2008). Física estadística . Física de Manchester (2.ª ed.). Chichester: John Wiley & Sons. ISBN 978-0471915331.
  2. ^ Young, Hugh D.; Friedman, Roger A.; Ford, Albert Lewis; Sears, Francis Weston; Zemansky, Mark Waldo. Física universitaria de Sears y Zemansky: con física moderna (12.ª ed.). San Francisco: Pearson, Addison-Wesley. ISBN 978-0-321-50130-1.
  3. ^ Enciclopedia de física (segunda edición), RG Lerner , GL Trigg, editores VHC, 1991, ISBN 3-527-26954-1 (Verlagsgesellschaft), ISBN 0-89573-752-3 (VHC Inc.)  
  4. ^ NA Krall y AW Trivelpiece, Principles of Plasma Physics, San Francisco Press, Inc., 1986, entre muchos otros textos sobre física básica del plasma.
  5. ^ ab Maxwell, JC (1860 A): Ilustraciones de la teoría dinámica de los gases. Parte I. Sobre los movimientos y colisiones de esferas perfectamente elásticas. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science , 4.ª serie, vol. 19, págs. 19-32. [1]
  6. ^ ab Maxwell, JC (1860 B): Ilustraciones de la teoría dinámica de los gases. Parte II. Sobre el proceso de difusión de dos o más tipos de partículas en movimiento entre sí. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science , 4.ª serie, vol. 20, págs. 21-37. [2]
  7. ^ Müller-Kirsten, HJW (2013). "2". Fundamentos de física estadística (2.ª ed.). World Scientific . ISBN 978-981-4449-53-3.OCLC 822895930  .
  8. ^ Serway, Raymond A.; Faughn, Jerry S. y Vuille, Chris (2011). Física universitaria, volumen 1 (novena edición). Cengage Learning. pág. 352. ISBN 9780840068484.
  9. ^ Gyenis, Balazs (2017). "Maxwell y la distribución normal: una historia coloreada de probabilidad, independencia y tendencia hacia el equilibrio". Estudios de historia y filosofía de la física moderna . 57 : 53–65. arXiv : 1702.01411 . Bibcode :2017SHPMP..57...53G. doi :10.1016/j.shpsb.2017.01.001. S2CID  38272381.
  10. ^ Boltzmann, L., "Weitere studien über das Wärmegleichgewicht unter Gasmolekülen". Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften en Viena, mathematisch-naturwissenschaftliche Classe , 66 , 1872, págs.
  11. ^ Boltzmann, L., "Über die Beziehung zwischen dem zweiten Hauptsatz der mechanischen Wärmetheorie und der Wahrscheinlichkeitsrechnung respektive den Sätzen über das Wärmegleichgewicht". Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften en Viena, Mathematisch-Naturwissenschaftliche Classe . Abt. II, 76 , 1877, págs. 373–435. Reimpreso en Wissenschaftliche Abhandlungen , vol. II, págs. 164–223, Leipzig: Barth, 1909. Traducción disponible en : http://crystal.med.upenn.edu/sharp-lab-pdfs/2015SharpMatschinsky_Boltz1877_Entropy17.pdf Archivado el 5 de marzo de 2021 en Wayback Machine.
  12. ^ Parker, Sybil P. (1993). Enciclopedia de Física McGraw-Hill (2.ª edición). McGraw-Hill. ISBN 978-0-07-051400-3.
  13. ^ Laurendeau, Normand M. (2005). Termodinámica estadística: fundamentos y aplicaciones. Cambridge University Press. pág. 434. ISBN 0-521-84635-8.

Lectura adicional

  • Tipler, Paul Allen; Mosca, Gene (2008). Física para científicos e ingenieros: con física moderna (6.ª ed.). Nueva York: WH Freeman. ISBN 978-0-7167-8964-2.
  • Shavit, Arthur; Gutfinger, Chaim (2009). Termodinámica: de los conceptos a las aplicaciones (2.ª ed.). CRC Press. ISBN 978-1-4200-7368-3.OCLC 244177312  .
  • Ives, David JG (1971). Termodinámica química . Química universitaria. Universidad Macdonald Technical and Scientific. ISBN 0-356-03736-3.
  • Nash, Leonard K. (1974). Elementos de termodinámica estadística . Principios de química (2.ª ed.). Addison-Wesley. ISBN 978-0-201-05229-9.
  • Ward, CA; Fang, G. (1999). "Expresión para predecir el flujo de evaporación de líquidos: enfoque de la teoría de la tasa estadística". Physical Review E . 59 (1): 429–440. doi :10.1103/physreve.59.429. ISSN  1063-651X.
  • Rahimi, P; Ward, CA (2005). "Cinética de la evaporación: enfoque de la teoría de la velocidad estadística". Revista internacional de termodinámica . 8 (9): 1–14.
  • "La distribución de velocidad de Maxwell" del Proyecto de demostraciones de Wolfram en Mathworld
Retrieved from "https://en.wikipedia.org/w/index.php?title=Maxwell–Boltzmann_distribution&oldid=1242061983"