Вириал

Вириал $G$ для множества $N$ точечных частиц в механике определяется как:

$G = \sum_{k=1}^N\mathbf{p}_k\cdot\mathbf{r}_k,$

где $\mathbf{r}_k$ и $\mathbf{p}_k$ — пространственные вектора координат и импульсов для $k$-й частицы.

Выражение «вириал» происходит от латинских слов «vis», «viris» — «сила» или «энергия». Оно было введено Клаузиусом в 1870 году.

Теорема о вириале

Для стабильной системы, связанной потенциальными силами, справедлива теорема о вириале:

$2\langle T\rangle=-\sum_{k=1}^N\langle\mathbf{F}_k\cdot\mathbf{r}_k\rangle,$

где $\langle T\rangle$ представляет среднюю полную кинетическую энергию и $\mathbf{F}_k$ — сила, действующая на $k$-ю частицу.

В частном случае, когда соответствующая силе потенциальная энергия взаимодействия $V(r)$ пропорциональна $n$-й степени расстояния между частицами $r$, вириальная теорема принимает простую форму

$2\langle T\rangle=n\langle U\rangle.$

Другими словами, удвоенная средняя полная кинетическая энергия $T$ равна $n$-кратной средней полной потенциальной энергии $U$.

Значение теоремы о вириале состоит в том, что она позволяет вычислить среднюю полную кинетическую энергию даже для очень сложных систем, что бросает вызов точным решениям, которые рассматривает, например, статистическая механика. Например, теорему о вириале можно использовать, чтобы вывести эквипарциальную теорему (теорема о равномерности распределении энергии по степеням свободы) или вычислить предел Чандрасекара для устойчивости белого карлика.

Производная по времени и усреднение

Производную по времени от вириала можно записать

$\frac{dG}{dt}=\sum_{k=1}^N\frac{d\mathbf{p}_k}{dt}\cdot\mathbf{r}_k+\sum_{k=1}^N\mathbf{p}_k\cdot\frac{d\mathbf{r}_k}{dt}=$

$=\sum_{k=1}^N\mathbf{F}_k\cdot\mathbf{r}_k+\sum_{k=1}^N m_k\frac{d\mathbf{r}_k}{dt}\cdot\frac{d\mathbf{r}_k}{dt}$

или в более простой форме

$\frac{dG}{dt}=2T+\sum_{k=1}^N\mathbf{F}_k\cdot\mathbf{r}_k.$

Здесь $m_k$ масса $k$-й частицы, $\mathbf{F}_k=\frac{d\mathbf{p}_k}{dt}$ — полная сила, действующая на частицу, а $T$ — полная кинетическая энергия системы

$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}.$

Усреднение этой производной за время $\tau$ определяется следующим образом:

$\left\langle\frac{dG}{dt}\right\rangle_\tau=\frac{1}{\tau}\int\limits_0^\tau\frac{dG}{dt}\,dt=\frac{1}{\tau}\int\limits_0^\tau dG=\frac{G(\tau)-G(0)}{\tau},$

откуда мы получим точное решение

$\left\langle\frac{dG}{dt}\right\rangle_\tau=2\langle T\rangle_\tau+\sum_{k=1}^N\langle\mathbf{F}_k\cdot\mathbf{r}_k\rangle_\tau.$

Вириальная теорема

Вириальная теорема утверждает:

Если $\left\langle\frac{dG}{dt}\right\rangle_\tau=0$, то

$2\langle T\rangle_\tau=-\sum_{k=1}^N\langle\mathbf{F}_k\cdot\mathbf{r}_k\rangle_\tau.$

Имеется несколько причин того, почему усреднение производной по времени исчезает, то есть $\left\langle\frac{dG}{dt}\right\rangle_\tau=0$. Одна часто цитируемая причина апеллирует к связанным системам, то есть системам, которые остаются ограниченными в пространстве. В этом случае вириал $G^{\mathrm{bound}}$ обычно ограничен двумя пределами, $G_\min$ и $G_\max$, и среднее стремится к нулю в пределе очень долгих времен $\tau$:

$\lim_{\tau\to\infty}\left|\left\langle\frac{dG^{\mathrm{bound}}}{dt}\right\rangle_\tau\right|=\lim_{\tau\to\infty}\left|\frac{G(\tau)-G(0)}{\tau}\right|\leqslant\lim_{\tau\to\infty}\frac{G_\max-G_\min}{\tau}=0.$

Если среднее значение производной по времени $\left\langle\frac{dG}{dt}\right\rangle_\tau\approx 0$, вириальная теорема имеет ту же степень приближения.

Соотношение с потенциальной энергией

Полная сила $\mathbf{F}_k$, действующая на частицу $k$, есть сумма всех сил действующих со стороны других частиц $j$ в системе

$\mathbf{F}_k=\sum_{j=1}^N\mathbf{F}_{jk},$

где $\mathbf{F}_{jk}$ — сила, действующая на частицу $j$ со стороны частицы $k$. Отсюда, слагаемое в производной по времени от вириала, содержащее силу, можно переписать в виде:

$\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.$

Поскольку отсутствует самодействие (то есть $\mathbf{F}_{jk}=0$, где $j=k$), мы получим:

$\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(\mathbf{r}_k-\mathbf{r}_j),$^[1]

где мы предположим, что выполняется третий закон Ньютона, то есть $\mathbf{F}_{jk}=-\mathbf{F}_{kj}$ (равны по модулю и противоположны по направлению).

