Черновики физика
Роман Парпалак

Простейшая модель электрона: электромагнитная масса и гиромагнитное отношение

16 марта 2011 года, 16:53

Продолжим развлекаться с классической электродинамикой. В прошлый раз мы подсчитали момент импульса в системе, состоящей из магнитного момента (маленькой катушки) в центре большой равномерно заряженной сферы

$$\vec{L}={2\over 3}{Q\vec{\mathfrak{m}}\over ca}.$$(1)

Но магнитное поле вокруг сферы можно создать без дополнительной катушки, просто закрутив ее. Получается, что часть момента импульса вращающейся сферы запасается в ее электромагнитном поле. Определим эту добавку.

Магнитный момент вращающейся сферы 

Мы собираемся показать, что магнитное поле вне равномерно заряженной сферы радиуса a и заряда Q, вращающейся с угловой скоростью ω, точно совпадает с полем магнитного диполя.

Начинаем с векторного потенциала

$$\vec{A}={1\over c}\int\!{\vec{j}\over r}\,dV={Q\over 4\pi ca^2}\int\!{\vec{v}\over r}\,dS={Q\over 4\pi ca^2}\int\!{\vec{\omega}\times a\vec{n}\over r}\,dS,$$

где $$r=|\vec{R}-a\vec{n}|$$ — расстояние от элемента поверхности dS, задаваемого радиус-вектором $$a\vec{n}$$, до точки $$\vec{R}$$, в которой вычисляется векторный потенциал. Учитывая линейность векторного произведения, получаем

$$\vec{A}={Q\over 4\pi ca}\,\vec{\omega}\times\int\!{\vec{n}\over \sqrt{R^2 + a^2-2Ra\cos\chi}}\,dS.$$

Интеграл $$\vec{I}$$ в последнем выражении — это усреднение единичного вектора $$\vec{n}$$ по направлениям с весом 1/r. Исходя из соображений симметрии ясно, что в результате интегрирования мы получим вектор, параллельный вектору $$\vec{R}$$. Подтвердим это вычислением.

В сферических координатах $$\inline \vec{n}=(\sin\theta\cos\varphi, \sin\theta\sin\varphi, \cos\theta)$$. Направим вектор $$\vec{R}=(0,0,R)$$ вдоль оси z. Тогда угол χ между векторами $$\inline \vec{n}$$ и $$\inline \vec{R}$$ совпадает со сферической координатой θ. Z-компонента интеграла

$$I_z=a^2\int\!\cos\theta\,{1\over \sqrt{R^2 + a^2-2Ra\cos\theta}}\,d(\cos\theta)\,d\varphi.$$

Известно, что функцию 1/r можно разложить по многочленам Лежандра (нас интересует область> a)

$${1\over \sqrt{R^2 + a^2-2Ra\cos\theta}}=\sum\limits_{l=0}^\infty{a^l\over R^{l+1}}P_l(\cos\theta).$$

Подынтегральное выражение состоит из произведения этого ряда на многочлен Лежандра $$P_1(\cos\theta)=\cos\theta$$. Тогда из условия ортогональности

$$\int\limits_{-1}^{1} P_k(x)P_l(x)\,dx={2\over 2k+1}\delta_{kl}$$

сразу следует ответ

$$I_z={4\pi\over 3}{a^3\over R^2}.$$

Компоненты Ix и Iy пропорциональны интегралам от периодических функций sin φ и cos φ, и поэтому равны нулю. Переходя от проекций к вектору, для векторного потенциала получаем

$$\vec{A}={Qa^2\over 3c}\,\vec{\omega}\times{\vec{R}\over R^3}.$$

Таким образом, магнитное поле вне вращающейся сферы совпадает с полем магнитного диполя

$$\vec{\mathfrak{m}}={Qa^2\vec{\omega}\over 3c}.$$(2)

Для полноты отметим, что аналогичными вычислениями легко показать однородность магнитного поля внутри вращающейся сферы.

Электромагнитная масса

Коэффициент в (1) можно переписать по-другому, если ввести понятие электромагнитной массы. Оно подробно разбирается в главе 28 выпуска 6 фейнмановских лекций (ниже мы воспроизводим некоторые вычисления оттуда и фактически на протяжении двух постов разбираем задачу 2 к этой главе).

Вслед за Фейнманом мы будем называть электромагнитной массой коэффициент пропорциональности между скоростью равномерного движения сферы и импульсом электромагнитного поля. В нерелятивистском случае

$$\vec{E}=Q{\vec{R}\over R^3},\quad\vec{H}={Q\over c}{\vec{v}\times\vec{R}\over R^3}.$$

Тогда импульс электромагнитного поля дается интегралом

$$\vec{P}={1\over 4\pi c}\int\!\vec{E}\times\vec{H}\,dV={Q^2\over 4\pi c^2}\int\!{1\over R^4}\,\vec{n}\times\left[\vec{v}\times\vec{n}\right]\,dV.$$

Такой интеграл (с точностью до коэффициента) мы вычисляли в прошлый раз, поэтому сейчас просто выпишем ответ

$$\vec{P}={2Q^2 \over 3 ac^2}\,\vec{v}=m_e \vec{v}.$$

Гиромагнитное отношение и модель электрона

Теперь мы можем переписать (1) в таком виде

$$\vec{L}={m_ec\over Q}\,\vec{\mathfrak{m}}.$$

