Вириал

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск

Вириал 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)\delta_{ik}-\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.