Итак берем гармоническое решение в вакууме.


И привязываем фоновую метрику в галилеевых координатах:

Вы спросите 
Erleker , почему именно так. Во-первых, все проверки ОТО в слабых полях подтверждаются именно при таком разбиении.
Удаленный наблюдатель находится в хорошем приближении именно в плоском пространстве-времени. 
Доказательство равенства инертной и гравитационных масс также проводится в таком разбиении. 
Можно считать, что инструмент астрономов - удаленные "неподвижные " звезды -  находятся именно в Минковском, введенным при таком разбиении.
Расчеты сближения двух  звезд в системе двойной звезды и потери на гравитационное излучение ведутся именно в гармонических координатах при введении плоской метрики (2a) и хорошо согласуются с наблюдениями.
Теперь переводит обе метрики в сферические координаты.

и плоскую:

Далее, чтобы соблюсти все формальности , заменяем 

 на 

, чтобы (3a)  была в стандартной форме (менять значки не буду).


Далее переводим (5a) к виду Пенлеве по формуле:

Здесь 

 - уже координата времени в Пенлеве . Получаем Пенлеве в сферических, как и вначале:

И плоскую метрику:

(вроде не накосячил, хотя проверьте). 
Теперь смотрите, у нас плоская фоновая метрика получилась в таком виде 
(8a), которая содержит сингулярность на горизонте.
Эта особенность говорит о применимости привязывания фоновой метрики к гармоническому решению. 
Вы не можете данную фоновую метрику распространить под горизонт.
И здесь возникает вторая проблема - неоднозначность. Поскольку 
Erleker, вы сами ссылаетесь на псевдотензор гравитационного поля
по Вайнбергу (7.6.4). Туда входят две метрики - 

 и 

. 
Значит , если вы откажитесь от метрики 
(8a) и замените ее на привычную плоскую в галилеевом виде, привязанном уже к метрике Пенлеве, то получите другой результат при вычислении энергии гравитационных волн (и в других вычислениях, где  используется псевдотензор).
-- 21.07.2017, 11:38 --(некоторые опечатки нашел, исправил).