Мы получили интересный результат: гиромагнитное отношение для вращающейся безмассовой заряженной сферы совпадает с гиромагнитным отношением электрона.

Обычно в литературе по квантовой механике утверждается, что нельзя представлять себе спин электрона как его вращение, так как в такой модели скорость точек на его поверхности будет больше скорости света. Сейчас мы в этом убедимся.

Подставим в (1) магнитный момент (2) и спин электрона ½:

$${\hbar\over 2}={2\over 3}\,{Q\over ca}\,{Qa^2\omega\over 3c},$$

откуда отношение экваториальной скорости к скорости света есть

$${a\omega\over c}={9\over 4}\,{\hbar c\over Q^2}={9\over 4\alpha}\approx308.$$

Несостоятельность простейшей модели электрона проявляется и в том, что полная энергия электрического поля

$$U={1\over 8\pi}\int\!{E^2}\,dV={1\over 8\pi}\int\limits_a^{\infty}{Q^2\over R^4}\,4\pi R^2\,dR={Q^2 \over 2a}={4\over 3}\,m_ec^2$$

отличается от ожидаемой величины mec2.

Ключевые слова: электродинамика | Комментарии (5)
Роман Парпалак

Парадокс Фейнмана, или потоки энергии в постоянных электромагнитных полях

14 марта 2011 года, 18:13

Я предлагаю отдохнуть от актуальных проблем теоретической физики и заняться более приземленными вопросами классической электродинамики.

В фейнмановских лекциях по физике (выпуск 6, глава 17) есть описание следующего парадокса.

Представим, что мы конструируем прибор, в котором имеется тонкий круглый пластмассовый диск, укрепленный концентрически на оси с хорошими подшипниками, так что он совершенно свободно вращается. На диске имеется катушка из проволоки — короткий соленоид, концентричный по отношению к оси вращения. Через этот соленоид проходит постоянный ток I от маленькой батареи, также укрепленной на диске. Вблизи края диска по окружности на равном расстоянии размещены маленькие металлические шарики, изолированные друг от друга и от соленоида пластмассовым материалом диска. Каждый из этих проводящих шариков заряжен одинаковым зарядом Q. Вся картина стационарна, и диск неподвижен.

Предположим, что случайно, а может и намеренно, ток в соленоиде прекратился, но, разумеется, без какого-либо вмешательства извне. Пока через соленоид шел ток, более или менее параллельно оси диска проходил магнитный поток. После того как ток прервался, поток этот должен уменьшиться до нуля. Поэтому должно возникать индуцированное электрическое поле, которое будет циркулировать по окружностям с центром на оси диска. Заряженные шарики на периферии диска будут все испытывать действие электрического поля, касательного к внешней окружности диска. Эта электрическая сила направлена для всех зарядов одинаково и, следовательно, вызовет у диска вращающий момент. Из этих соображений можно ожидать, что, когда ток в соленоиде исчезнет, диск начнет вращаться. Если нам известны момент инерции диска, ток в соленоиде и заряд шариков, то можно вычислить результирующую угловую скорость.

Но можно рассуждать и по-другому. Используя закон сохранения момента количества движения, мы могли бы сказать, что момент диска со всеми его пристройками вначале равен нулю, поэтому момент всей системы должен оставаться нулевым. Никакого вращения при остановке тока быть не должно. Какое из доказательств правильно? Повернется ли диск или нет? Мы предлагаем вам подумать над этим вопросом.

Решение парадокса заключается в том, что в присутствии электрических и магнитных полей имеются потоки энергии, описываемые вектором Пойнтинга $$\vec{S} \sim \vec{E} \times \vec{B}$$. В предложенной Фейнманом конфигурации эти потоки замкнуты. Поскольку поток энергии однозначно связан с плотностью импульса, наличие замкнутых потоков энергии свидетельствует о присутствии ненулевого момента импульса. Таким образом, в системе изначально был запас момента импульса, который после исчезновения магнитного поля был передан диску.

Однако есть люди, которые убеждены в том, что выражение для вектора потока энергии через векторное произведение полей годится только для переменных электромагнитных полей (действительно, перенос энергии в этом случае можно увидеть непосредственно). Им не нравится, что, согласно выражению для вектора Пойнтинга, энергия течет от источника постоянного тока к нагрузке не «по проводам», а снаружи. Основной аргумент сводится к тому, что на поток энергии легко влияют манипуляции с проводом, наличие примесей, разрывов; в то время как никакими телами снаружи, например, дополнительными зарядами или магнитами, остановить поток энергии к нагрузке не удается.

Противники применения вектора Пойнтинга в статике придумывают в этом случае другие выражения для плотности энергии, например, $$\varphi \vec{j}$$. Но такое выражение тоже не лишено недостатков, приписываемых вектору Пойнтинга. Даже если отбросить требование единообразного описания явлений и попытаться применить для зарядов и катушки выражение $$\varphi \vec{j}$$, мы сразу столкнемся с тем, что потенциал φ в точках катушки (а, значит, и и момент импульса системы) легко изменить, поместив систему внутрь большого проводящего заряженного ящика и меняя его заряд. Это плохо согласуется с гипотезой, по которой энергия, как и ток, течет по катушке.

