Articulo de referencia

distribución de Maxwell-Boltzmann

a>0 "},"support":{"wt":" x\\in (0;\\infty) "},"pdf":{"wt":" \\sqrt{\\frac{2}{\\pi}} \\, \\frac{x^2}{a^3} \\, \\exp\\left(\\frac{-x^2}{2a^2} \\right) \n(where {{math|exp}} is the...

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

Se definió y utilizó por primera vez para describir las velocidades de las 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 dichas partículas siguen lo que se conoce como estadística de Maxwell-Boltzmann , y la distribución estadística de velocidades se deriva igualando 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 (las 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 deT/metro{\displaystyle T/m}(la relación entre la temperatura y la masa de la partícula). [ 2 ]

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 fundamentales de los gases, incluyendo la presión y la difusión . [ 3 ] La distribución de Maxwell-Boltzmann se aplica fundamentalmente a las velocidades de las partículas en tres dimensiones, pero resulta que depende solo de la rapidez (la magnitud de la velocidad) de las partículas. Una distribución de probabilidad de la velocidad de las partículas indica qué velocidades son más probables: una partícula elegida al azar tendrá una velocidad seleccionada al azar de la distribución, y es más probable que se encuentre dentro de un rango de velocidades que dentro de 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, existen 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 de forma muy similar a un gas ideal y la distribución de velocidad de Maxwell es una excelente aproximación para dichos 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 por razones heurísticas . [ 5 ] [ 6 ] Posteriormente, 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 basándose en que maximiza la entropía del sistema. A continuación se presenta una lista de derivaciones:

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

Función de distribución

Para un sistema que contiene un gran número de partículas clásicas idénticas, no interactuantes y no relativistas en equilibrio termodinámico, la fracción de partículas dentro de un elemento infinitesimal del espacio de velocidad tridimensional d 3 v , centrado en un vector de velocidadv{\displaystyle \mathbf {v} }de magnitudv{\displaystyle v}, se da por F(v) d3v=[metro2πkBT]3/2exp(metrov22kBT) d3v,{\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} ,} dónde:

  • m es la masa de la partícula;
  • k B es la constante de Boltzmann ;
  • T es la temperatura termodinámica ;
  • F(v){\displaystyle f(\mathbf {v} )}es una función de distribución de probabilidad, debidamente normalizada de modo queF(v)d3v{\textstyle \int f(\mathbf {v} )\,d^{3}\mathbf {v} }En todas las velocidades es la unidad.
Funciones de densidad de probabilidad de velocidad 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 velocidades comod3v=dvincógnitadvydvz{\displaystyle d^{3}\mathbf {v} =dv_{x}\,dv_{y}\,dv_{z}}, para velocidades en un sistema de coordenadas cartesianas estándar, o comod3v=v2dvdΩ{\displaystyle d^{3}\mathbf {v} =v^{2}\,dv\,d\Omega }en un sistema de coordenadas esféricas estándar, dondedΩ=pecadoθ dϕ dθ{\displaystyle d\Omega =\sin {\theta }~d\phi ~d\theta }es un elemento de ángulo sólido yv2=|v|2=vincógnita2+vy2+vz2{\textstyle v^{2}=|\mathbf {v} |^{2}=v_{x}^{2}+v_{y}^{2}+v_{z}^{2}}.

Alternativamente, la función de distribución también puede escribirse en el espacio de momentos comoF(pag)d3pag=[12πmetrokBT]3/2exp(pag22metrokBT)d3pag{\displaystyle f(\mathbf {p} )\,d^{3}\mathbf {p} =\left[{\frac {1}{2\pi mk_{\text{B}}T}}\right]^{3/2}\exp \left(-{\frac {p^{2}}{2mk_{\text{B}}T}}\right)\,d^{3}\mathbf {p} }dóndepag=metrov{\displaystyle \mathbf {p} =m\mathbf {v} }es el vector de momento.

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 una distribución normal con una desviación estándar dekBT/metro{\textstyle {\sqrt {k_{\text{B}}T/m}}}: F(vincógnita) dvincógnita=metro2πkBTexp(metrovincógnita22kBT) dvincó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},} que se puede obtener integrando la forma tridimensional dada anteriormente sobre v y y v z .

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

