Ниже я разберу статью Babichev, Emond, Ramazanov, “Shrouded black holes in Einstein–Gauss–Bonnet gravity” как материал для доклада: сначала физическая идея, затем модель, после этого — уравнения и их вывод, аналитические режимы, численные решения, космология и физический смысл результатов. Авторы рассматривают скалярное поле с квадратичной связью с инвариантом Гаусса–Бонне и показывают, что около чёрной дыры может возникать тонкая скалярная «стенка», отделяющая две фазы поля. (2207.03944v2.pdf)

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

---

1. В чём главная идея статьи

Обычная ОТО описывает чёрную дыру Шварцшильда без какого-либо дополнительного скалярного поля вокруг неё. Авторы задают вопрос:

Что произойдёт, если в теорию гравитации добавить скалярное поле $\phi$, которое имеет собственный потенциал со спонтанным нарушением $Z_2$-симметрии и одновременно чувствует сильную кривизну через инвариант Гаусса–Бонне?

Ответ необычный.

Далеко от чёрной дыры кривизна мала, и скаляр находится в одном из минимумов

$$\phi=+v \quad\text{или}\quad \phi=-v.$$

Авторы выбирают $+v$.

Но около горизонта кривизна огромна. Связь

$$\ell^2\phi^2\mathcal G$$

с инвариантом Гаусса–Бонне $\mathcal G$ изменяет эффективный потенциал так, что около чёрной дыры минимум может переместиться к

$$\phi=0.$$

Следовательно,

$$\phi\simeq0 \qquad \text{возле горизонта},$$

а

$$\phi\rightarrow v \qquad \text{далеко от чёрной дыры}.$$

Поле обязано перейти между этими двумя состояниями. Этот переход происходит в оболочке — scalar wall, скалярной стенке.

Именно поэтому в названии статьи чёрная дыра называется shrouded, то есть буквально «окутанной» скалярной оболочкой. Авторы подчёркивают, что это отличается от обычной спонтанной скаляризации: здесь $Z_2$-симметрия уже нарушена далеко от объекта и восстанавливается в сильном гравитационном поле. (2207.03944v2.pdf)

---

2. Зачем вообще нужна такая модель

Мотивация состоит из нескольких частей.

2.1. Проверка ОТО в сильном поле

В Солнечной системе ОТО проверена очень хорошо. Поэтому особенно интересны модифицированные теории, которые почти совпадают с ОТО при слабой гравитации, но дают новый эффект рядом с чёрными дырами.

Наблюдения LIGO/Virgo/KAGRA как раз позволяют исследовать такой режим. Это одна из основных мотиваций авторов. (2207.03944v2.pdf)

2.2. Скаляризация

До этой работы было известно, что связь

$$\phi^2\mathcal G$$

может сделать решение $\phi=0$ неустойчивым рядом с компактным объектом. Тогда чёрная дыра самопроизвольно приобретает скалярные «волосы».

Но здесь авторы делают практически противоположную конструкцию.

Они вводят потенциал Хиггсовского типа, который сам по себе хочет

$$\phi=\pm v,$$

а кривизна около чёрной дыры, наоборот, хочет

$$\phi=0.$$

Поэтому возникает переход между двумя фазами.

2.3. Космологическая проблема обычной скаляризации

Есть ещё одна важная причина выбора знака связи.

При знаке, обычно используемом для скаляризации статических чёрных дыр, та же связь с Гауссом–Бонне может создавать огромную тахионную неустойчивость во время инфляции.

Авторы выбирают противоположный знак. Тогда инфляция становится устойчивой, хотя новая неустойчивость появляется позже, в радиационную эпоху. (2207.03944v2.pdf)

---

3. Действие теории: уравнение (1)

Модель задаётся действием

$$S= \int d^4x\,\sqrt{-g} \left[ \frac12 R -\frac12\nabla_\mu\phi\nabla^\mu\phi -V(\phi) -\ell^2\phi^2\mathcal G \right], \tag{1}$$

где

$$\mathcal G \equiv R^2-4R_{\mu\nu}R^{\mu\nu} +R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}.$$

В статье этот объект обозначен $R^2_{\rm GB}$. Это не просто квадрат кривизны $R$, а специальная комбинация кривизн — инвариант Гаусса–Бонне. (2207.03944v2.pdf)

Разберём действие по частям.

Первый член

$$\frac12R$$

— стандартное действие Эйнштейна–Гильберта.

Используются единицы

$$8\pi G_N=1.$$

Второй член

$$-\frac12(\nabla\phi)^2$$

— кинетическая энергия скаляра.

Третий

$$-V(\phi)$$

— собственный потенциал поля.

Последний

$$-\ell^2\phi^2\mathcal G$$

— связь поля с кривизной.

$\ell$ имеет размерность длины и задаёт масштаб, при котором модификация гравитации становится важной.

---

4. Потенциал: уравнение (2)

Авторы берут

$$V(\phi) = \frac{h^2}{4} \left(\phi^2-v^2\right)^2, \tag{2}$$

причём

$$v=\frac{\mu}{h}.$$

Это стандартный потенциал типа «мексиканская шляпа» для одного вещественного скаляра.

Раскроем скобки:

$$V= \frac{h^2}{4} \left( \phi^4-2v^2\phi^2+v^4 \right).$$

Его производная:

$$V_{,\phi} = h^2\phi(\phi^2-v^2).$$

Вторая производная:

$$V_{,\phi\phi} = h^2(3\phi^2-v^2).$$

Теперь в точке $\phi=0$

$$V_{,\phi\phi}(0) =-h^2v^2.$$

Поскольку

$$hv=\mu,$$

получаем

$$V_{,\phi\phi}(0)=-\mu^2.$$

То есть $\phi=0$ имеет отрицательную массу в квадрате и потому неустойчиво.

Минимумы удовлетворяют

$$V_{,\phi}=0,$$

то есть

$$\phi=0,\qquad\phi=\pm v.$$