В связи с такой критикой вектора Пойнтинга можно заняться интересным упражнением — непосредственно вычислить момент импульса электромагнитного поля через вектор Пойнтинга и сравнить его с моментом импульса, передаваемым диску. Даже если критиков эти вычисления не убедят, мы еще раз увидим красоту теории.

$$\begin{tikzpicture}[x=(-30:3cm),y=(30:3cm),z=(90:3cm),>=stealth] \def\R{1} \draw[->] (0,0,0) -- (0,0,0.3\R) node[right] {$\vec{\mathfrak{m}}$}; \draw[thin] (0,0,0) -- (.7\R,.7\R,0) node[below,pos=0.5] {$a$}; \node at (0,\R,0.5*\R) {$Q$}; \def\r{0.075} \draw[->,thick] (\r,0,0) arc (0:360:\r); \foreach \i in {0, 30,...,150} \draw [dotted] plot [domain=-180:180, samples=30, variable=\j] (\R*cos \i*sin \j,\R*sin \i*sin \j, \R*cos \j); \foreach \j in {0, 15,...,90} \draw [dotted] plot [domain=0:360, samples=60, variable=\i] (\R*cos \i*sin \j,\R*sin \i*sin \j, \R*cos \j); \end{tikzpicture}$$

Для упрощения вычислений изменим систему, предложенную Фейнманом. Будем рассматривать не заряды на краю диска, а равномерно заряженную сферу радиуса a, в центре которой находится небольшая катушка, обладающая магнитным моментом $$\vec{\mathfrak{m}}$$. Магнитный момент, находящийся в начале координат, создает в точке $$\vec{R}$$ векторный потенциал

$$\vec{A} = {\vec{\mathfrak{m}} \times \vec{R} \over R^3}$$(1)

и магнитное поле

$$\vec{H} = {3\vec{n}\,(\vec{\mathfrak{m}} \cdot \vec{n})-\vec{\mathfrak{m}} \over R^3}, \quad \vec{n} = {\vec{R}\over R}.$$

Вектор Поинтинга есть

$$\vec{S} = {c\over 4\pi}\vec{E}\times\vec{H}.$$

Он связан с плотностью импульса

$$\vec{P} = {\vec{S} \over c^2}.$$

Тогда, собирая вместе, получаем

$$\vec{P} = {1\over 4\pi c}{Q \over R^2} {1 \over R^3}\left[ \vec{n}\times \left(3\vec{n}\,(\vec{\mathfrak{m}} \cdot \vec{n}) -\vec{\mathfrak{m}} \right) \right] = {Q\over 4\pi c R^5}\left[\vec{\mathfrak{m}} \times \vec{n}\right].$$

Момент импульса

$$\vec{L}=\int\!\vec{R}\times\vec{P}\,dV=\int\! R\,{Q\over 4\pi c R^5}\,\vec{n}\times\left[\vec{\mathfrak{m}} \times \vec{n}\right]dV.$$

Введем сферические координаты так, что $$\vec{\mathfrak{m}}$$ направлен вдоль оси z, а угол θ есть угол между осью z и направлением $$\vec{R}$$. Тогда

$$\vec{L}=\int\limits_a^{\infty}\!{Q\over 4\pi c R^4}R^2dR \int\!\vec{n}\times\left[\vec{\mathfrak{m}} \times \vec{n}\right]d(\cos \theta)\,d\varphi.$$(2)

Первый интеграл берется элементарно. Второй есть усреднение двойного векторного произведения по всем направлениям $$\vec{n}$$, он кроме как от $$\vec{\mathfrak{m}}$$ ни от чего не зависит. Учитывая линейность интегрирования и векторного произведения, зависимость должна быть прямой пропорциональностью. Действительно, вычисления показывают, что второй интеграл равен

$$\vec{\mathfrak{m}}\int\!\sin^2\theta\,d(\cos \theta)\,d\varphi={8\pi\vec{\mathfrak{m}}\over 3}.$$

Окончательно получаем

$$\vec{L}={2\over 3}{Q\vec{\mathfrak{m}}\over ca}.$$

Теперь посмотрим, какой момент импульса будет передан сфере после исчезновения магнитного поля. Начнем с момента электрических сил

$$\vec{M}=\int\!\vec{R}\times\rho\vec{E}\,dV={Q\over 4\pi}\int\!\vec{R}\times\vec{E}\,d(\cos\theta)\,d\varphi.$$

Учитывая, что вихревое электрическое поле определяется формулой

$$\vec{E}=-{1\over c}{\partial \vec{A}\over\partial t},$$

для момента импульса имеем

$$\vec{L}=\int\!\vec{M}\,dt={Q\over 4\pi c}\int\!d(\cos\theta)\,d\varphi\,\vec{R}\times\!\int\!-{\partial \vec{A}\over\partial t}dt.$$

Векторный потенциал уменьшается от начального значения $$\vec{A}$$, задаваемого формулой (1), до нуля. Поэтому последний интеграл просто равен $$\vec{A}$$. Таким образом,

$$\vec{L}=\cfrac{Q}{4\pi c a}\int\!d(\cos\theta)\,d\varphi\,\vec{n}\times \left[\vec{\mathfrak{m}} \times \vec{n}\right].$$

Это выражение совпадает с (2).

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

Ключевые слова: электродинамика | Комментарии (25)
Роман Парпалак

Ультрабуст в пространстве (2+1)

6 марта 2011 года, 10:34