F(v)=[metro2πkBT]3/24πv2exp(metrov22kBT).{\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 da 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óna=kBT/metro.{\textstyle a={\sqrt {k_{\text{B}}T/m}}\,.} La distribución de Maxwell-Boltzmann es equivalente a la distribución chi con tres grados de libertad y parámetro de escala.a=kBT/metro.{\textstyle a={\sqrt {k_{\text{B}}T/m}}\,.}

La ecuación diferencial ordinaria más simple que satisface la distribución es: 0=kBTvF(v)+F(v)(metrov22kBT),F(1)=2π[metrokBT]3/2exp(metro2kBT);{\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 :0=a2incógnitaF(incógnita)+(incógnita22a2)F(incógnita),F(1)=1a32πexp(12a2).{\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}}} Mediante el método de valores medios de Darwin-Fowler , se obtiene la distribución de Maxwell-Boltzmann como resultado exacto.

Simulación de un gas bidimensional que se relaja hacia una distribución de velocidades de Maxwell-Boltzmann.

Relajación a la distribución de Maxwell-Boltzmann bidimensional

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

PAG(s<|v|<s+ds)=metroskBTexp(metros22kBT)ds{\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 se rige 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 se muestra una simulación de dinámica molecular (DM) en la que 900 partículas esféricas rígidas están restringidas a moverse en un rectángulo. Interactúan mediante colisiones perfectamente elásticas . El sistema se inicializa fuera de equilibrio, pero la distribución de velocidad (en azul) converge rápidamente a la distribución de Maxwell-Boltzmann bidimensional (en naranja). 

Velocidades típicas

Distribución de Maxwell-Boltzmann de la atmósfera solar.
La 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  kg1 Da  , y la temperatura es la temperatura efectiva de la fotosfera del Sol , T = 5800 K.V~{\displaystyle {\tilde {V}}},V¯{\displaystyle {\bar {V}}}, y V rms marcan las velocidades más probable, media y cuadrática media, respectivamente. Sus valores sonV~{\displaystyle {\tilde {V}}}9,79  km/s ,V¯{\displaystyle {\bar {V}}}11,05  km/s y V rms12,00  km/s .

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

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 calorífica (mayor energía interna a la misma temperatura) debido a su mayor número de grados de libertad , su energía cinética de traslación (y por lo tanto su velocidad) permanece inalterada. [ 8 ]

  • La velocidad más probable, v p , es la velocidad que tiene más probabilidades de poseer cualquier molécula (de la misma masa m ) en el sistema y corresponde al valor máximo o la moda de f ( v ) . Para hallarla, calculamos la derivada dFdv,{\displaystyle {\tfrac {df}{dv}},} Ponlo a cero y resuelve para v :dF(v)dv=8π[metro2πkBT]3/2v[metrov22kBT1]exp(metrov22kBT)=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}con la solución:metrovpag22kBT=1;vpag=2kBTmetro=2RTMETRO{\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}}}}dónde:
    • R es la constante de los gases ;
    • 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 :METRO=metronorteA.{\displaystyle M=mN_{\mathrm {A} }.}

    Para el nitrógeno diatómico ( N2 , el componente principal del aire ) [ nota 1 ] a temperatura ambiente (300  K ), esto da

    vpag28.31Jmetrool1K1 300K0,028kgramometrool1422metro/s.{\displaystyle v_{\text{p}}\approx {\sqrt {\frac {2\cdot 8.31\,\mathrm {J{\cdot }{mol}^{-1}K^{-1}} \ 300\,\mathrm {K} }{0.028\,\mathrm {{kg}{\cdot }{mol}^{-1}} }}}\approx 422\,\mathrm {m/s} .}
  • La velocidad media es el valor esperado de la distribución de velocidad, estableciendob=12a2=metro2kBT{\textstyle b={\frac {1}{2a^{2}}}={\frac {m}{2k_{\text{B}}T}}}:v=0vF(v)dv=4π[bπ]3/20v3mibv2dv=4π[bπ]3/212b2=4πb=8kBTπmetro=8RTπMETRO=2πvpag{\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}}}\\[1.4ex]&={\sqrt {\frac {4}{\pi b}}}={\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 mediav2{\displaystyle \langle v^{2}\rangle }es el momento bruto de segundo orden de la distribución de velocidad. La "velocidad cuadrática media"vrms{\displaystyle v_{\text{rms}}}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 , estableciendob=12a2=metro2kBT{\textstyle b={\frac {1}{2a^{2}}}={\frac {m}{2k_{\text{B}}T}}}:vrms=v2=[0v2F(v)dv]1/2=[4π(bπ)3/20v4mibv2dv]1/2=[4π(bπ)3/238(πb5)1/2]1/2=32b=3kBTmetro=3RTMETRO=32vpag{\displaystyle {\begin{aligned}v_{\text{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: vpag88.6% v<v<108.5% vvrms.{\displaystyle v_{\text{p}}\approx 88.6\%\ \langle v\rangle <\langle v\rangle <108.5\%\ \langle v\rangle \approx v_{\text{rms}}.}

La velocidad cuadrática media está directamente relacionada con la velocidad del sonido c en el gas, por do=γ3 vrmetros=F+23F vrmetros=F+22F vpag,{\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}},} dóndeγ=1+2F{\textstyle \gamma =1+{\frac {2}{f}}}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, nitrógeno diatómico (que se aproxima al aire ) a300  K ,F=5{\displaystyle f=5}[ nota 2 ] y do=715vrmetros68% vrmetros84% vpag353 metro/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} ,} El valor real del aire se puede aproximar utilizando el peso molar promedio del aire (29  g/mol ), obteniendo347  m/s a300  K (las correcciones por humedad variable son del orden del 0,1% al 0,6%).

La velocidad relativa promedio vrel|v1v2|=d3v1d3v2|v1v2|F(v1)F(v2)=4πkBTmetro=2v{\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}}} donde la distribución de velocidad tridimensional es F(v)[2πkBTmetro]3/2exp(12metrov2kBT).{\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 a coordenadas.=v1v2{\displaystyle \mathbf {u} =\mathbf {v} _{1}-\mathbf {v} _{2}}yU=12(v1+v2).{\textstyle \mathbf {U} ={\tfrac {1}{2}}(\mathbf {v} _{1}+\mathbf {v} _{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, queTmetrodo2kB{\displaystyle T\ll {\frac {mc^{2}}{k_{\text{B}}}}}. Para los electrones, la temperatura de los electrones debe serTmi5.93×109 K{\displaystyle T_{e}\ll 5.93\times 10^{9}~\mathrm {K} }Para la distribución de velocidades de partículas relativistas, véase la distribución de Maxwell-Jüttner .

estadística 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 en fundamentos mecánicos 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 están en las líneas de la derivación de Boltzmann de 1877, comenzando con un resultado conocido como estadística de Maxwell-Boltzmann (de la termodinámica estadística). La estadística de Maxwell-Boltzmann da el número promedio de partículas encontradas en un microestado de partícula única 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: hay constantesk{\displaystyle k}ydo{\displaystyle C}de tal manera que, para todosi{\displaystyle i}, registro(norteinorte)=1kmiiT+do.{\displaystyle -\log \left({\frac {N_{i}}{N}}\right)={\frac {1}{k}}\cdot {\frac {E_{i}}{T}}+C.} 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 demás partículas. Además, se supone que las partículas están en equilibrio térmico. [ 1 ] [ 12 ]

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

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 de normalización de modo que las razonesnortei:norte{\displaystyle N_{i}:N} 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).

Dado que la velocidad y la rapidez están relacionadas con la energía, la ecuación ( 1 ) puede utilizarse para derivar relaciones entre la temperatura y las velocidades de las partículas de gas. Basta con determinar la densidad de microestados en función de la energía, la cual se obtiene dividiendo el espacio de momentos en regiones de igual tamaño.

Distribución para el vector de momento

Se considera que la energía potencial es 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:

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:

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:

La constante de normalización se puede determinar reconociendo que la probabilidad de que una molécula tenga algún momento debe ser 1. Integrando la exponencial en la ecuación 4 sobre todos los p x , p y , y p z se obtiene un factor de +exp(pagincógnita2+pagy2+pagz22metrokBT)dpagincógnitadpagydpagz=[π2metrokBT]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:

Fpag(pagincógnita,pagy,pagz)=[12πmetrokBT]3/2exp(pagincógnita2+pagy2+pagz22metrokBT){\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 con distribución normal.pagincógnita{\displaystyle p_{x}},pagy{\displaystyle p_{y}}, ypagz{\displaystyle p_{z}}, con varianzametrokBT{\displaystyle mk_{\text{B}}T}Además, se puede observar que la magnitud del momento se distribuirá según una distribución de Maxwell-Boltzmann, cona=metrokBT{\textstyle a={\sqrt {mk_{\text{B}}T}}}La distribución de Maxwell-Boltzmann para el momento (o igualmente para las velocidades) se puede obtener de forma más fundamental utilizando el teorema H en equilibrio dentro del marco de la teoría cinética de los gases .

Distribución de la energía

La distribución de energía se encuentra imponiendo

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

Usando entonces ( 8 ) en ( 7 ), y expresando todo en términos de la energía E , obtenemos Fmi(mi)dmi=[12πmetrokBT]3/2exp(mikBT)4πmetro2metromi dmi=2miπ[1kBT]3/2exp(mikBT)dmi{\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}}} y finalmente

Fmi(mi)=2miπ[1kBT]3/2exp(mikBT){\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 de momento distribuidos normalmente, esta distribución de energía se puede escribir de forma equivalente como una distribución gamma , utilizando un parámetro de forma.kforma=3/2{\displaystyle k_{\text{shape}}=3/2}y un parámetro de escala,θescala=kBT.{\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 dividirFmi(mi)dmi{\displaystyle f_{E}(E)\,dE}en un conjunto de distribuciones chi-cuadrado , donde la energía por grado de libertad, ε, se distribuye como una distribución chi-cuadrado con un grado de libertad, [ 13 ]Fε(ε)dε=1πεkBT exp(εkBT)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 mantiene para cualquier número de grados de libertad. Por ejemplo, si las partículas son dipolos de masa rígidos con 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á según la distribución chi-cuadrado mencionada anteriormente con un grado de libertad, y la energía total se distribuirá según una distribución 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

Fvd3v=Fpag(dpagdv)3d3v{\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

Fv(vincógnita,vy,vz)=[metro2πkBT]3/2exp(metro(vincógnita2+vy2+vz2)2kBT){\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 ] alrededor de la velocidad v = [ v x , v y , v z ] es

Fv(vincógnita,vy,vz)dvincógnitadvydvz.{\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 con distribución normal.vincógnita{\displaystyle v_{x}},vy{\displaystyle v_{y}}, yvz{\displaystyle v_{z}}pero con variaciónkBT/metro{\textstyle k_{\text{B}}T/m}También se puede observar que la distribución de velocidad de Maxwell-Boltzmann para el vector velocidad [ v x , v y , v z ] es el producto de las distribuciones para cada una de las tres direcciones: Fv(vincógnita,vy,vz)=Fv(vincógnita)Fv(vy)Fv(vz){\displaystyle f_{\mathbf {v} }{\left(v_{x},v_{y},v_{z}\right)}=f_{v}(v_{x})f_{v}(v_{y})f_{v}(v_{z})} donde la distribución para una sola dirección es Fv(vi)=metro2πkBTexp(metrovi22kBT).{\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μvincógnita=μvy=μvz=0{\displaystyle \mu _{v_{x}}=\mu _{v_{y}}=\mu _{v_{z}}=0}y desviación estándarσvincógnita=σvy=σvz=kBT/metro{\textstyle \sigma _{v_{x}}=\sigma _{v_{y}}=\sigma _{v_{z}}={\sqrt {k_{\text{B}}T/m}}}, por lo que el vector tiene una distribución normal tridimensional, un tipo particular de distribución normal multivariada , con mediaμv=0{\displaystyle \mu _{\mathbf {v} }=\mathbf {0} }y covarianzaΣv=(kBTmetro)I{\textstyle \Sigma _{\mathbf {v} }=\left({\frac {k_{\text{B}}T}{m}}\right)I}, dóndeI{\displaystyle I}es la matriz identidad de 3 × 3 .

Una propiedad notable de la distribución del vector de velocidad es la independencia de la dirección, lo que significa que los componentes de la velocidad se distribuyen normalmente en cualquier dirección seleccionada, no solo en las tres direcciones base.incógnita{\displaystyle x},y{\displaystyle y}, yz{\displaystyle z}. [ 14 ]

Distribución para la velocidad

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

F(v)=2π[metrokBT]3/2v2exp(metrov22kBT).{\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).}

En el espacio n -dimensional

En un espacio n -dimensional, la distribución de Maxwell-Boltzmann se convierte en: F(v) dnortev=[metro2πkBT]norte/2exp(metro|v|22kBT) dnortev{\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: F(v) dv=Aexp(metrov22kBT)vnorte1 dv{\displaystyle f(v)~dv=A\exp \left(-{\frac {mv^{2}}{2k_{\text{B}}T}}\right)v^{n-1}~dv} dóndeA{\displaystyle A}es una constante de normalización.

El siguiente resultado integral resulta útil: 0vaexp(metrov22kBT)dv=[2kBTmetro]a+120miincógnitaincógnitaa/2dincógnita1/2=[2kBTmetro]a+120miincógnitaincógnitaa/2incógnita1/22dincógnita=[2kBTmetro]a+12Γ(a+12)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}}} dóndeΓ(z){\displaystyle \Gamma (z)}es la función Gamma . Este resultado se puede utilizar para calcular los momentos de la función de distribución de velocidad: v=0vvnorte1exp(metrov22kBT)dv0vnorte1exp(metrov22kBT)dv=2kBTmetro  Γ(norte+12)Γ(norte2){\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)}}}} que es la velocidad media mismavavgramo=v=2kBTmetro Γ(norte+12)Γ(norte2).{\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)}}.}

