En mecánica clásica, el teorema del virial es una ecuación general que relaciona la energía cinética total promedio
de un sistema con su energía potencial promedio
, donde los paréntesis angulares representan el promedio temporal de la magnitud contenida entre ellos. Matemáticamente, el teorema del virial establece que:
Donde Fk representa la fuerza sobre la partícula k-ésima, que está ubicada en la posición rk.
Aplicaciones
El teorema del virial permite calcular la energía cinética total promedio aun para sistemas muy complejos en los que es muy difícil obtener una solución exacta, tales como los relacionados en mecánica estadística; esta energía cinética total promedio se relaciona con la temperatura del sistema a través del teorema de equipartición. Un ejemplo de sus muchas aplicaciones es el uso del teorema del virial para calcular el límite de Chandrasekhar para la estabilidad de las estrellas enanas blancas. La palabra «virial» tiene su origen en vis, la palabra en Latín para «fuerza» o «energía», y Clausius en 1870 le dio su acepción técnica.[1]
Si la fuerza entre dos partículas cualesquiera del sistema es producida por una energía potencial V(r)=αr n que es proporcional a alguna potencia n de la distancia entre las partículas r, el teorema del virial adopta la forma:
En Termodinámica, el teorema del virial nos permite escribir un modelo que se aproxime a un gas real, que se encuentre en la Naturaleza. Para ello, se usa un desarrollo en potencias de 1/v, y se obtiene (en magnitudes molares):
Donde B(T), C(T), ..., son el segundo coeficiente del virial, tercer coeficiente del virial respectivamente. A este desarrollo también se le conoce con el nombre de desarrollo de Kammerlingh Onnes. Como ejemplo, el gas de van der Waals puede escribirse usando el desarrollo de Kammerlingh Onnes como (de nuevo, en magnitudes molares):
Por lo tanto, dos veces la energía cinética total
es igual a n veces la energía potencial total promedio
. Donde V(r) representa la energía potencial entre dos partículas, VTOT representa la energía potencial total del sistema, o sea la suma de la energía potencial V(r) sobre todos los pares de partículas en el sistema. Un ejemplo común de este sistema es una estrella que se existe gracias a su propia fuerza de gravedad, donde n es -1.
Aunque el teorema del virial depende de promediar la energía cinética total y la energía potencial total, esta presentación deja para un paso próximo el promediar.
Definiciones del virial y su derivada temporal
Para un grupo de
partículas puntuales, el momento de inercia escalar I con respecto al origen queda definido por la ecuación
![{\displaystyle I=\sum _{k=1}^{N}m_{k}\mathbf {r} _{k}^{2}=\sum _{k=1}^{N}m_{k}r_{k}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ea3405351653d8f0dc822795c912336d61536cc)
donde mk y rk representan la masa y la posición de la partícula késima. El virial escalar G queda definido por la ecuación
![{\displaystyle G=\sum _{k=1}^{N}\mathbf {p} _{k}\cdot \mathbf {r} _{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/439a1d79fa5ad1a5a3ff089335de5a48ac45395a)
donde pk es el vector momento de la partícula késima. Suponiendo que las masas son constantes, el virial G es la derivada temporal de este momento de inercia
![{\displaystyle G={\frac {1}{2}}{\frac {dI}{dt}}=\sum _{k=1}^{N}m_{k}{\frac {d\mathbf {r} _{k}}{dt}}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\mathbf {p} _{k}\cdot \mathbf {r} _{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e49006a02749006d00fa8c4b1cf93099b7cf1731)
A su vez, la derivada temporal del virial G es
![{\displaystyle {\frac {dG}{dt}}=\sum _{k=1}^{N}\mathbf {p} _{k}\cdot {\frac {d\mathbf {r} _{k}}{dt}}+\sum _{k=1}^{N}{\frac {d\mathbf {p} _{k}}{dt}}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}m_{k}{\frac {d\mathbf {r} _{k}}{dt}}\cdot {\frac {d\mathbf {r} _{k}}{dt}}+\sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dcbb3aa4a561082ebd1dadf74fb5bc106dfb1914)
o, en forma más simple,
![{\displaystyle {\frac {dG}{dt}}=2T+\sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d4c6d6bdf4a4d2645d793dc874aaac3d74351af)
Aquí
es la masa de la partícula
,
es la fuerza neta sobre la partícula y
es la energía cinética total del sistema
![{\displaystyle T={\frac {1}{2}}\sum _{k=1}^{N}m_{k}v_{k}^{2}={\frac {1}{2}}\sum _{k=1}^{N}m_{k}{\frac {d\mathbf {r} _{k}}{dt}}\cdot {\frac {d\mathbf {r} _{k}}{dt}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8df05d49ea38d5d0d6b109a02849a18d885a568a)
Conexión con la energía potencial entre partículas
La fuerza total
sobre la partícula
es la suma de todas las fuerzas que ejercen todas las otras partículas
en el sistema
![{\displaystyle \mathbf {F} _{k}=\sum _{j=1}^{N}\mathbf {F} _{jk}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37d348319d8c0966df13a29d95935f4dc0ede304)
donde
es la fuerza aplicada por la partícula
sobre la partícula
. Por lo tanto, el término de fuerza de la derivada temporal del virial resulta ser
![{\displaystyle \sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\sum _{j=1}^{N}\mathbf {F} _{jk}\cdot \mathbf {r} _{k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e407523424d60d5f567db12a0163a40a631429c0)
Dado que ninguna partícula actúa sobre sí misma (o sea,
siempre que
), se tiene
![{\displaystyle \sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\sum _{j<k}\mathbf {F} _{jk}\cdot \mathbf {r} _{k}+\sum _{k=1}^{N}\sum _{j>k}\mathbf {F} _{jk}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\sum _{j<k}\mathbf {F} _{jk}\cdot \left(\mathbf {r} _{k}-\mathbf {r} _{j}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d156ccaca488ba8db2e849d805f41348ba5e435c)
donde se ha supuesto que vale la tercera ley del movimiento de Newton, o sea,
(una reacción igual y opuesta).
A menudo sucede que las fuerzas son producto de una energía potencial
que es solo función de la distancia
entre las partículas
y
. Dado que la fuerza es el gradiente de la energía potencial, entonces resulta que
![{\displaystyle \mathbf {F} _{jk}=-\nabla _{\mathbf {r} _{k}}V=-{\frac {dV}{dr}}\left({\frac {\mathbf {r} _{k}-\mathbf {r} _{j}}{r_{jk}}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cae4e7b352fb369c1a1005f8cf8134d539413b96)
lo cual es igual y opuesto a
, la fuerza aplicada por la partícula
sobre la partícula
, lo que se puede confirmar mediante un cálculo explícito. Por lo tanto, el término fuerza de la derivada temporal del virial es
![{\displaystyle \sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\sum _{j<k}\mathbf {F} _{jk}\cdot \left(\mathbf {r} _{k}-\mathbf {r} _{j}\right)=-\sum _{k=1}^{N}\sum _{j<k}{\frac {dV}{dr}}{\frac {\left(\mathbf {r} _{k}-\mathbf {r} _{j}\right)^{2}}{r_{jk}}}=-\sum _{k=1}^{N}\sum _{j<k}{\frac {dV}{dr}}r_{jk}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d9fbb5949b4d502b73d4cf1d35942b20128020d4)
Por lo tanto, se tiene
![{\displaystyle {\frac {dG}{dt}}=2T+\sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=2T-\sum _{k=1}^{N}\sum _{j<k}{\frac {dV}{dr}}r_{jk}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/39041e4385fcdf7135453f944735f03fcf6d4cca)
Caso especial de fuerzas dependientes de potencias
Un caso especial común es aquel en el cual la energía potencial
entre dos partículas es proporcional a una potencia n de la distancia que las separa r
![{\displaystyle V(r_{jk})=\alpha r_{jk}^{n},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8312973ec44a0192254f9d068c79f62420401a46)
donde el coeficiente α y el exponente n son constantes. En estos casos, el término de fuerza de la derivada temporal del virial se expresa por la ecuación
![{\displaystyle -\sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=\sum _{k=1}^{N}\sum _{j<k}{\frac {dV}{dr}}r_{jk}=\sum _{k=1}^{N}\sum _{j<k}nV(r_{jk})=nV_{TOT}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f1e75c83a64b3bd8f8dc0036c5915dae36c235c)
donde
es la energía potencial total del sistema
![{\displaystyle V_{TOT}=\sum _{k=1}^{N}\sum _{j<k}V(r_{jk}).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a744a3b6493f42bfbe926e4f7c0406aa2e05dbc)
Por lo tanto, se tiene
![{\displaystyle {\frac {dG}{dt}}=2T+\sum _{k=1}^{N}\mathbf {F} _{k}\cdot \mathbf {r} _{k}=2T-nV_{TOT}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/58306b33e6b63fe658a91498813642e4e44606d2)
Para sistemas gravitatorios y para sistemas electrostáticos, el exponente n es -1, resultando la identidad de Lagrange
![{\displaystyle {\frac {dG}{dt}}={\frac {1}{2}}{\frac {d^{2}I}{dt^{2}}}=2T+V_{TOT}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6deb71e4e0cb2b15515e2ffcfc1bfe2bfd65f4a)
lo cual fue descubierto por Lagrange y posteriormente extendido por Jacobi.
Promedio temporal y el teorema del virial
El promedio de esta derivada en un lapso de tiempo
se define como
![{\displaystyle \left\langle {\frac {dG}{dt}}\right\rangle _{\tau }={\frac {1}{\tau }}\int _{0}^{\tau }{\frac {dG}{dt}}\,dt={\frac {1}{\tau }}\int _{0}^{\tau }dG={\frac {G(\tau )-G(0)}{\tau }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/165aa2fc01e3aeddcb0bcc198eda7af396a2f51d)
de donde se obtiene la siguiente ecuación exacta
![{\displaystyle \left\langle {\frac {dG}{dt}}\right\rangle _{\tau }=2\left\langle T\right\rangle _{\tau }+\sum _{k=1}^{N}\left\langle \mathbf {F} _{k}\cdot \mathbf {r} _{k}\right\rangle _{\tau }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b4392b126bd5c1987babf8e0cad98597c7793334)
El teorema del virial afirma que, si
, entonces
![{\displaystyle 2\left\langle T\right\rangle _{\tau }=-\sum _{k=1}^{N}\left\langle \mathbf {F} _{k}\cdot \mathbf {r} _{k}\right\rangle _{\tau }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a70087853eef9a4d1b4fe557393e25fae42eb199)
Existen numerosas razones por las cuales el promedio de la derivada temporal se puede anular, o sea
. Una razón que se menciona se aplica a sistemas constreñidos, o sea sistemas que se encuentran limitados a permanecer juntos por siempre. En este caso, el virial
por lo general queda acotado entre dos extremos,
y
, y el promedio tiende a cero en el límite de tiempos muy largos
![{\displaystyle \lim _{\tau \rightarrow \infty }\left|\left\langle {\frac {dG^{\mathrm {bound} }}{dt}}\right\rangle _{\tau }\right|=\lim _{\tau \rightarrow \infty }\left|{\frac {G(\tau )-G(0)}{\tau }}\right|\leq \lim _{\tau \rightarrow \infty }{\frac {G_{\max }-G_{\min }}{\tau }}=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e13faaed2cf087a09e6620672d3b729c1d6b7bc5)
Aun si el promedio de la derivada temporal
es aproximadamente cero, el teorema del virial vale con el mismo grado de aproximación.
Para fuerzas que obedecen a una ley de potencia con un exponente n, la ecuación general establece que
![{\displaystyle \langle T\rangle _{\tau }=-{\frac {1}{2}}\sum _{k=1}^{N}\langle \mathbf {F} _{k}\cdot \mathbf {r} _{k}\rangle _{\tau }={\frac {n}{2}}\langle V_{TOT}\rangle _{\tau }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7655766a9d42db042c2b4b97938ae61179ff10b7)
Para el caso de atracción gravitatoria, n es igual a -1 y la energía cinética promedio es igual a un medio de la energía potencial promedio negativa
![{\displaystyle \langle T\rangle _{\tau }=-{\frac {1}{2}}\langle V_{TOT}\rangle _{\tau }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e926b63fbe98e1b13bb0c5e21e562c8d2098435)
Este resultado es útil para sistemas gravitatorios complejos tales como sistemas solares o galaxias.
No es preciso que el promedio sea en el tiempo; se puede realizar un promedio de colectivo, con resultados equivalentes.
Si bien ha sido desarrollado para la mecánica clásica, el teorema del virial es también válido en el ámbito de la mecánica cuántica.
Extensiones del teorema del virial
En 1900, Lord Rayleigh publicó una generalización del teorema del virial.[2] Henri Poincaré utilizó una forma del teorema del virial en 1911 para el problema de determinar la estabilidad cosmológica.[3] En 1945 Ledoux desarrolló una forma variacional del teorema del virial.[4] Parker[5] Chandrasekhar[6] y Fermi[7] a su vez desarrollaron formas tensoriales del teorema del virial,
Inclusión de campos electromagnéticos
Es posible extender el teorema del virial para incluir campos eléctricos y magnéticos. El resultado es[8]
![{\displaystyle {\frac {1}{2}}{\frac {d^{2}}{dt^{2}}}I+\int _{V}x_{k}{\frac {\partial G_{k}}{\partial t}}d^{3}r=2(T+U)+W^{E}+W^{M}-\int x_{k}(p_{ik}+T_{ik})dS_{i},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69af3cd9783eee34518ae5b333bb7801b96cc632)
donde I es el momento de inercia, G es la densidad de momento del campo electromagnético, T es la energía cinética del "fluido", U es la energía «térmica» aleatoria de las partículas, WE y WM son la energía eléctrica y magnética contenidas en el volumen bajo consideración. Finalmente, pik es el tensor de presión del fluido expresdo en el sistema de coordenadas móvil local
,
y Tik es el tensor electromagnético de tensiones,
![{\displaystyle T_{ik}=\left({\frac {\varepsilon _{0}E^{2}}{2}}+{\frac {B^{2}}{2\mu _{0}}}\right)-\left(\varepsilon _{0}E_{i}E_{k}+{\frac {B_{i}B_{k}}{\mu _{0}}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54bd8565739836d23e3f71835c7fa7f9f750bac1)
Demostración
Empleando el formalismo lagrangiano definimos la siguiente magnitud:
Siendo
las coordenadas generalizadas y
los momentos generalizados.
A continuación calculamos
Suponiendo que en el sistema dado, las coordenadas y momentos generalizados están acotados, concluimos que:
Además, puesto que:
Obtenemos finalmente:
Véase también
Referencias
- ↑ Clausius, RJE (1870). «On a Mechanical Theorem Applicable to Heat». Philosophical Magazine, Ser. 4 40: 122–127.
- ↑ Lord Rayleigh (1903). Unknown.
- ↑ Poincare, H. Lectures on Cosmological Theories. Paris: Hermann.
- ↑ Ledoux, P. (1945). «On the Radial Pulsation of Gaseous Stars». Ap. J. 102: 143-153.
- ↑ Parker, E.N. (1954). «Tensor Virial Equations». Physical Review 96 (6): 1686-1689.
- ↑ Chandrasekhar, S; Lebovitz NR (1962). «Unknown». Ap. J. 136: 1037-1047.
- ↑ Chandrasekhar, S; Fermi E (1953). «Unknown». Ap. J. 118: 116.
- ↑ George Schmidt, Physics of High Temperature Plasmas (Second edition), Academic Press (1979), p.72
Bibliografía
- Goldstein, H (1980). Classical Mechanics (2nd. ed edición). Addison-Wesley. ISBN 0-201-02918-9.
- Collins, GW (1978). The Virial Theorem in Stellar Astrophysics. Pachart Press.
- Callen, Herbert B. (1985). Thermodynamics and an introduction to Thermostatistics.
Datos: Q620602