В прошлый раз, говоря о (2+1)-мерной гравитации, мы нашли метрику точечной массы:

$$ds^2=-c^2dt^2+d\rho^2+\alpha^2\rho^2d\varphi^2,\quad\alpha=1-4Gm/c^2.$$(1)

Теперь попробуем вывести метрику ультрарелятивистской частицы, пролетающей мимо наблюдателя с околосветовой скоростью. Сначала мы найдем метрику движущейся материальной точки, а затем применим операцию, называемую ультрабустом. Она заключается в одновременном устремлении скорости частицы к скорости света и массы к нулю, чтобы энергия оставалась конечной.

Лоренцев буст

Временно будем считать, что скорость света с = 1.

Для плоской метрики пространства Минковского переход в другую ИСО выполняется при помощи матрицы Лоренца Λ (такое преобразование оставляет метрику инвариантной):

$$g(v)=\Lambda^Tg\Lambda.$$(2)

Эту формулу можно применять не только к плоскому пространству, но и, например, к решению Шварцшильда. Действительно, метрика Шварцшильда асимптотически плоская, поэтому, формально проделав над ней такое преобразование, мы получим на бесконечности преобразования Лоренца из СТО. Следовательно, (2) описывает переход из системы отсчета покоя в систему, движущуюся относительно тела с постоянной скоростью.

Однако метрика (1) не является асимптотически плоской, поэтому обоснование возможности применить (2) в случае пространства (2+1) сложнее.

Для начала заменим переменные:

$$\left\{ \begin{array}{rl} x\!\!\!\!&=r\cos\varphi,\\ y\!\!\!\!&=r\sin\varphi,\\ \varepsilon\!\!\!\!&=1-\alpha^2. \end{array} $$(3)

В этих обозначениях метрика (1) принимает вид

$$g=\left[ \begin{array}{ccc} -1&0&0 \\ 0&1-{\dfrac {\varepsilon\,{y}^{2}}{{x}^{2}+{y}^{2}}}&{\dfrac {\varepsilon\,xy}{{x}^{2}+{y}^{2}}} \\ 0&{\dfrac {\varepsilon\,xy}{{x}^{2}+{y}^{2}}}&1-{\dfrac {\varepsilon\,{x}^{2}}{{x}^{2}+{y}^{2}}} \end{array} \right] .$$

Отсюда и из (1) видно, что ε → 0, когда → 0. Когда мы уменьшаем массу точки до нуля, пространство переходит в плоское. Результат цепочки преобразований: уменьшение массы до нуля — лоренцев буст — увеличение массы до первоначального значения совпадет с результатом от простого применения формулы (2). Этим и обосновывается возможность ее применения.

Матрица Лоренца есть

$$\Lambda=\left[ \begin {array}{ccc} {\cfrac {1}{\sqrt {1-v^{2}}}}&-{\cfrac {v}{\sqrt {1-v^{2}}}}&0 \\ -{\cfrac{v}{\sqrt {1-{v}^{2}}}}&{\cfrac {1}{\sqrt {1-v^{2}}}}&0 \\ 0&0&1\end {array} \right] .$$

Несложные, но объемные вычисления, которые лучше всего поручить компьютеру, дают

$$\Lambda^Tg\Lambda=\left[\begin{array}{ccc} -1-{v^2\over 1-v^2}{\varepsilon y^2\over x^2+y^2}&{v\over 1-v^2}{\varepsilon y^2\over x^2+y^2}&-{v\over\sqrt{1-v^2}}{\varepsilon xy \over x^2+y^2}}\\ \noalign{\medskip}{v\over 1-v^2}{\varepsilon y^2\over x^2+y^2}&1-{1\over 1-v^2}{\varepsilon y^2\over x^2+y^2}&{1\over\sqrt{1-v^2}}{\varepsilon xy \over x^2+y^2}}\\ \noalign{\medskip}-{v\over\sqrt{1-v^2}}{\varepsilon xy \over x^2+y^2}}&{1\over\sqrt{1-v^2}}{\varepsilon xy \over x^2+y^2}}&1-{\varepsilon x^2 \over x^2+y^2}\end {array}\right].$$

В этом выражении x — координата в сопутствующей системе отсчета. Она связана с координатой x1 в нашей системе отсчета преобразованием Лоренца

$$x = {x_1-vt \over \sqrt{1-v^2}}.$$

Подставим ее в предыдущее выражение и получим

$$g(v)=\left[\begin{array}{ccc}-1-{v^2\varepsilon y^2\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}&{v\varepsilon y^2\over\left(x_1-vt\right)^2+y^2\left( 1-v^2 \right)}&-{v\varepsilon\left(x_1-vt\right)y\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}\\ \noalign{\medskip}{v\varepsilon y^2\over \left(x_1-vt\right)^2+y^2\left( 1-v^2 \right)}&1-{\varepsilon y^2\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}&{\varepsilon\left( x_1-vt\right)y\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}\\ \noalign{\medskip}-{v\varepsilon\left( x_1-vt \right) y\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}&{\varepsilon\left(x_1-vt\right)y\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}&1-{\varepsilon\left(x_1-vt\right)^2\over\left(x_1-vt\right)^2+y^2\left(1-v^2\right)}\end{array}\right].$$

Ультрабуст

Используем результат, полученный выше, для вывода метрики материальной точки с нулевой массой, движущейся со скоростью света.

