Коллапсирующие сферы
Статья написана по мотивам обсуждений на разных форумах. Принципиально новых результатов здесь нет, идёт только суммирование и (очень частный) обзор известных результатов.
В основном, статья может представлять, равзе что, методический интерес.
Статья незакончена и, возможно, и не будет закончена. На текущий момент остаются следующие вопросы:
- досчитать сферу и частицу
- внешний вид сферы
- случай двух сфер
- разобраться с "запредельной" сферой
Введение
В статье рассматриваются массивные сферы с позиций ОТО.
В основном, используются Шварцшильдовские координаты, что, вообще говоря, требует отдельного рассмотрения продолжения решений под горизонт.
Хотя в отдельных случаях возможно задание единой системы координат во всём ПВ, в общем случае это либо невозможно, либо приводит к "катастрофическому" усложнению формул. Поэтому, будем придерживаться формализма, изложенного в главе 2 статьи , в и в .
Будут использоваться именно стандартные координаты, так как это позволяет (ввиду совпадения длин окружностей) ввести "единые" пространственные координаты.
Краткое изложение формализма в общем случае
Предполагается, что в ПВ имеется времениподобная поверхность \(\Sigma\), разделяющая ПВ на две области \(M_+\) и \(M_-\) в каждой из которых заданы различные системы координат \(X_+^i\) и \(X_-^i\). Требованием, связывающим эти СК, является совпадение метрик, индуцированных из каждой области на \(\Sigma\).
Считаем, что поверхность \(\Sigma\) параметризована внутренними координатами \(\zeta^\alpha\) и, соответственно, \[X_\pm^i=X_\pm^i(\zeta^\alpha)\]
Тогда, касательные вектора задаются соотношениями \[{e_\pm}_\alpha^i=\dpfs{X_\pm^i}{\zeta^\alpha}\label{covect}\]
Скалярное произведение этих векторов определяет индуцированную на поверхности метрику (первая фундаментальная форма) \[{g_\pm}_{\alpha\beta}={g_\pm}_{ij}{e_\pm}_\alpha^i{e_\pm}_\beta^j\label{metc}\]
Обе метрики должны совпадать в каждой точке поверхности \(\Sigma\): \({g_+}_{\alpha\beta}(\zeta)={g_-}_{\alpha\beta}(\zeta)=g_{\alpha\beta}(\zeta)\)
Тензор внешней кривизны (вторая фундаментальная форма) определяется через единичный ковектор \(n\), нормальный к поверхности, следующим образом (заметим, что направление этого ковектора должно быть одинаковым с обеих сторон поверхности) \[{K_\pm}_{\alpha\beta}={{n_\pm}_i}_{;j}{e_\pm}_\alpha^i{e_\pm}_\beta^j=-{n_\pm}_i{{e_\pm}_\alpha^i}_{;j}{e_\pm}_\beta^j\]
Тензор энергии-импульса поверхности (трёхмерный!) определяется через скачок тензора кривизны поверхности \([K_{\alpha\beta}]={K_+}_{\alpha\beta}-{K_-}_{\alpha\beta}\) и \([K]=g^{\alpha\beta}[K_{\alpha\beta}]\) уравнением Ланцоша (Lanczos) \[S_{\alpha\beta}=\frac{1}{8\pi}\left([K_{\alpha\beta}]-g_{\alpha\beta}[K]\right)\]
Последнее соотношение позволяет придать физически осмысленный вид задаче. В частности, "просто поверхность" должна иметь нулевой \(S_{\alpha\beta}\) и, соответственно, нулевой \([K_{\alpha\beta}]\), то есть совпадающие тензоры внешней кривизны.
Пример гравитирующей плоскости
Рассмотрим в декартовых координатах в области \(x\gt0\) метрику \[\dd{s}^2=(1+ax)^2\dd{t}^2-\dd{x}^2-\dd{y}^2-\dd{z}^2\]
Пусть при \(x=0\) имеется поверхность, параметризованная координатами \(t,y,z\). Очевидно, касательные вектора будут равны \[e^i_a=\Bmatrix{\pmatrix{1\\0\\0\\0}&\pmatrix{0\\0\\1\\0}&\pmatrix{0\\0\\0\\1}}\] а индуцированный метрический тензор \[g_{\alpha\beta}=\pmatrix{1&0&0\\0&-1&0\\0&0&-1}\]
Выберем направление ковектора нормали \(n_i\) наружу из нашей области. Тогда этот ковектор будет равен \[n_i=\pmatrix{0&-1&0&0}\] Его производная будет \[n_{i;j}=n_{i,j}-n_k\Gamma^k_{ij}=\pmatrix{a&0&0&0\\0&0&0&0\\0&0&0&0\\0&0&0&0}\] Откуда тензор внешней кривизны будет \[K_{\alpha\beta}=\pmatrix{a&0&0\\0&0&0\\0&0&0}\]
По другую сторону поверхности мы можем взять "симметричные" координаты и учесть, что направление нормального ковектора должно быть обратное. Или продолжить имеющиеся координаты, но в метрике изменить знак \(a\). Результат будет одинаковый. При этом очевидно, что индуцированная метрика будет та же самая. А для тензора кривизны будем иметь \[[K_{\alpha\beta}]=\pmatrix{2a&0&0\\0&0&0\\0&0&0}\\ [K]=2a\] Окончательно, \[S_{\alpha\beta}=\frac{a}{4\pi}\pmatrix{0&0&0\\0&1&0\\0&0&1}\] что означает нулевую плотность и пропорциональное ускорению давление.
Второй возможный вариант плоской статической вакуумной метрики есть \[\dd{s}^2=(1+ax)^{-2/3}\dd{t}^2-\dd{x}^2-(1+ax)^{4/3}(\dd{y}^2+\dd{z}^2)\]
Соответственно изменится производная \[n_{i;j}=\pmatrix{-\frac{a}{3}&0&0&0\\0&0&0&0\\0&0&-\frac{2a}{3}&0\\0&0&0&-\frac{2a}{3}}\] и \[K_{\alpha\beta}=\pmatrix{-\frac{a}{3}&0&0\\0&-\frac{2a}{3}&0\\0&0&-\frac{2a}{3}}\] Соответственно \[S_{alpha\beta}=\frac{1}{12\pi}\pmatrix{-4&0&0\\0&1&0\\0&0&1}\]
Движение массивной (заряженной) сферы в сферически симметричном поле
История движения сферы в ПВ как раз представляет собой пример поверхности \(\Sigma\), описанной в предыдущем разделе. Поскольку в данном случае мы имеем дело только с двумя областями, сохраним для их обозначения знаки \(\pm/+/-\). Величины без этих знаков будут либо относится к самой поверхности, либо быть тождественно равными для обеих областей.
Метрика вне сферы даётся теоремой Биркгофа.
При анализе будут использоваться сферические координаты. Это позволяет использовать угловые координаты общие для всего ПВ (и не использвать для них метки области).
Рассмотрим метрику вне сферы в виде \[\dd{s}^2=f_\pm\dd{t_\pm}^2-f_\pm^{-1}\dd{r_\pm}^2-r_\pm^2\dd{\mathbf\Omega}^2,\] где \[f_\pm(r)=1-\frac{2M_\pm}{r}+\frac{Q_\pm^2}{r^2}-H^2r^2,\\ f'_\pm(r)=\frac{2}{r^3}\left(M_\pm r-Q_\pm^2-H^2r^4\right)\] а \(M_\pm\) - параметр массы для соответствующей области, \(Q_\pm\) - заряд, а \(H\) - постоянная Хаббла (\(\sqrt{\Lambda/3}\)).
Параметризуем сферу координатами \((\tau,\theta,\varphi)\), где \(\tau\) - собственное время частиц на сфере. Точкой будет обозначаться производная по \(\tau\). Во внешних областях координаты сферы будут задаваться соотношениями \[r_\pm=R(\tau),\;t_\pm=t_\pm(\tau),\;\theta_\pm=\theta,\;\varphi_\pm=\varphi\] (в стандартных координатах с учётом условия \eqref{metc} координата \(r\) будет непрерывна на сфере)
Тогда, из \eqref{covect} имеем \begin{align} {e_\pm}_\tau^i &= (\dot{t_\pm},\dot{R},0,0)={u_\pm}^i\label{equ}\\ {e_\pm}_\theta^i &= (0,0,1,0)\\ {e_\pm}_\varphi^i &= (0,0,0,1) \end{align}
и нормальный единичный ковектор \[{n_\pm}_i=(\dot{R},-\dot{t_\pm},0,0)\label{eqn}\]
Заметим, что вектора \(u\) и \(n\) имеют довольно простой вид в силу того, что для параметризации поверхности было выбрано собственное время, и, следовательно, \[\dot{t_\pm}^2F_\pm-\dot{R}^2/F_\pm=1,\label{tau}\] где введено обозначение \(F_\pm(\tau)=f_\pm(R(\tau))\). Выполнение этих уравнений автоматически приводит к совпадению индуцированных метрик.
Введём полезную для дальнейшего ф-цию \(\beta_\pm(\tau)=\sqrt{\dot{R}^2+F_\pm}\). Тогда, можно переписать уравнения \eqref{tau} в виде \[\dot{t_\pm}F_\pm=\beta_\pm.\] Знак в этих уравнениях выбран таким образом, что бы соответствовать "естественным" доколлапсовым начальным условиям.
На всякий случай, явно выпишем \[\dot\beta_\pm=\frac{\dot{R}}{2\sqrt{\dot{R}^2+F_\pm}}\left(2\ddot{R}+{f_\pm}'\right)\]
Вычислим тензоры кривизны. \begin{alignat}{1} {K_\pm}^\tau_\tau\; &= -{n_\pm}_i{u_\pm}_{;j}^i{u_\pm}^j\; &= -{n_\pm}_i(\dot{{u_\pm}^i}+{\Gamma_\pm}^i_{kj}{u_\pm}^k{u_\pm}^j)\; &= \dot{\beta_\pm}/\dot{R}\\ {K_\pm}_\theta^\theta\; &= -{{n_\pm}_\theta}_{;\theta}/R^2\; &= {\Gamma_\pm}^r_{\theta\theta}n_r/R^2\; &= \beta_\pm/R\\ {K_\pm}_\varphi^\varphi\; &= -{{n_\pm}_\varphi}_{;\varphi}/R^2/cos^2\theta\; &= {\Gamma_\pm}^r_{\varphi\varphi}n_r/R^2/cos^2\theta\; &= \beta_\pm/R \end{alignat} (причём заметим, что эти уравнения справедливы для ф-ции \(f(r)\) любого вида).
Считая сферу идеальной жидкостью (\(S^\alpha_\beta=(p+\sigma)u^\alpha u_\beta-p \delta^\alpha_\beta\)) получим следующие соотношения \begin{alignat}{1} S^\tau_\tau\; &= \sigma\; &= -\frac{\beta_+-\beta_-}{4\pi R}\label{eqs}\\ S^\theta_\theta = S^\varphi_\varphi\; &= -p\; &= \frac{\dot{\beta_+}-\dot{\beta_-}}{8\pi \dot{R}} + \frac{\beta_+-\beta_-}{8\pi R}\label{eqp} \end{alignat} где \(\sigma\) - плотность сферы, а \(p\) - давление в ней.
Отметим, что в статическом случае будет \[\frac{\dot{\beta_\pm}}{\dot{R}}=\frac{{f_\pm}'}{2\sqrt{F_\pm}}=\frac{1}{R^2}\frac{M_\pm R-Q_\pm^2-H^2R^4}{\sqrt{R^2-2M_\pm R+Q_\pm^2-H^2R^4}}\label{stat_beta}\]
В частном случае пылевой сферы (\(p=0\)) из уравнения \eqref{eqp} получаем \((\beta_+-\beta_-)R = const\) и подставляя это выражение в \eqref{eqs} получаем "закон сохранения массы" \[4\pi R^2\sigma=m=const\label{eqm}\] и соответственно \[\beta_+=\beta_--\frac{m}{R}\]
Возводя в квадрат получим \[F_+-F_-=-2\sqrt{\dot{R}^2+F_-}\frac{m}{R}+\frac{m^2}{R^2}\] \[M_+-M_-=m\sqrt{\dot{R}^2+1-\frac{2M_-}{R}+\frac{Q_-^2}{R^2}-H^2R^2}-\frac{m^2}{2R}+\frac{Q_+^2-Q_-^2}{2R}\label{eqM}\] \[\dot{R}^2=\frac{(F_+-F_-)^2R^2}{4m^2}-\frac{(F_++F_-)}{2}+\frac{m^2}{4R^2}\] \[\dot{t_+}=\frac{1}{F_+}\sqrt{\dot{R}^2+F_+}=\frac{R^2(F_--F_+)-m^2}{2mRF_+}\] \[\dot{t_-}=\frac{1}{F_-}\sqrt{\dot{R}^2+F_-}=\frac{R^2(F_--F_+)+m^2}{2mRF_-}\] \[\dot{R}^2=\frac{(M_+-M_-)^2}{m^2}+\frac{(Q_+^2-Q_-^2)^2}{4m^2R^2}+\frac{m^2}{4R^2}+\frac{M_++M_-}{R}-\frac{Q_+^2+Q_-^2}{2R^2}-1+H^2R^2-\frac{(M_+-M_-)(Q_+^2-Q_-^2)}{m^2R}\] и далее для случая \(H=Q_-=Q_+=0\) \begin{alignat}{1} \dot{R}\; &=\pm\sqrt{\frac{(M_+-M_-)^2}{m^2}+\frac{M_++M_-}{R}+\frac{m^2}{4R^2}-1}\; &=\frac{\pm1}{m R}\sqrt{((M_+-M_-)^2-m^2)R^2+m^2(M_++M_-)R+m^4/4}\label{eqDu}\\ \dot{t_+}\; &=\frac{R}{R-2M_+}\sqrt{\dot{R}^2+1-\frac{2M_+}{R}}\; &=\frac{(M_+-M_-)R-m^2/2}{m(R-2M_+)}\\ \dot{t_-}\; &=\frac{R}{R-2M_-}\sqrt{\dot{R}^2+1-\frac{2M_-}{R}}\; &=\frac{(M_+-M_-)R+m^2/2}{m(R-2M_-)} \end{alignat}
При \(M_-=Q_-=H=0\) формула \eqref{eqM} имеет "наглядный" физический смысл: гравитационная масса равна кинетическая энергия (включая массу покоя) минус энергия связи плюс энергия электрического поля. При этом (\(Q_+=0\)), при "параболическом" движении (покой на бесконечности) гравитационная масса равна массе покоя, при "гиперболическом" - больше, а при "эллиптическом" соответственно меньше. Аналогично можно классифицировать движение и в общем случае.
Детали случая сжатия пылевой (незаряженной) сферы
Рассмотрим в подробностях случай \(p=0,\;M_-=Q_-=0\). Для удобства обозначим \(M_+=M\), \(q=\frac{Q}{m}\) и \(\gamma=\frac{M}{m}\), тогда \(\gamma\gt 1\) будет означать гиперболический случай и \(\gamma\lt 1\) - эллиптический.
Выпишем \[F_-=1-H^2R^2\] \[F_+=1-\frac{2M}{R}+\frac{q^2M^2}{\gamma^2R^2}-H^2R^2\] \[R^2(F_--F_+)=2MR-\frac{q^2M^2}{\gamma^2}\]
Тогда динамика сферы описывается следующими уравнениями: \begin{alignat}{1} \dot{R} &=-\sqrt{\gamma^2+\frac{M^2(1-q^2)^2}{4\gamma^2R^2}+\frac{M(1-q^2)}{R}-1+H^2R^2}\; &= -\sqrt{\left(\gamma+\frac{M(1-q^2)}{2\gamma R}\right)^2-(1-H^2R^2)}\label{r1}\\ \dot{t_-} &= \frac{1}{1-H^2R^2}\left(\gamma+\frac{M(1-q^2)}{2\gamma R}\right)\label{tm1}\\ \dot{t_+} &= \frac{1}{1-\frac{2M}{R}+\frac{q^2M^2}{\gamma^2R^2}-H^2R^2}\left(\gamma-\frac{M(1+q^2)}{2\gamma R}\right)\label{tp1} \end{alignat}
В гиперболическом случае, при \(H=0\), на бесконечности имеем \[\dot{R}=-\sqrt{\gamma^2-1}\] \[\dot{t_\pm}=\gamma\] и соответственно \[v_\infty=\sqrt{1-1/\gamma^2}\label{gamma}\]
В эллиптическом случае, опять же при \(H=0\), максимально достижимый радиус легко находится из \eqref{r1}: \[R_0=\frac{M(1-q^2)}{2\gamma(1-\gamma)}\] Легко видно, что при \(\gamma=(1+q^2)/2\) сфера не существует вне горизонта. Назовём этот случай предельным, так как при \(\gamma\lt (1+q^2)/2\) корректность нашего исходного построения оказывается под вопросом.
Уравнение \eqref{r1} легче всего интегрируется в параболическом случае при \(H=q=0\), но даже в этом случае мы получаем кубическое уравнение для \(R(\tau)\): \[(M-2R)\sqrt{M+4R}=6\sqrt{M}(\tau_0\pm\tau)\notag\]
Во всех допредельных случаях с точки зрения внешнего наблюдателя скорость сферы стремится к 0 по мере того как сфера приближается к горизонту. С точки зрения внутреннего наблюдателя прохождение сферой горизонта ничем особенным не выделено, и при \(R\) стремящемся к 0, скорость стремится к 1.
Далее будем полагать \(H=q=0\).
Рассмотрим чуть подробнее движение сферы вблизи горизонта.
Уравнение \eqref{r1} не имеет никаких особенностей при \(R=2M\), поэтому сделаем замену \(R=M(2+r)\) и проинтегрируем в предположении малости \(r\): \[-\dd{\tau}=\frac{M(2+r)\dd{r}}{\sqrt{(\gamma^2-1)(2+r)^2+2+r+\frac{1}{4\gamma^2}}}=M\left(\frac{4\gamma}{4\gamma^2-1}+\frac{2\gamma(4\gamma^2+1)}{(4\gamma^2-1)^3}r+O(r^2)\right)\dd{r}\approx\frac{4\gamma M}{4\gamma^2-1}\dd{r}\] откуда \[R/M=2+r\approx2-\frac{4\gamma^2-1}{4\gamma M}\tau\label{r2}\] и \begin{alignat}{1} \dot{t_-} &\approx\frac{4\gamma^2+1}{4\gamma} &\implies t_-\approx\frac{4\gamma^2+1}{4\gamma}\tau\\ \dot{t_+} &\approx\frac{2(4\gamma^2-1)}{4\gamma r}\approx\frac{-2M}{\tau} &\implies t_+\approx\left.-2M\ln(-\tau/M)\right|_{\tau\lt0}\label{tp2} \end{alignat}
Находя \(\tau\) из \eqref{tp2} и подставляя в \eqref{r2} получим \[R/M\approx2+\frac{4\gamma^2-1}{4\gamma}e^{-t_+/(2M)}\] Т.е., с точки зрения внешнего наблюдателя, сфера экспоненциально быстро приближается к горизонту.
Рассмотрим двух покоящихся наблюдателей непосредственно внутри и снаружи сферы.
С точки зрения этих наблюдателей движение сферы описывается следующими формулами: \begin{alignat}{1} \dd{l_-} &= \dd{r_-} &= \dot{R}\dd{\tau}\\ \dd{\tau_-} &= \dd{t_-} &= \dot{t_-}\dd{\tau}\\ v_- &= \dot{R}/\dot{t_-} &= -\frac{\sqrt{(\gamma^2-1)\gamma^2R^2+\gamma^2M R+M^2/4}}{\gamma^2R+M/2}\\ \dd{l_+} &= 1/\sqrt{F_+}\dd{r_+} &= \dot{R}/\sqrt{F_+}\dd{\tau}\\ \dd{\tau_+} &= \sqrt{F_+}\dd{t_+} &= \sqrt{F_+}\dot{t_+}\dd{\tau}\\ v_+ &= \dot{R}/(\dot{t_+}F_+) &= -\frac{\sqrt{(\gamma^2-1)\gamma^2R^2+\gamma^2M R+M^2/4}}{\gamma^2R-M/2} \end{alignat} Как нетрудно видеть, скорость сферы снаружи больше чем изнутри \[\frac{v_+}{v_-}=\frac{\gamma^2R+M/2}{\gamma^2R-M/2}\] а (радиальное) расстояние между парой неподвижных наблюдателей увеличивается как \[\ddfs{l_+}{l_-}=\frac{1}{\sqrt{F_+}}=\sqrt{\frac{R}{R-2M}}\]
ПВД
Рассмотрим \(\gamma=1\) и сделаем замену \(Q=\frac{1}{3}\sqrt{1+\frac{4R}{M}}\) \begin{alignat}{1} R\; &=M\frac{9Q^2-1}{4}\\ \dot{R}\; &=\frac{9MQ}{2}\dot{Q}\;&=-\frac{6Q}{9Q^2-1}\\ \dot{Q}\; &=-\frac{4}{3M(9Q^2-1)}\\ \dot{t_+}\; &=1+\frac{2}{3(Q^2-1)}\;&=\frac{3Q^2-1}{3(Q^2-1)}\;&=-\frac{M}{4}\left(27Q^2+15+\frac{16}{Q^2-1}\right)\dot{Q}\\ \dot{t_-}\; &=1+\frac{2}{9Q^2-1}\;&=\frac{9Q^2+1}{9Q^2-1}\;&=-\frac{3M}{4}(9Q^2+1)\dot{Q}\\ \end{alignat} \begin{alignat}{1} \tau\; &=\frac{3MQ}{4}(1-3Q^2)-\frac{3}{2}M\; &=\frac{1}{6}(M-2R)\sqrt{1+\frac{4R}{M}}-\frac{3}{2}M\\ t_+\; &=4M\acoth(Q)-\frac{3MQ}{4}(3Q^2+5)-M\; &=4M\acoth\left(\frac{1}{3}\sqrt{1+\frac{4R}{M}}\right)-\frac{4M+R}{3}\sqrt{1+\frac{4R}{M}}-M\\ t_-\; &=-\frac{3MQ}{4}(1+3Q^2)\; &=-\frac{1}{3}(M+R)\sqrt{1+\frac{4R}{M}} \end{alignat}
Графики зависимостей
Диаграмма Шварцшильда
ГС будущего
ГС прошлого
Последний сигнал
На этой диаграмме изображено ПВ внутри и снаружи коллапсирующей пылевой сферы (из покоя на бесконечности). Белая линия обозначает саму сферу. Область правее - внутренность сферы, левее - наружность. Горизонтальная ось соответствует "стандартной Шварцшильдовой" радиальной координате (точнее \(r/M\)). По вертикали события отмасштабированы таким образом, что бы во внутренней части "время шло равномерно".
Жёлтая линия (с "жирной точкой" в месте излома) соответствует горизонту событий (будущего).
Синяя линия (с "жирной точкой" в месте излома) соответствует последнему сигналу, на который удалённый наблюдатель может получить ответ.
Несколько графиков скоростей в зависимости от радиуса сферы
Приведены три семейства по три кривых, соответствующих гиперболическому (h), параболическому (p) и эллиптическому (e) случаям. Первое семейство (v-) описывает скорость (физическая совпадает с координатной) сферы с точки зрения внутреннего наблюдателя. Второе (v+) - "физическая" скорость сферы (измеренная наблюдателем, покоящимся непосредственно снаружи сферы). Третье (V+) - координатная скорость во внешней системе координат.
Статическая сфера
Скомбинируем уравнения \eqref{eqs}, \eqref{eqp} и \eqref{stat_beta} \[p=\frac{\sigma}{2}-\frac{1}{8\pi}\left(\frac{f_+'}{2\sqrt{F_+}}-\frac{f_-'}{2\sqrt{F_-}}\right)=\frac{\sigma}{2}-\frac{1}{8\pi R^2}\left(\frac{M_+R-Q_+^2-H^2R^4}{\sqrt{R^2-2M_+R+Q_+^2-H^2R^4}}-\frac{M_-R-Q_-^2-H^2R^4}{\sqrt{R^2-2M_-R+Q_-^2-H^2R^4}}\right)\] Кроме того, само уранение \eqref{eqs} даёт \[4\pi R\sigma=\sqrt{F_-}-\sqrt{F_+}=\sqrt{1-\frac{2M_-}{R}+\frac{Q_-^2}{R^2}-H^2R^2}-\sqrt{1-\frac{2M_+}{R}+\frac{Q_+^2}{R^2}-H^2R^2}\]
Считая внутренность сферы пустой (\(M_-=Q_-=0\)) и используя \eqref{eqm} и обозначения \(M_+=M\), \(q=\frac{Q_+}{m}\), \(\gamma=\frac{M}{m}\), \(\eta^2=1-H^2R^2\) и \(\mu=\frac{m}{R}\), получим условие равновесия пылевой сферы \[\mu+\frac{\eta^2-1}{\eta}=\frac{\gamma\mu-q^2\mu^2+\eta^2-1}{\sqrt{\eta^2-2\gamma\mu+q^2\mu^2}}\] и \[\mu=\eta-\sqrt{\eta^2-2\gamma\mu+q^2\mu^2}\] Выразим \(\gamma\) из последнего уравнения \[\gamma=\eta-\frac{1-q^2}{2}\mu\] и подставим в предыдущее \[\frac{1-q^2}{2}\mu=\frac{1-\eta^2}{\eta}\]
При нулевой лямбде получаем единственный возможный вариант пылевой сфере оставаться в статике это иметь удельный заряд \(q^2=1\) (что соответствует критической ЧД Рейснера-Нордстрёма).
\(\Pi=2\pi Rp\) \[\frac{4\Pi}{\mu}=1-\frac{(1-\mu/2)}{1-\mu}\]
Взаимодействие сферы с пробными частицами
Будем полагать, что частица движется в экваториальной плоскости.
Обозначим собственное время частицы как \(T\), а 4-скорость частицы через \(U_\pm^i\). Считая, что частица взаимодействует со сферой только гравитационно, получим следующие условия (??) на скорость частицы при пересечении сферы \[U_+^i{n_+}_i=U_-^i{n_-}_i\label{Un}\] и \[{g_+}_{ij}{U_+}^i{e_+}_\alpha^j={g_-}_{ij}{U_-}^i{e_-}_\alpha^j\label{Ue}\] Эти уравнения означают, что с точки зрения наблюдателя на сфере частица пролетает сквозь сферу без "искажения" траектории.
Из \eqref{Un} и \eqref{eqn} получаем \[U_+^t\dot{R}-U_+^r\dot{t_+}=U_-^t\dot{R}-U_-^r\dot{t_-}\] Из \eqref{Ue} и \eqref{equ} получим \[U_+^tF_+\dot{t_+}-U_+^r\dot{R}/F_+=U_-^tF_-\dot{t_-}-U_-^r\dot{R}/F_-\]
Домножая, вычитая и учитывая \eqref{tau} получим \[U_+^t=U_-^t(F_-\dot{t_-}\dot{t_+}-\dot{R}^2/F_+)+U_-^r\dot{R}(\dot{t_-}/F_+-\dot{t_+}/F_-)\] \[U_+^r=U_-^t\dot{R}(F_-\dot{t_-}-F_+\dot{t_+})+U_-^r(F_+\dot{t_+}\dot{t_-}-\dot{R}^2/F_-)\] Обратные соотношения получаются, очевидно, сменой знаков + и - (а детерминант матрицы перехода равен 1).
Рассмотрим для примера случай покоящейся частицы внутри пылевой сферы (\(M_-=0,\; p=0,\; M_+=1,\; U_-^t=1,\; U_-^r=0\)). Тогда, после пролёта сферы (считаем, что это происходит над горизонтом) 4-скорость частицы окажется \[(U_+^t,\; U_+^r)=\left(\frac{R^2-R-m^2/2}{R(R-2)},\; -\frac{\sqrt{(1-m^2)R^2+m^2R+m^4/4}}{R^2}\right)\] Или, скорость частицы по отношению к скорости (сжатия) сферы будет \[\frac{R-m^2/2}{R^2-R-m^2/2}\notag\]
В более общем случае пылевой сферы (\(p=0\)) будет \[U_+^t=U_-^t\frac{R^2-(M_++M_-)R-m^2/2}{R(R-2M_+)}+U_-^r\frac{\sqrt{((M_+-M_-)^2-m^2)R^2+m^2(M_++M_-)R+m^4/4}}{(R-2M_+)(R-2M_-)}\] \[U_+^r=U_-^t\frac{\sqrt{((M_+-M_-)^2-m^2)R^2+m^2(M_++M_-)R+m^4/4}}{R^2}+U_-^r\frac{R^2-(M_++M_-)R-m^2/2}{R(R-2M_-)}\] \[U_+^\theta=U_-^\theta,\quad U_+^\varphi=U_-^\varphi\] или, полагая \(M_-=0,\; M_+=1\), \[U_+^t=U_-^t\frac{R^2-R-m^2/2}{R(R-2)}+U_-^r\frac{\sqrt{(1-m^2)R^2+m^2R+m^4/4}}{R(R-2)}\] \[U_+^r=U_-^t\frac{\sqrt{(1-m^2)R^2+m^2R+m^4/4}}{R^2}+U_-^r\frac{R^2-R-m^2/2}{R^2}\] \[U_+^\theta=U_-^\theta,\quad U_+^\varphi=U_-^\varphi\]
Рассмотрим частицу со скоростью \(U_-^i=(1,V\cos\alpha,0,V/r\sin\alpha)/\sqrt{1-V^2}\) - т.е. движущуюся под углом \(\alpha\) к "радиусу". Тогда, скорость после пролёта сферы будет \[U_+^t=\frac{R^2-R-m^2/2+V\cos\alpha\sqrt{(1-m^2)R^2+m^2R+m^4/4}}{R(R-2)\sqrt{1-V^2}}\] \[U_+^r=\frac{\sqrt{(1-m^2)R^2+m^2R+m^4/4}+V\cos\alpha(R^2-R-m^2/2)}{R^2\sqrt{1-V^2}}\] \[U_+^\varphi=\frac{V\sin\alpha}{R\sqrt{1-V^2}}\]
Случай двух пылевых сфер
Во всей этой главе подразумевается \(H=Q=0\) во всех областях пространства и \(p=0\) на всех сферах.
Для определённости будем считать, что сферы "коллапсируют" (\(\dot{R}<0\)).
Кроме того, введём обозначение "производной по радиусу": \(X'=\dot{X}/\dot{R}\). В частности, \(\tau'=1/\dot{R}\) и \(X'=\dot{X}\tau'\).
Движение пылевой сферы с массой внутри
Повторим уравнение \eqref{eqDu} \[\dot{R}=-\frac{1}{R}\sqrt{((M_+-M_-)^2/m^2-1)R^2+(M_++M_-)R+m^2/4}\]
Обозначим \(\gamma=\frac{M_+-M_-}{m}\) и \(\kappa=\frac{M_++M_-}{m}\) и введём \(\rho=R/m\). Тогда \[m\rho\dot{\rho}=-\sqrt{(\gamma^2-1)\rho^2+\kappa\rho+1/4}\] или \[\tau'=-\frac{2\rho}{\sqrt{4(\gamma^2-1)\rho^2+4\kappa\rho+1}}\]
Уравнения для временных координат будут в этих обозначениях: \[\dot{t_+}=\frac{\gamma\rho-1/2}{\rho-(\kappa+\gamma)}\\ \dot{t_-}=\frac{\gamma\rho+1/2}{\rho-(\kappa-\gamma)}\] или \[t_+'=-\frac{2\gamma\rho-1}{\rho-(\kappa+\gamma)}\frac{\rho}{\sqrt{4(\gamma^2-1)\rho^2+4\kappa\rho+1}}\\ t_-'=-\frac{2\gamma\rho+1}{\rho-(\kappa-\gamma)}\frac{\rho}{\sqrt{4(\gamma^2-1)\rho^2+4\kappa\rho+1}}\]
Параболический случай \(\gamma=1\)
\[\tau'=-\frac{2\rho}{\sqrt{4\kappa\rho+1}}\]
Введём \(X=\sqrt{1+4\kappa\rho}\). \[\rho=\frac{X^2-1}{4\kappa}\\ X'=\frac{2\kappa}{mX}\] \[\frac{\tau'}{m}=-\frac{X^2-1}{4\kappa^2}X'\] \[6\kappa^2\frac{\tau}{m}=(1-2\kappa\rho)\sqrt{1+4\kappa\rho}-1\] (константа интегрирования выбрана так, что бы было \(R(0)=0\)).
\[t_+'=-\frac{2\rho-1}{\rho-(\kappa+1)}\frac{\rho}{\sqrt{4\kappa\rho+1}}=-m\frac{X^2-1-2\kappa}{X^2-1-4\kappa^2-4\kappa}\frac{X^2-1}{4\kappa^2}X'\\ t_-'=-\frac{2\rho+1}{\rho-(\kappa-1)}\frac{\rho}{\sqrt{4\kappa\rho+1}}=-m\frac{X^2-1+2\kappa}{X^2-1-4\kappa^2+4\kappa}\frac{X^2-1}{4\kappa^2}X'\] \[6\kappa^2\frac{t_+-{t_+}_0}{m}=12\kappa^2(\kappa+1)\acoth\frac{\sqrt{1+4\kappa\rho}}{2\kappa+1}-(2\kappa\rho+6\kappa^2+3\kappa-1)\sqrt{1+4\kappa\rho}\\ 6\kappa^2\frac{t_--{t_-}_0}{m}=12\kappa^2(\kappa-1)\acoth\frac{\sqrt{1+4\kappa\rho}}{2\kappa-1}-(2\kappa\rho+6\kappa^2-3\kappa-1)\sqrt{1+4\kappa\rho}\]
Гиперболический случай \(\gamma\gt1\)
\[m\rho\dot{\rho}=-\sqrt{(\gamma^2-1)\rho^2+\kappa\rho+1/4}=-\sqrt{\gamma^2-1}\sqrt{\rho^2+\frac{\kappa}{\gamma^2-1}\rho+\frac{1}{4(\gamma^2-1)}}\] \[m\rho\dot{\rho}=-\sqrt{\gamma^2-1}\sqrt{(\rho+\frac{\kappa}{2(\gamma^2-1)})^2+\frac{\gamma^2-1-\kappa^2}{4(\gamma^2-1)^2}}\]
Делаем замену \(x=\rho+\frac{\kappa}{2(\gamma^2-1)}\): \[2x\dot{x}-\frac{\kappa}{(\gamma^2-1)}\dot{x}=-\frac{2\sqrt{\gamma^2-1}}{m}\sqrt{x^2-\frac{\kappa^2-\gamma^2+1}{4(\gamma^2-1)^2}}\] \[\frac{\dot{x^2}}{\sqrt{x^2-\frac{\kappa^2-\gamma^2+1}{4(\gamma^2-1)^2}}}-\frac{\frac{\kappa}{\gamma^2-1}\dot{x}}{\sqrt{x^2-\frac{\kappa^2-\gamma^2+1}{4(\gamma^2-1)^2}}}=-2\frac{\sqrt{\gamma^2-1}}{m}\] \[\sqrt{\rho^2+\frac{\kappa}{\gamma^2-1}\rho+\frac{1}{4(\gamma^2-1)}}-\frac{\kappa}{\gamma^2-1}\ln\left(x\frac{2(\gamma^2-1)}{\sqrt{\kappa^2-\gamma^2+1}}+\sqrt{(x\frac{2(\gamma^2-1)}{\sqrt{\kappa^2-\gamma^2+1}})^2-1}\right)=C-2\frac{\sqrt{\gamma^2-1}\tau}{m}\] \[\sqrt{4(\gamma^2-1)\rho^2+4\kappa\rho+1}-\frac{2\kappa}{\sqrt{\gamma^2-1}}\ln\left(2\rho\sqrt{\gamma^2-1}+\frac{\kappa}{\sqrt{\gamma^2-1}}+\sqrt{4(\gamma^2-1)\rho^2+4\kappa\rho+1}\right)=1-\frac{2\kappa}{\sqrt{\gamma^2-1}}\ln\left(\frac{\kappa}{\sqrt{\gamma^2-1}}+1\right)-4(\gamma^2-1)\frac{\tau}{m}\]
\[m\dot{\rho}|_{\rho=\infty}=-\sqrt{\gamma^2-1}\]
Эллиптический случай \(\gamma\lt1\)
Решается аналогично предыдущему случаю, только вместо \(\acosh\) будет \(\arccos\). Поэтому выпишем сокращённо.
\[m\rho\dot{\rho}=-\sqrt{\kappa\rho+1/4-(1-\gamma^2)\rho^2}=-\sqrt{1-\gamma^2}\sqrt{\frac{\kappa}{1-\gamma^2}\rho+\frac{1}{4(1-\gamma^2)}-\rho^2}\] \[m\rho\dot{\rho}=-\sqrt{1-\gamma^2}\sqrt{\frac{1-\gamma^2+\kappa^2}{4(1-\gamma^2)^2}-(\rho-\frac{\kappa}{2(1-\gamma^2)}\rho)^2}\]
Делаем ту же самую замену \(x=\rho-\frac{\kappa}{2(1-\gamma^2)}\): \[\frac{\dot{x^2}}{\sqrt{\frac{1-\gamma^2+\kappa^2}{4(1-\gamma^2)^2}-x^2}}+\frac{\kappa}{1-\gamma^2}\frac{\dot{x}}{\sqrt{\frac{1-\gamma^2+\kappa^2}{4(1-\gamma^2)^2}-x^2}}=-\frac{2\sqrt{1-\gamma^2}}{m}\] \[\sqrt{4\kappa\rho+1-4(1-\gamma^2)\rho^2}+\frac{2\kappa}{\sqrt{1-\gamma^2}}\arccos\left(\frac{2(1-\gamma^2)\rho-\kappa}{\sqrt{1-\gamma^2+\kappa^2}}\right)=1+\frac{2\kappa}{\sqrt{1-\gamma^2}}\arccos\left(\frac{-\kappa}{\sqrt{1-\gamma^2+\kappa^2}}\right)+4(1-\gamma^2)\frac{\tau}{m}\]
\[\rho_{max}=\frac{\kappa+\sqrt{1-\gamma^2+\kappa^2}}{2(1-\gamma^2)}\]
Параболическое движение двух сфер
Для простоты рассмотрения ограничимся случаем параболического движения обеих сфер.
В данном случае мы имеем две поверхности, 3 области и 4 стороны поверхностей. Введём кодировку: всё внутреннее будем обозначать индексом "-", всё внешнее - "+", промежуточное - "*". Стороны сфер будем обозначать двумя индексами, где первый соответствует сфере. И обозначим \(k=\frac{m_-}{m_+}\).
То есть, \[M_-=0,\ M_*=\gamma_-m_-=m_-,\ M_+=\gamma_+m_++\gamma_-m_-=m_++m_-\\ \kappa_-=\gamma_-=1,\ \kappa_+=\gamma_++2\gamma_-k=1+2k\] \[6\frac{\tau_-}{m_-}=(1-2\rho_-)\sqrt{1+4\rho_-}-1\\ 6(1+2k)^2\frac{\tau_+}{m_+}=(1-2(1+2k)\rho_+)\sqrt{1+4(1+2k)\rho_+}-1\] \[3\frac{t_{-+}-{t_{-+}}_0}{m_-}=12\acoth\frac{\sqrt{1+4\rho_-}}{3}-(\rho_-+4)\sqrt{1+4\rho_-}\\ 3\frac{t_{--}-{t_{--}}_0}{m_-}=-(\rho_-+1)\sqrt{1+4\rho_-}\\ 3(1+2k)^2\frac{t_{++}-{t_{++}}_0}{m_+}=6(1+2k)^2(2+2k)\acoth\frac{\sqrt{1+4(1+2k)\rho_+}}{3+4k}-((1+2k)\rho_++4+15k+12k^2)\sqrt{1+4(1+2k)\rho_+}\\ 3(1+2k)^2\frac{t_{+-}-{t_{+-}}_0}{m_+}=12k(1+2k)^2\acoth\frac{\sqrt{1+4(1+2k)\rho_+}}{1+4k}-((1+2k)\rho_++1+9k+12k^2)\sqrt{1+4(1+2k)\rho_+}\]
Для дальнейшего упрощения примем \(k=1\) (то есть, сферы имеют одинаковую массу, которую будем обозначать \(m\)). \[6\frac{\tau_-}{m}=(1-2\rho_-)\sqrt{1+4\rho_-}-1\\ 54\frac{\tau_+}{m}=(1-6\rho_+)\sqrt{1+12\rho_+}-1\] \[3\frac{t_{-+}-{t_{-+}}_0}{m}=12\acoth\frac{\sqrt{1+4\rho_-}}{3}-(\rho_-+4)\sqrt{1+4\rho_-}\\ 3\frac{t_{--}-{t_{--}}_0}{m}=-(\rho_-+1)\sqrt{1+4\rho_-}\\ 27\frac{t_{++}-{t_{++}}_0}{m}=216\acoth\frac{\sqrt{1+12\rho_+}}{7}-(3\rho_++31)\sqrt{1+12\rho_+}\\ 27\frac{t_{+-}-{t_{+-}}_0}{m}=108\acoth\frac{\sqrt{1+12\rho_+}}{5}-(3\rho_++22)\sqrt{1+12\rho_+}\]
Единственный вопрос, который надо решить, состоит в том, что \(t_{+-}\) и \(t_{-+}\) относятся к одной области, и соответствующие "константы интегрирования" должны быть согласованы. Но, на самом деле, это вопрос "начальных данных" - каковы радиусы сфер в некий "нулевой момент" в промежуточной области.