$\phi=0$ — максимум, а

$$\phi=\pm v$$

— минимумы.

На бесконечности авторы выбирают

$$\phi_\infty=v. \tag{3}$$

При этом

$$V(v)=0,$$

поэтому этот вакуум не создаёт дополнительной космологической постоянной. (2207.03944v2.pdf)

---

5. Самая важная физика: эффективная масса

Теперь посмотрим, что делает член Гаусса–Бонне.

При фиксированном фоне можно формально определить

$$V_{\rm eff}(\phi) = V(\phi) +\ell^2\mathcal G\,\phi^2.$$

Почему плюс? Потому что в лагранжиане стоит

$$-V-\ell^2\phi^2\mathcal G = -V_{\rm eff}.$$

Около $\phi=0$

$$V(\phi) \simeq V(0)-\frac12\mu^2\phi^2+\ldots$$

и поэтому

$$V_{\rm eff} \simeq V(0) + \frac12 \left( -\mu^2+2\ell^2\mathcal G \right)\phi^2.$$

Значит,

$$\boxed{ m_{\rm eff}^2 = -\mu^2+2\ell^2\mathcal G }.$$

Вот центральная формула всей статьи.

Далеко от чёрной дыры

$$\mathcal G\rightarrow0,$$

поэтому

$$m_{\rm eff}^2\rightarrow-\mu^2<0.$$

Точка $\phi=0$ неустойчива, и поле выбирает

$$\phi=\pm v.$$

Рядом с чёрной дырой $\mathcal G$ возрастает. Если

$$2\ell^2\mathcal G>\mu^2,$$

то

$$m_{\rm eff}^2>0,$$

и теперь $\phi=0$ становится устойчивым.

Это и есть восстановление симметрии сильной кривизной.

---

6. Радиус перехода $r_{\rm cross}$

Для чёрной дыры Шварцшильда

$$R_{\mu\nu}=0, \qquad R=0.$$

Поэтому инвариант Гаусса–Бонне совпадает с квадратом тензора Римана:

$$\mathcal G = R_{\mu\nu\rho\sigma} R^{\mu\nu\rho\sigma}.$$

Для Шварцшильда

$$\mathcal G = \frac{12r_S^2}{r^6}.$$

Авторы определяют $r_{\rm cross}$ условием, когда два вклада в эффективную массу сравниваются:

$$\mu^2 \simeq 2\ell^2\mathcal G(r_{\rm cross}). \tag{5}$$

Подставляем:

$$\mu^2 = 2\ell^2 \frac{12r_S^2}{r_{\rm cross}^6}.$$

То есть

$$\mu^2 = \frac{24\ell^2r_S^2}{r_{\rm cross}^6}.$$

Отсюда

$$r_{\rm cross}^6 = \frac{24\ell^2r_S^2}{\mu^2}.$$

Берём корень шестой степени:

$$r_{\rm cross} = 24^{1/6} \left( \frac{\ell r_S}{\mu} \right)^{1/3}.$$

Поскольку

$$24^{1/6}\simeq1.7,$$

получаем

$$\boxed{ r_{\rm cross} \simeq 1.7 \left( \frac{\ell r_S}{\mu} \right)^{1/3} }. \tag{6}$$

Это один из важнейших результатов.

Чтобы переход происходил снаружи горизонта,

$$r_{\rm cross}\gtrsim r_S.$$

Подставляем оценку:

$$1.7 \left( \frac{\ell r_S}{\mu} \right)^{1/3} \gtrsim r_S.$$

Возводим в куб:

$$1.7^3\frac{\ell r_S}{\mu} \gtrsim r_S^3.$$

Следовательно,

$$\mu \lesssim 1.7^3 \frac{\ell}{r_S^2}.$$

Так как

$$1.7^3\simeq4.9,$$

получаем

$$\boxed{ \mu\lesssim 5\frac{\ell}{r_S^2} }. \tag{7}$$

Именно это условие говорит: масса поля не должна быть слишком большой, иначе обычный потенциал удерживает $\phi$ возле $v$ вплоть до горизонта и никакой стенки не возникает. (2207.03944v2.pdf)

---

7. Приближение пробного поля

Сначала авторы игнорируют обратное влияние $\phi$ на метрику.

То есть считают, что пространство-время точно шварцшильдовское, а поле $\phi$ распространяется на этом фиксированном фоне.

Такое приближение называется test field approximation. Оно особенно хорошо работает при

$$\ell\ll r_S$$

или

$$v\ll1.$$

(2207.03944v2.pdf)

Метрика записывается как

$$ds^2 = -e^\nu dt^2 + e^\lambda dr^2 +r^2d\Omega^2, \tag{9}$$

а для Шварцшильда

$$e^\nu=e^{-\lambda} = 1-\frac{r_S}{r}. \tag{10}$$

Обозначим

$$f(r)=1-\frac{r_S}{r}.$$

Тогда

$$e^\nu=f, \qquad e^\lambda=\frac1f.$$

---

8. Вывод скалярного уравнения

Варьируем действие по $\phi$.

Кинетический член даёт

$$\delta \left[ -\frac12(\nabla\phi)^2 \right] \longrightarrow \Box\phi.$$

Потенциал даёт

$$-V_{,\phi}.$$

Гаусс–Бонне-связь

$$-\ell^2\phi^2\mathcal G$$

даёт

$$-2\ell^2\phi\mathcal G.$$

Поэтому ковариантное уравнение:

$$\boxed{ \Box\phi - V_{,\phi} - 2\ell^2\mathcal G\,\phi =0 }.$$

Для статического $\phi(r)$