Будем переходить к пределу  c, → 0 таким образом, чтобы энергия

$$E={mc^2 \over \sqrt{1-\dfrac{v^2}{c^2}}}$$

материальной точки оставалась конечной. С учетом (1) и (3) легко найти, что

$$\varepsilon={8GE \over {c^4}}\sqrt {1-\dfrac{v^2}{c^2}}-{16G^2E^2 \over c^8}\left( {1-\dfrac{v^2}{c^2}} \right).$$

Кроме того, нам понадобится известное соотношение

$$\lim \limits_{a\to +0} \dfrac{a}{z^2+a^2}=\pi \delta (z).$$

В конечном итоге мы получим

$$g(c ) =\left[\begin{array}{ccc} {-c^2-\dfrac{8\pi G}{c^2}E\left| y \right|\delta \left( {x_1 -ct} \right)}& {\dfrac{8\pi G}{c^3}E\left| y \right|\delta \left( {x_1 -ct} \right)}& 0\\{\dfrac{8\pi G}{c^3}E\left| y \right|\delta \left( {x_1 -ct} \right)}& {1-\dfrac{8\pi G}{c^4}E\left| y \right|\delta \left( {x_1 -ct} \right)}& 0 \\ 0& 0& 1\\ \end{array}\right] .$$

В переменных u = x1 − ct, x1 + ct интервал принимает вид

$$ds^2=-du\,dv+dy^2-\dfrac{8\pi G}{c^4}E\,\delta (u) \left| y \right|du^2.$$

Полезно сравнить этот ответ с ультрабустом в (3+1) [1]:

$$ds^2=-du\,dv+dr^2+r^2 d\theta ^2-8GE\,\delta (u) \ln r du^2.$$

Результаты отличаются функцией от перпендикулярного расстояния, входящей в коэффициент перед du2. В двумерном случае это линейная функция, а в трехмерном — логарифм. Эти функции являются функциями Грина уравнения Лапласа размерности на 2 меньше, чем размерность соответствующего пространства-времени. Как видим, результаты получились похожими.

Ссылки

[1] P. C. Aichelburg, R. U. Sexl (Vienna U.). «On the Gravitational field of a massless particle». May 1970. Published in Gen.Rel.Grav.2:303-312,1971.

[2] S. Deser, Alan R. Steif (Brandeis U.). «Gravity theories with lightlike sources in D = 3.» Published in Class.Quant.Grav.9:L153-L160,1992.

Ключевые слова: гравитация | Комментарии (9)
Роман Парпалак

Гравитация в пространстве (2+1)

2 марта 2011 года, 23:21

Мы рассмотрим некоторые особенности гравитации в трехмерном пространстве-времени (две пространственных координаты плюс время). Оказывается, такое пространство в присутствии масс локально не искривляется, однако в нем появляются глобальные топологические эффекты. Например, длина окружности, охватывающей материальную точку, меньше, чем 2πR.

Тензор Римана в (2+1)

Как известно, тензор Римана обладает следующими симметриями:

$$R_{abcd}=R_{cdab},\quad R_{abcd}=-R_{bacd}=-R_{abdc},$$

где каждый индекс пробегает значения 0, 1, 2. Непосредственным перебором легко проверить, что у тензора Римана остается 6 независимых ненулевых компонент: R0101, R0112, R0120, R1212, R1220, R2020. Число независимых компонент тензора Риччи тоже 6. Оказывается, в трехмерном пространстве-времени тензор Римана можно выразить через тензор Риччи:

$$R_{abcd}=g_{ac}Q_{bd} + g_{bd}Q_{ac}-g_{ad}Q_{bc}-g_{bc}Q_{ad},$$(1)

где QabRab − ¼gabR, RRaa.

Точечная масса

Вычислим статическую аксиально-симметричную метрику в присутствии точечной массы. Как обычно, начинаем с уравнений Эйнштейна

$$R_{ab} -{R \over 2} g_{ab} = {8 \pi G \over c^4} T_{ab}.$$

Взяв след обеих частей, получаем

$$R = -{16 \pi G \over c^4} T.$$(2)

В пустом пространстве вокруг тела Tab ≡ 0, следовательно, из (2), уравнения Эйнштейна и (1) получаем R ≡ 0, Rab ≡ 0, Rabcd ≡ 0. Таким образом, пустое пространство (2+1) должно быть локально плоским. Это, в частности, означает, что в рассматриваемом случае нет гравитационных волн.