Часто случается, что силы могут быть получены из потенциальной энергии $V$, которая является функцией только расстояния $r_{jk}$ между точечными частицами $j$ и $k$. Поскольку сила — это градиент потенциальной энергии с обратным знаком, мы имеем в этом случае

$\mathbf{F}_{jk} = -\nabla_{\mathbf{r}_k}V=-\frac{dV}{dr}\frac{\mathbf{r}_k-\mathbf{r}_j}{r_{jk}},$

который равен по модулю и противоположен по направлению вектору $\mathbf{F}_{kj} = -\nabla_{\mathbf{r}_j}V$ — силе, которая действует со стороны частицы $k$ на частицу $j$, как можно показать простыми вычислениями. Отсюда силовое слагаемое в производной от вириала по времени равно

$\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-\mathbf{r}_j)=-\sum_{k=1}^N\sum_{j<k} \frac{dV}{dr}\frac{(\mathbf{r}_k-\mathbf{r}_j)^2}{r_{jk}}=-\sum_{k=1}^N\sum_{j<k}\frac{dV}{dr}r_{jk}.$

Применение к силам, зависящим от расстояния степенным образом

Часто оказывается, что потенциальная энергия $V$ имеет вид степенной функции

$V(r_{jk})=\alpha r_{jk}^n,$

где коэффициент $\alpha$ и показатель $n$ — константы. В таком случае, силовое слагаемое в производной от вириала по времени задаётся следующими уравнениями

$-\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})=nU,$

где $U$ — полная потенциальная энергия системы:

$U=\sum_{k=1}^N\sum_{j<k}V(r_{jk}).$

В таких случаях, когда среднее от производной по времени от вириала $\left\langle\frac{dG}{dt}\right\rangle_\tau=0$, выполняется уравнение

$\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 U\rangle_\tau.$

Обычно приводимый пример — гравитационное притяжение, для которого $n=-1$. В том случае, средняя кинетическая энергия — половина средней отрицательной потенциальной энергии

$\langle T\rangle_\tau=-\frac{1}{2}\langle U\rangle_\tau.$

Этот результат является замечательно полезным для сложных гравитационных систем, типа солнечная система или галактика, и также выполняется для электростатической системы, для которой $n=-1$ также.

Хотя это выражение получено для классической механики, вириальная теорема также верна для квантовой механики.

Учёт электромагнитных полей

Вириальную теорему можно обобщить на случай электрических и магнитных полей. Результат:^[2]

$\frac{1}{2}\frac{d^2}{dt^2}I+\int\limits_V x_k\frac{\partial G_k}{\partial t}\,d^3r=2(T+U)+W^E+W^M-\int x_k(p_{ik}+T_{ik})\,dS_i,$

где $I$ — момент инерции, $G$ — вектор Пойнтинга, $T$ — кинетическая энергия «жидкости», $U$ — случайная тепловая энергия частиц, $W^E$ и $W^M$ — энергия электрического и магнитного поля в рассматриваемом объёме системы, $p_{ik}$ — тензор давления жидкости выраженный в локальной движущейся системе координат, сопутствующей жидкости:

$p_{ik}=\Sigma n^\sigma m^\sigma\langle v_iv_k\rangle^\sigma-V_iV_k\Sigma m^\sigma n^\sigma$

и $T_{ik}$ — тензор энергии-импульса электромагнитного поля:

$T_{ik}=\left(\frac{\varepsilon_0E^2}{2}+\frac{B^2}{2\mu_0}\right)-\left(\varepsilon_0E_iE_k+\frac{B_iB_k}{\mu_0}\right).$

Плазмоид — ограниченная конфигурация магнитных полей и плазмы. С помощью вириальной теоремы легко показать, что любая такая конфигурация расширяется, если не сдерживается внешними силами. В конечной конфигурации поверхностный интеграл исчезнет без оказывающих давление стен или магнитных катушек. Так как все другие слагаемые справа положительные, ускорение момента инерции также будет положительно. Легко оценить время расширения $\tau$. Если полная масса $M$ ограничена в пределах радиуса $R$, то момент инерции — примерно $MR^2$, и левая сторона в вириальной теореме — $MR^2/\tau^2$. Слагаемые справа составляют в целом величину порядка $pR^3$, где $p$ — большее из плазменного давления или магнитного давления. Приравнивая эти два члена и учитывая, что $~M= m_i n V \sim m_i n R^3$, $~p \sim n kT$, $~ c^2_s \sim \frac {kT} {m_i}$, где $m_i$ есть масса иона, $n$ – концентрация ионов, $~ V \sim R^3$ – объём плазмоида, $k$ – постоянная Больцмана, $T$ – температура, для $\tau$ находим:

$~\tau\sim R/c_s,$

где $c_s$ является скоростью ионной акустической волны (или волны Альфена, если магнитное давление выше, чем плазменное давление). Таким образом, время жизни плазмоида, как ожидают, будет равняться по порядку величины акустическому (альфеновскому) времени прохождения.

См. также

• Вириальное разложение

Примечания

1. ↑ Доказательство этого равенства
2. ↑ Schmidt G. Physics of High Temperature Plasmas. — Second edition. — Academic Press, 1979. — p. 72.

Литература

• Goldstein H. Classical Mechanics. — 2nd. ed. — Addison-Wesley, 1980. — ISBN 0-201-02918-9.