$$\Box\phi = \frac1{\sqrt{-g}} \partial_r \left( \sqrt{-g}g^{rr}\phi' \right).$$

Для шварцшильдовской метрики

$$\sqrt{-g}=r^2\sin\theta, \qquad g^{rr}=f.$$

Поэтому

$$\Box\phi = \frac1{r^2} \frac{d}{dr} \left( r^2f\phi' \right).$$

Раскрываем производную:

$$\frac1{r^2} \left[ r^2f\phi'' + (2rf+r^2f')\phi' \right].$$

То есть

$$\Box\phi = f\phi'' + \left( \frac{2f}{r}+f' \right)\phi'.$$

Уравнение становится

$$f\phi'' + \left( \frac{2f}{r}+f' \right)\phi' - V_{,\phi} - 2\ell^2\mathcal G\phi =0.$$

Делим на $f$:

$$\phi'' + \left( \frac{f'}f+\frac2r \right)\phi' - \frac1fV_{,\phi} - \frac{2\ell^2\mathcal G}{f}\phi =0.$$

Поскольку

$$\nu'=\frac{f'}f,$$

получаем форму уравнения (11).

Можно также непосредственно подставить

$$\mathcal G=\frac{12r_S^2}{r^6},$$

и тогда

$$\boxed{ \phi'' + \left( \nu'+\frac2r \right)\phi' - e^{-\nu}V_{,\phi} - \frac{24\ell^2r_S^2}{r^6} e^{-\nu}\phi =0. }$$

Это физически гораздо прозрачнее громоздкой записи (11) в статье. (2207.03944v2.pdf)

---

9. Откуда берётся граничное условие на горизонте

На горизонте

$$r=r_S$$

имеем

$$f(r_S)=0.$$

Поэтому в уравнении появляются члены $\sim1/f$, которые потенциально расходятся.

Положим

$$r=r_S+\epsilon.$$

При малом $\epsilon$

$$f = 1-\frac{r_S}{r_S+\epsilon} \simeq \frac{\epsilon}{r_S}.$$

Следовательно,

$$\frac1f \simeq \frac{r_S}{\epsilon}$$

и

$$\frac{f'}f \simeq \frac1\epsilon.$$

Берём только сингулярные части уравнения:

$$\frac{\phi'_S}{\epsilon} - \frac{r_SV_{,\phi}(\phi_S)}{\epsilon} - \frac{24\ell^2\phi_S}{r_S^3\epsilon} \simeq0.$$

Для регулярного решения коэффициент перед $1/\epsilon$ должен обратиться в ноль:

$$\boxed{ \phi'_S - \frac{24\ell^2}{r_S^3}\phi_S - r_SV_{,\phi}(\phi_S) =0. } \tag{13}$$

Это не дополнительный произвольный физический постулат. Это просто условие отсутствия сингулярности скалярного профиля на горизонте.

На бесконечности второе условие:

$$\phi(\infty)=v.$$

---

10. Почему в статье три режима

Авторы выделяют

$$\begin{cases} \ell\ll r_S,&\text{Case I},\\[2mm] \ell\gg r_S,\quad \mu\ll1/\sqrt{\ell r_S}, &\text{Case II},\\[2mm] \ell\gg r_S,\quad \mu\gg1/\sqrt{\ell r_S}, &\text{Case III}. \end{cases} \tag{15}$$

(2207.03944v2.pdf)

Характерная длина

$$\boxed{L_{\rm GB}\sim\sqrt{\ell r_S}}$$

появляется естественным образом из скалярного уравнения.

Вторая длина —

$$L_\mu\sim\mu^{-1},$$

комптоновская длина скаляра.

Поэтому Cases II и III фактически отличаются тем, какая длина больше:

Case II

$$\mu^{-1}\gg\sqrt{\ell r_S}.$$

Поле лёгкое, потенциал меняется относительно медленно, стенка широкая.

Case III

$$\mu^{-1}\ll\sqrt{\ell r_S}.$$

Собственный потенциал сильнее и заставляет поле быстрее перейти в $v$, поэтому стенка становится гораздо тоньше.

---

11. Case I: слабая связь

Если

$$\ell\ll r_S,$$

поле почти везде близко к

$$\phi=v.$$

Пишем

$$\phi=v+\delta\phi, \qquad |\delta\phi|\ll v.$$

Интересно, почему квадратичная связь становится фактически линейной.

Раскроем

$$\phi^2 = (v+\delta\phi)^2 = v^2+2v\delta\phi+\delta\phi^2.$$

Следовательно,

$$-\ell^2\phi^2\mathcal G = -\ell^2v^2\mathcal G -2\ell^2v\delta\phi\,\mathcal G +\ldots$$

Первый член — постоянный коэффициент перед $\mathcal G$. В четырёх измерениях интеграл чистого $\mathcal G$ является топологическим и локальные уравнения движения не меняет.

Поэтому ведущий динамический эффект:

$$-2\ell^2v\delta\phi\,\mathcal G,$$

то есть фактически линейная связь скаляра с Гауссом–Бонне.

Авторы получают вдали

$$\phi\simeq v \left[ 1- \frac{10A\ell^2}{r_Sr} e^{-\sqrt2\mu r} \right], \tag{16}$$

где $A\sim1$.

Главное здесь не число $10A$, а зависимость

$$\delta\phi\propto\ell^2.$$

То есть эффект действительно пертурбативный.

---

12. Аналитический вывод в приложении

Теперь самая математически интересная часть статьи.

Авторы делят пространство на три области:

1. непосредственно около горизонта;
2. промежуточная область;
3. дальняя область.

(2207.03944v2.pdf)

---

13. Область I: непосредственно около горизонта

При

$$r-r_S\ll r_S \tag{59}$$

скалярное уравнение упрощается до

$$\phi'' + \frac{\phi'}{r-r_S} - \frac{24\ell^2} {r_S^3(r-r_S)} \phi =0. \tag{60}$$

Введём

$$y=r-r_S.$$

Тогда

$$\phi_{yy} + \frac1y\phi_y - \frac{24\ell^2}{r_S^3y}\phi =0.$$

Теперь зададим

$$z=2\sqrt{ \frac{24\ell^2}{r_S^3}y }.$$

То есть

$$z= \sqrt{96}\, \frac{\ell}{r_S} \sqrt{\frac r{r_S}-1}.$$

При такой замене уравнение превращается в

$$\frac{d^2\phi}{dz^2} + \frac1z \frac{d\phi}{dz} -\phi=0.$$

Это модифицированное уравнение Бесселя порядка ноль.

Общее решение:

$$\phi = A I_0(z)+B K_0(z).$$

Но

$$K_0(z)$$

логарифмически расходится при $z\to0$, то есть на горизонте.

Поэтому регулярность требует

$$B=0.$$

И получаем

$$\boxed{ \phi= \phi_S I_0 \left[ \sqrt{96} \frac{\ell}{r_S} \sqrt{\frac r{r_S}-1} \right]. } \tag{61}$$

Это важный результат: при большой $\ell/r_S$ функция $I_0$ растёт чрезвычайно быстро по мере удаления от горизонта. Поэтому, чтобы далеко от горизонта поле не стало огромным, исходное значение

$$\phi_S$$

должно быть экспоненциально маленьким.

Так математически появляется восстановление симметрии:

$$\phi_S\approx0.$$

(2207.03944v2.pdf)

---

14. Промежуточная область

При

$$r_S\ll r\ll r_{\rm cross} \tag{62}$$

можно использовать ньютоновское приближение и пренебречь собственным потенциалом по сравнению с членом Гаусса–Бонне.

Получаем

$$\boxed{ \phi'' +\frac2r\phi' -\frac{24\ell^2r_S^2}{r^6}\phi =0. } \tag{63}$$

Теперь покажу, откуда возникает странный бесселевский индекс $1/4$.

Введём

$$x= \frac{\sqrt6\,\ell r_S}{r^2}.$$

Тогда

$$\frac{dx}{dr} = -\frac{2x}{r}.$$

Для любой функции $\phi(x)$

$$\phi' = -\frac{2x}{r}\phi_x.$$

Ещё одна производная:

$$\phi'' = \frac{4x^2}{r^2}\phi_{xx} + \frac{6x}{r^2}\phi_x.$$

Следовательно,

$$\phi''+\frac2r\phi' = \frac{4x^2}{r^2}\phi_{xx} + \frac{2x}{r^2}\phi_x.$$

Кроме того,

$$x^2 = \frac{6\ell^2r_S^2}{r^4}.$$

Поэтому

$$\frac{24\ell^2r_S^2}{r^6} = \frac{4x^2}{r^2}.$$

Уравнение (63) становится

$$4x^2\phi_{xx} + 2x\phi_x - 4x^2\phi=0.$$

Делим на $4$:

$$x^2\phi_{xx} + \frac12x\phi_x -x^2\phi=0.$$

Теперь положим

$$\phi=x^{1/4}y(x).$$

После подстановки получаем

$$x^2y'' +xy' - \left( x^2+\frac1{16} \right)y=0.$$

Это модифицированное уравнение Бесселя

$$x^2y''+xy'-(x^2+\nu^2)y=0$$

с

$$\nu=\frac14.$$

Поэтому

$$y=C_1K_{1/4}(x)+C_2I_{1/4}(x).$$

И именно отсюда появляется формула (64) статьи. (2207.03944v2.pdf)

---

15. Далеко от чёрной дыры

При

$$r\gg r_{\rm cross}$$

эффект Гаусса–Бонне слаб.

Пишем

$$\phi=v+\tilde\phi, \qquad |\tilde\phi|\ll v.$$

Нужно линеаризовать потенциал.

Мы уже получили

$$V_{,\phi\phi}(v)=2\mu^2.$$

Следовательно,

$$V_{,\phi} \simeq 2\mu^2\tilde\phi.$$

И уравнение становится

$$\tilde\phi'' +\frac2r\tilde\phi' - 2\mu^2\tilde\phi=0. \tag{66}$$

Это обычное массивное радиальное уравнение.

Общее решение:

$$\tilde\phi = \frac{ A e^{+\sqrt2\mu r} + D e^{-\sqrt2\mu r} }{r}.$$

Первый член растёт на бесконечности, поэтому

$$A=0.$$

Получаем

$$\boxed{ \phi= v+ \frac{D}{r} e^{-\sqrt2\mu r}. } \tag{67}$$

То есть поле возвращается к вакууму $v$ по закону Юкавы. (2207.03944v2.pdf)

---

16. Case II: большая $\ell$, маленькая $\mu$

В этом случае около горизонта аргумент бесселевской функции большой:

$$x= \frac{\sqrt6\ell r_S}{r^2}\gg1.$$

Используем

$$K_\nu(x) \simeq \sqrt{\frac{\pi}{2x}}e^{-x},$$ $$I_\nu(x) \simeq \frac{e^x}{\sqrt{2\pi x}}.$$

Регулярное глобальное решение требует сильного подавления растущей $I_{1/4}$-компоненты.

В итоге поле имеет характерную зависимость

$$\phi \propto x^{-1/4}e^{-x}.$$

Так как

$$x= \frac{\sqrt6\ell r_S}{r^2},$$

получаем

$$\boxed{ \phi \sim v \left( \frac{r^2}{\sqrt6\ell r_S} \right)^{1/4} \exp \left[ -\frac{\sqrt6\ell r_S}{r^2} \right]. } \tag{17}$$

Поэтому около горизонта поле экспоненциально мало.

Затем при

$$r\sim\sqrt{\ell r_S}$$

аргумент $x\sim1$, экспоненциальное подавление исчезает и поле начинает переходить к $v$.

Полный профиль удобно записать как

$$\boxed{ \phi\simeq \frac v2 \left( \frac{\sqrt6\ell r_S}{r^2} \right)^{1/4} K_{1/4} \left( \frac{\sqrt6\ell r_S}{r^2} \right). } \tag{19}$$

Отсюда сразу виден характерный радиус:

$$\boxed{ r_{\rm wall}\sim\sqrt{\ell r_S}. }$$

И поскольку в этой формуле другой длины нет,

$$\boxed{ \delta_{\rm wall}\sim\sqrt{\ell r_S}. }$$

То есть в Case II стенка сравнительно широкая. (2207.03944v2.pdf)

---

17. Что показывают рисунки 1 и 2

На рис. 1, стр. 6 изображено

$$\frac{\phi(r)}{v}.$$

В слабом режиме оранжевая кривая почти всё время близка к единице.

В сильном режиме поле около горизонта близко к нулю, а потом резко поднимается к единице.

Это непосредственно и есть стенка.

На рис. 2, стр. 7 показана производная

$$\frac{\phi'}v.$$

Там стенка превращается в пик.

Поэтому очень удобно определять:

• положение стенки — положение максимума $\phi'$;
• толщину стенки — ширину этого пика.

Авторы показывают, что при увеличении $\mu$ стенка смещается к горизонту и становится уже. (2207.03944v2.pdf)

---

18. Case III: большая $\mu$

Теперь

$$\mu\sqrt{\ell r_S}\gg1.$$

Здесь собственный потенциал становится существенным ещё до достижения масштаба $\sqrt{\ell r_S}$.

Поэтому главный радиус стенки теперь определяется не условием $x\sim1$, а именно сменой знака эффективной массы:

$$m_{\rm eff}^2=0.$$

Значит,

$$\boxed{ r_{\rm wall} \sim r_{\rm cross} \sim 1.7 \left( \frac{\ell r_S}{\mu} \right)^{1/3}. } \tag{22}$$

При росте $\mu$

$$r_{\rm cross}\propto\mu^{-1/3}$$

уменьшается.

Поэтому стенка приближается к горизонту. (2207.03944v2.pdf)

---

19. Откуда берётся необычная ширина стенки, уравнение (23)

Это одна из наиболее красивых частей статьи.

Вблизи стенки уравнение приблизительно имеет вид

$$\phi'' + \left( \mu^2- \frac{24\ell^2r_S^2}{r^6} \right)\phi =0. \tag{86}$$

Обозначим

$$q(r) = \mu^2- \frac{24\ell^2r_S^2}{r^6}.$$

По определению $r_{\rm cross}$

$$q(r_{\rm cross})=0.$$

Разложим около этой точки:

$$q(r) \simeq q'(r_{\rm cross}) (r-r_{\rm cross}).$$

Производная:

$$q'(r) = 144 \frac{\ell^2r_S^2}{r^7}.$$

Поэтому

$$q(r) \simeq \frac{144\ell^2r_S^2} {r_{\rm cross}^7} (r-r_{\rm cross}). \tag{87}$$

Уравнение:

$$\phi'' + \kappa (r-r_{\rm cross})\phi=0,$$

где

$$\kappa= \frac{144\ell^2r_S^2}{r_{\rm cross}^7}.$$

Вводим безразмерную координату

$$x= \kappa^{1/3} (r-r_{\rm cross}).$$

Тогда

$$\frac{d}{dr} = \kappa^{1/3}\frac{d}{dx},$$ $$\frac{d^2}{dr^2} = \kappa^{2/3}\frac{d^2}{dx^2}.$$

Уравнение принимает вид

$$\frac{d^2\phi}{dx^2} +x\phi=0. \tag{88}$$

Это уравнение Эйри.

Характерная длина перехода соответствует

$$\Delta x\sim1.$$

Значит,

$$\delta_{\rm wall} \sim \kappa^{-1/3}.$$

То есть

$$\delta_{\rm wall} \sim \left( \frac{r_{\rm cross}^7} {144\ell^2r_S^2} \right)^{1/3}. \tag{90}$$

Теперь

$$r_{\rm cross} \propto \left( \frac{\ell r_S}{\mu} \right)^{1/3}.$$

Следовательно,

$$r_{\rm cross}^7 \propto \ell^{7/3} r_S^{7/3} \mu^{-7/3}.$$

Делим на $\ell^2r_S^2$:

$$\propto \ell^{1/3} r_S^{1/3} \mu^{-7/3}.$$

И берём кубический корень:

$$\boxed{ \delta_{\rm wall} \sim \frac{(\ell r_S)^{1/9}} {\mu^{7/9}}. } \tag{23}$$

Вот почему для большой $\mu$ стенка может быть намного тоньше радиуса Шварцшильда. (2207.03944v2.pdf)

---

20. Теперь включаем обратное влияние на геометрию

До этого метрика была фиксированной.

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

В разделе 4 авторы решают одновременно уравнения для

$$\nu(r),\qquad \lambda(r),\qquad \phi(r).$$

Получаются уравнения (24), (25), (26) и (27). (2207.03944v2.pdf)

Четыре написанных уравнения не являются независимыми: из-за тождеств Бианки достаточно трёх.

Удобно использовать:

• $tt$-компоненту Эйнштейновых уравнений;
• $rr$-компоненту;
• скалярное уравнение.

---

21. Граничные условия полной системы

На горизонте требуется

$$e^\nu\rightarrow0,$$ $$e^{-\lambda}\rightarrow0. \tag{28}$$

Авторы раскладывают

$$e^\nu = a_1(r-r_S)+\ldots,$$ $$e^{-\lambda} = b_1(r-r_S)+\ldots,$$ $$\phi = \phi_S + \phi'_S(r-r_S)+\ldots. \tag{29}$$

Это стандартное разложение регулярного невырожденного горизонта.

---

22. Вывод формулы для $b_1$

Возьмём уравнение (24) и оставим наиболее сингулярные части.

Поскольку

$$e^{-\lambda} \simeq b_1(r-r_S),$$

то

$$e^\lambda \simeq \frac1{b_1(r-r_S)}.$$

А

$$\lambda' \simeq -\frac1{r-r_S}.$$

После подстановки коэффициент перед

$$\frac1{r-r_S}$$

должен обращаться в ноль.

Это даёт

$$-b_1 \left( r_S-8\ell^2\phi_S\phi'_S \right) + 1-r_S^2V(\phi_S) =0.$$

Следовательно,

$$\boxed{ b_1 = \frac{ 1-r_S^2V(\phi_S) }{ r_S-8\ell^2\phi_S\phi'_S }. } \tag{30}$$

Именно так условие регулярного горизонта связывает наклон метрики с полем на горизонте.

Авторы затем аналогичным способом используют остальные уравнения и получают граничное условие (32) для $\phi'_S$. (2207.03944v2.pdf)

---

23. Как решается численная задача

Это двухточечная краевая задача:

на горизонте задана регулярность, а на бесконечности нужно

$$\phi\rightarrow v.$$

Авторы используют два метода.

Shooting method

Стартуют около горизонта и выбирают некоторое

$$\phi_S.$$

Интегрируют наружу.

Если на бесконечности получилось не $v$, меняют $\phi_S$.

Повторяют, пока

$$\phi(\infty)=v.$$

Relaxation method

Задают условия сразу на двух границах, а алгоритм ищет весь профиль целиком.

Авторы используют оба метода и получают одинаковые решения. (2207.03944v2.pdf)

---

24. Что показывает полный численный расчёт

Рис. 3 демонстрирует главное:

стенка не является артефактом приближения пробного поля.

Когда метрика становится динамической, тот же эффект остаётся.

Причём стенка может заметно изменить геометрию.

На рис. 4 и 5 показано отклонение метрических функций от шварцшильдовских.

Чем стенка тоньше, тем сильнее локальное отклонение метрики от ОТО при фиксированных $\ell$ и $v$. Авторы подчёркивают это и в заключении. (2207.03944v2.pdf)

Это физически естественно: тонкая стенка означает большой градиент

$$\phi',$$

а энергия градиента содержит

$$(\phi')^2.$$

Поэтому тонкая стенка может сильно гравитировать даже при не очень большом $v$.

---

25. Но решения существуют не для любых параметров

Это важное ограничение работы.

Если обратное влияние скаляра становится слишком сильным, дифференциальные уравнения теряют регулярность.

После исключения $\nu''$ скалярное уравнение имеет структуру

$$\mathcal D(r)\phi'' +\ldots=0.$$

Если

$$\mathcal D(r)=0,$$

коэффициент перед старшей производной исчезает.

Тогда обычная краевая задача перестаёт быть хорошо определённой.

Авторы получают приблизительное условие

$$\frac{ 128\ell^4\phi^2r_S^2 }{ r^5(r+16\ell^2\phi\phi') } \lesssim1. \tag{34}$$

В наиболее интересной области около $r_{\rm cross}$

$$\phi\simeq v, \qquad \phi'\simeq0.$$

Поэтому

$$\frac{ 128\ell^4v^2r_S^2 }{ r_{\rm cross}^6 } \lesssim1.$$

Но из определения $r_{\rm cross}$

$$r_{\rm cross}^6 \simeq \frac{24\ell^2r_S^2}{\mu^2}.$$

Подставляем:

$$\frac{128}{24} \ell^2\mu^2v^2 \lesssim1.$$

То есть

$$v \lesssim \frac1{\sqrt{128/24}} \frac1{\mu\ell}.$$

Численно

$$\frac1{\sqrt{128/24}}\simeq0.43.$$

И статья записывает

$$\boxed{ v\lesssim \frac{0.4}{\mu\ell}. } \tag{35}$$

Это ограничивает максимально возможную обратную реакцию поля. (2207.03944v2.pdf)

---

26. Космология: почему знак связи так важен

Для плоской FLRW-Вселенной

$$ds^2=-dt^2+a^2(t)d\mathbf x^2$$

инвариант Гаусса–Бонне равен

$$\boxed{ \mathcal G = 24H^2\frac{\ddot a}{a} = 24H^2(H^2+\dot H). } \tag{36}$$

(2207.03944v2.pdf)

Во время инфляции

$$H\simeq\text{const},$$ $$\frac{\ddot a}{a}\simeq H^2,$$

поэтому

$$\mathcal G\simeq24H^4.$$

Вклад в эффективную массу:

$$2\ell^2\mathcal G = 48\ell^2H^4.$$

То есть

$$\boxed{ m_{\rm infl}^2 = 48\ell^2H_{\rm infl}^4. } \tag{37}$$

Если

$$\ell^2<0,$$

масса отрицательная — огромная тахионная неустойчивость.

Если же

$$\ell^2>0,$$

как в этой статье,

$$m_{\rm infl}^2>0.$$

Поле быстро прижимается к

$$\phi\simeq0.$$

Таким образом, выбранный авторами знак хорош для инфляции.

---

27. Почему возникает проблема в радиационную эпоху

В радиационную эпоху

$$a(t)\propto t^{1/2}.$$

Поэтому

$$H=\frac1{2t}.$$

Кроме того,

$$\frac{\ddot a}{a} = -\frac1{4t^2}.$$

Тогда

$$\mathcal G = 24 \frac1{4t^2} \left( -\frac1{4t^2} \right) = -\frac{3}{2t^4}.$$

Вклад в массу:

$$2\ell^2\mathcal G = -\frac{3\ell^2}{t^4}.$$

Теперь он отрицательный.

Однородное скалярное уравнение

$$\ddot\phi+3H\dot\phi+m_{\rm eff}^2\phi=0$$

даёт

$$\boxed{ \ddot\phi + \frac{3}{2t}\dot\phi - \frac{3\ell^2}{t^4}\phi =0. } \tag{39}$$

То есть та самая связь, которая стабилизировала поле при инфляции, становится тахионной в замедляющейся Вселенной. (2207.03944v2.pdf)

---

28. Решение космологического уравнения

Уравнение (39) снова сводится к модифицированным функциям Бесселя.

Вводится переменная приблизительно

$$z=\frac{\sqrt3\,\ell}{t}$$

и выносится множитель

$$\left(\frac{\ell}{t}\right)^{1/4}.$$

Получается

$$\boxed{ \phi(t) = C_1 \left(\frac{\ell}{t}\right)^{1/4} I_{1/4} \left( \frac{\sqrt3\ell}{t} \right) + C_2 \left(\frac{\ell}{t}\right)^{1/4} K_{1/4} \left( \frac{\sqrt3\ell}{t} \right). } \tag{40}$$

На первый взгляд $K$ при большом аргументе мал,

$$K_\nu(z)\sim e^{-z}.$$

Но важно, что коэффициент $C_2$, получаемый из начальных условий при $t=t_{\rm reh}$, содержит обратный экспоненциальный множитель

$$e^{+\sqrt3\ell/t_{\rm reh}}.$$

В результате поле после reheating может получить огромный рост.

Это и есть проблема радиационной эпохи. (2207.03944v2.pdf)

---

29. Ограничение на температуру reheating

Чтобы экспонента не стала катастрофической, авторы требуют примерно

$$\frac{\ell}{t_{\rm reh}} \lesssim100.$$

В радиационную эпоху

$$H=\frac1{2t}.$$

А стандартная связь температуры и Хаббла:

$$H(T) = \sqrt{ \frac{\pi^2g_*(T)}{90} } \frac{T^2}{M_{\rm Pl}}. \tag{43}$$

Комбинируя это с условием выше, авторы получают

$$\boxed{ T_{\rm reh} \lesssim 10 \left[ \frac{90}{\pi^2g_*} \right]^{1/4} \sqrt{\frac{M_{\rm Pl}}{\ell}}. } \tag{44}$$

Для

$$\ell\sim r_{S,\odot}$$

получается примерно

$$\boxed{ T_{\rm reh}\lesssim3\,{\rm GeV}. } \tag{45}$$

В абстракте и заключении это округляется как ограничение порядка $1$ GeV; речь идёт именно об оценке по порядку величины, а не о точном пороге. (2207.03944v2.pdf)

---

30. Космологические domain walls

Есть ещё одна проблема.

Во время инфляции поле загоняется к

$$\phi\simeq0.$$

После этого, когда симметрия снова нарушается, разные причинно несвязанные области могут выбрать

$$+v$$

или

$$-v.$$

Между ними возникают обычные космологические доменные стенки.

Для потенциала

$$V=\frac{h^2}{4}(\phi^2-v^2)^2$$

их натяжение равно

$$\boxed{ \sigma = \frac{2\sqrt2}{3}hv^3. } \tag{46}$$

Это можно вывести напрямую:

для статической kink-стенки

$$\sigma = \int_{-v}^{v} d\phi\,\sqrt{2V(\phi)}.$$

Имеем

$$\sqrt{2V} = \frac h{\sqrt2}(v^2-\phi^2)$$

между $-v$ и $v$.

Поэтому

$$\sigma = \frac h{\sqrt2} \int_{-v}^{v} (v^2-\phi^2)d\phi.$$

Интеграл:

$$\left[ v^2\phi-\frac{\phi^3}{3} \right]_{-v}^{v} = \frac43v^3.$$

Значит,

$$\sigma = \frac{4}{3\sqrt2}hv^3 = \frac{2\sqrt2}{3}hv^3.$$

Совпадает с (46).

Космологические ограничения на такие стены получаются очень жёсткими. Авторы отмечают, что проблему можно убрать небольшим явным нарушением $Z_2$, тогда два вакуума перестают быть точно вырожденными и стены со временем исчезают. (2207.03944v2.pdf)

---

31. Возможное перепроизводство тёмной материи

После перехода в нарушенную фазу поле может осциллировать около

$$\phi=v.$$

Такие когерентные осцилляции массивного скаляра ведут себя как холодная материя.

Осцилляции начинаются, когда

$$H_*\sim\mu.$$

Используя

$$H(T)\propto\frac{T^2}{M_{\rm Pl}},$$

получаем

$$\boxed{ T_* \simeq \left( \frac{90}{\pi^2g_*} \right)^{1/4} \sqrt{\mu M_{\rm Pl}}. } \tag{48}$$

Начальная амплитуда порядка

$$v=\frac{\mu}{h}.$$

Энергетическая плотность осцилляций порядка

$$\rho_{\rm osc} \sim \mu^2v^2.$$

Так как

$$v^2=\frac{\mu^2}{h^2},$$

то

$$\rho_{\rm osc} \sim \frac{\mu^4}{h^2}.$$

После начала осцилляций

$$\rho_{\rm osc}\propto a^{-3},$$

поэтому

$$\boxed{ \rho_{\rm osc}(t) \simeq \frac{\mu^4}{h^2} \left( \frac{a_*}{a} \right)^3. } \tag{49}$$

Требование не перепроизвести наблюдаемую материю даёт сильное ограничение

$$\boxed{ \mu \lesssim 2 \left( \frac{M_{\rm Pl}}v \right)^4 10^{-27}\ {\rm eV}. } \tag{51}$$

Авторы предлагают простой выход: разрешить $\phi$ распадаться в более лёгкие частицы. (2207.03944v2.pdf)

---

32. Скорость гравитационных волн

После GW170817 известно, что скорость гравитации практически совпадает со скоростью света.

Авторы проверяют это явно.

В их модели

$$c_T = 1+ 16\ell^2 \frac{ H\phi\dot\phi -\phi\ddot\phi -\dot\phi^2 }{ M_{\rm Pl}^2 - 16\ell^2H\phi\dot\phi }. \tag{52}$$

В поздней Вселенной поле находится почти в минимуме эффективного потенциала.

Условие минимума:

$$V_{,\phi} + 2\ell^2\mathcal G\phi =0.$$

Для $\phi\neq0$

$$h^2(\phi^2-v^2) + 2\ell^2\mathcal G =0.$$

Отсюда

$$\phi^2 = v^2 - \frac{2\ell^2\mathcal G}{h^2}.$$

Но

$$h^2=\frac{\mu^2}{v^2}.$$

Поэтому

$$\phi^2 = v^2 \left( 1- \frac{2\ell^2\mathcal G}{\mu^2} \right).$$

При малой поправке

$$\sqrt{1-\epsilon} \simeq 1-\frac\epsilon2.$$

Следовательно,

$$\boxed{ \phi \simeq v \left( 1- \frac{\ell^2\mathcal G}{\mu^2} \right). } \tag{53}$$

Это уравнение статьи мы фактически только что вывели полностью.

Поскольку

$$\mathcal G\sim H^4,$$

то

$$\dot{\mathcal G}\sim H^5.$$

Следовательно,

$$|\dot\phi| \sim \frac{\ell^2v}{\mu^2}H^5,$$

с численным коэффициентом порядка $10^2$:

$$|\dot\phi| \sim 100 \frac{\ell^2H^5v}{\mu^2}. \tag{54}$$

Отсюда авторы оценивают

$$\boxed{ |c_T-1| \sim 1000 \frac{\ell^4H^6}{\mu^2} \frac{v^2}{M_{\rm Pl}^2}. } \tag{55}$$

Если взять характерный пример

$$\ell\sim\mu^{-1}\sim3\,{\rm km},$$ $$v\sim0.1M_{\rm Pl},$$ $$H=H_0,$$

получается

$$\boxed{ |c_T-1|\sim10^{-135}. } \tag{56}$$

То есть ограничение GW170817 выполняется с колоссальным запасом. (2207.03944v2.pdf)

---

33. Какова главная новизна статьи

Теперь можно сформулировать очень компактно.

В стандартной спонтанной скаляризации схема примерно такая:

$$\phi=0 \quad\text{далеко}$$

и

$$\phi\neq0 \quad\text{около ЧД}.$$

Здесь всё наоборот:

$$\phi=v \quad\text{далеко}$$

и

$$\phi\simeq0 \quad\text{около ЧД}.$$

Поэтому появляется поверхность перехода

$$0\longrightarrow v,$$

то есть скалярная стенка.

---

34. Физический смысл трёх режимов в одной таблице

Режим	Что происходит
$\ell\ll r_S$	Небольшая поправка к $\phi=v$, стенки практически нет
$\ell\gg r_S,\;\mu\sqrt{\ell r_S}\ll1$	Сильное подавление $\phi$ около горизонта, широкая стенка $r_{\rm wall}\sim\delta_{\rm wall}\sim\sqrt{\ell r_S}$
$\ell\gg r_S,\;\mu\sqrt{\ell r_S}\gg1$	Стенка около $r_{\rm cross}$, становится очень тонкой: $\delta_{\rm wall}\propto\mu^{-7/9}$

---

35. Какие результаты статьи действительно важны

Я бы выделил пять.

Первый: существует новый тип статических чёрных дыр с пространственно разделёнными фазами скалярного поля:

$$\phi\simeq0 \quad\text{внутри}, \qquad \phi\simeq v \quad\text{снаружи}.$$

Второй: положение стенки регулируется кривизной и параметрами

$$r_{\rm wall}\sim \left( \frac{\ell r_S}{\mu} \right)^{1/3}$$

в режиме Case III.

Третий: толщина стенки может быть существенно меньше $r_S$:

$$\delta_{\rm wall} \sim \frac{(\ell r_S)^{1/9}}{\mu^{7/9}}.$$

Четвёртый: стенка может давать заметные локальные отклонения метрики от Шварцшильда, и эффект сохраняется после учёта обратной реакции.

Пятый: модель имеет нетривиальные космологические ограничения. Она лучше обычной скаляризации во время инфляции, но сталкивается с проблемами во время радиационной эпохи, с космологическими domain walls и возможным перепроизводством тёмной материи. Авторы прямо подчёркивают, что для полной феноменологической состоятельности сценарий, вероятно, потребует дополнительных ингредиентов. (2207.03944v2.pdf)

---

36. Что статья НЕ доказывает

Это тоже важно сказать на докладе.

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

Они прямо пишут, что анализ устойчивости — отдельная будущая задача. (2207.03944v2.pdf)

То есть логика такая:

$$\text{регулярное статическое решение существует}$$

не означает автоматически

$$\text{оно устойчиво к любым возмущениям}.$$

Также в работе не моделируется слияние двух таких чёрных дыр и не исследуется вращающийся аналог во всей полноте.

---

37. Как я бы построил ваш доклад

Если выступление примерно на 20–30 минут, логика может быть такой:

1. Проблема: зачем искать отклонения от ОТО около чёрных дыр.
2. Предыдущая идея: спонтанная скаляризация.
3. Новая идея статьи: Higgs-like потенциал + $\phi^2\mathcal G$.
4. Главная физическая формула $$m_{\rm eff}^2=-\mu^2+2\ell^2\mathcal G.$$
5. Crossover radius $$r_{\rm cross}\sim(\ell r_S/\mu)^{1/3}.$$
6. Показать рис. 1 и 2 и объяснить возникновение стенки.
7. Объяснить Cases I, II, III.
8. Для математики показать вывод $$\phi\sim x^{1/4}K_{1/4}(x)$$ и вывод Airy-оценки $$\delta_{\rm wall}\sim(\ell r_S)^{1/9}\mu^{-7/9}.$$
9. Показать рис. 3–5: обратная реакция не убивает эффект.
10. Космология: инфляция хороша, radiation domination проблематична.
11. GW170817: $c_T$ практически равна единице.
12. В заключении подчеркнуть: это не просто «hair», а фазовая оболочка вокруг чёрной дыры.

Самая полезная фраза для устного объяснения всей статьи может звучать так:

Кривизна играет роль среды, которая меняет эффективную массу скалярного поля: далеко от чёрной дыры симметрия нарушена и $\phi=v$, а около горизонта сильная кривизна восстанавливает симметрию и делает $\phi=0$ устойчивым. Между двумя областями возникает скалярная доменная стенка.

Это фактически вся физика статьи в одном предложении.