$$\begin{tikzpicture}[line width=0.21mm,scale=1.0545] \def\re{1.2} \def\rta{0.32} \def\rtb{0.5} \def\a{310} \def\b{110} \def\c{1} \def\d{0.3} \def\lx{0.6} \def\ly{3.8} \draw[white!30!blue!70!green] (\re,0) arc (0:\a:\re) coordinate(y); \draw[thin] (y) -- (0,0) -- (\re,0); \draw[->,very thin] (\rta,0) arc (0:\a:\rta) node[pos=0.83,below] {$2\pi\alpha$}; \draw[->,very thin] (0,0) -- (\b:\re) (\rtb,0) arc (0:\b:\rtb) node[pos=0.75,above] {$\varphi'$}; \draw[white!30!blue!70!green] (\lx,-\ly) arc (-105:105:0.3 and \c) coordinate(x); \draw[thin] (x) -- +(-1.4,-\c) -- (\lx,-\ly); \draw[thin,dashed,white!30!blue!70!green] (\lx,-\ly) arc (255:105:0.3 and \c); \end{tikzpicture}$$

Метрика такого пространства с учетом симметрий имеет следующий вид:

$$ds^2=-c^2dt^2+d\rho^2+\alpha^2\rho^2d\varphi^2,$$

где $$\alpha^2$$ — некоторая константа, а φ меняется от 0 до 2π. Введем новую угловую координату φ′ = αφ. Константа α и будет определять глобальные гравитационные эффекты. Заглядывая вперед, скажем, что α < 1. Таким образом, длина окружности, охватывающей точечную массу, будет меньше, чем 2πR.

Время входит в интервал с постоянным коэффициентом и не смешивается с другими координатами, поэтому скаляр кривизны R определяется только пространственной частью метрики.

Конический дефект

Двумерное пространство, описываемое такой метрикой, соответствует конической поверхности. В каждой точке поверхности, за исключением вершины, тензор Риччи равен нулю. В вершине имеется расходимость, связанная с расходимостью плотности энергии точечного тела.

Обычно связь между массой (которая пропорциональна коэффициенту в кривизне, задаваемой дельта-функцией) и величиной дефекта α определяется через эйлерову характеристику поверхности (которая равна сумме интеграла от кривизны по поверхности и интеграла от «внешней кривизны» границы по самой границе). Мы же пойдем более наглядным и методически более простым путем, «размывая» точечную массу.

$$\begin{tikzpicture}[line width=0.21mm,scale=1.0545] \def\c{1} \def\cc{0.5} \def\beta{6} \def\g{18} \draw (0,0) arc (-90-\beta:90+\beta:0.2 and \c) -- +(-1,-\c+\cc) arc (90+\beta:-90-\beta:0.1 and \cc) -- cycle; \draw[thin,dashed] (0,0) arc (270-\beta:90+\beta:0.2 and \c) ++(-1,-\c+\cc) arc (90+\beta:270-\beta:0.1 and \cc); \draw (-1,\c+0.489) arc (90+\g:270-\g:0.3 and 0.523); \end{tikzpicture}$$

Можно было бы рассмотреть некоторое распределение массы в ограниченной области и выяснить, как будет меняться α, когда размер области стремится к нулю. Однако мы будем действовать противоположным образом, что избавит нас от необходимости решать уравнения Эйнштейна. Деформируем коническую поверхность, «сгладив» вершину в сферический сегмент радиуса r, и выясним, какому распределению масс соответствует такое пространство.

Известно, что скаляр кривизны двумерной поверхности равен удвоенной гауссовой кривизне:

$$R = {2 K} = {2 \over r_1r_2},$$

где r1 и r2 — главные радиусы кривизны поверхности. Для сферической поверхности радиуса r скаляр кривизны R = 2/r2.

$$\begin{tikzpicture}[line width=0.4mm,scale=1.0545]\small \def\r{1.8} \coordinate[label=above left:$A$] (A) at (-0.5*\r,0.866*\r); \coordinate[label=below:$B$] (B) at (-0.5*\r,0); \coordinate[label=below left:$C$] (C1) at (-\r,0); \coordinate[label=below:$O$] (O) at (0,0); \coordinate[label=below:$P$] (P) at (-2*\r,0); \coordinate (A3) at (0,1.333*0.866*\r); \draw[thin] (0,\r) arc (90:120:\r) -- node[above] {$\rho$} (P) -- node[pos=0.625,above] {$h$} node[pos=0.19,above,inner sep=1] {$\beta$} (O) -- node[right] {$r$} (A) |- (B); \draw[] (A3) -- (A) arc (120:190:\r); \draw[line width=0.21mm,opacity=0] (-2*\r-0.2,-0.4) rectangle (0.2,2.1) \end{tikzpicture}$$

Установим некоторые геометрические соотношения. Длина окружности в основании конуса равна $$2\pi|AB|=2\pi\rho\sin\beta$$, длина той же линии на развертке конуса равна $$2\pi\alpha|AP|=2\pi\alpha\rho$$, откуда $$\alpha=\sin\beta$$. Тогда площадь $$S=2\pi r\cdot |BC|$$ сферического сегмента высоты $$|BC|=h=r(1-\sin\beta)$$ выражается как $$2\pi r^2(1-\alpha)$$.

Проинтегрируем (2):

$$\int T \sqrt{-g}\, d^2x =- {c^4 \over 16\pi G}\int R \sqrt{-g}\, d^2 x=-{c^4 \over 16\pi G}RS = -{c^4 \over 4G} (1- \alpha).$$

Учитывая аддитивность массы (отсутствие гравитационного дефекта массы следует из того, что тензор Римана в пустом пространстве нулевой и пространство локально галилеево; мы обсудим это ниже), получаем

$$mc^2 =\int {T_{00}} \sqrt{-g}\, d^2x=-\int T \sqrt{-g}\, d^2x={c^4 \over 4G} (1-\alpha) ,$$

$$\alpha =1-{4Gm \over c^2}.$$

Интеграл от скаляра кривизны по поверхности не зависит от размера области «сглаживания», поэтому полученный результат справедлив и для точечной массы, когда r → 0. Как видим, для обычных тел с положительной массой α < 1. Видно также, что масса рассматриваемого тела не может превышать величину c2/4G.

Геодезические

$$ \begin{tikzpicture}[line width=0.21mm,scale=1.0545] \def\re{1.4} \def\ri{0.7} \def\ll{1.3} \def\l{0.7} \def\a{310} \def\x{118} \draw (0,0) -- (\re,0) arc (0:\a:\re) -- cycle; \draw (0:\ri) -- +(\x:\ll) (\a:\ri) -- +(180+\x+\a:\l); \end{tikzpicture} $$

Решая систему уравнений для геодезических

$$\frac{d^2x^\lambda }{d s^2} + \Gamma^{\lambda}_{~\mu \nu }\frac{dx^\mu }{d s}\frac{dx^\nu }{d s} = 0,$$

которая в наших координатах принимает вид

$$\left\{ \begin{array}{l}\ddot{\varphi} + \dfrac{2}{r}\dot{r}\dot{\varphi}=0, \\ \ddot{r}-r\dot{\varphi}^2\alpha^2=0, \\ \end{array} \right. $$

легко убедиться в том, что геодезические для пространства с точечной массой — это прямые $$r\sim 1/\sin(\alpha\varphi)$$ на развертке конуса.

Можно спроецировать геодезические на плоскость, перпендикулярную оси конуса. Тогда получится, что точечная масса «искривляет» пространство, или «отклоняет» движущиеся частицы. Однако легко видеть, что угол «отклонения» не зависит от прицельного параметра. Таким образом, уравнения Эйнштейна в (2+1) не содержат ньютоновское тяготение как предельный случай.

Более того, как показывает изучение геодезических, в определенном смысле гравитационное взаимодействие в (2+1) отсутствует. Как видно из вышеприведенной системы, для любых r0 и φ0 «кривая» покоя r = r0, φ = φ0 является геодезической. В обычном четырехмерном случае это не так. Например, для шварцшильдовского решения не существует такой системы координат, в которой тело в любой точке оставалось бы в покое.

Обобщения

Мы убедились в том, что точечной массе в трехмерном пространстве-времени соответствует плоскость с вырезанным углом (в вершине которого сама материальная точка) и отождествленными точками на противоположных сторонах разреза. Подобная картина характерна и для нескольких точек. С каждым телом связан вырезанный угол, величина которого пропорциональна массе тела. Сумма масс всех тел также не может превышать c2/4G.

Помимо бесконечной двумерной поверхности, можно представить еще и замкнутую поверхность, топологически эквивалентную сфере. Очевидно, масса такой замкнутой Вселенной есть c2/2G. Сфера будет реализовываться в случае равномерного распределения вещества. Дискретные массы будут образовывать многогранники. Например, восемь одинаковых материальных точек могут дать куб.

Ключевые слова: гравитация | Комментарии (8)
Роман Парпалак

Эффект Унру

28 февраля 2011 года, 15:37

Суть эффекта Унру заключается в том, что равноускоренный наблюдатель начинает видеть вокруг себя равновесное тепловое излучение, в то время как наблюдатель в инерциальной системе отсчета не видит ничего. В работе «Is there Unruh radiation?» авторов G. W. Ford и R. F. O'Connell есть вывод формулы для температуры Унру. Проследим за этим выводом.

Модель

Рассмотрим струну в пространстве (1 + 1) с лагранжианом

$$L=\int dy \left[ {\frac{\sigma }{2}}\left({\frac{\partial u}{\partial t}}\right)^{2}-{\frac{\tau }{2}}\left({\frac{\partial u}{\partial y}}\right)^{2}\right].$$

Для случая скалярного поля нужно взять σ = 1/4π, τ = с2/4π. Несложно получить уравнение движения

$${\frac{\partial ^{2}u}{\partial t^{2}}}-c^{2}{\frac{\partial^{2}u}{\partial y^{2}}}=0.$$

Его решение выписывается через ряд Фурье

$$u(y,t)=\sum_{k}\sqrt{{\frac{\hbar }{2\sigma L\omega }}}\left(a_{k}e^{i(ky- \omega t)}+a_{k}^{\dag }e^{-i(ky-\omega t)}\right),$$

где L — длина струны; частота и импульс связаны дисперсионным соотношением ω c|k|; сумма по импульсам пробегает значения, кратные 2π/L.

Квантование

Теперь проквантуем эту систему, потребовав выполнения коммутационных соотношений

$$\lbrack a_{k},a_{k^{\prime }}^{\dag }]=\delta _{k^{\prime}k},\quad\lbrack a_{k},a_{k^{\prime }}]=0.$$