v2=0v2vnorte1exp(metrov22kBT)dv0vnorte1exp(metrov22kBT)dv=[2kBTmetro]Γ(norte+22)Γ(norte2)=[2kBTmetro]norte2=nortekBTmetro{\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}}} lo que proporciona una velocidad cuadrática media vrms=v2=nortekBTmetro.{\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: dF(v)dv=Aexp(metrov22kBT)[metrovkBTvnorte1+(norte1)vnorte2]=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 arroja la velocidad más probable ( modo ).vpag=(norte1)kBT/metro.{\textstyle v_{\text{p}}={\sqrt {\left(n-1\right)k_{\text{B}}T/m}}.}

Extensión a gases reales

Las deducciones muestran que la validez de la distribución de velocidad de Maxwell-Boltzmann se limita a gases ideales. Se conoce una generalización de la fórmula a todos los gases (ideales y reales por igual), cuya deducción parte del hecho de que las propiedades de los gases ideales y reales deben ser independientes de la dirección. La fórmula obtenida contienepagVmetro{\displaystyle pV_{\text{m}}}términos en lugar deRT{\displaystyle RT}, dóndepag{\displaystyle p}es la presión,Vmetro{\displaystyle V_{\text{m}}}es el volumen molar de la muestra de gas: [ 15 ]

F(v)=[METRO2πpagVmetro]3/24πv2exp(METROv22pagVmetro).{\displaystyle f(v)={\biggl [}{\frac {M}{2\pi pV_{\text{m}}}}{\biggr ]}^{{3}/{2}}\,4\pi v^{2}\exp \left(-{\frac {Mv^{2}}{2pV_{\text{m}}}}\right).}

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 ,3RTMETROmetro{\displaystyle {\frac {3RT}{M_{\text{m}}}}}sigue siendo la energía cinética de traslación media . El hecho de que el nitrógeno sea diatómico solo afecta al valor de la masa molar 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. 1 2 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 (2008). 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 (2.ª edición), RG Lerner , GL Trigg, VHC publishers, 1991, ISBN 3-527-26954-1(Verlagsgesellschaft), ISBN 0-89573-752-3(VHC Inc.)
  4. NA Krall y AW Trivelpiece, Principios de física de plasmas, San Francisco Press, Inc., 1986, entre muchos otros textos sobre física básica de plasmas
  5. 1 2 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, pp. 19-32.
  6. 1 2 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, pp. 21-37.
  7. Müller-Kirsten, HJW (2013). "2". Fundamentos de física estadística (2.ª ed.). World Scientific . ISBN  978-981-4449-53-3OCLC 822895930 
  8. Serway, Raymond A.; Faughn, Jerry S. y Vuille, Chris (2011). Física universitaria, volumen 1 (9.ª ed.). 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 en 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 de McGraw-Hill (2.ª ed.). 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.
  14. Rapp-Kindner, I.; Ősz, K.; Lente, G. (2025). "La ley de los gases ideales: derivaciones y antecedentes intelectuales" . ChemTexts . 11 (1): 1. doi : 10.1007/s40828-024-00198-9 .
  15. Lente, G. (2025). "La independencia de la dirección como propiedad clave para derivar una distribución de velocidad de partículas en gases reales" . Journal of Mathematical Chemistry . 63. doi : 10.1007/s10910-025-01742-9 .

Lecturas adicionales

  • 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-3OCLC 244177312 
  • Ives, David JG (1971). Termodinámica química . Química universitaria. 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 estadística de la velocidad". Revista Internacional de Termodinámica . 8 (9): 1– 14.
  • "La distribución de velocidad de Maxwell" del proyecto de demostraciones de Wolfram en Mathworld.