Как утверждается, для струны в тепловом равновесии при температуре T можно определить следующие вакуумные средние:

$$\left\langle a_{k}a_{k^{\prime }}^{\dag }+a_{k^{\prime }}^{\dag}a_{k}\right\rangle =\,\mbox{cth}\, \frac{\hbar \omega }{2kT}\,\delta_{k^{\prime}k}, \quad\left\langle a_{k}a_{k^{\prime }}+a_{k^{\prime }}a_{k}\right\rangle =0.$$

Действительно, в первом выражении легко распознать $$2\bar{n}_k + 1$$. Среднее число квантов $$\bar{n}_k$$ дается статистикой Бозе — Эйнштейна, откуда и получается гиперболический котангенс.

Термодинамическое равновесие

Мы будем изучать поведение корреляционной функции поля

$$C(\Delta y,\Delta t)\equiv {\frac{1}{2}}\left\langle u(y_{1},t_{1}) u(y_{2},t_{2}) + u(y_{2},t_{2}) u(y_{1},t_{1}) \right\rangle ,$$

где Δy = y1 − y2, Δt = t1 − t2. Рассматривать спектральную плотность и пространственное распределение излучения было бы нагляднее. Но можно заниматься и корреляционной функцией, ведь она связана со спектральной плотностью (и, видимо, пространственным распределением) преобразованием Фурье.

После выполнения вычислений и перехода к бесконечной длине (→ ∞), связанного с заменой суммирования интегрированием, получаем

$$C(\Delta y,\Delta t)=\frac{\hbar }{4\pi \sigma }\int\limits_{-\infty }^{\infty }dk \frac{1}{\omega }\,\mbox{cth}\, \frac{\hbar \omega }{2kT}\cos \left( k\Delta y-\omega \Delta t\right).$$

Это выражение, вообще-то, расходится в области больших длин волн (или малых k). Но его можно, как обычно, разделить на сумму конечной части, зависящей от Δy и Δt, и бесконечной, не зависящей от этих переменных:

$$C(\Delta y,\Delta t) = const -\frac{\hbar }{4\pi \sigma c}\left(\ln \mbox{sh} \frac{\pi kT}{\hbar }\left(\Delta t-\frac{\Delta y}{c}\right)+ \ln \mbox{sh} \frac{\pi kT}{\hbar }\left(\Delta t+\frac{\Delta y}{c}\right)\right).$$

Корреляционная функция в фиксированной точке (Δy = 0) принимает вид

$$C(0,\Delta t)=const-\frac{\hbar }{2\pi \sigma c}\ln \mbox{sh}\, \frac{\pi kT\Delta t}{\hbar }.$$(1)

Еще нам понадобится корреляционная функция при нулевой температуре

$$C_{0}(\Delta y,\Delta t)\equiv {\frac{\hbar }{4\pi \sigma }}\int\limits_{-\infty}^{\infty }\frac{dk}{\omega }\cos (k\Delta y-\omega \Delta t).$$

Здесь нужно выделять конечную часть по-другому. После преобразований получается

$$C_{0}(\Delta y,\Delta t)=const-{\frac{\hbar }{4\pi \sigma c}}\ln \left\vert\Delta t^{2}-\frac{\Delta y^{2}}{c^{2}}\right\vert.$$(2)

Равноускоренное движение

Движение под действием постоянной силы F описывается в СТО известным уравнением

$${\frac{d}{dt}}{\frac{mv}{\sqrt{1-v^{2}/c^{2}}}}=F,$$

где под скоростью v понимается dy/dt. Его решение легко найти:

$$y={\frac{mc^{2}}{F}} \, \mbox{ch}\, \frac{F\tau }{mc},\quad t={\frac{mc}{F}}\,\mbox{sh}\, \frac{F\tau }{mc},\quad -\infty <\tau <\infty .$$

Параметр τ совпадает с собственным временем $$\inline\int\!dt\sqrt{1-v^{2}/c^{2}}}$$. Отсюда для двух точек на мировой линии можно получить, что

$$\sqrt{\Delta t^{2}-\Delta y^{2}/c^{2}}={\frac{2mc}{F}}\,\mbox{sh}\, \frac{F\Delta\tau }{2mc},$$(3)

где Δτ = τ1 − τ2.

Температура Унру

Из (2) и (3) получаем, что функция корреляции поля с нулевой температурой для точек вдоль мировой линии равноускоренного наблюдателя зависит только от Δτ:

$$C_{0}(\Delta y,\Delta t)=const-{\frac{\hbar }{2\pi \sigma c}}\ln \mbox{sh}\, \frac{F\Delta \tau }{2mc}.$$

Сравнение последнего выражения с (1) показывает, что поле с нулевой температурой будет выглядеть для движущегося равноускоренно наблюдателя так, как будто обладает температурой Унру

$$kT={\frac{\hbar F}{2\pi mc}}.$$

Выводы

Отметим, что температура Унру очень мала. Так, для ускорения, совпадающего с ускорением свободного падения, температура Унру равна 4·10−20 К.

Эффект Унру в некотором смысле аналогичен излучению Хокинга. Действительно, для равноускоренного наблюдателя существует так называемый риндлеровский горизонт, аналогичный горизонту событий черной дыры.

Были предложения проверить эффект Унру, наблюдая дополнительное излучение за счет тепловых флуктуаций ускоренно движущейся частиц, например, электронов, освещенных мощными лазерами. Однако ряд авторов опровергает наличие дополнительного излучения, заявляя о компенсации возможного испускания поглощением энергии вакуумных (уже теплых!) полей. Например, далее в упомянутой статье разбирается пример осциллятора, связанного со скалярным полем, и прямым вычислением показывается отсутствие излучения.

Несмотря на малую величину, эффект Унру имеет важное философское значение. Действительно, этот эффект позволяет в принципе определить абсолютное ускорение системы отсчета. Таким образом опровергается принцип Маха в формулировке, утверждающей, что «имеет значение только ускорение относительно неподвижных звезд».

Ключевые слова: равновесное излучение, квантовая теория поля, гравитация | Комментарии (20)

← сюда